Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earthquakes”

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 451 records · Page 25

Post-hazard condition assessment of nuclear piping-equipment systems: Novel approach to feature extraction and deep learning

Over the past decade, the use of artificial intelligence techniques in the field of health-monitoring has gained significant interest, especially for structures such as building and bridges. However, applications to industrial systems such as equipment-piping systems in nuclear plants have not been explored. Here, in this paper, it is shown that the existing techniques developed for buildings and bridges cannot be extended directly to equipment-piping systems as the response of such systems is governed by multiple localized modes unlike that in buildings and bridges. This paper proposes a new approach that consists of three key aspects: (i) a novel vector of degradation-sensitive features extracted from measured data, (ii) using a deep Artificial Neural Network (ANN) for diagnosis of degradation location and degradation severity, and (iii) consideration of uncertainty in degradation severity when training the ANN. Degradation in piping-equipment systems can occur due to flow-accelerated erosion and corrosion. These locations can potentially exhibit damage such as localized yielding or initiation of cracking due to an external event such as an earthquake. Moreover, such locations can at times go undetected by current inspection techniques. Therefore, a robust framework is needed for detection of degradation after a seismic event. This manuscript proposes a proof-of-concept framework, which utilizes data collected from sensors to generate a deep ANN database for predicting degraded locations and severity in a piping-equipment system. Degradation severity is classified as minor, moderate, and severe. In the suggested methodology, a novel vector of degradation-sensitive features is extracted from the sensor data to train the ANN. A simple piping-equipment system is selected to demonstrate feature extraction as a means to simplify pattern recognition, explore the design and parameters of an ANN, and develop a sensor placement strategy. The effectiveness of the proposed framework is demonstrated on a realistic primary safety system of a two-loop nuclear reactor. It is shown that the proposed post-hazard condition assessment framework is able to detect degraded locations along with the severity levels, including minor degradation, with considerably higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗

Seismic tomography of compressional wave velocity and attenuation structure for Makushin Volcano, Alaska

Located in the central portion of the Alaska-Aleutian subduction zone, Makushin Volcano is among the most active volcanoes in the United States and has been classified as a very high threat by the National Volcano Early Warning System based on its eruptive history and proximity to the City of Unalaska and international air routes. In 2015, five standalone seismic stations and three mini seismic arrays of 15 stations each were deployed on Unalaska Island to supplement the Alaska Volcano Observatory permanent seismic network. This temporary array was operational for one year. Taking advantage of the increased azimuthal coverage and the array's increased earthquake detection capability, we developed body-wave Vp and Vp/Vs seismic images of the velocity structure beneath the volcano. Our results show a complex structure with the upper 4 km of the crust. The shallow high-Vp features possibly delineate remnant magma pathways or conduits. Low-Vp regions are found east of the caldera at approximately 4–8 km depth. This is in agreement with previous seismic tomographic work and InSAR-based geodetic models, which had identified this region as a possible long-term source of magma. A three-dimensional compressional wave attenuation model complemented the velocity models. In general, we observe regions of corresponding low P-wave velocity and high attenuation east of the summit caldera between 5 and 7 km depth, supporting the presence of a melt-rich zone. The absolute Vp values (5.5–6.0 km/s) in this region are consistent with the presence of a magma body with a mafic signature. Here, we suggest that a triggering mechanism that involves mafic injections into a more evolved shallow magma chamber is representative of the current magmatic system below the volcano, similarly to what petrological studies have inferred for past eruptions.

58 GEOSCIENCES↗

Recent seismicity at Ceboruco Volcano (Mexico)

