scholarly journals 3D Linked Subduction, Dynamic Rupture, Tsunami, and Inundation Modeling: Dynamic Effects of Supershear and Tsunami Earthquakes, Hypocenter Location, and Shallow Fault Slip

2021 ◽  
Vol 9 ◽  
Author(s):  
Sara Aniko Wirp ◽  
Alice-Agnes Gabriel ◽  
Maximilian Schmeller ◽  
Elizabeth H. Madden ◽  
Iris van Zelst ◽  
...  

Physics-based dynamic rupture models capture the variability of earthquake slip in space and time and can account for the structural complexity inherent to subduction zones. Here we link tsunami generation, propagation, and coastal inundation with 3D earthquake dynamic rupture (DR) models initialized using a 2D seismo-thermo-mechanical geodynamic (SC) model simulating both subduction dynamics and seismic cycles. We analyze a total of 15 subduction-initialized 3D dynamic rupture-tsunami scenarios in which the tsunami source arises from the time-dependent co-seismic seafloor displacements with flat bathymetry and inundation on a linearly sloping beach. We first vary the location of the hypocenter to generate 12 distinct unilateral and bilateral propagating earthquake scenarios. Large-scale fault topography leads to localized up- or downdip propagating supershear rupture depending on hypocentral depth. Albeit dynamic earthquakes differ (rupture speed, peak slip-rate, fault slip, bimaterial effects), the effects of hypocentral depth (25–40 km) on tsunami dynamics are negligible. Lateral hypocenter variations lead to small effects such as delayed wave arrival of up to 100 s and differences in tsunami amplitude of up to 0.4 m at the coast. We next analyse inundation on a coastline with complex topo-bathymetry which increases tsunami wave amplitudes up to ≈1.5 m compared to a linearly sloping beach. Motivated by structural heterogeneity in subduction zones, we analyse a scenario with increased Poisson's ratio of ν = 0.3 which results in close to double the amount of shallow fault slip, ≈1.5 m higher vertical seafloor displacement, and a difference of up to ≈1.5 m in coastal tsunami amplitudes. Lastly, we model a dynamic rupture “tsunami earthquake” with low rupture velocity and low peak slip rates but twice as high tsunami potential energy. We triple fracture energy which again doubles the amount of shallow fault slip, but also causes a 2 m higher vertical seafloor uplift and the highest coastal tsunami amplitude (≈7.5 m) and inundation area compared to all other scenarios. Our mechanically consistent analysis for a generic megathrust setting can provide building blocks toward using physics-based dynamic rupture modeling in Probabilistic Tsunami Hazard Analysis.

2021 ◽  
Author(s):  
Sara Aniko Wirp ◽  
Alice-Agnes Gabriel ◽  
Elizabeth H. Madden ◽  
Maximilian Schmeller ◽  
Iris van Zelst ◽  
...  

<p>Earthquake rupture dynamic models capture the variability of slip in space and time while accounting for the structural complexity which is characteristic for subduction zones. The use of a geodynamic subduction and seismic cycling (SC) model to initialize dynamic rupture (DR) ensures that initial conditions are self-consistent and reflect long-term behavior. We extend the 2D geodynamical subduction and SC model of van Zelst et al. (2019) and use it as input for realistic 3-dimensional DR megathrust earthquake models. We follow the subduction to tsunami run-up linking approach described in Madden et al. (2020), including a complex subduction setup along with their resulting tsunamis. The distinct variation of shear traction and friction coefficients with depth lead to realistic average rupture speeds and dynamic stress drop as well as efficient tsunami generation. </p><p>We here analyze a total of 14 subduction-initialized 3D dynamic rupture-tsunami scenarios. By varying the hypocentral location along arc and depth, we generate 12 distinct unilateral and bilateral earthquakes with depth-variable slip distribution and directivity, bimaterial, and geometrical effects in the dynamic slip evolutions. While depth variations of the hypocenters barely influence the tsunami behavior, lateral varying nucleation locations lead to a shift in the on-fault slip which causes time delays of the wave arrival at the coast. The fault geometry of our DR model that arises during the SC model is non-planar and includes large-scale roughness. These features (topographic highs) trigger supershear rupture propagation in up-dip or down-dip direction, depending on the hypocentral depth.</p><p>In two additional scenarios, we analyze variations in the energy budget of the DR scenarios. In the SC model, an incompressible medium is assumed (ν=0.5) which is valid only for small changes in pressure and temperature. Unlike in the DR model where the material is compressible and a different Poisson’s ratio (ν=0.25) has to be assigned. Poisson’s ratios between 0.1 and 0.4 stand for various compressible materials. To achieve a lower shear strength of all material on and off the megathrust fault and to facilitate slip, we increase the Poisson ratio in the DR model to ν=0.3 which is consistent with basaltic rocks. As a result, larger fault slip is concentrated at shallower depths and generates higher vertical seafloor displacement and earthquake moment magnitude respectively. Even though the tsunami amplitudes are much higher, the same dynamic behavior as in the twelve hypocenter-variable models can be observed. Lastly, we increase fracture energy by changing the critical slip distance in the linear slip-weakening frictional parameterization. This generates a tsunami earthquake (Kanamori, 1972) characterized by low rupture velocity (on average half the amount of s-wave speed) and low peak slip rate, but at the same time large shallow fault slip and moment magnitude. The shallow fault slip of this event causes the highest vertical seafloor uplift compared to all other simulations. This leads to the highest tsunami amplitude and inundation area while the wavefront hits the coast delayed compared to the other scenarios.</p>


