Numerical assessment of a dynamical relaxation data assimilation scheme for a catchment hydrological model

Ruud Hurkmans, Claudio Paniconi, Peter A Troch

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

A dynamical relaxation scheme for assimilating observation data into a three-dimensional Richards equation-based numerical model was presented in Paniconi et al. (2003. Advances in Water Resources 26: 161-178). The technique, known as Newtonian relaxation or nudging, relies on a forcing term to reduce the difference between computed and observed values of a state variable such as soil moisture content. The forcing term contains data quality and nudging influence factors, and, importantly, spatio-temporal weighting functions that determine the manner and extent of spreading of a state variable beyond its measurement points and times. In this paper, a series of numerical experiments is run for a small catchment in southern Belgium to investigate the performance of the nudging algorithm. In a first set of runs for a short (10-day) simulation period, the model's sensitivity to influence radii and forcing strength parameters is examined. Here, we find that prediction errors decrease and numerical cost increases for increasing values of the parameters, except for the vertical radius of influence, where intermediate values produced the lowest prediction errors. Based on a compromise between numerical and physical results for this 10-day experiment, 'optimal' values were selected for these parameters, and a second, much longer (8-month) set of simulations is used to explore the effect of the observation frequency on the quality of the assimilation. Here, it is found that prediction errors are lowest at an intermediate frequency, whereas serious numerical difficulties occur at a high frequency. As could be expected, computational costs are highest when frequent observations are assimilated. However, since CPU times increase only very slowly as more observation datasets are added, we can conclude that the nudging method is computationally efficient, albeit quiet sensitive in terms of numerical performance to its parameter settings, and thus requires further study to devise more robust formulations of the algorithm.

Original languageEnglish (US)
Pages (from-to)549-563
Number of pages15
JournalHydrological Processes
Volume20
Issue number3
DOIs
StatePublished - Feb 28 2006
Externally publishedYes

Fingerprint

data assimilation
catchment
prediction
Richards equation
data quality
simulation
moisture content
soil moisture
experiment
water resource
parameter
cost

Keywords

  • Data assimilation
  • Groundwater modelling
  • Newtonian nudging
  • Richards equation
  • Subsurface flow
  • Unsaturated zone

ASJC Scopus subject areas

  • Water Science and Technology

Cite this

Numerical assessment of a dynamical relaxation data assimilation scheme for a catchment hydrological model. / Hurkmans, Ruud; Paniconi, Claudio; Troch, Peter A.

In: Hydrological Processes, Vol. 20, No. 3, 28.02.2006, p. 549-563.

Research output: Contribution to journalArticle

@article{d43b17256a5a42279e0d1f3f0841ccd8,
title = "Numerical assessment of a dynamical relaxation data assimilation scheme for a catchment hydrological model",
abstract = "A dynamical relaxation scheme for assimilating observation data into a three-dimensional Richards equation-based numerical model was presented in Paniconi et al. (2003. Advances in Water Resources 26: 161-178). The technique, known as Newtonian relaxation or nudging, relies on a forcing term to reduce the difference between computed and observed values of a state variable such as soil moisture content. The forcing term contains data quality and nudging influence factors, and, importantly, spatio-temporal weighting functions that determine the manner and extent of spreading of a state variable beyond its measurement points and times. In this paper, a series of numerical experiments is run for a small catchment in southern Belgium to investigate the performance of the nudging algorithm. In a first set of runs for a short (10-day) simulation period, the model's sensitivity to influence radii and forcing strength parameters is examined. Here, we find that prediction errors decrease and numerical cost increases for increasing values of the parameters, except for the vertical radius of influence, where intermediate values produced the lowest prediction errors. Based on a compromise between numerical and physical results for this 10-day experiment, 'optimal' values were selected for these parameters, and a second, much longer (8-month) set of simulations is used to explore the effect of the observation frequency on the quality of the assimilation. Here, it is found that prediction errors are lowest at an intermediate frequency, whereas serious numerical difficulties occur at a high frequency. As could be expected, computational costs are highest when frequent observations are assimilated. However, since CPU times increase only very slowly as more observation datasets are added, we can conclude that the nudging method is computationally efficient, albeit quiet sensitive in terms of numerical performance to its parameter settings, and thus requires further study to devise more robust formulations of the algorithm.",
keywords = "Data assimilation, Groundwater modelling, Newtonian nudging, Richards equation, Subsurface flow, Unsaturated zone",
author = "Ruud Hurkmans and Claudio Paniconi and Troch, {Peter A}",
year = "2006",
month = "2",
day = "28",
doi = "10.1002/hyp.5921",
language = "English (US)",
volume = "20",
pages = "549--563",
journal = "Hydrological Processes",
issn = "0885-6087",
publisher = "John Wiley and Sons Ltd",
number = "3",

}