Magma and related hydrothermal fluid movement, and their interaction with solid rock, in active volcanic regions, generate a wide variety of seismic waves whose characterization can mitigate the risk of a potential eruption. Located in the western region of the Trans-Mexican Volcanic Belt, Ceboruco Volcano, whose last eruptive period was 1870–1875, is considered to be one of the most hazardous volcanoes in Mexico. We have conducted a detailed study of the seismicity in the surroundings of Ceboruco's volcanic edifice to assess the current state of this volcano. A dense temporary seismic network with 25 seismic stations in an area of 16 km × 16 km was deployed between November 2016 and July 2017, as part of the P-24 project of the CeMIEGeo consortium; this effort has allowed the detection of 81 earthquakes concentrated beneath the crater with depths between 4 and 8 km. In this study, we observe that the recorded seismicity occurs in swarms, and we specifically identify four sequences that we characterize in detail via the first focal mechanisms available for this volcano. Our results suggest a change in the local seismicity distribution compared to earlier observations, which reported seismic activity near the volcano edifice associated with fluid migration along zones of weakness related to the extensional stresses of the Tepic-Zacoalco rift. The changes in seismic patterns and obtained focal mechanisms are consistent with observed fluid effects at many geothermal sites worldwide, but also could suggest resumption of activity at this currently dormant volcano.

58 GEOSCIENCES↗

Hydrocarbon production reduces natural methane seeps in the Santa Barbara channel

Abundant natural hydrocarbon seepage in the Santa Barbara Channel includes large amounts of gas, mainly methane. Numerous seeps are focused along a faulted anticlinal axis where in 1982 ARCO placed two large (900 m 2 each) steel gas capture tents on the seabed adjacent to Platform Holly. The seep rate into the tent at high tide is about 5 percent less compared with low tide (2m tidal range). An earthquake in 2013 (M = 4.8) did not affect seepage rates into the steel tents, indicating that gas was not released by fault dilation. The CO 2 and H 2 S content of the captured gas are anomalously low compared to typical Holly gas due to bacterial and sea water-gas interaction at the seep tent. Platform Holly production had an overall effect of reducing gas seepage into the tents over time. Wells beneath the tents at 1-km depth had a much greater influence on the seepage rate than more distal wells. The nearby wells had a direct effect on seepage rate when they were put on production (seepage decreased) versus shut in (seepage increased). The response time was within hours, indicating a direct hydraulic connection between the reservoir at 1 km subsea and the seabed. Permeability of the approximate 1-km fault pathway between the reservoir and the seep tent is calculated to be about 3 × 10 -14 m 2 (30 millidarcys). Gradually, the seep rate declined from peak flow in 1990 (4.5×10 4 m 3 dy-1) to less than 7.1×10 3 m 3 dy -1 in 2013. In June of 2013, a horizontal well was completed high on the structural anticline between Platform Holly and directly beneath the seep tents. Within three weeks of turning on this new well, seepage into the tents completely ceased, ending more than 30 years of natural methane seepage captured by the tents.

58 GEOSCIENCES↗

Structural modeling and dynamic analysis of condensate storage tanks in nuclear power plants

Condensate storage tanks (CSTs) in nuclear power plants (NPPs) are classified as critical equipment capable of surviving strong shaking in a design basis earthquake to assure the ability to subsequently provide cooling water. Structural modeling and dynamic analysis of CSTs are complicated due to fluid-structure interaction (FSI) and coupling issues, which could have significant effects on seismic response. The capability of a number of approximate models to capture seismic response is investigated using dynamic analysis results from detailed finite element models. Modal analyses and time history analyses are carried out using both 2D and 3D models to predict the dynamic behavior of CSTs as a function of ground motion intensity. Although 2D simplified models can be used to quickly evaluate the dynamic response of CSTs when there is substantial margin to failure, a detailed 3D model is required when it is necessary to examine a limit state associated with the failure mode of the tank. The findings show good agreement between the natural frequencies of convective modes determined from the 2D and 3D CST models. The difference between the frequencies of impulsive modes were found to be equal to or less than 25%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Seismic isolation: A pathway to standardized advanced nuclear reactors

