Inversion of multicomponent seismic time shifts for reservoir pressure and length: a feasibility study

Pressure drops associated with reservoir production generate excess stress and strain that cause travel‐time shifts of reflected waves. Here, we invert time shifts of P‐, S‐, and PS‐waves measured between baseline and monitor surveys for pressure reduction and reservoir length. The inversion results can be used to estimate compaction‐induced stress and strain changes around the reservoir. We implement a hybrid inversion algorithm that incorporates elements of gradient, global/genetic, and nearest neighbour methods and permits exploration of the parameter space while simultaneously following local misfit gradients. Our synthetic examples indicate that optimal estimates of reservoir pressure from P‐wave data can be obtained using the reflections from the reservoir top. For S‐waves, time shifts from the top of the reservoir can be accurately inverted for pressure if the noise level is low. However, if noise contamination is significant, it is preferable to use S‐wave data (or combined shifts of all three modes) from reflectors beneath the reservoir. Joint wave type inversions demonstrate improvements over any single pure mode. Reservoir length can be estimated using the time shifts of any mode from the reservoir top or deeper reflectors. We also evaluate the differences between the actual strain field and those corresponding to the best‐case inversion results obtained using P‐ and S‐wave data. Another series of tests addresses the inversion of the time shifts for the pressure drops in two‐compartment reservoirs, as well as for the associated strain field. Numerical testing shows that a potentially serious source of error in the inversion is a distortion in the strain‐sensitivity coefficients, which govern the magnitude of stiffness changes. This feasibility study suggests which wave types and reflector locations may provide the most accurate estimates of reservoir parameters from compaction‐induced time shifts.     

The role of the hyporheic zone across stream networks

Many hyporheic papers state that the hyporheic zone is a critical component of stream ecosystems, and many of these papers focus on the biogeochemical effects of the hyporheic zone on stream solute loads. However, efforts to show such relationships have proven elusive, prompting several questions: Are the effects of the hyporheic zone on stream ecosystems so highly variable in place and time (or among streams) that a consistent relationship should not be expected? Or, is the hyporheic zone less important in stream ecosystems than is commonly expected? These questions were examined using data from existing groundwater modelling studies of hyporheic exchange flow at five sites in a fifth‐order, mountainous stream network. The size of exchange flows, relative to stream discharge (Q_HEF:Q), was large only in very small streams at low discharge (area ≈ 100 ha; Q < 10 l/s). At higher flows (flow exceedance probability > 0·7) and in all larger streams, Q_HEF:Q was small. These data show that biogeochemical processes in the hyporheic zone of small streams can substantially influence the stream's solute load, but these processes become hydrologically constrained at high discharge or in larger streams and rivers. The hyporheic zone may influence stream ecosystems in many ways, however, not just through biogeochemical processes that alter stream solute loads. For example, the hyporheic zone represents a unique habitat for some organisms, with patterns and amounts of upwelling and downwelling water determining the underlying physiochemical environment of the hyporheic zone. Similarly, hyporheic exchange creates distinct patches of downwelling and upwelling. Upwelling environments are of special interest, because upwelling water has the potential to be thermally or chemically distinct from stream water. Consequently, micro‐environmental patches created by hyporheic exchange flows are likely to be important to biological and ecosystem processes, even if their impact on stream solute loads is small. Published in 2011 by John Wiley & Sons, Ltd.     

Investigation of Small‐Scale Preferential Flow with a Forced‐Gradient Tracer Test