2021 ◽  
Vol 22 (11) ◽  
pp. 5793
Author(s):  
Brianna M. Quinville ◽  
Natalie M. Deschenes ◽  
Alex E. Ryckman ◽  
Jagdeep S. Walia

Sphingolipids are a specialized group of lipids essential to the composition of the plasma membrane of many cell types; however, they are primarily localized within the nervous system. The amphipathic properties of sphingolipids enable their participation in a variety of intricate metabolic pathways. Sphingoid bases are the building blocks for all sphingolipid derivatives, comprising a complex class of lipids. The biosynthesis and catabolism of these lipids play an integral role in small- and large-scale body functions, including participation in membrane domains and signalling; cell proliferation, death, migration, and invasiveness; inflammation; and central nervous system development. Recently, sphingolipids have become the focus of several fields of research in the medical and biological sciences, as these bioactive lipids have been identified as potent signalling and messenger molecules. Sphingolipids are now being exploited as therapeutic targets for several pathologies. Here we present a comprehensive review of the structure and metabolism of sphingolipids and their many functional roles within the cell. In addition, we highlight the role of sphingolipids in several pathologies, including inflammatory disease, cystic fibrosis, cancer, Alzheimer’s and Parkinson’s disease, and lysosomal storage disorders.


2021 ◽  
Author(s):  
Kor de Jong ◽  
Marc van Kreveld ◽  
Debabrata Panja ◽  
Oliver Schmitz ◽  
Derek Karssenberg

<p>Data availability at global scale is increasing exponentially. Although considerable challenges remain regarding the identification of model structure and parameters of continental scale hydrological models, we will soon reach the situation that global scale models could be defined at very high resolutions close to 100 m or less. One of the key challenges is how to make simulations of these ultra-high resolution models tractable ([1]).</p><p>Our research contributes by the development of a model building framework that is specifically designed to distribute calculations over multiple cluster nodes. This framework enables domain experts like hydrologists to develop their own large scale models, using a scripting language like Python, without the need to acquire the skills to develop low-level computer code for parallel and distributed computing.</p><p>We present the design and implementation of this software framework and illustrate its use with a prototype 100 m, 1 h continental scale hydrological model. Our modelling framework ensures that any model built with it is parallelized. This is made possible by providing the model builder with a set of building blocks of models, which are coded in such a manner that parallelization of calculations occurs within and across these building blocks, for any combination of building blocks. There is thus full flexibility on the side of the modeller, without losing performance.</p><p>This breakthrough is made possible by applying a novel approach to the implementation of the model building framework, called asynchronous many-tasks, provided by the HPX C++ software library ([3]). The code in the model building framework expresses spatial operations as large collections of interdependent tasks that can be executed efficiently on individual laptops as well as computer clusters ([2]). Our framework currently includes the most essential operations for building large scale hydrological models, including those for simulating transport of material through a flow direction network. By combining these operations, we rebuilt an existing 100 m, 1 h resolution model, thus far used for simulations of small catchments, requiring limited coding as we only had to replace the computational back end of the existing model. Runs at continental scale on a computer cluster show acceptable strong and weak scaling providing a strong indication that global simulations at this resolution will soon be possible, technically speaking.</p><p>Future work will focus on extending the set of modelling operations and adding scalable I/O, after which existing models that are currently limited in their ability to use the computational resources available to them can be ported to this new environment.</p><p>More information about our modelling framework is at https://lue.computationalgeography.org.</p><p><strong>References</strong></p><p>[1] M. Bierkens. Global hydrology 2015: State, trends, and directions. Water Resources Research, 51(7):4923–4947, 2015.<br>[2] K. de Jong, et al. An environmental modelling framework based on asynchronous many-tasks: scalability and usability. Submitted.<br>[3] H. Kaiser, et al. HPX - The C++ standard library for parallelism and concurrency. Journal of Open Source Software, 5(53):2352, 2020.</p>


