Finite-fault Model Research Articles

The 2017 Ischia Earthquake (Southern Italy): Source Mechanism and Rupture Model From the Inversion of a Near-Source Strong Motion Record

With the aim to investigate the rupture complexity and the radiated wave field of the 2017, Mw 3.9, Ischia earthquake, south-west of Naples (Italy), we used finite-fault modeling to invert the near-source (< 1-km epicentral distance) horizontal-component velocity records of the accelerometric station (IOCA) and searched for the best-fit kinematic rupture parameters. This analysis showed that the rupture nucleated at about 600 m west of IOCA and 1.1-km depth, along a 1 km, NW-SE striking fault (i.e., thrust with right-lateral component), with a rupture velocity of about 0.7 km/s. The retrieved rupture model coupled with multipath reverberations effects related to a thin, low-velocity near-surface volcanic sedimentary layer, well explains the observed long ground motion duration and the large amplitudes recorded all over the island. Finally, the apparent source time function (STF), obtained from inverse modeling using a theoretical Green’ function (GF), is validated by implementing an empirical GF (EGF) analysis.

Earthquake Rupture and Tsunami Generation of the 2015 Mw 5.9 Bonin Event Revealed by In Situ Pressure Gauge Array Observations and Integrated Seismic and Tsunami Wave Simulation

AbstractOn September 1, 2015, an Mw 5.9 interplate earthquake occurred near the Bonin Trench. An array of in situ ocean‐bottom absolute pressure gauges (APGs) observed its tsunami generation field consisting of static and dynamic pressure changes due to tsunami and crustal deformation and due to seismic motion, respectively, under much higher station density than ever reported. We propose an approach to synthesize the pressure change inside the focal area, which reproduces the APG waveforms well. We further successfully estimate the finite fault model of the earthquake and constrain the rupture duration only from the APG data. The relatively low stress drop and seismic wave radiation from the fault model may suggest the nature of a tsunami earthquake. The in situ APGs have a large potential to reveal the broadband source process of the earthquake, which is essential to understand the seismotectonics in the subduction zone.

Coseismic and Early Postseismic Deformation Due to the 2021 M7.4 Maduo (China) Earthquake

AbstractThe 2021 Maduo earthquake ruptured a 150 km‐long left‐lateral fault in the northeast Tibet. We used Synthetic Aperture Radar data collected by the Sentinel‐1A/B satellites within days of the earthquake to derive a finite fault model and investigate the details of slip distribution with depth. We generated coseismic interferograms and pixel offsets from different look directions corresponding to the ascending and descending satellite orbits. At the eastern end the rupture bifurcated into two sub‐parallel strands, with larger slip on the northern strand. Inversions of coseismic displacements show maximum slip to the east of the epicenter. The averaged coseismic slip has a peak at depth of 3–4 km, similar to slip distributions of a number of shallow strike‐slip earthquakes. Postseismic observations over several weeks following the Maduo earthquake reveal surface slip with amplitude up to 0.1 m that at least partially eliminated the coseismic slip deficit in the uppermost crust.

The March 2021 Tyrnavos, central Greece, doublet (Μw6.3 and Mw6.0): Aftershock relocation, faulting details, coseismic slip and deformation

On 3 March 2021, the Mw6.3 Tyrnavos earthquake shook much of the Thessalia region, leading to extensive damage in many small towns and villages in the activated area. The first main shock was followed in the next day, on 4th of March 2021, by an “equivalent” main shock with Mw6.0 in the adjacent fault segment. These are the largest earthquakes to strike the northeastern part of Thessalia since the M6.3, 1941 Larissa earthquake. The main shocks triggered extensive liquefaction mainly along the banks of the Titarisios tributary where alluvial flood deposits most probably amplified the ground motions. Our seismic monitoring efforts, with the use of recordings of the regional seismological network along with a dense local network that was installed three days after the seismic excitation initiation, led to the improved understanding the geometry and kinematics of the activated faults. The aftershocks form a north–northwest–trending, east–northeast–dipping, ~40 km long distribution, encompassing the two main ruptures along with minor activated structures, consistent with the rupture length estimated from analysis of regional waveform data and InSAR modeling. The first rupture was expanded bilaterally, the second main shock nucleated at its northern tip, where from this second rupture propagated unilaterally to the north–northwest. The focal mechanisms of the two main shocks support an almost pure normal faulting, similar to the aftershocks fault plane solution determined in this study. The strong ground motion of the March 3 main shock was computed with a stochastic simulation of finite fault model. Coseismic displacements that were detected using a dense GPS / GNSS network of five permanent stations located the Thessaly region, have shown an NNE–SSW extension as expected from the nature and location of the causative fault. Coulomb stress changes due to the coseismic slip of the first main shock, revealed that the hypocentral region of the second main shock was brought closer to failure by more than 10 bars.