A new tracer experiment (referred to as MADE‐5) was conducted at the well‐known Macrodispersion Experiment (MADE) site to investigate the influence of small‐scale mass‐transfer and dispersion processes on well‐to‐well transport. The test was performed under dipole forced‐gradient flow conditions and concentrations were monitored in an extraction well and in two multilevel sampler (MLS) wells located at 6, 1.5, and 3.75 m from the source, respectively. The shape of the breakthrough curve (BTC) measured at the extraction well is strongly asymmetric showing a rapidly arriving peak and an extensive late‐time tail. The BTCs measured at seven different depths in the two MLSs are radically different from one another in terms of shape, arrival times, and magnitude of the concentration peaks. All of these characteristics indicate the presence of a complex network of preferential flow pathways controlling solute transport at the test site. Field‐experimental data were also used to evaluate two transport models: a stochastic advection‐dispersion model (ADM) based on conditional multivariate Gaussian realizations of the hydraulic conductivity field and a dual‐domain single‐rate (DDSR) mass‐transfer model based on a deterministic reconstruction of the aquifer heterogeneity. Unlike the stochastic ADM realizations, the DDSR accurately predicted the magnitude of the concentration peak and its arrival time (within a 1.5% error). For the multilevel BTCs between the injection and extraction wells, neither model reproduced the observed values, indicating that a high‐resolution characterization of the aquifer heterogeneity at the subdecimeter scale would be needed to fully capture 3D transport details.     

Lessons Learned from 25 Years of Research at the MADE Site

Field studies at well‐instrumented research sites have provided extensive data sets and important insights essential for development and testing of transport theories and mathematical models. This paper provides an overview of over 25 years of research and lessons learned at one of such field research sites on the Columbus Air Force Base in Mississippi, commonly known as the Macrodispersion Experiment (MADE) site. Since the mid‐1980s, field data from the MADE site have been used extensively by researchers around the world to explore complex contaminant transport phenomena in highly heterogeneous porous media. Results from field investigations and modeling analyses suggested that connected networks of small‐scale preferential flow paths and relative flow barriers exert dominant control on solute transport processes. The classical advection‐dispersion model was shown to inadequately represent plume‐scale transport, while the dual‐domain mass transfer model was found to reproduce the primary observed plume characteristics. The MADE site has served as a valuable natural observatory for contaminant transport studies where new observations have led to better understanding and improved models have sprung out analysis of new data.     

Principal Slip Zones in Limestone: Microstructural Characterization and Implications for the Seismic Cycle (Tre Monti Fault, Central Apennines, Italy)

Earthquakes in central Italy, and in other areas worldwide, often nucleate within and rupture through carbonates in the upper crust. During individual earthquake ruptures, most fault displacement is thought to be accommodated by thin principal slip zones. This study presents detailed microstructural observations of the slip zones of the seismically active Tre Monti normal fault zone. All of the slip zones cut limestone, and geological constraints indicate exhumation from <2?km depth, where ambient temperatures are ?100°C. Scanning electron microscope observations suggest that the slip zones are composed of 100% calcite. The slip zones of secondary faults in the damage zone contain protocataclastic and cataclastic fabrics that are cross-cut by systematic fracture networks and stylolite dissolution surfaces. The slip zone of the principal fault has much more microstructural complexity, and contains a 2?C10?mm thick ultracataclasite that lies immediately beneath the principal slip surface. The ultracataclasite itself is internally zoned; 200?C300???m-thick ultracataclastic sub-layers record extreme localization of slip. Syn-tectonic calcite vein networks spatially associated with the sub-layers suggest fluid involvement in faulting. The ultracataclastic sub-layers preserve compelling microstructural evidence of fluidization, and also contain peculiar rounded grains consisting of a central (often angular) clast wrapped by a laminated outer cortex of ultra-fine-grained calcite. These ??clast-cortex grains?? closely resemble those produced during layer fluidization in other settings, including the basal detachments of catastrophic landslides and saturated high-velocity friction experiments on clay-bearing gouges. An overprinting foliation is present in the slip zone of the principal fault, and electron backscatter diffraction analyses indicate the presence of a weak calcite crystallographic preferred orientation (CPO) in the fine-grained matrix. The calcite c-axes are systematically inclined in the direction of shear. We suggest that fluidization of ultracataclastic sub-layers and formation of clast-cortex grains within the principal slip zone occurred at high strain rates during propagation of seismic ruptures whereas development of an overprinting CPO occurred by intergranular pressure solution during post-seismic creep. Further work is required to document the range of microstructures in localized slip zones that cross-cut different lithologies, and to compare natural slip zone microstructures with those produced in controlled deformation experiments.     

