Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

Trace explosive residue detection of HMX and RDX in post-detonation dust from an open-air environment

Explosives are often used in industry, geology, mining, and other applications, but it is not always clear what remains after a detonation or the fate and transport of any residual material. The goal of this study was to determine to what extent intact molecules of high explosive (HE) compounds are detectable and quantifiable from post-detonation dust and particulates in a field experiment with varied topography. We focused on HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane), which is less studied in field detonation literature, as the primary explosive material and RDX (1,3,5-Trinitroperhydro-1,3,5-triazine) as the secondary material. The experiment was conducted at Site 300, Lawrence Livermore National Laboratory’s Experimental Test Site, in California, USA. Two 20.4 kg and one 40.8 kg above ground explosions (primarily comprised of LX-14, an HMX-based polymer-bonded high explosive) were detonated on an open-air firing area on separate days. The complex terrain of the firing area (e.g., buildings, berm, low-height obstacles) was advantageous to study HE deposition in relation to plume dynamics. Three types of samples were collected up to 100 m away from each shot: surface swipes of aluminum plates, surface swipes of fixed objects, and filters from air samples. We used atmospheric flow tube-mass spectrometry (AFT-MS) to quantify picogram levels of molecular residue of HE material in the post-detonation dust. An aliquot of sample extract in methanol (e.g., 1 µL of 0.5 mL) was placed onto a resistive material and then thermally desorbed into the AFT-MS. We successfully detected and quantified both HMX and RDX in many of the samples. Based on mass (pg) detected and solution dilution, we back-calculated the mass collected on the swipe or filter (ng per sample). The aerial distribution of molecular residue was consistent with the path of the plume, which was strongly determined by wind speed and direction at the time of each shot. The quantity of material detected appeared to correlate more with distance from the shot and the wind conditions than with shot size. This study demonstrates that the picogram detection levels of AFT-MS are well-suited for quantification of analytes (e.g., HMX and RDX) in environmental samples.

atmospheric flow tube-mass spectrometry (AFT-MS), ↗

Gaining Perspective on Unconventional Well Design Choices through Play-level Application of Machine Learning Modeling

The recent development of unconventional oil and gas (O&G) reservoirs has led to an abundant hydrocarbon supply, both domestically and globally. However, there is a continued push to develop new and innovative approaches to improve exploration and extraction efficiencies and overall well productivity moving forward. Substantial improvements in unconventional O&G development are expected through optimized well completion and stimulation strategies aimed at maximizing well productivity. Optimizing well designs will require tailoring to the distinctive geologic conditions present for any newly placed well. To better evaluate the impact of well design attributes and their associated interactions on productivity in a major unconventional play, multivariate machine learning-based models that use empirical datasets were developed. A gradient boosted regression tree (GBRT) algorithm was applied. GBRT has been narrowly investigated for O&G applications but enables straightforward parametric importance and influence evaluation, as well as assessment of parameter interaction effects. Models were trained on well design and locational parameters that serve as a proxy for variable geologic conditions to estimate two types of productivity indicator response variables strongly correlated to estimated ultimate recovery (EUR). The dataset utilized consists of over 7,000 well observations that cover the majority of the productive region of the Marcellus Shale. Model performance was evaluated and algorithm parameters tuned by analyzing the goodness-of-fit for simulated results against observed data in a cross-validation approach. Models were found capable of 73–79 percent prediction accuracy on held out testing data of gas equivalent production and can be used to inform future well design and placement decisions for increasing EUR per well and improving overall field-level recovery. Study results indicate that Marcellus well performance improves most with upscaling perforated interval lengths and water and proppant volumes per foot; but relative productivity improvements are spatially dependent across the play. Finally, optimal combinations of water and proppant on well performance were found to vary depending on well location, emphasizing the utility of data-driven models capable of broad application across a play of interest for informing tailored well design approaches prior to their field deployment.

04 OIL SHALES AND TAR SANDS↗

Transportability of exogenous microbial community correlates with interwell connectivity in deep aquifers

