ALBERT

Description: Soil temperature has been recognized as a property that strongly influences a myriad of hydro-biogeochemical processes and reflects how various properties modulate the soil thermal flux. In spite of its importance, our ability to acquire soil temperature data with high spatial and temporal resolution and coverage is limited because of the high cost of equipment, the difficulties of deployment, and the complexities of data management. Here we propose a new strategy that we call distributed temperature profiling (DTP) for improving the characterization and monitoring near-surface thermal properties through the use of an unprecedented number of laterally and vertically distributed temperature measurements. We developed a prototype DTP system, which consists of inexpensive, low-impact, low-power, and vertically resolved temperature probes that independently and autonomously record soil temperature. The DTP system concept was tested by moving sequentially the system across the landscape to identify near-surface permafrost distribution in a discontinuous permafrost environment near Nome, Alaska, during the summertime. Results show that the DTP system enabled successful acquisition of vertically resolved profiles of summer soil temperature over the top 0.8 m at numerous locations. DTP also enabled high-resolution identification and lateral delineation of near-surface permafrost locations from surrounding zones with no permafrost or deep permafrost table locations overlain by a perennially thawed layer. The DTP strategy overcomes some of the limitations associated with – and complements the strengths of – borehole-based soil temperature sensing as well as fiber-optic distributed temperature sensing (FO-DTS) approaches. Combining DTP data with co-located topographic and vegetation maps obtained using unmanned aerial vehicle (UAV) and electrical resistivity tomography (ERT) data allowed us to identify correspondences between surface and subsurface property distribution and in particular between topography, vegetation, shallow soil properties, and near-surface permafrost. Finally, the results highlight the considerable value of the newly developed DTP strategy for investigating the significant variability in and complexity of subsurface thermal and hydrological regimes in discontinuous permafrost regions.

Description: Improving our ability to estimate the parameters that control water and heat fluxes in the shallow subsurface is particularly important due to their strong control on recharge, evaporation and biogeochemical processes. The objectives of this study are to develop and test a new inversion scheme to simultaneously estimate subsurface hydrological, thermal and petrophysical parameters using hydrological, thermal and electrical resistivity tomography (ERT) data. The inversion scheme – which is based on a nonisothermal, multiphase hydrological model – provides the desired subsurface property estimates in high spatiotemporal resolution. A particularly novel aspect of the inversion scheme is the explicit incorporation of the dependence of the subsurface electrical resistivity on both moisture and temperature. The scheme was applied to synthetic case studies, as well as to real datasets that were autonomously collected at a biogeochemical field study site in Rifle, Colorado. At the Rifle site, the coupled hydrological-thermal-geophysical inversion approach well predicted the matric potential, temperature and apparent resistivity with the Nash–Sutcliffe efficiency criterion greater than 0.92. Synthetic studies found that neglecting the subsurface temperature variability, and its effect on the electrical resistivity in the hydrogeophysical inversion, may lead to an incorrect estimation of the hydrological parameters. The approach is expected to be especially useful for the increasing number of studies that are taking advantage of autonomously collected ERT and soil measurements to explore complex terrestrial system dynamics.

Description: This paper compares and integrates different strategies to characterize the variability of end-of-winter snow depth and its relationship to topography in ice-wedge polygon tundra of Arctic Alaska. Snow depth was measured using in situ snow depth probes and estimated using ground-penetrating radar (GPR) surveys and the photogrammetric detection and ranging (phodar) technique with an unmanned aerial system (UAS). We found that GPR data provided high-precision estimates of snow depth (RMSE  =  2.9 cm), with a spatial sampling of 10 cm along transects. Phodar-based approaches provided snow depth estimates in a less laborious manner compared to GPR and probing, while yielding a high precision (RMSE  =  6.0 cm) and a fine spatial sampling (4 cm × 4 cm). We then investigated the spatial variability of snow depth and its correlation to micro- and macrotopography using the snow-free lidar digital elevation map (DEM) and the wavelet approach. We found that the end-of-winter snow depth was highly variable over short (several meter) distances, and the variability was correlated with microtopography. Microtopographic lows (i.e., troughs and centers of low-centered polygons) were filled in with snow, which resulted in a smooth and even snow surface following macrotopography. We developed and implemented a Bayesian approach to integrate the snow-free lidar DEM and multiscale measurements (probe and GPR) as well as the topographic correlation for estimating snow depth over the landscape. Our approach led to high-precision estimates of snow depth (RMSE  =  6.0 cm), at 0.5 m resolution and over the lidar domain (750 m × 700 m).