Sensitivity analysis of CDOM spectral slope in artificial and natural samples: an application in the central eastern Mediterranean Basin

In the past two decades, optical properties of chromophoric dissolved organic matter (CDOM) in marine environments have been extensively studied. Many of these studies report CDOM properties for the offshore environment where this complex mixture of optically active compounds is strongly diluted. Nevertheless, autochthonous and allochthonous sources have been identified and sinks related to photodegradation and bacterial activity have been demonstrated. The calculation of the spectral slope of the CDOM absorption curve has been proven to be useful and is often reported. However, a rigorous uncertainty analysis of the slope calculation is rarely reported. In this paper, we propose a method to evaluate the uncertainty of CDOM spectral slope calculated between 270 and 400 nm, using both naturally sampled and artificial solutions. We use these results to study the ultra-oligotrophic waters of the Mediterranean Sea (central eastern basin), where little is known about CDOM spatial distribution. We show that dilutions of both artificial and natural samples produce a Gaussian distribution of spectral slopes, indicating that consistent values may be determined, with a typical uncertainty of ±0.0004 nm⁻¹ when absorption at 300 nm was greater then 0.1 m⁻¹ (0.1 m pathlength). Comparing the distribution of spectral slopes from central eastern basin samples to a Gaussian distribution, we show differences between measurements that were significantly different. These values allow us to distinguish possible sources (algal derived CDOM), sinks (e.g. photo-bleaching) at different depths. We propose a subdivision of CDOM compounds into refractory and semilabile/refractory pools and evaluate the CDOM spectral slope of algal derived CDOM released at or near deep chlorophyll maximum.     

Export and degassing of terrestrial carbon through watercourses draining a temperate podzolized catchment

We measured spatial and temporal variations in carbon concentrations, isotopic compositions and exports during a complete hydrological cycle in nine watercourses draining a lowland forested podzolized catchment, flowing into the Arcachon lagoon (France). In addition, integrated fluxes of CO₂ across the water-atmosphere interface were estimated to assess the relative importance of CO₂ evasion versus lateral carbon transport at the catchment scale. Watercourse similarities and specificities linked to the local catchment characteristics are discussed and compared with other riverine systems. Low concentrations of suspended particulate matter and particulate organic carbon (POC) were generally measured in all the watercourses (8.4 ± 3.4 and 1.6 ± 0.6 mg L^?1, respectively), reflecting limited mechanical soil erosion. The generally high POC content in the suspended matter (20 %), low Chl a concentrations (1.3 ± 1.4 μg L^?1) and the relatively constant δ¹³C-POC value (near ?28 ‰) throughout the year reveal this POC originates from terrestrial C₃ plant and soil detritus. The presence of podzols leads to high levels of dissolved organic carbon (DOC; 6.6 ± 2.2 mg L^?1). Similarly, high dissolved inorganic carbon (DIC) concentrations were measured in the Arcachon lagoon catchment (5.9 ± 2.2 mg L^?1). The δ¹³C-DIC value around ?20 ‰ throughout the year in many small watercourses reveals the predominance of terrestrial carbon mineralisation and silicate rock weathering in soils as the major DIC source. With pCO₂ between 1,000 and 10,000 ppmv, all watercourses were a source of CO₂ to the atmosphere, particularly during the low river stage. Organic carbon parameters remained relatively stable throughout the year, whereas DIC parameters showed strong seasonal contrasts closely linked to the hydrological regime and hyporheic flows. In total, the carbon export from the Arcachon watershed was estimated at 15,870 t C year^?1 or 6 t C km^?2 year^?1, mostly exported to the lagoon as DOC (35 %), DIC (24 %) and lost as CO₂ degassing to the atmosphere (34 %).     