2021 ◽  
Author(s):  
Jagdish Chandra Vyas ◽  
Martin Galis ◽  
Paul Martin Mai

<p>Geological observations show variations in fault-surface topography not only at large scale (segmentation) but also at small scale (roughness). These geometrical complexities strongly affect the stress distribution and frictional strength of the fault, and therefore control the earthquake rupture process and resulting ground-shaking. Previous studies examined fault-segmentation effects on ground-shaking, but our understanding of fault-roughness effects on seismic wavefield radiation and earthquake ground-motion is still limited.  </p><p>In this study we examine the effects of fault roughness on ground-shaking variability as a function of distance based on 3D dynamic rupture simulations. We consider linear slip-weakening friction, variations of fault-roughness parametrizations, and alternative nucleation positions (unilateral and bilateral ruptures). We use generalized finite difference method to compute synthetic waveforms (max. resolved frequency 5.75 Hz) at numerous surface sites  to carry out statistical analysis.  </p><p>Our simulations reveal that ground-motion variability from unilateral ruptures is almost independent of  distance from the fault, with comparable or higher values than estimates from ground-motion prediction equations (e.g., Boore and Atkinson, 2008; Campbell and Bozornia, 2008). However, ground-motion variability from bilateral ruptures decreases with increasing distance, in contrast to previous studies (e.g., Imtiaz et. al., 2015) who observe an increasing trend with distance. Ground-shaking variability from unilateral ruptures is higher than for bilateral ruptures, a feature due to intricate seismic radiation patterns related to fault roughness and hypocenter location. Moreover, ground-shaking variability for rougher faults is lower than for smoother faults. As fault roughness increases the difference in ground-shaking variabilities between unilateral and bilateral ruptures increases. In summary, our simulations help develop a fundamental understanding of ground-motion variability at high frequencies (~ 6 Hz) due small-scale geometrical fault-surface variations.</p>


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Steve Noone ◽  
Alison Branch ◽  
Melissa Sherring

Purpose Positive behavioural support (PBS) as a framework for delivering quality services is recognised in important policy documents (CQC, 2020; NICE, 2018), yet there is an absence in the literature on how this could be implemented on a large scale. The purpose of this paper is to describe a recent implementation of a workforce strategy to develop PBS across social care and health staff and family carers, within the footprint of a large integrated care system. Design/methodology/approach A logic model describes how an initial scoping exercise led to the production of a regional workforce strategy based on the PBS Competence Framework (2015). It shows how the creation of a regional steering group was able to coordinate important developmental stages and integrate multiple agencies into a single strategy to implement teaching and education in PBS. It describes the number of people who received teaching and education in PBS and the regional impact of the project in promoting cultural change within services. Findings This paper demonstrates a proof of concept that it is possible to translate the PBS Competency Framework (2015) into accredited courses. Initial scoping work highlighted the ineffectiveness of traditional training in PBS. Using blended learning and competency-based supervision and assessment, it was possible to create a new way to promote large-scale service developments in PBS supported by the governance of a new organisational structure. This also included family training delivered by family trainers. This builds on the ideas by Denne et al. (2020) that many of the necessary building blocks of implementation already exist within a system. Social implications A co-ordinated teaching and education strategy in PBS may help a wide range of carers to become more effective in supporting the people they care for. Originality/value This is the first attempt to describe the implementation of a framework for PBS within a defined geographical location. It describes the collaboration of health and social care planners and a local university to create a suite of courses built around the PBS coalition competency framework.