Description: Quantitative characterization of soil organic carbon (OC) content is essential due to its significant impacts on surface–subsurface hydrological–thermal processes and microbial decomposition of OC, which both in turn are important for predicting carbon–climate feedbacks. While such quantification is particularly important in the vulnerable organic-rich Arctic region, it is challenging to achieve due to the general limitations of conventional core sampling and analysis methods, and to the extremely dynamic nature of hydrological–thermal processes associated with annual freeze–thaw events. In this study, we develop and test an inversion scheme that can flexibly use single or multiple datasets – including soil liquid water content, temperature and electrical resistivity tomography (ERT) data – to estimate the vertical distribution of OC content. Our approach relies on the fact that OC content strongly influences soil hydrological–thermal parameters and, therefore, indirectly controls the spatiotemporal dynamics of soil liquid water content, temperature and their correlated electrical resistivity. We employ the Community Land Model to simulate nonisothermal surface–subsurface hydrological dynamics from the bedrock to the top of canopy, with consideration of land surface processes (e.g., solar radiation balance, evapotranspiration, snow accumulation and melting) and ice–liquid water phase transitions. For inversion, we combine a deterministic and an adaptive Markov chain Monte Carlo (MCMC) optimization algorithm to estimate a posteriori distributions of desired model parameters. For hydrological–thermal-to-geophysical variable transformation, the simulated subsurface temperature, liquid water content and ice content are explicitly linked to soil electrical resistivity via petrophysical and geophysical models. We validate the developed scheme using different numerical experiments and evaluate the influence of measurement errors and benefit of joint inversion on the estimation of OC and other parameters. We also quantify the propagation of uncertainty from the estimated parameters to prediction of hydrological–thermal responses. We find that, compared to inversion of single dataset (temperature, liquid water content or apparent resistivity), joint inversion of these datasets significantly reduces parameter uncertainty. We find that the joint inversion approach is able to estimate OC and sand content within the shallow active layer (top 0.3 m of soil) with high reliability. Due to the small variations of temperature and moisture within the shallow permafrost (here at about 0.6 m depth), the approach is unable to estimate OC with confidence. However, if the soil porosity is functionally related to the OC and mineral content, which is often observed in organic-rich Arctic soil, the uncertainty of OC estimate at this depth remarkably decreases. Our study documents the value of the new surface–subsurface, deterministic–stochastic inversion approach, as well as the benefit of including multiple types of data to estimate OC and associated hydrological–thermal dynamics.

Description: Quantitative characterization of soil organic carbon (OC) content is essential due to its significant impacts on surface–subsurface hydrological-thermal processes and microbial decomposition of OC, which both in turn are important for predicting carbon-climate feedbacks. While such quantification is particularly important in the vulnerable organic-rich Arctic region, it is challenging to achieve due to the general limitations of conventional core sampling and analysis methods, and to the extremely dynamic nature of hydrological-thermal processes associated with annual freeze-thaw events. In this study, we develop and test an inversion scheme that can flexibly use single or multiple datasets, including soil water liquid, temperature and electrical resistivity data (ERT), to estimate the vertical distribution of OC content. We subsequently explore the control of OC on hydrological-thermal behavior. We employ the Community Land Model to simulate nonisothermal surface-subsurface hydrological dynamics from the bedrock to the top of canopy, with consideration of land surface processes and ice/liquid water phase transitions. For inversion, we combine a deterministic and an adaptive Markov chain Monte Carlo (MCMC) optimization algorithm to estimate posterior distributions of desired model parameters. For hydrological-thermal to geophysical variable transformation, the simulated subsurface temperature, liquid and ice water content are explicitly linked to the soil apparent resistivity via petrophysical and geophysical models. We validate the developed scheme using different numerical experiments and evaluate the influence of measurement errors and benefit of joint inversion on the estimation of OC and other parameters. We also quantified the propagation of uncertainty from the estimated parameters to prediction of hydrological-thermal responses. We find that compared to inversion of single dataset (either temperature or liquid or apparent resistivity), joint inversion of these datasets significantly reduces parameter uncertainty. We find that the joint inversion approach is able to estimate OC and sand content within the shallow active layer (0.3 m) with high reliability. Due to the small variations of temperature and moisture within the shallow permafrost (0.6 m), the approach is unable to estimate OC with confidence. However, if the soil porosity is functionally related to the OC and mineral content, the uncertainty of OC estimate at this depth remarkably decreases. Our study documents the value of the new surface-subsurface, deterministic-stochastic inversion approach, as well as the benefit of including multiple types of data to estimate OC and associated hydrological-thermal dynamics.

