NIOM Analysis of Earthquake Observation Records in a Rockfill Dam

ABSTRACT We applied the normalized input–output minimization method (a method developed for the analysis of propagation times in vertical array records) to long-term earthquake observation records from Aratozawa Dam (in Kurihara, Miyagi prefecture, Japan), spanning the period from July 1992 to December 2019 to determine the propagation velocity of seismic waves in the embankment, and investigated changes in soil properties. As a result, we showed that (1) the velocities of S and P waves in the upper section were 449 and 993 m/s, respectively, prior to the strong earthquake motions derived from earthquake records from January 1997 through October 2001, whereas 608 and 1538, respectively, in the lower section, (2) in the Iwate–Miyagi Nairiku earthquake, the S-wave velocity in the upper section decreased to 158 m/s in the principal shock, and (3) in subsequent minor earthquakes the propagation velocity increased more or less in proportion with the logarithm of the number of elapsed days, requiring three years or longer to return to the initial value, (4) although similar changes were observed in the Great East Japan earthquake of 2011, the reduction in propagation velocity that remained after the principal shock was smaller than in the case of the Iwate–Miyagi Nairiku earthquake, and it was judged that there were no large effects on the dam body such as those that occurred in the Iwate–Miyagi Nairiku earthquake, and furthermore (5) in the principal shock of the Iwate–Miyagi Nairiku earthquake, the shear modulus in the upper part of the dam body decreased from 400 to 50 MPa (with a maximum shear strain of 10−3), resulting in more pronounced changes than in the lower section, whereas the damping ratio increased by at least 10% in the lower section during the principal shock of the Iwate–Miyagi Nairiku earthquake, resulting in much greater changes than in the upper section.

Estimation of crustal deformation parameters and strain build-up in Northwest Himalaya using GNSS data measurements

GPS measurements have proved extremely useful in quantifying strain accumulation rate and assessing seismic hazard in a region. Continuous GPS measurements provide estimates of secular motion used to understand the earthquake and other geodynamic processes. GNSS stations extending from the South of India to the Higher Himalayan region have been used to quantify the strain build-up rate in Central India and the Himalayan region to assess the seismic hazard potential in this realm. Velocity solution has been determined after the application of Markov noise estimated from GPS time series data. The recorded GPS data are processed along with the closest International GNSS stations data for estimation of daily basis precise positioning. The baseline method has been used for the estimation of the linear strain rate between the two stations. Whereas the principal strain axes, maximum shear strain, rotation rate, and crustal shortening rate has been calculated through the site velocity using an independent approach; least-square inversion approach-based triangulation method. The strain rate analysis estimated by the triangulation approach exhibits a mean value of extension rate of 26.08 nano-strain/yr towards N131°, the compression rate of –25.38 nano-strain/yr towards N41°, maximum shear strain rate of 51.47 nano-strain/yr, dilation of –37.57 nano-strain/yr and rotation rate of 0.7°/Ma towards anti-clockwise. The computed strain rate from the Baseline method and the Triangulation method reports an extensive compression rate that gradually increases from the Indo-Gangetic Plain in South to Higher Himalaya in North. The slip deficit rate between India and Eurasia Plate in Kumaun Garhwal Himalaya has been computed as 18±1.5 mm/yr based on elastic dislocation theory. Thus, in this study, present-day surface deformation rate and interseismic strain accumulation rate in the Himalayan region and the Central Indian region have been estimated for seismic hazard analysis using continuous GPS measurements.

Shear induced deformation twinning evolution in thermoelectric InSb

Twin boundary (TB) engineering has been widely applied to enhance the strength and plasticity of metals and alloys, but is rarely adopted in thermoelectric (TE) semiconductors. Our previous first-principles results showed that nanotwins can strengthen TE Indium Antimony (InSb) through In–Sb covalent bond rearrangement at the TBs. Herein, we further show that shear-induced deformation twinning enhances plasticity of InSb. We demonstrate this by employing large-scale molecular dynamics (MD) to follow the shear stress response of flawless single-crystal InSb along various slip systems. We observed that the maximum shear strain for the $$(111)[11\bar 2]$$ ( 111 ) [ 11 2 ¯ ] slip system can be up to 0.85 due to shear-induced deformation twinning. We attribute this deformation twinning to the “catching bond” involving breaking and re-formation of In–Sb bond in InSb. This finding opens up a strategy to increase the plasticity of TE InSb by deformation twinning, which is expected to be implemented in other isotypic III–V semiconductors with zinc blende structure.