Subsurface resource engineering operations often utilize continuous injection of externally-sourced water into geological reservoirs for formation pressure maintenance, resource recovery or energy/waste storage. Such injected water generally contains naturally occurring microbes. Little is known, however, about how the injectate microbes transport through geological media as a community, how such transportability is affected by injector-producer connectivity, and whether such knowledge can be utilized for flowpath characterization. In this study, we analyzed daily-to-weekly timeseries microbial community data from the injected- and produced-fluids of a ten-month flow test at a deep, well-characterized engineered aquifer. We found that the injectate microbial community was distinct from the indigenous community at the amplicon sequence variant (ASV) level, and that the transportability of injectate community towards a given producer, quantified by an “nASV-Overlap” metric we propose, had strong and significant positive correlation with known injector-producer connectivities at our site. This suggests that the better the connectivity, the higher the probability for more injectate species to flow through the interwell region and arrive at a producer. Because interwell connectivity is an important yet usually unknown parameter in subsurface resource engineering, such correlation in turn points to nASV-Overlap as a useful indicator of interwell connectivity for aquifer characterization and long-term monitoring. Based on our findings, an nASV-Overlap-based microbial tracing approach was developed for characterizing and monitoring the relative connectivities across multiple producers with a given injector. A side-by-side comparison between the new nASV-Overlap approach and traditional artificial tracer methods is presented, and their respective strengths and limitations are discussed.

Deep biosphere↗

Real-Time Atomic-Scale Structural Analysis Resolves the Amorphous to Crystalline CaCO 3 Mechanism Controversy

Amorphous calcium carbonate (ACC) occurs as a precursor to geological and biogenic calcium carbonate (CaCO 3 ), yet its transformation pathways and reaction mechanisms remain inconsistent and controversial. In this study, we investigated the transformation of ACC to calcite under both solution and dry conditions, in the presence and absence of impurity ions, utilizing operando time-resolved synchrotron X-ray diffraction (TRXRD) and reactive transport modeling. Results demonstrate that TRXRD techniques allow us to differentiate dissolution-reprecipitation versus solid-state transformation mechanisms for amorphous to crystalline phase transitions. Specifically, we observe that in environments with abundant water, ACC transforms to calcite through a dissolution-reprecipitation mechanism. This features an activation energy of 63 ± 2 kJ/mol and unit cell volume contraction during calcite crystal growth. Conversely, under water-limited conditions, ACC to calcite transformation proceeds via a solid-state transformation mechanism, with an activation energy of 210 ± 2 kJ/mol, three times greater than the dissolution-reprecipitation route, and a unit cell expansion during crystalline calcite growth. Further, to illustrate the magnitude of these effects, the rates of calcite growth were similar during dissolution-reprecipitation at 3 °C [0.00207(35) s –1 ] and solid-state transformation at 280 °C [0.00134(11) s –1 ]. Moreover, the incorporation of an impurity, strontium, significantly retards the rate of calcite growth while expanding its unit cell but whose incorporation is history dependent. Reactive transport modeling of the dissolution–precipitation kinetics suggests that ACC must be dissolving as compact aggregates. These various transformation mechanisms drive diverse geological and biological carbonate formations, impacting their use as paleoenvironmental markers and functional materials synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recent Advances in Corrosion Science Applicable To Disposal of High-Level Nuclear Waste

High level radioactive waste is accumulating at temporary storage locations around the world and will eventually be placed in deep geological repositories. Many different types of geological formation are under consideration, but all will eventually allow water to contact the metallic waste containers considering that the period of performance required to allow sufficient decay of dangerous radionuclides is on the order of 105 – 106 years. Corrosion of the containers and waste forms in the aqueous repository environment is therefore a concern. This review describes the recent advances of the field of materials corrosion based on the work of the Center for Performance and Design of Nuclear Waste Forms and Containers (WastePD) to address the issues associated with the long-term performance assessment and the design of materials with improved performance, where performance is defined as resistance to aqueous corrosion. Glass, crystalline ceramics, and metals are discussed separately, and then the near field interactions of these different materials classes are addressed. Recommendations for future directions are provided.

Frankel, Gerald↗

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗

Self-Assembly of Silica Nanoparticles at Water–Hydrocarbon Interfaces: Insights from In Operando Small-Angle X-ray Scattering Measurements and Molecular Dynamics Simulations

