Improving understanding of Arctic ecosystem climate feedback and parameterization of models that simulate freeze-Thaw dynamics require advances in quantifying soil and snow properties. Due to the significant spatiotemporal variability of soil properties and the limited information provided by point-scale measurements (e.g., cores), geophysical methods hold potential for improving soil and permafrost characterization. In this study, we evaluate the use of a ground-penetrating radar (GPR) to estimate thaw layer thickness, snow depth, and ice-wedge characteristics in an ice-wedge-dominated tundra region near Barrow, AK, USA. To this end, we analyze GPR and point-scale measurements collected along several parallel transects at the end of the growing season and the end of frozen season. In addition, we compare the structural information extracted from the GPR data with electrical resistivity tomography (ERT) information about ice-wedge characteristics. Our study generally highlights the value of GPR data collected in the frozen season, when conditions lead to the improved GPR signal-To-noise ratio, facilitate data acquisition, and reduce acquisition-related ecosystem disturbance relative to growing season. We document for the first time that GPR data collected during the frozen season can provide reliable estimates of active layer thickness and geometry of ice wedges. We find that the ice-wedge geometry extracted from GPR data collected during the frozen season is consistent with ERT data, and that GPR data can be used to constrain the ERT inversion. Consistent with recent studies, we also find that GPR data collected during the frozen season can provide good estimates of snow thickness, and that GPR data collected during the growing season can provide reliable estimate thaw depth. Our quantification of the value of the GPR and ERT data collected during growing and frozen seasons paves the way for coupled inversion of the datasets to improve understanding of permafrost variability.

Permafrost physicochemical parameters play a key role in controlling the response of permafrost carbon to climate change. We studied the physicochemical parameters of permafrost in an Arctic tundra region to evaluate (1) how soil parameters vary with depth and whether and how they are interrelated, (2) whether and how permafrost soil differs from its overlaying active layer, and (3) whether soil property-depth relationships are different across geomorphic features (e.g., low, flat, and high centered polygons). We also explored the possible biogeochemical processes that led to these soil characteristics and how they may affect biogeochemical reactions upon permafrost thaw. We observed (1) consistent relationships between soil property and depth and between major parameters, (2) large contrasts of key soil parameters between active layer and permafrost, indicative of potentially different response of the permafrost carbon to warming when compared to the active layer, and (3) a correlation between soil hydraulic conductivity and topographic features that impacts soil hydrologic processes. Our analysis suggests that the permafrost has a marine-derived chemical signature that differs from the active layer and shapes the physicochemical fingerprints of the different geomorphic features. Specifically, we revealed the unique signatures of the high center polygons, indicative of possible microbial activity at depth (>1 m). Our study suggested consistent key soil parameter-depth correlations while demonstrating complex lateral and vertical variabilities. These results are valuable for identifying approaches to upscale point-based measurements and for improving model parameterization to predict permafrost carbon behavior and feedback under future climate.

Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiagvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

Understanding the interactions among agricultural processes, soil, and plants is necessary for optimizing crop yield and productivity. This study focuses on developing effective monitoring and analysis methodologies that estimate key soil and plant properties. These methodologies include data acquisition and processing approaches that use unmanned aerial vehicles (UAVs) and surface geophysical techniques. In particular, we applied these approaches to a soybean farm in Arkansas to characterize the soil-plant coupled spatial and temporal heterogeneity, as well as to identify key environmental factors that influence plant growth and yield. UAV-based multitemporal acquisition of high-resolution RGB (red-green-blue) imagery and direct measurements were used to monitor plant height and photosynthetic activity. We present an algorithm that efficiently exploits the high-resolution UAV images to estimate plant spatial abundance and plant vigor throughout the growing season. Such plant characterization is extremely important for the identification of anomalous areas, providing easily interpretable information that can be used to guide near-real-time farming decisions. Additionally, high-resolution multitemporal surface geophysical measurements of apparent soil electrical conductivity were used to estimate the spatial heterogeneity of soil texture. By integrating the multiscale multitype soil and plant datasets, we identified the spatiotemporal co-variance between soil properties and plant development and yield. Our novel approach for early season monitoring of plant spatial abundance identified areas of low productivity controlled by soil clay content, while temporal analysis of geophysical data showed the impact of soil moisture and irrigation practice (controlled by topography) on plant dynamics. Our study demonstrates the effective coupling of UAV data products with geophysical data to extract critical information for farm management.