Group characterization of impact-induced, in vivo human brain kinematics

Brain movement during an impact can elicit a traumatic brain injury, but tissue kinematics vary from person to person and knowledge regarding this variability is limited. This study examines spatio-temporal brain–skull displacement and brain tissue deformation across groups of subjects during a mild impact in vivo . The heads of two groups of participants were imaged while subjected to a mild (less than 350 rad s −2 ) impact during neck extension (NE, n = 10) and neck rotation (NR, n = 9). A kinematic atlas of displacement and strain fields averaged across all participants was constructed and compared against individual participant data. The atlas-derived mean displacement magnitude was 0.26 ± 0.13 mm for NE and 0.40 ± 0.26 mm for NR, which is comparable to the displacement magnitudes from individual participants. The strain tensor from the atlas displacement field exhibited maximum shear strain (MSS) of 0.011 ± 0.006 for NE and 0.017 ± 0.009 for NR and was lower than the individual MSS averaged across participants. The atlas illustrates common patterns, containing some blurring but visible relationships between anatomy and kinematics. Conversely, the direction of the impact, brain size, and fluid motion appear to underlie kinematic variability. These findings demonstrate the biomechanical roles of key anatomical features and illustrate common features of brain response for model evaluation.

Kinematics of Crustal Deformation along the Central Himalaya

Using an updated set of GPS surface velocities, the present study provides fault locking behavior and slip rate distribution of the Main Himalayan Thrust (MHT) along the central Himalaya. The two-dimensional velocity field is inverted through Bayesian inversion to estimate fault geometry and kinematic parameters of the MHT along the central Himalaya. The modeling results reveal that: (1) MHT is fully locked in the upper flat (0-9 km), partially locked along the mid-crustal ramp (15-21 km), and it is creeping in the deeper flat (> 21 km); (2) there is an insignificant slip rate of MHT along the locked-to-creeping transition zone, indicating its partially coupled/locked behavior; (3) along the deeper flat of the MHT, the estimated creeping rate is ~16.3 mm/yr, ~14.7 mm/yr, and ~14.3 mm/yr along western, central, and eastern Nepal, respectively; and (4) along the MHT on the upper crust, the modeled locking width turns out to be 97 km, 106 km, and 129 km in the western, central, and eastern Nepal, respectively. In addition, the posterior probability distribution of the locking width exhibits a bimodal Gaussian distribution coinciding with the two ramp geometry of the MHT along the western Nepal. Along the foothills of the Higher Himalaya, the inferred locking line is also aligned to the estimated maximum shear strain concentration and observed seismicity along the central Himalaya. With a general agreement to the previous geodetic results, geological estimates, and background seismicity, our findings provide a promising avenue of the contemporary crustal deformation along the Nepal Himalaya. The estimated inversion results in a Bayesian framework exhibit updated fault kinematics of the MHT and hence provides valuable inputs for seismic hazard assessment along the central Himalaya.

Analytical Method for Preliminary Seismic Design of Tunnels

Buried structures are categorized based on their shape, size and location. These main categories are near surface structures (e.g., pipes and other facilities), large section structures (e.g., tunnels, subways, etc.), and vertical underground structures (e.g., shafts and ducts). Seismic assessments of these structures are important in areas close to severe seismic sources. Seismic design of tunnels requires calculation of the deformation in surrounding geological formations. The seismic hazard on a site is usually expressed as a function of amplitude parameters of free-field motion. Therefore, simplified relations between depth and parameters of ground motion are necessary for preliminary designs. The objective of this chapter is to study and review the main analytical seismic methods which are used to develop a simple relationship between maximum shear strain, maximum shear stress and other seismic parameters.

Air leaks: stapling affects porcine lungs biomechanics