Simulation of the 1980 eruption of Mount St. Helens using the ash-tracking model PUFF

The dispersal of volcanic ash from the May 18, 1980 eruption of Mount St. Helens (MSH) has been simulated using the Lagrangian ash-tracking model PUFF. Previous applications of the model were limited to smaller, short-lived eruptions with ash dispersal occurring mainly within the troposphere. Two high-resolution atmospheric reanalysis datasets (ERA-40 and NCEP/NCAR-40) allowed MSH ash cloud dispersal to be simulated up to 30 km elevation. The 1980 eruption was divided into two distinct eruptive phases, (1) an initial, relatively short-lived blast/surge phase that injected ash up to 30 km and (2) a subsequent nine-hour plinian phase that maintained an average eruption column height of 16 km. Using PUFF, the two phases of the MSH eruption were modeled separately based on a range of individual input parameters and then combined to produce an integrated simulation of the entire eruption. The trajectory and areal extent of the modeled atmospheric ash cloud best match the actual distribution of MSH ash when input parameters are set to values inferred from satellite and radar data collected on May 18, 1980. The prevailing wind field exerts the strongest control on the advection and ultimate position of the modeled ash cloud, making the maximum column height and the vertical distribution of ash the most sensitive of the PUFF input parameters for this event. The results indicate that the PUFF model works well at simulating the dispersal of ash injected well into the lower stratosphere from a moderate, relatively long-lived eruption, such as MSH. However, attempts to use PUFF to recreate some granulometric aspects of the MSH fallout deposit, such as the maximum particle size as a function of distance from source, were not successful. PUFF consistently predicts much greater fallout distances for small ash particles (< 500 µm) than actually observed in the MSH deposit. The effective settling velocities used by the PUFF model appear to be too slow to accurately predict fallout distances of small ash particles. As a consequence the PUFF model may overestimate the duration of ash loading in the atmosphere associated with the distal fine ash component of explosive eruptions.     

The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability

The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport.     

Strong slope-aspect control of regolith thickness by bedrock foliation

The porous near-surface layer of the Earth's crust – the critical zone – constitutes a vital reservoir of water for ecosystems, provides baseflow to streams, guides recharge to deep aquifers, filters contaminants from groundwater, and regulates the long-term evolution of landscapes. Recent work suggests that the controls on regolith thickness include climate, tectonics, lithology, and vegetation. However, the relative paucity of observations of regolith structure and properties at landscape scales means that theoretical models of critical zone structure are incompletely tested. Here we present seismic refraction and electrical resistivity surveys that thoroughly characterize subsurface structure in a small catchment in the Santa Catalina Mountains, Arizona, USA, where slope-aspect effects on regolith structure are expected based on differences in vegetation. Our results show a stark contrast in physical properties and inferred regolith thickness on opposing slopes, but in the opposite sense of that expected from environmental models and observed vegetation patterns. Although vegetation (as expressed by normalized difference vegetation index [NDVI]) is denser on the north-facing slope, regolith on the south-facing slope is four times thicker (as indicated by lower seismic velocities and resistivities). This contrast cannot be explained by variations in topographic stress or conventional hillslope morphology models. Instead, regolith thickness appears to be controlled by metamorphic foliation: regolith is thicker where foliation dips into the topography, and thinner where foliation is nearly parallel to the surface. We hypothesize that, in this catchment, hydraulic conductivity and infiltration capacity control weathering: infiltration is hindered and regolith is thin where foliation is parallel to the surface topography, whereas water infiltrates deeper and regolith is thicker where foliation intersects topography at a substantial angle. These results suggest that bedrock foliation, and perhaps by extension sedimentary layering, can control regolith thickness and must be accounted for in models of critical zone development. © 2020 John Wiley & Sons, Ltd.