The EGS Collab project performed well-monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory to inform challenges in implementing enhanced geothermal system (EGS) technology. This project, supported by the US Department of Energy, gathered data and observations from the field tests and compared these to simulation results to understand processes and to build confidence in numerical modeling of the processes. The project consisted of 3 Experiments, each comprising test and testbed design, many individual tests, numerical simulation, and analysis. The Experiments were performed in two deep underground testbeds at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. Field experiments are now complete, significant data sets have been collected and analyzed, and some analysis continues. Experiments using underground test facilities have many advantages in that they allow: • Three-dimensional characterization of the stimulated volume by complementary geophysical methods surrounding the experiment • Using techniques that are currently not applicable under geothermal condition to provide processes insight • Comprehensive tracer testing and detailed characterization of complex fluid movements • Understanding the geometry of the stimulated network at the meso-scale and its implications for effective fracture surface area, rock block size, and heat exchange. Underground testing has its own set of complications, however, which affect the ability to perform tests as desired and affect the experiment results. Included here are the inability to flow at desired injection rates due to stress gradients caused by drift cooling, and the need to strongly limit induced seismicity because the distance to people and equipment. Experiment 1 examined hydraulic fracturing at a depth of 1.5 km in a well-characterized phyllite. Eight subhorizontal boreholes were used in this Experiment. Geophysical monitoring instrumentation was deployed in six boreholes to monitor stimulation events and flow tests. The other two boreholes were used to perform and carefully measure water injection and production. More than a dozen stimulations and nearly one year of flow tests in the testbed were performed. Detailed observations of processes occurring during stimulation and dynamic flow tests were collected and analyzed. Flow tests using ambient-temperature and chilled water were performed with intermittent tracer tests to examine system behavior. We achieved adaptive control of the tests using close monitoring of rapidly disseminated data and near-real-time simulation. Numerical simulation was critical in answering key experimental design questions, forecasting fracture behavior, and analyzing results. We were successful in performing many simulations in near-real-time in conjunction with the field experiments, with more detailed simulations performed later. Experiment 2 was intended to examine hydraulic shearing of natural fractures at a depth of 1.25 km in amphibolite. The stresses, rock type, and fracture conditions are different than in Experiment 1. The testbed consists of 9 boreholes, in addition to 2 exploratory characterization boreholes. Four boreholes drilled as two fans of 2 monitoring holes contained grouted-in monitoring sensors. The remaining five open boreholes drilled in a five-spot pattern were adaptively used for injection, production, and monitoring. Approximately five fracture set orientations were encountered in the testbed along with a low-stress rhyolite sill at 35 m below the access drift in exploratory well TV4100. The testbed was designed to optimize the potential for shear stimulation while also avoiding the low-stress rhyolite. Experiment 2 focused on stimulating a fracture in the most likely orientation to shear, however shear stimulation did not occur probably due to cementation from natural secondary mineralization. Other fracture sets encountered were also cemented and had orientations less likely to shear. Experiment 3 consisted of several stimulations in the same testbed as Experiment 2 allowing different stimulation approaches including ramped-rate injection, rapid injection, and oscillating-pressure injection. Ultimately these methods created hydraulic fractures, one of which was used for a medium-duration cold water injection test. The major findings of the EGS Collab Project include: 1. Significant shear stimulation did not occur during our stimulation attempts. Shear stimulation may occur but under a limited set of conditions not encountered. 2. Our stimulations resulted in hydraulic fractures that required hydraulic propping. Pumping at pressures exceeding the minimum principal stress may not be feasible in an enhanced geothermal system. 3. The systems we generated were complex hydraulic fracture/natural fracture systems, and these systems changed over time in response to applied pressures and flowrates and to unknown stimuli. 4. The project attempted alternative stimulation methods, which did not provide significant flow improvement. 5. Thermal breakthrough was not achieved as designed, most likely because flow to production boreholes was not adequate. 6. The combination of geophysical tools used provided excellent understanding of many important processes. 7. Microearthquakes (MEQs) didn’t necessarily identify flow paths. 8. Engineering tools bounding expected seismicity are needed. This report provides a summary of tests and analyses performed for EGS Collab Experiment 2 (Shear Stimulation in Testbed 2) and Experiment 3 (Alternative Stimulation methods). Much of the EGS Collab work has been published in journals and conference papers, presented in conferences, included in written reports, and submitted in data sets to the Geothermal Data Repository (GDR). The entirety of these written works is included as an appendix to this report, and this report serves as a summary and framework pointing to these published papers, presentations, and reports.

We conducted a time-lapse seismic experiment utilizing automated active seismic source and sensor arrays to monitor a reactivated fault within the Opalinus clay formation at the Mont Terri Rock Laboratory (Switzerland), an analog caprock for geologic carbon storage. A series of six brine injections were conducted into the so-called Main Fault to reactivate it. Seismic instrumentation in five monitoring boreholes on either side of the fault was used to continuously probe changes in P-wave travel-times associated with fault displacement and leakage. We performed time-lapse travel-time tomography on five hundred sequential data sets; this revealed a zone of decreased P-wave velocity, up to 16 m/s, during each injection cycle, followed by a velocity increase during shut-in. These results demonstrate varying elastic property perturbations, both spatially and temporally, along the fault plane during reactivation. We then interpreted these velocity changes in terms of fault dilation induced by pressurized fluids along the fault.