Indirect boundary element method solution to the seismic ground motion of near-fault sedimentary valley

In this paper, the indirect boundary element method (IBEM) is developed to solve the ground motion effect of the near-fault complex site. By establishing a two-dimensional kinematic finite fault model, the amplification effect of seismic response of sedimentary valley under the continuous dislocation of reverse faults is quantitatively analyzed. Meanwhile, the influence of the parameters such as the buried depth of the fault top, the fault dip-angle, the fault distance of sedimentary valley, the fault-plane shape, and the rupture velocity of fault on the seismic response are studied. The results show that the sedimentary valley has a pronounced amplification effect on the seismic response induced by near-fault, and the peak value of the acceleration response spectrum of the analyzed model can be amplified by 99.3%. In the sedimentary valley, the duration of earthquake motion is significantly prolonged, and the long-period velocity pulse appears, and the first velocity pulse can be amplified by 68.7%. The near-fault ground motion damages have typical concentration-effect and hanging wall effect. This study provides a new and effective method, which is of great significance for seismic zoning of complex sites near-fault. It will be of guiding significance for regional planning and seismic fortification of near-fault sites.

Inferring Critical Slip-Weakening Distance from Near-Fault Accelerogram of the 2014 Mw 6.2 Ludian Earthquake

Abstract To better assess potential earthquake hazards requires a better understanding of fault friction and rupture dynamics. Critical slip-weakening distance (Dc) as one of the key friction parameters, however, is hard to determine on natural faults. For strike-slip earthquakes, we may directly estimate the Dc from Dc″—the double near-fault ground displacement at the time of the peak velocity (Fukuyama and Mikumo, 2007). Yet near-fault observations are very few, and, thus, there were only limited earthquakes with such Dc″ estimation. In 2014, an Mw 6.2 strike-slip event—the Ludian earthquake—occurred in southwest China. The strong-motion station (LLT) that is ∼0.45 km from the fault recorded the earthquake and enabled us to estimate Dc″ from the accelerograms. We inspect the polarity of the accelerometers and compare the integrated velocities with waveforms of nearby broadband stations. We also analyze the particle motion at the LLT station and retrieve the earthquake initiation at the intersection of the conjugated faults. We then apply the baseline correction to the seismograms, recover the ground velocities and displacements, and obtain the value of Dc″=0.1 m at the station. The recovered final displacements are compared with the predicted ground displacements of a finite-fault model. The discrepancy of fault-parallel displacements might imply limited underestimation of Dc″, and the estimated upper limit is 0.3 m. Comparison between the Dc″ and final slip on the fault patch follows the scaling of previous larger earthquakes. Analysis of the near-fault accelerometer data enhances our understanding on the earthquake source of the Ludian earthquake. This case extends the lower magnitude boundary of the Dc″ values obtained from natural faults and opens a window into the friction property in the seismically active region.

Complex Source Behaviors and Spatiotemporal Evolution of Seismicity During the 2015–2016 Earthquake Sequence in Cushing, Oklahoma