Description: This paper presents a time-lapse application of electrical methods (Electrical Resistivity Tomography – ERT – and Mise-à-la-Masse – MALM) for monitoring plant roots and their activity (root water uptake) during a controlled infiltration experiment. The use of non-invasive geophysical monitoring is of increasing interest as these techniques provide time-lapse imaging of processes that otherwise can only be measured at few specific spatial locations. The experiment here described was conducted in a vineyard in Bordeaux (France) and was focused on the behaviour of two neighbouring grapevines. The joint application of ERT and MALM has several advantages. While ERT in time-lapse mode is sensitive to changes in soil electrical resistivity and thus to the factors controlling it (mainly soil water content, in this context), MALM uses DC current injected in a tree stem to image where the plant-root system is in effective electrical contact with the soil at locations that are likely to be the same where root water uptake (RWU) takes place. Thus ERT and MALM provide complementary information about the root structure and activity. The experiment shows that the region of likely electrical current sources produced by MALM does not change significantly during the infiltration study time in spite of the strong changes of electrical resistivity caused by changes in soil water content. This fact, together with the evidence that current injection in the soil produces totally different patterns, corroborates the idea that this application of MALM highlights the active root density in the soil. When considering the electrical resistivity changes (as measured by ERT) inside the stationary volume of active roots delineated by MALM, the overall tendency is towards a resistivity increase, which can be linked to a decrease in soil water content caused by root water uptake. On the contrary, when considering the soil volume outside the MALM-derived root water uptake region, the electrical resistivity tends to decrease as an effect of soil water content increase caused by the infiltration. The results are particularly promising, and the method can be applied to a variety of scales including the laboratory scale where direct evidence of roots structure and root water uptake can help corroborate the approach. Once fully validated, the joint use of MALM and ERT can be used as a valuable tool to study the activity of roots under a wide variety of field conditions.

Description: The investigation of plant roots is inherently difficult and often neglected. Being out of sight, roots are often out of mind. Still, roots play a key role in the exchange of mass and energy between soil and the atmosphere, let alone the many practical applications in agriculture. In this paper, we propose a method for roots imaging based on the joint use of two electrical non-invasive methods, Electrical Resistivity Tomography (ERT) and Mise-a-la-Masse (MALM). The approach is based on the key assumption that the plant root system acts as an electrically conductive body, so that injecting electrical current in the plant stem will ultimately result in the injection of current in the subsoil through the root system, and particularly through the root terminations via hair roots. Evidence from field data, showing that voltage distribution is very different whether current is injected in the tree stem or in the ground, strongly supports this hypothesis. The proposed procedure involves a stepwise inversion of both ERT and MALM data that ultimately leads to the identification of electrical resistivity distribution, and of the current-injection root distribution in the three-dimensional soil space. This, in turn, is a proxy to the active (hair) root density in the ground. We tested the proposed procedure on synthetic data and, more importantly, on field data collected in vineyard, where the estimated depth of the root zone proved to be in agreement with literature on similar crops. The proposed non-invasive approach is a step forward towards a better quantification of roots structure and functioning.

Description: The investigation of plant roots is inherently difficult and often neglected. Being out of sight, roots are often out of mind. Nevertheless, roots play a key role in the exchange of mass and energy between soil and the atmosphere, in addition to the many practical applications in agriculture. In this paper, we propose a method for roots imaging based on the joint use of two electrical noninvasive methods: electrical resistivity tomography (ERT) and mise-à-la-masse (MALM). The approach is based on the key assumption that the plant root system acts as an electrically conductive body, so that injecting electrical current into the plant stem will ultimately result in the injection of current into the subsoil through the root system, and particularly through the root terminations via hair roots. Evidence from field data, showing that voltage distribution is very different whether current is injected into the tree stem or in the ground, strongly supports this hypothesis. The proposed procedure involves a stepwise inversion of both ERT and MALM data that ultimately leads to the identification of electrical resistivity (ER) distribution and of the current injection root distribution in the three-dimensional soil space. This, in turn, is a proxy to the active (hair) root density in the ground. We tested the proposed procedure on synthetic data and, more importantly, on field data collected in a vineyard, where the estimated depth of the root zone proved to be in agreement with literature on similar crops. The proposed noninvasive approach is a step forward towards a better quantification of root structure and functioning.