Geosphere ◽  
2020 ◽  
Author(s):  
Katherine A. Guns ◽  
Richard A Bennett ◽  
Joshua C. Spinler ◽  
Sally F. McGill

Assessing fault-slip rates in diffuse plate boundary systems such as the San Andreas fault in southern California is critical both to characterize seis­mic hazards and to understand how different fault strands work together to accommodate plate boundary motion. In places such as San Gorgonio Pass, the geometric complexity of numerous fault strands interacting in a small area adds an extra obstacle to understanding the rupture potential and behavior of each individual fault. To better understand partitioning of fault-slip rates in this region, we build a new set of elastic fault-block models that test 16 different model fault geometries for the area. These models build on previ­ous studies by incorporating updated campaign GPS measurements from the San Bernardino Mountains and Eastern Transverse Ranges into a newly calculated GPS velocity field that has been removed of long- and short-term postseismic displacements from 12 past large-magnitude earthquakes to estimate model fault-slip rates. Using this postseismic-reduced GPS velocity field produces a best- fitting model geometry that resolves the long-standing geologic-geodetic slip-rate discrepancy in the Eastern California shear zone when off-fault deformation is taken into account, yielding a summed slip rate of 7.2 ± 2.8 mm/yr. Our models indicate that two active strands of the San Andreas system in San Gorgonio Pass are needed to produce sufficiently low geodetic dextral slip rates to match geologic observations. Lastly, results suggest that postseismic deformation may have more of a role to play in affecting the loading of faults in southern California than previously thought.


2020 ◽  
Author(s):  
Valère Lambert ◽  
Nadia Lapusta

Abstract. Substantial insight into earthquake source processes has resulted from considering frictional ruptures analogous to cohesive-zone shear cracks from fracture mechanics. This analogy holds for slip-weakening representations of fault friction that encapsulate the resistance to rupture propagation in the form of breakdown energy, analogous to fracture energy, prescribed in advance as if it were a material property of the fault interface. Here, we use numerical models of earthquake sequences with enhanced weakening due to thermal pressurization of pore fluids to show how accounting for thermo-hydro-mechanical processes during dynamic shear ruptures makes breakdown energy rupture-dependent. We find that local breakdown energy is neither a constant material property nor uniquely defined by the amount of slip attained during rupture, but depends on how that slip is achieved through the history of slip rate and dynamic stress changes during the rupture process. As a consequence, the frictional breakdown energy of the same location along the fault can vary significantly in different earthquake ruptures that pass through. These results suggest the need for re-examining the assumption of pre-determined frictional breakdown energy common in dynamic rupture modeling and for better understanding of the factors that control rupture dynamics in the presence of thermo-hydro-mechanical processes.


2019 ◽  
Vol 9 (2) ◽  
pp. 152
Author(s):  
Rahmat Setyo Yuliatmoko ◽  
Telly Kurniawan

The amount of stress released by an earthquake can be calculated with a stress drop, the stress ratio before and after an earthquake where the stress accumulated in a fault or a subduction zone is immediately released during an earthquake. The purpose of this research is to calculate the amount of stress drop in faults and subduction in Maluku and Halmahera and their variations and relate them to the geological conditions in the area so that the tectonic characteristics in the area can be identified. This research employed mathematical analysis and the Nelder Mead Simplex nonlinear inversion methods. The results show that Maluku and Halmahera are the area with complex tectonic conditions and large earthquake impacts. The Maluku sea earthquake generated a stress drop of 0.81 MPa with a reverse fault mechanism in the zone of subduction, while for the Halmahera earthquake the stress drop value was 52.72 MPa, a typical strike-slip mechanism in the fault zone. It can be concluded that there is a difference in the stress drop between the subduction and fault zones; the stress drop in the fault was greater than that in the subduction zone due to different rock structure and faulting mechanisms as well as differences in the move slip rate that plays a role in the process of holding out the stress on a rock. This information is very important to know the amount of pressure released from the earthquake which has a very large impact as part of disaster mitigation measures.


Sign in / Sign up

Export Citation Format

Share Document