the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Impact of parameter updates on soil moisture assimilation in a 3D heterogeneous hillslope model
Natascha Brandhorst
Insa Neuweiler
Abstract. Models of variably saturated subsurface flow require knowledge of the soil hydraulic parameters. However, the determination of these parameters in heterogeneous soils is not easily feasible and subject to large uncertainties. As the modeled soil moisture is very sensitive to these parameters, especially the saturated hydraulic conductivity, porosity and the parameters describing the retention and relative permeability functions, it is likewise highly uncertain. Data assimilation can be used to handle and reduce both, state and parameter uncertainty. In this work, we apply the ensemble Kalman filter (EnKF) to a three-dimensional heterogeneous hillslope model and investigate the influence of updating the different soil hydraulic parameters on the accuracy of the estimated soil moisture. We further examine the usage of a simplified layered soil structure instead of the fully resolved heterogeneous soil structure in the ensemble. It is shown that the best estimates are obtained when performing a joint update of porosity and van Genuchten parameters and optionally the saturated hydraulic conductivity. The usage of a simplified soil structure gave decent estimates of spatially averaged soil moisture in combination with parameter updates but led to a failure of the EnKF and very poor soil moisture estimates at non-observed locations.
Natascha Brandhorst and Insa Neuweiler
Status: closed
-
RC1: 'Comment on hess-2022-311', Anonymous Referee #1, 26 Oct 2022
The comment was uploaded in the form of a supplement: https://hess.copernicus.org/preprints/hess-2022-311/hess-2022-311-RC1-supplement.pdf
-
AC1: 'Reply on RC1', Natascha Brandhorst, 23 Nov 2022
The comment was uploaded in the form of a supplement: https://hess.copernicus.org/preprints/hess-2022-311/hess-2022-311-AC1-supplement.pdf
-
AC1: 'Reply on RC1', Natascha Brandhorst, 23 Nov 2022
-
RC2: 'Comment on hess-2022-311', Anonymous Referee #2, 01 Dec 2022
This is a very clearly written, well-constructed article that documents a numerical study of data assimilation using an integrated numerical model. I found this article scientifically interesting, as I also think the readership of HESS would. I recommend publication pending minor revisions. I have some general and specific comments listed below.
General comments:
The integrated model used here incorporates overland flow.  However, it's unclear what the role of overland flow was in the study as I don't believe it was used in the DA.  Did the presence of overland runoff modify the soil moisture in some way (beyond acting as a boundary condition at the bottom of the hill slope)? Would the results be essentially the same with a Richards’ only hill slope model?
The authors found that porosity was a particularly sensitive parameter in the DA. This makes sense, as they state (e.g. 615), as this limits the total amount of water available in the soil. However, this sensitivity is likely larger on the wet side of the soil moisture curve, on the dry side other parameters (processes) may play a larger role. The aridity of the simulations is driven by the meteorological forcing used in the experiment, did the authors consider the impact a different forcing dataset might have on these findings?
Specific comments:
Section 3.1.1: Are the random fields for ln(K), ln(alpha), phi and n spatially correlated? Â I realize they are correlated with each other, but are the fields disordered in space or correlated with some spatial correlation structure? Â There is a lot of evidence in the literature demonstrating the spatial correlation of random fields and this should be discussed in the manuscript. Â If the fields are correlated, I recommend including some discussion of the correlation model and associated parameters.
Section 3.1/Table 2: For the numerical simulations were the random fields constrained in some way to prevent non-physical parameter values? Â Or even parameter values that would be outside the range of solution (VG parameters such as n that result in eqs 2, 3 being nondifferentiable).
Citation: https://doi.org/10.5194/hess-2022-311-RC2 -
AC2: 'Reply on RC2', Natascha Brandhorst, 07 Dec 2022
The comment was uploaded in the form of a supplement: https://hess.copernicus.org/preprints/hess-2022-311/hess-2022-311-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Natascha Brandhorst, 07 Dec 2022
Status: closed
-
RC1: 'Comment on hess-2022-311', Anonymous Referee #1, 26 Oct 2022
The comment was uploaded in the form of a supplement: https://hess.copernicus.org/preprints/hess-2022-311/hess-2022-311-RC1-supplement.pdf
-
AC1: 'Reply on RC1', Natascha Brandhorst, 23 Nov 2022
The comment was uploaded in the form of a supplement: https://hess.copernicus.org/preprints/hess-2022-311/hess-2022-311-AC1-supplement.pdf
-
AC1: 'Reply on RC1', Natascha Brandhorst, 23 Nov 2022
-
RC2: 'Comment on hess-2022-311', Anonymous Referee #2, 01 Dec 2022
This is a very clearly written, well-constructed article that documents a numerical study of data assimilation using an integrated numerical model. I found this article scientifically interesting, as I also think the readership of HESS would. I recommend publication pending minor revisions. I have some general and specific comments listed below.
General comments:
The integrated model used here incorporates overland flow.  However, it's unclear what the role of overland flow was in the study as I don't believe it was used in the DA.  Did the presence of overland runoff modify the soil moisture in some way (beyond acting as a boundary condition at the bottom of the hill slope)? Would the results be essentially the same with a Richards’ only hill slope model?
The authors found that porosity was a particularly sensitive parameter in the DA. This makes sense, as they state (e.g. 615), as this limits the total amount of water available in the soil. However, this sensitivity is likely larger on the wet side of the soil moisture curve, on the dry side other parameters (processes) may play a larger role. The aridity of the simulations is driven by the meteorological forcing used in the experiment, did the authors consider the impact a different forcing dataset might have on these findings?
Specific comments:
Section 3.1.1: Are the random fields for ln(K), ln(alpha), phi and n spatially correlated? Â I realize they are correlated with each other, but are the fields disordered in space or correlated with some spatial correlation structure? Â There is a lot of evidence in the literature demonstrating the spatial correlation of random fields and this should be discussed in the manuscript. Â If the fields are correlated, I recommend including some discussion of the correlation model and associated parameters.
Section 3.1/Table 2: For the numerical simulations were the random fields constrained in some way to prevent non-physical parameter values? Â Or even parameter values that would be outside the range of solution (VG parameters such as n that result in eqs 2, 3 being nondifferentiable).
Citation: https://doi.org/10.5194/hess-2022-311-RC2 -
AC2: 'Reply on RC2', Natascha Brandhorst, 07 Dec 2022
The comment was uploaded in the form of a supplement: https://hess.copernicus.org/preprints/hess-2022-311/hess-2022-311-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Natascha Brandhorst, 07 Dec 2022
Natascha Brandhorst and Insa Neuweiler
Natascha Brandhorst and Insa Neuweiler
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
276 | 85 | 14 | 375 | 5 | 4 |
- HTML: 276
- PDF: 85
- XML: 14
- Total: 375
- BibTeX: 5
- EndNote: 4
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1