Search NASA⌕ Search

Engineering topics

Chiang, Andrea

Publications and source records attributed to Chiang, Andrea.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probabilistic Seismic Hazard Assessment for Azerbaijan

Probabilistic Seismic Hazard Assessments (PSHA) underpin the determination of seismic loads in most contemporary seismic provisions of building codes around the world. Modern building codes are migrating towards using the entire uniform hazard spectrum at a range of vibration periods (typically up to 4 s or 10 s) rather than a single peak value such as peak ground acceleration (PGA), thus requiring a larger range of PSHA outputs. The hazard maps in the current version of the building code of Azerbaijan (2011) are in terms of intensity and peak ground acceleration. Recognizing the need for an up-to-date seismic hazard assessment in the country, Seismic Cooperation Program (SCP) under the Lawrence Livermore National Laboratory (LLNL) undertook, in coordination with the Republican Seismic Survey Center of the Azerbaijan National Academy of Sciences and the Azerbaijan Scientific Research Institute of Construction and Architecture, a PSHA study that reflects new seismic data recorded locally and recent research conducted nationally and regionally since the last update to the building code. The PSHA framework for this project was designed to help develop a new earthquake catalogue, to incorporate a novel characterization of ground motions that specifically reflects the attenuation characteristics in the eastern Caucasus, and to generate hazard information in a form that is useful for an update of the current building code or the development of a new building code for Azerbaijan. The project also aimed to provide training and support for the local seismologists and engineers related to the seismic hazard models, probabilistic seismic hazard results and their use towards changes in the building code.

58 GEOSCIENCES↗
DOE OSTI2024

Relocation of the 8 September 2023 High Atlas, Morocco, Earthquake Aftershock Sequence

The earthquake that occurred on 8 September 2023, with a magnitude of 6.8, was the most destructive earthquake event in Morocco in the past decade. This earthquake took place in the Al Haouz region, located in the western part of the High Atlas Mountain range. To better understand what caused and triggered this earthquake, the earthquake catalogs including P and S arrival times were collected from the Moroccan seismic network and combined with regional data from the International Seismological Centre. The mainshock and aftershocks were relocated by using iLoc, a state-of-the-art single-event location algorithm, and then by the multiple event location double-difference algorithm, hypoDD. The improved earthquake relocations using iLoc and the double-difference methods provide sharper lineation of seismicity and agree well with tomographic images of the earthquake zone. Finally, the seismicity distribution and the focal mechanism of the mainshock indicate that the earthquake sequence has occurred along the South Atlas fault system.

58 GEOSCIENCES↗
DOE OSTI2024

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗
DOE OSTI2024

Report on the sensitivity kernel construction and updating the SPiRaL model with waveform data

The WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes research leading towards the prediction of entire seismic and acoustic waveforms produced by natural and manmade events, including explosions. A key aspect of this research is the development of Earth (seismic) models at multiple scales including crustal, regional, and global. The work described here pertains to the global-scale seismic tomography effort led by LLNL.

58 GEOSCIENCES↗
DOE OSTI2024

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗
DOE OSTI2024

Central Asia Seismic Hazard Assessment (CASHA): A Probabilistic Seismic Hazard Assessment for Kazakhstan, Kyrgyzstan and Tajikistan

Probabilistic seismic hazard assessments (PSHA) underpin the calculation of earthquake loads in most building codes around the world. In Central Asia, the building codes are slowly being updated to incorporate some of the contemporary concepts of seismic hazard representation. There is also a regional desire to coordinate hazard assessments and building code modernization. However, some challenges remain. Expertise in the region related to seismic hazard assessments is still largely compartmentalised, requiring a significant amount of training and capacity building in seismic hazard assessment related topics. In addition, there are vast amounts of seismic data (bulletin and waveforms), both from analogue and digital eras, that the region’s countries stored but until recently did not use or share among themselves or with the broader seismological community around the world. Finally, after the collapse of the Soviet Union in the 1990s, many of the countries’ seismic networks suffered a major setback with the lack of attention and budget to update existing equipment and installation of new instruments. In order to address these issues, the United States Department of Energy through Lawrence Livermore National Laboratory (LLNL) initiated a project in 2016 to engage and train local scientists in Central Asia to install new equipment, to enhance the quality of seismic monitoring and reporting, to improve and harmonise the regional earthquake catalogue, and to conduct national probabilistic seismic hazard assessments using the new and improved datasets. To achieve the seismic hazard assessment related goals, a series of workshops were held in Almaty, Kazakhstan; Bishkek, Kyrgyzstan; and Dushanbe, Tajikistan from 2016 until 2020. During the time that the COVID-19 pandemic restricted travel, workshops continued online (22 online workshops were hosted in two years). Finally, in May 2022, an in-person workshop in Istanbul, Turkey brought together all project participants along with civil engineers engaged with building code activities in their respective countries, providing a platform to discuss the implementation of the hazard models into updates of building codes in each country, as well as to discuss model parameters, sensitivity analyses and model results in terms of hazard maps, uniform hazard spectra and hazard deaggregations. The workshops were a combination of lectures and hands-on exercises, and included international participation as well as local scientists and engineers. The workshops served several purposes, including training, coordination of data collection, interactions between local earth scientists and engineers, and brainstorming and knowledge exchange among local and international experts. This report outlines the new earthquake catalogue compilation effort and the PSHA project undertaken in Kyrgyzstan, Tajikistan, and Kazakhstan as part of this initiative. The southern part of this region is tectonically active with moderate to high levels of both shallow crustal seismic activity and occurrence of deeper earthquakes under the Hindu Kush and Pamir mountain ranges. Deeper earthquakes also occur near southwestern Kazakhstan, under the eastern Greater Caucasus and Caspian Sea. Large portions of central and northern Kazakhstan, on the other hand, are in stable continental regions with low levels of seismic activity. This study systematically compiles and improves all available data on local seismicity, active faults, and ground motion attenuation characteristics of the region; and builds a framework to enable a contemporary PSHA to be carried out with the engagement of local scientists. While the project was regional, the seismic hazard assessments are primarily driven by the countries’ own national preferences and understanding of data collection, interpretation, and validation of results.