TY - JOUR

T1 - Numerical assessment of a dynamical relaxation data assimilation scheme for a catchment hydrological model

AU - Hurkmans, Ruud

AU - Paniconi, Claudio

AU - Troch, Peter A

PY - 2006/2/28

Y1 - 2006/2/28

N2 - A dynamical relaxation scheme for assimilating observation data into a three-dimensional Richards equation-based numerical model was presented in Paniconi et al. (2003. Advances in Water Resources 26: 161-178). The technique, known as Newtonian relaxation or nudging, relies on a forcing term to reduce the difference between computed and observed values of a state variable such as soil moisture content. The forcing term contains data quality and nudging influence factors, and, importantly, spatio-temporal weighting functions that determine the manner and extent of spreading of a state variable beyond its measurement points and times. In this paper, a series of numerical experiments is run for a small catchment in southern Belgium to investigate the performance of the nudging algorithm. In a first set of runs for a short (10-day) simulation period, the model's sensitivity to influence radii and forcing strength parameters is examined. Here, we find that prediction errors decrease and numerical cost increases for increasing values of the parameters, except for the vertical radius of influence, where intermediate values produced the lowest prediction errors. Based on a compromise between numerical and physical results for this 10-day experiment, 'optimal' values were selected for these parameters, and a second, much longer (8-month) set of simulations is used to explore the effect of the observation frequency on the quality of the assimilation. Here, it is found that prediction errors are lowest at an intermediate frequency, whereas serious numerical difficulties occur at a high frequency. As could be expected, computational costs are highest when frequent observations are assimilated. However, since CPU times increase only very slowly as more observation datasets are added, we can conclude that the nudging method is computationally efficient, albeit quiet sensitive in terms of numerical performance to its parameter settings, and thus requires further study to devise more robust formulations of the algorithm.

AB - A dynamical relaxation scheme for assimilating observation data into a three-dimensional Richards equation-based numerical model was presented in Paniconi et al. (2003. Advances in Water Resources 26: 161-178). The technique, known as Newtonian relaxation or nudging, relies on a forcing term to reduce the difference between computed and observed values of a state variable such as soil moisture content. The forcing term contains data quality and nudging influence factors, and, importantly, spatio-temporal weighting functions that determine the manner and extent of spreading of a state variable beyond its measurement points and times. In this paper, a series of numerical experiments is run for a small catchment in southern Belgium to investigate the performance of the nudging algorithm. In a first set of runs for a short (10-day) simulation period, the model's sensitivity to influence radii and forcing strength parameters is examined. Here, we find that prediction errors decrease and numerical cost increases for increasing values of the parameters, except for the vertical radius of influence, where intermediate values produced the lowest prediction errors. Based on a compromise between numerical and physical results for this 10-day experiment, 'optimal' values were selected for these parameters, and a second, much longer (8-month) set of simulations is used to explore the effect of the observation frequency on the quality of the assimilation. Here, it is found that prediction errors are lowest at an intermediate frequency, whereas serious numerical difficulties occur at a high frequency. As could be expected, computational costs are highest when frequent observations are assimilated. However, since CPU times increase only very slowly as more observation datasets are added, we can conclude that the nudging method is computationally efficient, albeit quiet sensitive in terms of numerical performance to its parameter settings, and thus requires further study to devise more robust formulations of the algorithm.

KW - Data assimilation

KW - Groundwater modelling

KW - Newtonian nudging

KW - Richards equation

KW - Subsurface flow

KW - Unsaturated zone

UR - http://www.scopus.com/inward/record.url?scp=33644773832&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33644773832&partnerID=8YFLogxK

U2 - 10.1002/hyp.5921

DO - 10.1002/hyp.5921

M3 - Article

AN - SCOPUS:33644773832

VL - 20

SP - 549

EP - 563

JO - Hydrological Processes

JF - Hydrological Processes

SN - 0885-6087

IS - 3

ER -