Standardizing advanced nuclear reactors is a pathway to substantially reducing their overnight capital cost and achieving parity with other power sources, including renewables and fossil fuels. The seismic load case has thwarted standardization of nuclear power plants because site-specific seismic hazard and local near-surface geology has triggered soil-structure-interaction analysis, design, equipment qualification, regulatory review, and licensing, ensuring that each build is different. To achieve standardized or site-independent certified advanced reactor designs, the impact of the seismic load case on the engineering and construction cost and time must be substantially mitigated. Seismic isolation is a mature technology that has been used for more than 30 years in non-nuclear sectors to substantially reduce earthquake demands in buildings and other infrastructure. In this paper, seismic isolation is used to enable standardization of advanced reactor designs, aimed at the complete re-use of a site-independent, certified design and repeated procurement of safety-class equipment. A pathway to standardized designs using seismic isolation is demonstrated for two fundamentally different advanced reactors: a molten salt reactor and a high temperature gas reactor. Each reactor building is equipped with three specialized pieces of safety-class equipment, namely, a reactor vessel, a steam generator, and a control rod drive mechanism housing that is attached to the reactor head. Analysis is performed per ASCE and ASME standards to design the buildings and the equipment for two base conditions: conventional (fixed base) and base isolated. The impact of the seismic load case is characterized for the reinforced concrete walls in the buildings and for the equipment, measured using vessel wall thickness and horizontal accelerations. Here the analysis results show that the fixed-base buildings, designed for a site of low seismic hazard (peak ground acceleration, PGA = 0.15 g) could be constructed at a site of much greater seismic hazard (PGA = 0.7 g) if seismic base isolation is employed. Importantly, the scope of the site-specific analysis, design, and qualification would be limited to the seismic isolators and the isolated substructure, drastically reducing plant-specific engineering, review, and licensing, and time to construction start. Regulatory challenges and opportunities with standardized reactor designs are identified.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mid-height seismic isolation of equipment in nuclear power plants: Numerical simulations and design recommendations

Seismic isolation can mitigate earthquake demands on safety-class equipment and is being considered for application to some advanced nuclear reactors. For tall, slender vessels that could represent an advanced reactor, a steam generator or a heat exchanger, mid-height seismic isolation has been shown to be beneficial and practical. Here, this paper focuses on numerical modeling of a mid-height isolated tall, slender vessel with two primary goals 1) investigate and quantify the benefits of mid-height seismic isolation for a range of support structure stiffnesses, isolation systems, and seismic inputs, and 2) provide analysis and design recommendations for isolation of safety-class equipment. Results of response-history analysis were in good agreement with experimental measurements and demonstrated that mid-height isolation can substantially reduce seismic demands on tall, slender vessels for a range of support structure stiffnesses, isolation systems, and seismic inputs. Importantly, the reductions in horizontal spectral accelerations in the mid-height isolated vessel from the non-isolated condition were not affected by the stiffness of the support structure. Recommendations are made for analysis and design of isolated equipment and testing of isolators used for equipment protection in nuclear facilities, which are also applicable to other industries.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Series of semi-Markov processes to model infrastructure resilience under multihazards

Civil infrastructure systems are subjected to multiple hazards, including natural and anthropogenic, that disrupt their function or the level of service offered. Estimating the function recovery of these systems (or how soon normalcy of operations will be restored) when subjected to repeated hazard events by considering the inter-event dependencies is an important problem in multihazard infrastructure resilience. However, this problem has been less addressed in the field. This paper proposes a series of semi-Markov processes model to capture the inter-event dependencies in infrastructure recovery when subjected to successive hazard events. Recovery after each new hazard event is represented by a unique semi-Markov process that models the reduced recovery rates and the increased recovery times caused by the system’s incomplete recovery from the preceding event. Two novel formulations of the inter-event dependency modeling, namely Maximal Effects Dependency (considers the worst impact of two successive hazard events) and Cumulative Effects Dependency (considers the aggregated impacts of two successive hazard events), are proposed and discussed. The model is demonstrated by considering the following applications: Three-state system subjected to deterministic and random occurrences of identical hazard events; and Multihazard resilience of a building in Charleston, SC, considering earthquake and hurricane hazards. Results indicate that considering inter-event dependencies in recovery modeling can lead to lesser-predicted resilience, thereby affecting resilience-based decision-making.

42 ENGINEERING↗

Resilience-based Recovery Scheduling of Transportation Network in Mixed Traffic Environment: A Deep-Ensemble-Assisted Active Learning Approach