To achieve tunable controls on the interactions of siliceous materials in subsurface environments bearing water and hydrocarbons, it is essential to determine the influence of hydrocarbon–water interfaces on the self-assembly of the silica particles. The hydrophilic silica particle have a tendency to aggregate on the water front of the interface. Self-assembly of silica nanoparticles proceeds via the migration of these nanoparticles to the water–hydrocarbon interface followed by aggregation at the interface. Fractal-like morphologies of assembled silica nanoparticles at water–hydrocarbon interfaces are observed. Rapid assembly of the hydrophilic silica nanoparticles at water–hydrocarbon interfaces corresponds to an overall reduction in the surface tension of water–toluene and water–heptane systems. These studies demonstrate the silica aggregation that is the precursor of silica polymerization, or nucleation and growth is influenced by the presence of hydrocarbons in subsurface geologic environments. These insights were derived from in operando ultrasmall and small-angle X-ray scattering (USAXS/SAXS) measurements, cryo-scanning electron microscopy (Cryo-SEM) imaging, and classical molecular dynamics simulations. As a result, these studies are intended to inform current and future efforts aimed at tuning silica reactivity in subsurface geologic environments to enhance permeability, the design and use of silica-based proppants to enhance fractures, and the development of effective strategies to control silica-based scaling behavior in subsurface reservoirs.

42 ENGINEERING↗

Alteration of Fractured Foamed Cement Exposed to CO2-Saturated Water: Implications for Well Integrity

Geologic CO2 storage (GCS) is a method to mitigate the adverse impact of global climate change. Potential leakage of CO2 from fractured cement at the wellbore poses a risk to the feasibility of GCS. Foamed cement is widely applied in deepwater wells where fragile geologic formations cannot support the weight of conventional cement. Thus, it is critical to know whether fractures in foamed cement self-seal in a similar manner as conventional cement systems. This study is the first to investigate the changes in physical and chemical attributes of foamed cement under dynamic flow conditions using CO2-saturated water. Self-sealing of fractures in the cement was observed at a solution flow rate of 0.1 mL/min and a pressure of 6.9 MPa. The formation of CaCO3 precipitates in pore spaces and fractures led to a decrease in permeability by 1 order of magnitude. The extents of self-sealing in foamed cement samples, specifically the 20 and 30% air volume formulations, were similar to that of conventional cements. We attribute this to the greater alteration depth in the foamed cement, which compensated for the reduced availability of Portlandite and higher initial porosity. The results can be used to evaluate the risk of leakage associated with foamed cement.

Geologcial Carbon Storage, Well Integrity in CO2 S↗

The Value of Reversible Carbon Storage in a Zero-Emissions World

Atmospheric carbon dioxide removal (CDR) is required to stabilize global temperature. CDR can be achieved via ecosystem-based approaches that are cost-effective but reversible (e.g., soil and forest management) or by more durable but expensive approaches (e.g., direct air capture coupled with geologic storage). Here, we examine trade-offs between these approaches, focusing on timing, climate impacts, and cost. We simulated reversible carbon accrual for a range of CDR contract structures using a general minimalist model of ecosystem carbon cycling, and parameterized it to simulate US agricultural soil management─specifically cover cropping─as a case study. We then quantified the resulting impact on atmospheric carbon and global temperature using a climate model emulator. We find that maintaining a patchwork of reversible CDR projects by replacing lapsed projects with new projects can reduce warming by 22–195 μ°C in 2100 and that the magnitude of this cooling effect depends on how effectively the patchwork is maintained. Long-term maintenance of reversible CDR projects requires institutional stability that cannot be guaranteed over multiple decades. Consequently, effective CDR ultimately requires replacing reversible projects with durable projects. To address this problem, we modeled the cost of replacing reversible agricultural soil CDR with geologic CDR. We found that using reversible CDR as a bridge to durable CDR is potentially more cost-effective as a global cooling strategy (0.20–0.81 billion USD per μ°C avoided) than perpetual maintenance of reversible CDR (0.32–1.31 billion USD per μ°C avoided) or an immediate transition to durable CDR (1.37–2.19 billion USD per μ°C avoided). However, we emphasize that institutional commitments to maintain reversible CDR projects cannot be guaranteed. Reliance on reversible CDR as a bridge to durable CDR therefore carries an unknown amount of risk and will only function if efforts to maintain reversible CDR are robust.

carbon capture and storage↗

Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

02 PETROLEUM↗