AbstractIn September and October 2015, three M4+ earthquakes occurred as a sequence along a fault northwest of the Cushing city, Oklahoma, followed by another M5 earthquake in November 2016. While previous studies have shown that moderate‐size earthquakes in Oklahoma are likely induced by wastewater injections, it is still not clear what controls the rupture process and spatiotemporal evolutions of seismicity during individual sequences. In this study, we investigated the rupture process of these four M4‐5 events in 2015–2016 with finite fault model (FFM) inversions, and computed the static stress changes during this sequence. We found that the rupture processes of four M4‐5 earthquakes were very complex, and each of them had several subevents with different rupture directivities. The 2016 M5 earthquake started near the region where three M4+ events initiated, but the majority of the slip occurred a few kilometers away in the northeast direction. In comparison, the 2015 M4.3 event mainly ruptured toward the southwest direction. Due to data limitation and inversion uncertainties, we were unable to constrain the rupture directivities for the other two M4+ events. The foreshocks 3 days before the first M4+ earthquake in 2015 occurred in a region of positive shear stress changes caused by previous earthquakes in 2014–2015 on unmapped faults several kilometers to the south. Our results suggest small‐scale heterogeneity in controlling complex seismicity and rupture patterns in the 2015–2016 Cushing sequence, and possible triggering of this sequence by a small stress perturbation on order of a few kilopascals.

Probing the Fault Complexity of the 2017 Ms 7.0 Jiuzhaigou Earthquake Based on the InSAR Data

On 8 August 2017, a surface wave magnitude (Ms) 7.0 earthquake occurred at the buried faults extending to the north of the Huya fault. Based on the coseismic deformation field obtained from interferometric synthetic aperture radar (InSAR) data and a series of finite fault model tests, we propose a brand-new two-fault model composed of a main fault and a secondary fault as the optimal model for the Jiuzhaigou earthquake, in which the secondary fault is at a wide obtuse angle to the northern end of the main fault plane. Results show that the dislocation distribution is dominated by sinistral slip, with a significant shallow slip deficit. The main fault consists of two asperities bounded by an aftershock gap, which may represent a barrier. In addition, most aftershocks are located in stress shadows and appear a complementary pattern with the coseismic high-slip regions. We propose that the aftershocks are attributable to the background tectonic stress, which may be related to the velocity-strengthening zones.

Seismic zoning of Tabriz area by stochastic finite fault model considering site-specific soil effects

Abstract Tabriz, as one of the most earthquake-prone cities in the Iran plateau, has experienced enormous earthquakes that have even destroyed the city altogether. Considering this seismological background and the vicinity of Tabriz's northwestern fault, reducing the possible earthquake losses can be highly useful by scrutinizing the strong ground motion resulting from the fault activation. To this end, a stochastic finite-fault ground motion simulation (EXSIM) method was applied as an important means for predicting the ground motion near the epicenter of the earthquake. EXSIM is an open-source stochastic finite-fault simulation algorithm that generates the time-series of the earthquake's ground motion. Based on the findings, the peak horizontal acceleration reached 0.83 g in the northern parts by creating artificial accelerograms and Tabriz's seismic zonation. In comparison, it reduced by 0.48 g by departing from the fault in the city's southern parts. Additionally, providing a seismic zonation map in Tabriz revealed that stopping the construction in the north parts while extending the settlement construction to the south part of the city are considered vital and unavoidable. Also, by applying the magnification and effects of the soil layers above the bedrock, it was further found that the existence of the loose layer with low strength and compaction intensify the application of seismic acceleration on the near-surface structures in the central, west, and southwest parts of the target area.

Rapid estimation of tsunami earthquake magnitudes at local distance

A tsunami earthquake is an earthquake event that generates abnormally high tsunami waves considering the amplitude of the seismic waves. These abnormally high waves relative to the seismic wave amplitude are related to the longer rupture duration of such earthquakes compared with typical events. Rapid magnitude estimation is essential for the timely issuance of effective tsunami warnings for tsunami earthquakes. For local events, event magnitude estimated from the observed displacement amplitudes of the seismic waves, which can be obtained before estimation of the seismic moment, is often used for the first tsunami warning. However, because the observed displacement amplitude is approximately proportional to the moment rate, conventional magnitudes of tsunami earthquakes estimated based on the seismic wave amplitude tend to underestimate the event size. To overcome this problem, we investigated several methods of magnitude estimation, including magnitudes based on long-period displacement, integrated displacement, and multiband amplitude distribution. We tested the methods using synthetic waveforms calculated from finite fault models of tsunami earthquakes. We found that methods based on observed amplitudes could not estimate magnitude properly, but the method based on the multiband amplitude distribution gave values close to the moment magnitude for many tsunami earthquakes. In this method, peak amplitudes of bandpass filtered waveforms are compared with those of synthetic records for an assumed source duration and fault mechanism. We applied the multiband amplitude distribution method to the records of events that occurred around the Japanese Islands and to those of tsunami earthquakes, and confirmed that this method could be used to estimate event magnitudes close to the moment magnitudes.