Devising effective post-hazard recovery strategies is critical in enhancing the resilience of transportation networks (TNs). However, existing work does not consider the multiclass users’ travel behavior in network functionality quantification and the metaheuristic solution procedures often suffer from extensive computational burden due to the exploration need in large solution space and the expensive functionality quantification. Here, this study develops a bilevel decision-making framework for the resilience-based recovery scheduling of the TN in a mixed traffic environment with connected and autonomous vehicles (CAVs) and human-driven vehicles (HDVs). The lower level quantifies the TN's functionality over time considering different travel behavior of CAV and HDV users arisen from their different levels of information perception. The upper level presents a novel deep-ensemble-assisted active learning approach to balance optimization performance and computational cost. This framework can help decision makers better quantify the TN's functionality to support effective recovery scheduling of TN with different mixed traffic scenarios ranging from HDV-only to future CAV-dominant traffic. The optimization approach bears the potential to be extended to solving general large-scale network recovery scheduling problems effectively and efficiently. The proposed methodology is demonstrated using a real-world traffic network in Southern California under earthquake considering deterministic and stochastic repair durations.

42 ENGINEERING↗

Time domain probabilistic seismic risk analysis using ground motion prediction equations of Fourier amplitude spectra

Modeling of Fourier amplitude spectra (FAS) of seismic motions has gained much attention in engineering seismology. In the past few years, several ground motion prediction equations (GMPEs) and inter-frequency correlation structure of FAS have been established. Due to many preferable characteristics of FAS, probabilistic seismic hazard/risk analysis is rapidly changing from ergodic, spectrum acceleration Sa(T 0 )-based approach to non-ergodic, site-specific, FAS-based approach. This paper presents time domain intrusive framework for probabilistic seismic risk analysis using GMPE of FAS. Herein, methodology for time domain stochastic ground motion modeling based on GMPEs of FAS is presented in some detail. The simulated uncertain motions are modeled as a random process and represented by polynomial chaos Karhunen-Loève expansion. The random process excitations are further propagated into the uncertain structural system using Galerkin stochastic finite element method (SFEM). Probabilistic evolution of structural response is solved, and such solution is used to develop seismic risk for any damage state. The presented framework is illustrated through seismic risk analysis of a four-story building subjected to possible earthquakes from two strike slip faults. The influences of the epistemic uncertainties in source stress drop Δσ and site attenuation κ0 on seismic risk are investigated. The need for non-ergodic seismic risk analysis with source-specific and site specific characterizations is emphasized.

58 GEOSCIENCES↗

A Markov framework for generalized post-event systems recovery modeling: From single to multihazards

State-dependent models can be used to represent the system recovery process as a series of stochastic transitions from lower to higher functional states. However, the applications of these models have been limited in scope and there is a lack of a generalized recovery modeling framework. A generalized framework would permit a robust forecasting of systems and system-of-systems recovery under multiple hazards, and more broadly, would contribute to community disaster preparedness. This paper develops a generalized post hazard-event recovery modeling framework based on state-dependent Markov-type processes. We then apply the proposed framework to solve a spectrum of problems that range from hind-casting single-system recovery following a single hazard event to forecasting post-event trajectories under multiple hazards and modeling the recovery of a system-of-systems. First, Markov chains are used to hind-cast the observed recovery for a portfolio of buildings affected by the 2014 South Napa, California, earthquake. Next, Markov processes are used to formulate a parametric post hazard-event recovery model, which can be updated using Bayesian statistics when relevant datasets become available. Semi-Markov processes are then used to develop a more general model of single hazard recovery, which accounts for the intensity of the loading and level of damage caused by the event. Semi-Markov processes with non-renewal features are then used to account for multihazard interactions in a post-event recovery model, and applied to a case study that involves a community in Charleston, South Carolina. Lastly, Markov-type processes are combined with Bayesian networks to model the recovery of residential, commercial, educational, and industrial buildings (system-of-systems) following a hazard event. Overall, these applications demonstrate the versatility of the Markov framework towards handling recovery problems with varying levels of complexity.

42 ENGINEERING↗

Seasonal and Long‐Term Groundwater Unloading in the Central Valley Modifies Crustal Stress