Porous Colloidal Nanoparticles as Injectable Multimodal Contrast Agents for Enhanced Geophysical Sensing

Injecting fluids into underground geologic structures is crucial for the development of long-term strategies for managing captured carbon and facilitating sustainable energy extraction operations. Here, we have previously reported that the injection of metal–organic frameworks (MOFs) into the subsurface can enhance seismic monitoring tools to track fluids and map complex structures, reduce risk, and verify containment in carbon storage reservoirs because of their absorption capacity of low-frequency seismic waves. Here, we demonstrate that water-based Cr/Zn/Zr MOF colloidal suspensions (nanofluids) are multimodal geophysical contrast agents that enhance near-wellbore logging tools. Based on experimental fluid-only measurements, MIL-101(Cr), ZIF-8, and UiO-66 nanofluids have distinct complex conductivity and/or low-field nuclear magnetic resonance (NMR) signatures that are relevant to field-deployed technologies, implying the potential to enhance near-wellbore monitoring of CO 2 injection and associated processes with downhole logging tools. Small- and wide-angle X-ray scattering characterization of ~0.5 wt % MIL-101(Cr) suspensions confirmed phase stability and provided insight into the fractal nature of colloidal nanoparticles. Finally, low-field (2 MHz) NMR measurements of MIL-101(Cr) nanofluid injection into a prototypical Berea sandstone demonstrate how paramagnetic high-surface area MOFs may dominate the relaxation times of hydrogen-bearing fluids in porous geologic matrices, enhancing the mapping of near-surface and near-wellbore transport pathways and advancing sustainable subsurface energy technologies.

NMR↗

Discriminating Underground Nuclear Explosions Leading To Late-Time Radionuclide Gas Seeps

Utilizing historical data from the U.S. nuclear test program and freely available barometric pressure data, we performed an analytical barometric-pumping efficiency analysis to determine factors resulting in late-time radionuclide gas seeps from underground nuclear explosions. In this work, we considered 16 underground nuclear explosions with similar geology and test setup, of which five resulted in the measurement of late-time radionuclide gas concentrations at the ground surface. Additionally, the factors we considered include barometric frequency and amplitude, depth of burial, air-filled porosity, intact-rock permeability, fracture aperture, and fracture spacing. The analysis indicates that the best discriminators of late-time radionuclide gas seeps for these explosions are barometric frequency and amplitude and air-filled porosity. While geologic information on fracture aperture and spacing is not available for these explosions, the sensitivity of barometric-pumping efficiency to fracture aperture indicates that it would likely also be a good discriminator.

58 GEOSCIENCES↗

Preliminary Analysis of Source Physics Experiment Explosion-Triggered Microseismicity Using the Back-Projection Method

A series of four chemical explosions were detonated in a deep borehole within the Yucca Flat Dry Alluvium Geology (DAG) at the Nevada National Security Site between 2018 and 2019. The two larger chemical explosions of 50 tons (DAG-2) and 10 tons (DAG-4) TNT equivalent yield triggered energetic aftershock sequences numbering 1392 and 347 microearthquakes, respectively, within the first 10 days. No significant aftershock activity was observed for the two smaller 1-ton explosions (DAG-1 and DAG-3). Here, we used a back-projection method based on travel-time migration and stacking of signal-to-noise ratio traces to detect, associate and locate aftershocks from a subset of 22-geophones within a larger 2 urn:x-wiley:21699313:media:jgrb54910:jgrb54910-math-0001 2 km seismic array surrounding the borehole. The aftershocks located within 300 m of the borehole and the depths were above the working points of 300 and 50 m depths of DAG-2 and DAG-4, respectively, ruling out triggering slip on geologic faults or disturbances beneath neighboring collapse craters. DAG-2 and DAG-4 aftershocks decayed at similar rates, with power-law exponents of p = 1.48 and p = 1.49, respectively. These decay rates are comparable to aftershocks sequences triggered by earthquakes and historical nuclear explosions at Yucca Flat. A smooth power-law aftershock decay within the first 10 days suggests a triggering mechanism from explosion generated stress relaxation due to the diffusion of high gas pressures in the cavity and radial fractures. A more random and episodic aftershock rate would be expected due to cavity collapse or falling rubble in chimney formation.

58 GEOSCIENCES↗

