VZJ
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 25 February 2008
Published in Vadose Zone J 7:238-248 (2008)
DOI: 10.2136/vzj2007.0087
© 2008 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Full Text Free
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Looms, M. C.
Right arrow Articles by Hansen, T. M.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Looms, M. C.
Right arrow Articles by Hansen, T. M.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Looms, M. C.
Right arrow Articles by Hansen, T. M.
Related Collections
Right arrow Downhole/Borehole Methods
Right arrow Ground Penetrating Radar, GPR
Right arrow Variably Saturated Fluid Flow

Identifying Unsaturated Hydraulic Parameters Using an Integrated Data Fusion Approach on Cross-Borehole Geophysical Data

Majken C. Loomsa,*, Andrew Binleyb, Karsten H. Jensena, Lars Nielsena and Thomas M. Hansenc

a Univ. of Copenhagen, Dep. of Geography and Geology, Øster Voldgade 10, DK-1350 Copenhagen K, Denmark
b Lancaster Univ., Dep. of Environmental Science, Lancaster, LA1 4YQ, UK
c Univ. of Copenhagen, Niels Bohr Institute, Juliane Maries Vej 28, DK-2100 Copenhagen Ø, Denmark


Figure 1
View larger version (29K):
[in this window]
[in a new window]

 
FIG. 1. Schematic drawing of the field site setup at Arrenæs, Denmark. The light gray area indicates the infiltration area.

 

Figure 2
View larger version (22K):
[in this window]
[in a new window]

 
FIG. 2. (a) Background moisture content profiles estimated using cross-borehole ground penetrating radar (GPR) and electrical resistivity tomography (ERT); (b) sediment samples from a nearby well (d10, d50 and d90 are the 10th, 50th and 90th percentiles of the grain size distribution; data provided by Copenhagen Energy), and (c) the five-layered model used in the hydrologic forward simulation. The boundaries of the different materials were deduced from (a) and (b).

 

Figure 3
View larger version (35K):
[in this window]
[in a new window]

 
FIG. 3. The methodology used for the integrated data fusion (GPR, ground penetrating radar; ERT, electrical resistivity tomography; {theta}r and {theta}s, residual and saturated moisture content; {alpha} and n, empirical parameters; Ks, saturated hydraulic conductivity).

 

Figure 4
View larger version (24K):
[in this window]
[in a new window]

 
FIG. 4. The misfit values of the one-layered models. For illustration purposes, only the 500 realizations with the lowest misfit values are included in the figure. The parameter ranges of the five realizations with the lowest misfit values are indicated with red; {theta}r is residual moisture content, {alpha} and n are empirical parameters, Ks is saturated hydraulic conductivity.

 

Figure 5
View larger version (32K):
[in this window]
[in a new window]

 
FIG. 5. The simulated zero-offset profile (ZOP) travel times resulting from the one-layered model for three selected days: (a) Day 0, (b) Day 1, and (c) Day 8. The five realizations with the lowest misfit values are highlighted with red. The measured travel time profile is highlighted with green.

 

Figure 6
View larger version (23K):
[in this window]
[in a new window]

 
FIG. 6. The misfit values of the one-layered models using only data collected from Day 8. For illustration purposes, only the 500 realizations with the lowest misfit values are included in the figure. The parameter ranges of the five realizations with the lowest misfit values are indicated with red; an arrow points out the realization with a higher saturated hydraulic conductivity value; {theta}r is residual moisture content, {alpha} and n are empirical parameters, Ks is saturated hydraulic conductivity.

 

Figure 7
View larger version (51K):
[in this window]
[in a new window]

 
FIG. 7. The simulated zero-offset profile (ZOP) travel times resulting from the one-layered model for two selected days: (a) Day 1 and (b) Day 8. Only data collected from Day 8 were used to calculate the misfit values. The five realizations with the lowest misfit values are highlighted with red. The measured travel time profile is highlighted with green. An arrow points out the realization with a higher saturated hydraulic conductivity value.

 

Figure 8
View larger version (54K):
[in this window]
[in a new window]

 
FIG. 8. Travel time profiles from the five-layered model for two selected days: (a) Day 1 and (b) Day 8. The five realizations with the lowest misfit values are highlighted with red. The measured travel time profile is highlighted with green.

 

Figure 9
View larger version (36K):
[in this window]
[in a new window]

 
FIG. 9. A synthetic evaluation of the five-layered model analysis using the two data types combined. For illustration purposes, only the 500 realizations with the lowest misfit values are included in the figure. The parameter ranges of the five realizations with the lowest misfit values are indicated with red and the "true" parameter value is shown with a red diamond.

 

Figure 10
View larger version (40K):
[in this window]
[in a new window]

 
FIG. 10. A synthetic evaluation of the five-layered model analysis using only ground penetrating radar (GPR) data for (a) Day 8 and (b) Days 1 to 8. For illustration purposes, only the 500 realizations with the lowest misfit values are included in the figure. The parameter ranges of the five realizations with the lowest misfit values are indicated with red and the "true" parameter value is shown with a red diamond; {theta}r is residual moisture content, {alpha} and n are empirical parameters, Ks is saturated hydraulic conductivity.

 

Figure 11
View larger version (31K):
[in this window]
[in a new window]

 
FIG. 11. The misfit values of the five-layered model using only electrical resistivity tomography (ERT) data, using only ground penetrating radar (GPR) data, and combining ERT and GPR data in the analysis. For illustration purposes, only the 500 realizations with the lowest misfit are included in the figure. Only the misfit plot of the top three layers of the saturated hydraulic conductivity, Ks, and the empirical parameter n are shown.

 

Figure 12
View larger version (26K):
[in this window]
[in a new window]

 
FIG. 12. Travel time profiles as they developed with time. The five best realizations of the combined analysis (red) are shown along with the measured profiles (green).

 

Figure 13
View larger version (38K):
[in this window]
[in a new window]

 
FIG. 13. Moisture content profiles of the five best realizations of the combined analysis as they developed with time. The background moisture content profile (Day 0) is illustrated with black in the subsequent seven measurement times, while the range of all the realizations of moisture content curves is indicated with gray.

 

Figure 14
View larger version (15K):
[in this window]
[in a new window]

 
FIG. 14. Relative solute concentration profiles of the five best realizations of the combined analysis as they developed with time. The range of all the realizations of solute curves is indicated with gray.

 





HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Soil Science Society of America Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome
Copyright © 2008 by the Soil Science Society of America.