Abstract Changes in terrestrial water content cause elastic deformation of the Earth's crust. This deformation is thought to play a role in modulating crustal stress and seismicity in regions where large water storage fluctuations occur. Groundwater is an important component of total water storage change in California, helping to drive annual water storage fluctuations and loss during periods of drought. Here we use direct estimates of groundwater volume loss during the 2007–2010 drought in California's Central Valley obtained from high resolution Interferometric Synthetic Aperture Radar‐based vertical land motion data to investigate the effect of groundwater volume change on the evolution of the stress field. We show that GPS‐derived elastic load models may not capture the contribution of groundwater to terrestrial water loading, resulting in an underestimation of nontectonic crustal stress change. We find that groundwater unloading during the drought causes Coulomb stress change of up to 5.5 kPa and seasonal fluctuations of up to 2.6 kPa at seismogenic depth. We find that faults near the Valley show the largest stress change and the San Andreas fault experiences only ~40 Pa of Coulomb stress change over the course of a year from groundwater storage change. Annual Coulomb stress change peaks dominantly in the fall, when the groundwater level is low; however, some faults experience peak stress in the spring when groundwater levels are higher. Additionally, we find that periods of increased stress correlate with higher than average seismic moment release but are not correlated with an increase in the number of earthquakes. This indicates groundwater loading likely contributes to nontectonic loading of faults, especially near the Valley edge, but is not a dominant factor in modulation of seismicity in California because the amplitude of stress change declines rapidly with distance from the Valley. By carefully quantifying and spatially locating groundwater fluctuations, we will improve our understanding of what drives nontectonic stress and forces that modulate seismicity in California.

58 GEOSCIENCES↗

Identifying Different Classes of Seismic Noise Signals Using Unsupervised Learning

Abstract Proper classification of nontectonic seismic signals is critical for detecting microearthquakes and developing an improved understanding of ongoing weak ground motions. We use unsupervised machine learning to label five classes of nonstationary seismic noise common in continuous waveforms. Temporal and spectral features describing the data are clustered to identify separable types of emergent and impulsive waveforms. The trained clustering model is used to classify every 1 s of continuous seismic records from a dense seismic array with 10–30 m station spacing. We show that dominate noise signals can be highly localized and vary on length scales of hundreds of meters. The methodology demonstrates the complexity of weak ground motions and improves the standard of analyzing seismic waveforms with a low signal‐to‐noise ratio. Application of this technique will improve the ability to detect genuine microseismic events in noisy environments where seismic sensors record earthquake‐like signals originating from nontectonic sources.

Johnson, Christopher W.↗

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↗

The Preservation of Climate‐Driven Landslide Dams in Western Oregon

Abstract Bedrock landsliding, including the formation of landslide dams, is a predominant geomorphic process in steep landscapes. Clarifying the importance of hydrologic and seismic mechanisms for triggering deep‐seated landslides remains an ongoing effort, and formulation of geomorphic metrics that predict dam preservation is crucial for quantifying secondary landslide hazards. Here, we identify >200 landslide‐dammed lakes in western Oregon and utilize dendrochronology and enhanced 14 C dating (“wiggle matching”) of “ghost forests” to establish slope failure timing at 20 sites. Our dated landslide dataset reveals bedrock landsliding has been common since the last Cascadia Subduction Zone earthquake in January 1700 AD. Our study does not reveal landslides that date to 1700 AD. Rather, we observe temporal clustering of at least four landslides in the winter of 1889/1890 AD, coincident with a series of atmospheric rivers that generated one of the largest regionally recorded floods. We use topographic and field analyses to assess the relation between dam preservation and topographic characteristics of the impounded valleys. In contrast to previous studies, we do not observe systematic scaling between dam size and upstream drainage area, though dam stability indices for our sites correspond with “stable” dams elsewhere. Notably, we observe that dams are preferentially preserved at drainage areas of ∼1.5 to 13 km 2 and valley widths of ∼25 to 80 m, which may reflect the reduced downstream influence of debris flows and the accumulation of mature conifer trees upstream from landslide‐dammed lake outlets. We suggest that wood accumulation upstream of landslide dams tempers large stream discharges, thus inhibiting dam incision.

Struble, William T.↗

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↗