During thoracic operations, surgical staplers resect cancerous tumors and seal the spared lung. However, post-operative air leaks are undesirable clinical consequences: staple legs wound lung tissue. Subsequent to this trauma, air leaks from lung tissue into the pleural space. This affects the lung's physiology and patients' recovery. The objective is to biomechanically and visually characterize porcine lung tissue with and without staples in order to gain knowledge on air leakage following pulmonary resection. Therefore, a syringe pump filled with air inflates and deflates eleven porcine lungs cyclically without exceeding 10 cmH2O of pressure. Cameras capture stereo-images of the deformed lung surface at regular intervals while a microcontroller simultaneously records the alveolar pressure and the volume of air pumped. The raw images are then used to compute tri-dimensional displacements and strains with the Digital Image Correlation method (DIC). Air bubbles originated at staple holes of inner row from exposed porcine lung tissue due to torn pleural on costal surface. Compared during inflation, left upper or lower lobe resections have similar compliance (slope of the pressure vs volume curve), which are 9% lower than healthy lung compliance. However, lower lobes statistically burst at lower pressures than upper lobes (p-value<0.046) in ex vivo conditions confirming previous clinical in vivo studies. In parallel, the lung deformed mostly in the vicinity of staple holes and presented maximum shear strain near the observed leak location. To conclude, a novel technique DIC provided concrete evidence of the post-operative air leaks biomechanics. Further studies could investigate causal relationships between the mechanical parameters and the development of an air leak.

Experimental Investigation on Shear Failure Mechanism of Rock Mass with Intermittent Joints

There are a large number of discontinuous weak planes distributed in the natural rock mass, which makes the sliding failure of rock mass along the intermittent structural plane very complex. To investigate the shear failure mechanism of rock mass with intermittent joints and study the influence of different joint heights on the shear failure mode of the rock mass, direct shear tests were carried out by presetting a series of jointed rock specimens with different undulating heights. During the shear loading, digital image correlation (DIC) technology was employed to monitor the surface strain field of the specimens in real time. The results show that the fluctuation height has a significant effect on the evolution of shear strain. With the increase of shear load, the maximum shear strain of the jointed specimens with different undulating heights first increases slowly and then increases rapidly. When the undulating height is 5 mm, the failure of the specimen is dominated by the rock sliding along prefabricated joints. When the undulating height is larger than 10 mm, the shear fracture of the rock becomes dominant. With the increase of the undulating height, more penetrating cracks perpendicular to the preexisting joints appear between the serrated surfaces, and the shear fracture phenomenon is more obvious.

Crustal strain-rate fields estimated from GNSS data with a Bayesian approach and its correlation to seismic activity

&lt;p&gt;We proposed a new Bayesian approach to estimate continuous crustal strain-rate fields from spatially discrete displacement-rate data, based on Global Navigation Satellite System (GNSS) observations, under the prior constraint on spatial flatness of the strain-rate fields. The optimal values of the hyperparameters in the model of strain-rate fields are determined by using Akaike's Bayesian Information Criterion. A methodological merit of this approach is that, by introducing a two-layer Delaunay tessellation technique, the time-consuming computation of strain rates can be omitted through the model estimation process. We applied the Bayesian approach to GNSS displacement-rate data in Mainland China and examined the correlation between the estimated strain-rate fields and seismic activity by using Molchan&amp;#8217;s Error Diagram. The results show that the increase rate of maximum shear strain is positively correlated with the occurrence of earthquakes, indicating the strain rate can be used to augment probability earthquake models for background seismicity forecasting.&lt;/p&gt;

Stress-independent parameters for stress-strain relationship and damping in nonlinear one-dimensional seismic site response analysis

The modulus reduction and damping curves represent the nonlinear behavior of soil under cyclic load. In the literature, those curves were produced from lab tests of soil at particular confining stresses. This study developed a set of parameters that can be used to normalize the modulus reduction and damping curves to be stress-independent. The proposed formulations for the stress-independent parameters were implemented in the finite element code SRAP and validated through producing shear modulus reduction and damping curves that match the existed ones. Nonlinear 1D seismic site response analyses were conducted for centrifuge experiments to verify the developed computer code. Comparisons of the analysis results between SRAP and another computer code were presented in terms of maximum and minimum displacement, peak ground acceleration, maximum shear strain profiles, and response spectra. Keywords: backbone curve; hysteretic damping; dynamic soil model; stress-independent parameters; finite element method; nonlinear 1D seismic site response analysis.