Simulation of the MW = 7.9 Gulf of Alaska earthquake on January 23, 2018, by the stochastic finite fault model

On the premise of constraining the seismogenic fault structure of a future large earthquake, we proposed using empirical equations to determine the length, width, seismic moment and slip distribution of a large seismogenic fault plane and using the stochastic finite fault model to predict future large earthquakes. The ground motion time histories and response spectra recorded by 12 seismic stations on bedrock during the MW = 7.9 Gulf of Alaska earthquake on January 23, 2018, were simulated. The simulation error determined by the average ratio of the simulated spectrum amplitude to the recorded spectrum amplitude varied between 1.08 and 0.92 in the period range of 0–10 s, and the standard deviation of the simulation error at different frequencies did not exceed 1; the 95% confidence interval also did not change significantly with the period. The above analyses show that our simulation results reflect the mean ground motion. To further discuss the reliability of predicting future large earthquakes by the stochastic finite fault model, we redistributed the initial rupture point and slip distribution on the seismogenic fault plane by the quasi-random method, and the simulation errors and simulation results of the redistribution model were similar to those of the previous model. Further research confirmed that our method for obtaining the seismic source parameters is viable and that the stochastic finite fault model for the prediction of future large earthquakes is reliable, especially for large far-field earthquakes. The seismic stations that we used are all situated on bedrock on one side of the fault and do not involve rupture directivity, i.e., the seismic wave pathways may be similar, so the simulation results are ideal. However, if the rupture directivity, different site conditions, surface topography and basin effects are considered, it will be necessary to amend the proposed method.

Present-day orogenic processes in the western Kalpin nappe explored by interseismic GNSS measurements and coseismic InSAR observations of the 2020Mw 6.1 Kalpin event

SUMMARYAs the largest and most active intracontinental orogenic belt on Earth, the Tien Shan (TS) is a natural laboratory for understanding the Cenozoic orogenic processes driven by the India–Asia collision. On 19 January 2020, a Mw 6.1 event stuck the Kalpin region, where the southern frontal TS interacts with the Tarim basin. To probe the local ongoing orogenic processes and potential seismic hazard in the Kalpin region, both interseismic and instantaneous deformation derived from geodetic observations are used in this study. With the constraint of interseismic global navigation satellite system (GNSS) velocities, we estimate the décollement plane parameters of the western Kalpin nappe based on a 2-D dislocation model, and the results suggest that the décollement plane is nearly subhorizontal with a dip of ∼3° at a depth of 24 km. Then, we collect both Sentinel-1 and ALOS-2 satellite images to capture the coseismic displacements caused by the 2020 Kalpin event, and the interferometric synthetic aperture radar (InSAR) images show a maximum displacement of 7 cm in the line of sight near the epicentral region. With these coseismic displacement measurements, we invert the source parameters of this event using a finite-fault model. We determine the optimal source mechanism in which the fault geometry is dominated by thrust faulting with an E–W strike of 275° and a northward dip of 11.2°, and the main rupture slip is concentrated within an area 28.0 km in length and 10.3 km in width, with a maximum slip of 0.3 m at a depth of 6–8 km. The total released moment of our preferred distributed slip model yields a geodetic moment of 1.59 × 1018 N$ \cdot $m, equivalent to Mw 6.1. The contrast of the décollement plane depth from interseismic GNSS and the rupture depth from coseismic InSAR suggests that a compression still exists in the Kalpin nappe forefront, which is prone to frequent moderate events and may be at risk of a much more dangerous earthquake.