58 GEOSCIENCES↗
DOE OSTI2024

Compilation of a Comprehensive Earthquake Catalog and Relocations in the Caucasus Region

Instrumental seismic monitoring has a long history in the Caucasus and started in 1899 when the first seismograph was installed in Tbilisi, Georgia. Much of the analog paper records from this time period are preserved in the Tbilisi archives because Georgia served as the regional data center. In the 1990s, due to the collapse of the Soviet Union and the political turmoil in the region, the analog networks and the communication between the newly formed national networks deteriorated. In Georgia, for the next 13 yr, the seismic network coverage was poor until the 2002 Tbilisi earthquake. Following this earthquake, the first permanent digital seismic station in Georgia was established in Tbilisi in 2003. The digital era progressively improved the ability to collect and archive data and today more than a hundred broadband seismic stations (including temporary arrays) are operating in the southern Caucasus. Until recently, the region lacked a coordinated effort to catalog all analog and digital era data collected by different countries into a single repository. As a result of collaboration between Lawrence Livermore National Laboratory, the Ilia State University, and the Republican Seismic Survey Center of Azerbaijan, a comprehensive earthquake catalog was compiled for the Caucasus and neighboring areas as part of a broader probabilistic seismic hazard assessment project. Here this project digitized Soviet-era paper bulletins, compiled a unified earthquake catalog from regional bulletins, developed 1D reference velocity model, and used it to relocate the events. The final catalog contains 16,963 events with magnitudes 3.7 and above, bringing together all the available data sets in the Caucasus region from 1900 to 2015, significantly improving locations, and generating the most complete earthquake catalog in the region, temporally and geographically.

58 GEOSCIENCES↗
DOE OSTI2023

CANVAS: An Adjoint Waveform Tomography Model of California and Nevada

Abstract We present the California‐Nevada Adjoint Simulations (CANVAS) model, an adjoint waveform tomography model of the crust and uppermost mantle of the states of California and Nevada. We used WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549 ) as the starting model and iteratively decreased the minimum period of CANVAS from 30 to 12 s. CANVAS was iterated in two distinct stages: the first stage with source mechanisms from the Global Centroid Moment Tensor (GCMT) catalog and the second stage with inverted moment tensors (MT) using the CANV_WUS model (Doody et al., 2023, https://doi.org/10.1029/2023jb026463 ). We show that updating the MTs with 3D Green's functions improved waveform fits and azimuthal coverage of windowed data used to calculate the gradients. As for the model itself, we improved waveform fits over WUS256, particularly in the dispersed surface waves. CANVAS resolved tectonic features seen in other models and accurately defined the depth to basement of major basins, including the Central Valley and the Ventura Basin. We propose CANVAS as a starting model for crustal tomography models on smaller scales.

58 GEOSCIENCES↗
DOE OSTI2023

Relocation of the Seismicity of the Caucasus Region

Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.

58 GEOSCIENCES↗
DOE OSTI2023

Bayesian inference for the seismic moment tensor using regional waveforms and a data-derived distribution of velocity models

The largest source of uncertainty in any source inversion is the velocity model used to construct the transfer function employed in the forward model that relates observed ground motion to the seismic moment tensor. However, standard inverse procedures often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. We attempt to incorporate this uncertainty into an estimation of the seismic moment tensor using a distribution of velocity models calculated in a prior effort based on different and complementary data sets. The posterior distribution of velocity models is then used to construct Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional waveform data. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and can be interpreted in terms of most-probable source-type.

58 GEOSCIENCES↗
DOE OSTI2023

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