Distributed Acoustic Sensing as a Distributed Hydraulic Sensor in Fractured Bedrock

Distributed acoustic sensing (DAS) was originally intended to measure oscillatory strain at frequencies of 1 Hz or more on a fiber optic cable. Recently, measurements at much lower frequencies have opened the possibility of using DAS as a dynamic strain sensor in boreholes. A fiber optic cable mechanically coupled to a geologic formation will strain in response to hydraulic stresses in pores and fractures. A DAS interrogator can measure dynamic strain in the borehole, which can be related to fluid pressure through the mechanical compliance properties of the formation. Because DAS makes distributed measurements, it is capable of both locating hydraulically active features and quantifying the fluid pressure in the formation. We present field experiments in which a fiber optic cable was mechanically coupled to two crystalline rock boreholes. The formation was stressed hydraulically at another well using alternating injection and pumping. The DAS instrument measured oscillating strain at the location of a fracture zone known to be hydraulically active. Rock displacements of less than 1 nm were measured. Laboratory experiments confirm that displacement is measured correctly. These results suggest that fiber optic cable embedded in geologic formations may be used to map hydraulic connections in three-dimensional fracture networks. A great advantage of this approach is that strain, an indirect measure of hydraulic stress, can be measured without beforehand knowledge of flowing fractures that intersect boreholes. The technology has obvious applications in water resources, geothermal energy, CO 2 sequestration, and remediation of groundwater in fractured bedrock.

54 ENVIRONMENTAL SCIENCES↗

Potential Chemical Impacts of Subsurface CO 2 : An Integrated Experimental and Numerical Assessment for a Case Study of the Ogallala Aquifer

Leakage from geologic CO 2 sequestration (GCS) reservoirs to overlying underground sources of drinking water (USDW) is a tangible risk. This study is an integrated assessment that combines column experiments and reactive transport simulations of sediments sampled from the Ogallala aquifer above an active commercial-scale GCS site (the Farnsworth Unit in northern Texas). Experimental and simulation results suggest that carbonate mineral (calcite and dolomite) dissolution is the most significant reaction following CO 2 intrusion, and is also the dominant source of trace metal release. Cation exchange is another key mechanism controlling trace metal release by cation interference. Most of the trace metals, including Ba, Sr, As, Pb, and Zn, show a short-term release and quickly drop to the baseline values, suggesting low risk to the overlying USDW quality. Other trace metals, such as Mn and U, exhibit a tangible increase of their concentrations in the beginning, and drop to a higher level compared to the baseline, which may become a potential concern of long-term USDW quality change with CO2introduction. This study provides a comprehensive example of combining laboratory experiments and simulations for assessment of CO 2 -sediment interactions with combined release mechanisms in shallow groundwater aquifers. Data presented here provides useful insights for quantitative risk assessment and effective public education regarding CO 2 geological sequestration, and trace metal reactive transport studies in shallow groundwater aquifers.

54 ENVIRONMENTAL SCIENCES↗

Preferential Solute Transport in Low Permeability Zones During Spontaneous Imbibition in Heterogeneous Porous Media

Multiphase flows in porous media|and the associated solute transport processes|are controlled by a combination of gravity, capillary, and viscous forces. Geologic heterogeneity influences flow and transport processes, yet gaps exist in our understanding of how heterogeneity impacts capillary-driven transport, such as during spontaneous imbibition. Here we use positron emission tomography, combined with a newly developed method for conducting spontaneous imbibition experiments, to observe solute transport into low permeability regions during imbibition. Using an experimentally parameterized 2D numerical model, we demonstrate how capillary-driven ow near the imbibition front can carry solutes into low permeability regions. This process displaces solutes from high permeability zones while the cumulative amount of solute in low permeability zones increases- opposite of what is observed under fully saturated solute transport conditions. These results highlight the complex flow and solute transport behavior that arise during multiphase displacements in heterogeneous geologic porous media.

58 GEOSCIENCES↗

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗

Strong Upper-Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

We use earthquake-based adjoint tomography to invert for three-dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land-based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L-BFGS optimization algorithm, improving data fit and introducing P- and S-wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. Here, the inversion reveals increased along-strike heterogeneity on the margin. In Cook Strait we observe a low-velocity zone interpreted as deep sedimentary basins. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone.

58 GEOSCIENCES↗