Across-fault ground motions and their effects on some bridges in the 1999 Chi-Chi earthquake

Simply-supported bridges are vulnerable to surface fault rupture as evidenced by several fault-crossing bridges in the 1999 Chi-Chi earthquake. To investigate the seismic collapse mechanism of simply-supported bridges crossing the fault, across-fault ground motions are firstly simulated in the present study. In particular, based on a previously developed fault model of the 1999 Chi-Chi earthquake, broadband across-fault ground motions at six fault-crossing bridges are simulated using the hybrid deterministic-stochastic method, in which the low- and high-frequency components are computed using the deterministic Green’s function method and the stochastic finite-fault modeling method, respectively. The simulation results indicate that the hybrid deterministic-stochastic method can give reasonable predictions to the across-fault ground motions. Furthermore, utilizing the explicit dynamic finite element (FE) code LS-DYNA, behaviors of a three-span simply-supported bridge under a selected pair of across-fault ground motions are numerically simulated. Numerical results indicate that the structural responses and collapse mechanisms are dominated by the low-frequency ground motions. The large differential static offset across the fault is the main reason for the collapse of the simply-supported bridges. This study contributes understandings for the across-fault ground motions and the collapse mechanism of some bridges in the 1999 Chi-Chi earthquake.

A New Ground Motion Model (GMM) for Northeast India (NEI) and Its Adjacent Countries for Interface Earthquakes Considering Both Strong Motion Records and Simulated Data

A ground motion model (GMM) for interface subduction zone earthquakes of Northeast India (NEI) and its adjacent countries is developed for the first time. Countries adjacent to NEI are Bangladesh, Bhutan, China, Myanmar and Nepal. High-magnitude earthquakes occur frequently in these regions due to buildup of high-stress parameters in the subduction zone of the Indian tectonic plate. Strong motion data are too few and sparse to develop a robust GMM for this region. We used both finite-fault simulations and a stochastic point-source model in developing our GMM. In our GMM, we used 50,000 ground motion samples which were stochastically simulated for different moment magnitudes (Mw) of 5.0–9.0 and hypocentral distances of 30–300 km using a point-source seismological stochastic model and finite fault model. In this study, we calculated stress drop (∆σ), quality factor Q(f) and all other region-specific seismic input parameters from the past strong motion records of interface subduction zone earthquakes of NEI and its adjacent countries. We used these seismic input parameters in ground motion simulation. Sensitivity analyses of the input parameters were also performed to check the bias of the present model. Our GMM was validated by comparing it with the existing NEI interface strong motion records. We compared our GMM with other GMMs developed for interface subduction zone earthquakes for different regions in the world. We also compared our GMM with point-source and finite-fault simulation models. Ground motion parameters estimated using the point-source model are comparatively higher than the finite-fault simulation model. Horizontal components of peak ground acceleration (PGA) and spectral acceleration (Sa) can be estimated for NEI using our GMM.

Historical earthquakes, tsunamis and real-time earthquake monitoring for tsunami advisory in the South China Sea region

The South China Sea Tsunami Advisory Center (SCSTAC) established by China, under the aegis of UNESCO's Intergovernmental Oceanographic Commission (UNESCO/IOC), had inaugurated commencing its full operation on November 5, 2019. This center is operating 24 × 7 h and round-the-clock shift to monitor tsunami hazard that pose a serious threat to countries which bordering the South China Sea (SCS) region. Prior to the official operation, SCSTAC had taken action for the last 10 years in upgrading their technology capability which are real-time earthquake monitoring and the processing system that is crucial to be able providing the international standard of tsunami warning services in the SCS and its adjacent areas. This paper briefly reviewed on the initiation steps and its developments of the South China Sea region Tsunami Warning and Mitigation Systems, tectonic setting as well as the characteristics of historical earthquakes and tsunamis in the region. In addition, we highlighted the structure and basic functions of the earthquake monitoring and processing system, earthquake location, source mechanism solution and finite fault model inversion using the real-time seismic waveform data from regional and global seismographic networks that will result in the rapid source parameters estimation for a larger earthquake in tsunami warning. Numerous simulations and hands-on events have shown that the preliminary earthquake parameters could be determined less than 8 min after earthquake. The W phase method is used and be able to produce rapid and reliable estimation of the moment magnitude and source mechanism for larger events within 10–15 min from earthquake origin time. A finite fault model can be acquired just after the earthquake event via computing teleseismic body-wave inversion program. The earthquake monitoring and processing system provide accurate and reliable information in contributing to tsunami warning services, thus promoting the development of tsunami warning technologies, which enhancing the tsunami warning capability and tsunami emergency responses. These high-end technology can be used in facilitating others such as marine disaster prevention, mitigation and its risk reduction.

