Geostatistical inference using crosshole ground-penetrating radar

High-resolution tomographic images obtained from crosshole geophysical measurements have the potential to provide valuable information about the geostatistical properties of unsaturated-zone hydrologic-state va riables such as moisture content. Under drained or quasi-steady-state conditions, the moisture content will reflect the variation of the physical properties of the subsurface, which determine the flow patterns in the unsaturated zone. Deterministic least-squares inversion of crosshole groundpenetrating-radar GPR traveltimes result in smooth, minimumvariance estimates of the subsurface radar wave velocity structure, which may diminish the utility of these images for geostatistical inference. We have used a linearized stochastic inversion technique to infer the geostatistical properties of the subsurface radar wave velocity distribution using crosshole GPR traveltimes directly. Expanding on a previous study, we have determined that it is possible to obtain estimates of global variance and mean velocity values of the subsurface as well as the correlation lengths describing the subsurface velocity structures. Accurate estimation of the global variance is crucial if stochastic realizations of the subsurface are used to evaluate the uncertainty of the inversion estimate. We have explored the full potential of the geostatistical inference method using several synthetic models of varying correlation structures and have tested the influence of different assumptions concerning the choice of covariance function and data noise level. In addition, we have tested the methodology on traveltime data collected at a field site in Denmark. There, inferred correlation structures indicate that structural differences exist between two areas located approximately 10 m apart, an observation confirmed by a GPR reflection profile. Furthermore, the inferred values of the subsurface global variance and the mean velocity have been corroborated with moisturecontent measurements, obtained gravimetrically from samples collected at the field site

Mapping Root-Zone Soil Moisture Using a Temperature–Vegetation Triangle Approach with an Unmanned Aerial System: Incorporating Surface Roughness from Structure from Motion

High resolution root-zone soil moisture (SM) maps are important for understanding the spatial variability of water availability in agriculture, ecosystems research and water resources management. Unmanned Aerial Systems (UAS) can flexibly monitor land surfaces with thermal and optical imagery at very high spatial resolution (meter level, VHR) for most weather conditions. We modified the temperature⁻vegetation triangle approach to transfer it from satellite to UAS remote sensing. To consider the effects of the limited coverage of UAS mapping, theoretical dry/wet edges were introduced. The new method was tested on a bioenergy willow short rotation coppice site during growing seasons of 2016 and 2017. We demonstrated that by incorporating surface roughness parameters from the structure-from-motion in the interpretation of the measured land surface-atmosphere temperature gradients, the estimates of SM significantly improved. The correlation coefficient between estimated and measured SM increased from not significant to 0.69 and the root mean square deviation decreased from 0.045 m3∙m−3 to 0.025 m3∙m−3 when considering temporal dynamics of surface roughness in the approach. The estimated SM correlated better with in-situ root-zone SM (15⁻30 cm) than with surface SM (0⁻5 cm) which is an important advantage over alternative remote sensing methods to estimate SM. The optimal spatial resolution of the triangle approach was found to be around 1.5 m, i.e. similar to the length scale of tree-crowns. This study highlights the importance of considering the 3-D fine scale canopy structure, when addressing the links between surface temperature and SM patterns via surface energy balances. Our methodology can be applied to operationally monitor VHR root-zone SM from UAS in agricultural and natural ecosystems

Advancing measurements and representations of subsurface heterogeneity and dynamic processes: towards 4D hydrogeology

Essentially all hydrogeological processes are strongly influenced by the subsurface spatial heterogeneity and the temporal variation of environmental conditions, hydraulic properties, and solute concentrations. This spatial and temporal variability generally leads to effective behaviors and emerging phenomena that cannot be predicted from conventional approaches based on homogeneous assumptions and models. However, it is not always clear when, why, how, and at what scale the 4D (3D + time) nature of the subsurface needs to be considered in hydrogeological monitoring, modeling, and applications. In this paper, we discuss the interest and potential for the monitoring and characterization of spatial and temporal variability, including 4D imaging, in a series of hydrogeological processes: (1) groundwater fluxes, (2) solute transport and reaction, (3) vadose zone dynamics, and (4) surface–subsurface water interactions. We first identify the main challenges related to the coupling of spatial and temporal fluctuations for these processes. We then highlight recent innovations that have led to significant breakthroughs in high-resolution space–time imaging and modeling the characterization, monitoring, and modeling of these spatial and temporal fluctuations. We finally propose a classification of processes and applications at different scales according to their need and potential for high-resolution space–time imaging. We thus advocate a more systematic characterization of the dynamic and 3D nature of the subsurface for a series of critical processes and emerging applications. This calls for the validation of 4D imaging techniques at highly instrumented observatories and the harmonization of open databases to share hydrogeological data sets in their 4D components