Revisiting the seismogenic characteristics of stable continental interiors: The case of three Indian events

Kinematic source process models using tele-seismic waveform inversion of three stable continental region (SCR) earthquakes from India (1993 Mw 6.2 Latur, 1997 Mw 5.8 Jabalpur and 2001 Mw 7.6 Bhuj) are presented in this paper to demonstrate the slip evolution and stress drop. These finite-fault models are methodologically unique to constrain the source dimensions, compared with assumptions and quasi-observations around blind faults. Our results show that these events do have a compact zone of singular asperity breakage within the Indian crust. Whereas the Bhuj and Jabalpur events have their rupture majorly restricted within the lower crustal regions. The Latur event ruptured a very shallow crustal asperity. The estimated rupture velocities are in the range of 2.6–3.2 km/s, Bhuj event the slowest and lengthiest. Our results do not favour an updip shallow component of asperity breakage for the Bhuj earthquake, as evident from lack of surface ruptures. Study also indicates that the 1993 Latur event occurred on a fault with high stress concentration which indicates stronger fault asperities or slip within a newer fault. Models on 1997 Jabalpur event imply higher frictional strength and very brittle nature at the lower crustal regions of the Indian slab, indicating a mechanically very strong lower crust. We conclude that an event like Latur can occur anywhere within continental interiors of Indian SCR, or elsewhere globally, and is an underestimated seismic hazard.

Source mechanism and triggered large aftershocks of the Mw 6.5 Ambon, Indonesia earthquake

On September 26, 2019, a significant, Mw 6.5, earthquake shook Ambon region, Indonesia, causing severe damage on the Island(s) of Ambon. Due to the complexity of the fault network and in-situ stress field it was, up to now, not possible to define the fault plane using data from the BMKG regional seismic network. In this study, we analyze the fault plane of the 2019 Ambon earthquake and the reactivation potential of the surrounding faults using local networks. Eleven stations were deployed to monitor the aftershocks from October 18th to December 15th, 2019 augmented with data of four regional stations. During the monitoring period, 1,778 events were identified comprised of 10,938 P- and 10,315 and S- wave arrival times. The locations of aftershock were determined in a stepwise approach, i.e. (i) initial location determination using a non-linear approach, (ii) updating the velocity model, and (iii) relative double-difference relocation. Slip inversion using teleseismic data was performed to infer of high strain relief of the mainshock and to compute its associated static stress transfer (ΔCFF). Based on aftershock distribution and finite fault modeling, we conclude that the Mw 6.5 Ambon earthquake occurred on a N-S oriented fault plane. Two clusters consisting of ~60% of total events are located at both tips of the plane. Another cluster ~30% was sharply aligned in a NE-SW trend, 10 km westward, starting by an Mw 5.2 event on November 2nd, 2019. The b-value of the NE-SW events is ~0.25 lower than the other clusters with a b-value of 0.85±0.14. ΔCFF imparted by the mainshock caused ~0.5 Bar stress increase on the NE-SW fault. We concluded that the NE-SW trend was the reactivation of a preexisting fault crossing Ambon Island. The triggered large aftershock caused further significant damages to already weakened infrastructure and, thus, had the largest mapped damage area.

Benchmarking the Optimal Time Alignment of Tsunami Waveforms in Nonlinear Joint Inversions for the Mw 8.8 2010 Maule (Chile) Earthquake

Finite-fault models for the 2010Mw8.8 Maule, Chile earthquake indicate bilateral rupture with large-slip patches located north and south of the epicenter. Previous studies also show that this event features significant slip in the shallow part of the megathrust, which is revealed through correction of the forward tsunami modeling scheme used in tsunami inversions. The presence of shallow slip is consistent with the coseismic seafloor deformation measured off the Maule region adjacent to the trench and confirms that tsunami observations are particularly important for constraining far-offshore slip. Here, we benchmark the method of Optimal Time Alignment (OTA) of the tsunami waveforms in the joint inversion of tsunami (DART and tide-gauges) and geodetic (GPS, InSAR, land-leveling) observations for this event. We test the application of OTA to the tsunami Green’s functions used in a previous inversion. Through a suite of synthetic tests we show that if the bias in the forward model is comprised only of delays in the tsunami signals, the OTA can correct them precisely, independently of the sensors (DART or coastal tide-gauges) and, to the first-order, of the bathymetric model used. The same suite of experiments is repeated for the real case of the 2010 Maule earthquake where, despite the results of the synthetic tests, DARTs are shown to outperform tide-gauges. This gives an indication of the relative weights to be assigned when jointly inverting the two types of data. Moreover, we show that using OTA is preferable to subjectively correcting possible time mismatch of the tsunami waveforms. The results for the source model of the Maule earthquake show that using just the first-order modeling correction introduced by OTA confirms the bilateral rupture pattern around the epicenter, and, most importantly, shifts the inferred northern patch of slip to a shallower position consistent with the slip models obtained by applying more complex physics-based corrections to the tsunami waveforms. This is confirmed by a slip model refined by inverting geodetic and tsunami data complemented with a denser distribution of GPS data nearby the source area. The models obtained with the OTA method are finally benchmarked against the observed seafloor deformation off the Maule region. We find that all of the models using the OTA well predict this offshore coseismic deformation, thus overall, this benchmarking of the OTA method can be considered successful.

Operational Capabilities during Crisis: The Chilean Seismographic Network

AbstractThe severe mobility restrictions imposed countrywide by authorities of Chile in response to the pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have impacted all areas of activities. Major difficulties began in March 2020, with partial quarantines in the capital city, which later extended to other cities in the country, and it soon transformed in total confinement. We examine the evolution of the pandemic and its consequences on the field and headquarters operations of the National Seismological Centre (CSN), as remote stations could not be visited to carry out maintenance work. Several indicators, as a function of time, such as station operability percentage, timely reports, and ability to deliver requested information, reveal some negative impact on the uptime of stations but not in the capability of the CSN to fulfill its mission, which is to deliver timely seismic information to emergency services. The largest event in the country within this period took place on 3 June 2020, activating in a timely manner, the new tools on finite-fault modeling being developed within the center.

Experimental and Numerical Investigation of Old Masonry Wall Using a Macro-Modeling Approach

Background: A 3-D finite element model of the internal masonry wall of a 103-year-Old Senate hall, Allahabad University, has been modeled using macro-modeling approaches. The masonry wall is an excellent example of Indo-Saracenic style architecture used by Britishers during the late 19th Century, which is a unification of the Mughal and Colonial architecture. Methods: Non-destructive Test (NDT) has been conducted to estimate is compressive strength and Young’s modulus of the wall. Compressive strength of the brick masonry and stone arch was estimated in the range of 10.5-12.5 MPa and 18.6-21.2 MPa, respectively, whereas Young’s Modulus was estimated in the range of 1800-5000 MPa and 5500-8000 MPa (outlier not considered). Finite Element model was prepared using the macro-modeling approach. Results: The gravity load analysis shows that the wall is stable, and its geometrical configuration is safe with maximum Von-Mises stress of 5.38 MPa and deformation of 2.27 mm. The results of the first six modes are presented. Further, in the absence of a recorded ground motion for the Prayag city, synthetic ground motion is simulated for 25th April 2015 Nepal earthquake (Mw) using a stochastic finite fault model. Conclusion: Evaluated behaviour of the internal masonry wall is shown in the form of acceleration, deformation and stress response.