Search

Interferometric synthetic aperture radar (InSAR) decorrelation that creates great challenges to phase unwrapping has been a critical issue for mapping large earthquake deformation. Some studies have proposed a 'remove-and-return model' solution to tackle this problem, but it has not been fully validated yet, and therefore has rarely been applied to real earthquake scenarios. In this study, we use the 2023 M w 7.8 and 7.6 earthquake doublet in Turkey and Syria as a case example to develop an iterative modelling method for InSAR-based coseismic mapping. We first derive surface deformation fields using Sentinel-1 offset tracking and Sentinel-2 optical image correlation, and invert them for an initial coseismic slip model, based on which we simulate InSAR coseismic phase measurements. We then remove the simulated phase from the actual Sentinel-1 phase and conduct unwrapping. The simulated phase is added back to the unwrapped phase to produce the final phase measurements. Comparing to the commonly used unwrapping method, our proposed approach can significantly improve coherence and reduce phase gradients, enabling accurate InSAR measurements. Combining InSAR, offset tracking and optical image correlation, we implement a joint inversion to obtain an optimal coseismic slip model. Our model shows that slip on the Çardak Fault is concentrated on a ∼100 km segment; to both ends, slip suddenly diminished. On the contrary, rupture on the East Anatolian Fault Zone propagated much longer as its geometry is fairly smooth. The iterative coseismic modelling method is proven efficient and can be easily applied to other continental earthquakes. [ABSTRACT FROM AUTHOR]

Earthquake source parameters can be estimated using seismological observations, but the identification of the fault responsible is often complicated by location uncertainties and the inherent ambiguity between nodal planes. Satellite Interferometric Synthetic Aperture Radar (InSAR) can be used to observe ground deformation and model fault geometry but is limited by climate conditions (water vapour) and ground coverage (dense vegetation). In the tropics, the atmosphere is dynamic and most regions are densely vegetated, making detecting coseismic deformation challenging. Here, we perform a systematic inspection of coseismic interferograms from Sentinel-1 SAR images, to assess their suitability for detecting coseismic deformation in Costa Rica. Using data from the seismological network, we target seven earthquakes between 2016 and 2020 with depths |$\le \, 20$| km and magnitudes M w 5.3–6.2. For each event, we use the seismic parameters to compute line-of-sight displacements for ascending and descending geometries and for both nodal planes and generate 12- and 24-d coseismic interferograms where available. We obtain interferograms with coseismic displacement signals for three of the seven earthquakes. We invert the geodetic data to retrieve the earthquake source parameters but the lack of offshore geodetic coverage causes trade-offs between parameters and large uncertainties. The Jacó and Golfito earthquakes likely occurred on the subduction interface and the geodetic locations were 6–9 km closer to the coast than previous seismic estimates. The Burica earthquake occurred on a shallow steeply dipping thrust fault in the outer forearc. For the other earthquakes, no coseismic deformation was detected due to atmospheric noise or poor coherence. These results demonstrate the suitability of 12-d Sentinel-1 interferograms for monitoring shallow earthquakes of magnitude > M w 5.7 in Central America. This approach can be used to begin a surface deformation catalogue for the region, which will ultimately help improve the understanding of active deformation processes and improve hazard maps. [ABSTRACT FROM AUTHOR]

Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson's ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth's structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M 9 event, or distributed deformation along plate boundaries as constrained from InSAR. [ABSTRACT FROM AUTHOR]

High quality Interferometric Synthetic Aperture Radar (InSAR) interferograms are essential for determining surface deformation from InSAR time-series. However, InSAR interferograms are usually polluted by spatially correlated errors (SCEs), especially the unmodelled atmospheric phase delays. To mitigate spatially correlated errors and improve the quality of InSAR interferograms, we propose a new approach to incorporate the Global Navigation Satellite System (GNSS) data from continuously operating reference stations for enhancing InSAR interferograms via modelling SCEs as signals and solving the signals together with the systematic parameters using least squares collocation (LSC), which is referred to as the LSC-GInSAR approach. Our improvement for the GInSAR method of Neely et al. can correct more SCEs. The Sentinel-1 data over the southern Central Valley of California, USA, are processed with our LSC-GInSAR approach, which is compared to the GInSAR approach. The performance of the LSC-GInSAR approach is evaluated by external GNSS displacements. The results show that the LSC-GInSAR approach can effectively mitigate medium-to-long-wavelength SCEs. The displacements resolved by LSC-GInSAR are more consistent with the cGNSS observations than those resolved by GInSAR, with an average root mean square improvement of 14.3 per cent. In addition, the LSC-GInSAR approach reduced the average standard deviations of all 276 InSAR interferograms from 14.2 to 11.0 mm compared to that of the GInSAR approach. [ABSTRACT FROM AUTHOR]

In the present study, the wear behavior of the Fe–14Cr–1Ti–1W–0.5L2O3 oxide dispersion strengthened (ODS) steels has been reported. The ball-on-disc wear tester with silicon carbide ball as counter surface was used to perform the wear tests. The effect of normal load (10–40 N with the interval of 10 N) on the wear behavior of the alloy has been investigated. Furthermore, to understand the operative wear mechanism of the alloy, detailed microstructural and compositional analyses have also been carried out using scanning electron microscopy. The highest width of the wear track was obtained as 896 ± 78 µm at 40 N normal load. The wear rate of the alloy was found to increase from 4.54  ±  0.12 × 10 - 6 mm3/m to 25.48  ±  1.15 × 10 - 6 mm3/m with the normal load from 10 to 40 N at constant sliding velocity and distance. Moreover, the microstructure analyses of the worn surfaces, wear debris and counter surface have shown the prominent wear mechanism to be a combination of adhesion, abrasion, oxidation and plastic deformation. The addition of lanthana was found to reduce the wear rate of ferritic steels. [ABSTRACT FROM AUTHOR]

Atmospheric rivers (ARs) are efficient mechanisms for transporting atmospheric moisture from low latitudes to the Antarctic continent. AR events induce intense snowfall episodes, which increase crustal deformation. Here, we used an AR detection algorithm, via a spatial matrix operation to quantify the contribution of AR-induced snowfall in the Antarctic continent, 2010–2019. Our results reveal that the AR snowfall contribution to Antarctica primarily ranged from 9.28 to 29.73 per cent from 2010 to 2019, and there was an evident increasing trend from 2015 to 2019 (20.66 per cent in 2015 to 29.30 per cent in 2019). AR-induced snowfall is one of the factors influencing the surface deformation of the Antarctic continent, based on the hourly AR snowfall deformation calculations for the Antarctic continent, both the average and maximum crustal deformation or displacement tends to be greatest near the coastline, while the displacement is less affected by AR further inland. [ABSTRACT FROM AUTHOR]

SUMMARY: Large continental earthquakes are infrequent, but they are important for understanding active tectonics. Many pre-1990, magnitude greater than 7 earthquakes have not been well documented due to a lack of data. The declassification of KeyHole (KH) imagery provides opportunities for investigating those historical events. In this study, we present new geodetic observations of surface deformation from the 1976 November 24 Ms 7.3 Chaldiran earthquake in eastern Turkey. We use various historical KH imagery and modern Pleiades stereo imagery to determine the coseismic deformation, providing new constraints on the 1976 rupture. Based on the image measurements, we resolve the coseismic slip at depth. The inverted slip from a simple elastic dislocation model, ∼4.0 m at shallow depth of the Hidirmentes segment, is larger than the surface displacement, around 3 m from previous field surveys, suggesting that the coseismic deformation is a combination of slip on the primary fault with 30 ± 16 per cent off-fault deformation. We also measure cumulative offsets at two sites with well-preserved features, and find a constant ratio of horizontal to vertical offsets throughout the past few earthquake cycles. Both the observed constant slip ratio in this study and offset clusters from previous field measurements support that the Chaldiran fault exhibits a characteristic slip behaviour in the late Quaternary. By combining all the geodetic and geological observations in the broader region, we also present a brief analysis of the slip rate of the Chaldiran fault. [ABSTRACT FROM AUTHOR]

SUMMARY: Unzen volcano, located on Shimabara Peninsula, Nagasaki, Japan, is an active volcano that has been intensively monitored since 1989, one year before the most recent eruption in 1990–1995. Previous earthquake and surface deformation studies have revealed that magma is transported obliquely from a magma reservoir beneath Tachibana Bay, to the west of Shimabara Peninsula. Here, we conduct broad-band magnetotelluric (MT) surveys at 99 sites around Shimabara Peninsula to investigate the crustal structure beneath Unzen volcano that is related to magma migration. A 3-D resistivity model that is constructed from 25 broad-band MT sites and 45 telluric sites shows a broad high-resistivity zone beneath Shimabara Peninsula and low-resistivity zones to the west and east of the peninsula. An unexpected observation is the spatial alignment of the high-resistivity zone with a seismic low-velocity zone (LVZ) at 3–15 km depth.  Quantitative analysis indicates this high-resistivity zone contains < 0.7 per cent melt under the assumption that the melt is stored in a good porosity network, while < 11 per cent melt in relatively poor pore network. We infer this high-resistivity, LVZ to be a highly crystallized mush zone (HCMZ) with low permeability. The hypocentres and pressure sources of the 1990–1995 eruption are distributed along the boundary between the high- and low-resistivity zones beneath the western part of the peninsula. We therefore conclude that the magma migrated along a structural boundary that possessed a relatively high permeability. Previous studies have suggested that eruptible magma is usually transported vertically upward through the centre of the mush zone, whereas the present results reveal that magma can be transported along the upper boundary of an HCMZ. [ABSTRACT FROM AUTHOR]

Monitoring and modeling ground surface deformation are crucial for the dynamic assessment of geothermal resources and sustainable exploitation in a geothermal field. In this study, we extract the deformation in the Yangbajing geothermal field by the small baseline-synthetic aperture radar interferometry (SBAS-InSAR) method using 141 Sentinel-1A images collected between March 2017 and November 2021. The InSAR result indicates both uplift and subsidence in the geothermal field. Subsequently, we use a dual-source model combining a dipping ellipsoid and a rectangular surface to model the shallow reservoir that contracted to cause the complex subsidence field in the north of the geothermal field. The shallow reservoir that expanded to cause the uplift in the south is modeled by an ellipsoid source. The parameters inversion is processed by the nonlinear Bayesian inversion method which has been applied to search the optimal parameters setting in a priori space and evaluate the uncertainties by the confidence intervals. To validate our inversion results, we collect the data from 31 wells including extraction and injection wells, detailed tectonic survey data and geothermal isotherms of the Yangbajing geothermal field. The modeled shallow reservoir in the north is within a 140 °C ground isotherm and a 160 °C underground isotherm. The thickness and depth of the reservoir slightly exceed the data from producing wells, probably due to the structural subsidence of the reservoir. The geometric structure is consistent with the channel faults that control the formation and development of the shallow reservoirs. The modeled reservoir under the rebound area in the south is also validated by comparing the data. Furthermore, we estimate the volume loss and recovery of the contracting and expanding reservoirs. We analyse the deformation mechanisms by considering the tectonic formation, the reservoir structure, and the extraction and rechargeability of the reservoirs. This research suggests the deformation pattern of shallow porous reservoirs that have similar layer composition and tectonic structure elsewhere in the world. Moreover, it provides a theoretical model to explore the parameters and volume change of geothermal reservoirs in plateau-embedded basins around the world. Our work is significant for the conservation and development of geothermal energy in geothermal fields that experience deformation, both subsidence and uplift. [ABSTRACT FROM AUTHOR]

The 2017 M w 5.5 Pohang earthquake occurred near an enhanced geothermal system site and generated thousands of aftershocks, the largest of which, a M w 4.6 earthquake, occurred 87 d after the mainshock. Redistribution of the groundwater pressure perturbed by the mainshock has been suggested as a cause of the post-seismic stress changes triggering several aftershocks, including the time-delayed event. However, to date, possible effects of variations in pore pressure on the aftershock occurrence have not been quantified in this region. Therefore, we conducted poroelastic modelling to evaluate this contribution to spatiotemporal distribution of the aftershocks, including the delayed event, using a fully coupled hydromechanical code. To construct a poroelastic model, a segmented fault geometry and a layered lithological structure were used. In addition, we utilized a kinematic slip model, a split-node algorithm and in-situ properties to simulate reliable coseismic and post-seismic behaviour. Our reference model successfully reproduced coseismic surface deformation in a line-of-sight direction, comparable to the corresponding observation from interferometric synthetic aperture radar, and was calibrated using groundwater measurement in a well. In addition to constructing the reference model, a series of numerical simulations were conducted to explore the effects and sensitivities of various hydraulic conductivities. Finally, the modelled Coulomb stress changes and spatiotemporal distribution of the aftershocks were analysed to elucidate the transient triggering mechanisms based on conditional statements to classify the mechanisms into several subsets. The classification showed that the poroelastic effect driven by depth/conductivity-dependent fluid diffusion is more critical to aftershock occurrence than the diffusion in the entire simulation time, and we propose that the delayed earthquake of Mw 4.6 could be correlated with poroelastic triggering rather than diffusion triggering. Furthermore, we inferred that this poroelastic effect could contribute to decay of aftershocks, particularly relatively small-magnitude aftershocks, as well as slow this decay in bedrocks. However, the proposed model does not explain all of the observed aftershocks, and other driving forces or triggering mechanisms need to be considered. [ABSTRACT FROM AUTHOR]

Saturn's rings are composed of icy grains, most in the mm to m size ranges, undergoing several collisions per orbit. Their collective behaviour generates a remarkable array of structures over many orders of magnitude, much of it not well understood. On the other hand, the collisional properties and parameters of individual ring particles are poorly constrained; usually, N -body simulations and kinetic theory employ hard-sphere models with a coefficient of restitution ϵ that is constant or a decreasing function of impact speed. Due to the plastic deformation of surface regolith, however, it is likely that ϵ will be more complicated, at the very least a non-monotonic function. We undertake N -body simulations with the REBOUND code with non-monotonic ϵ laws to approximate surfaces that are friable but not sticking. Our simulations reveal that such ring models can support two thermally stable steady states for the same (dynamical) optical depth: a cold and a warm state. If the ring breaks up into radial bands of one or the other, we find that warmer states tend to migrate into the colder states via a coherent travelling front. We also find stationary 'viscous' fronts, which connect states of different optical depths, but the same angular momentum flux. We discuss these preliminary results and speculate on their implications for structure formation in Saturn's B and C-rings, especially with respect to structures that appear in Cassini images but not in occultations. [ABSTRACT FROM AUTHOR]

We consider a large class of physical fields |$u$| written as double inverse Fourier transforms of some functions |$F$| of two complex variables. Such integrals occur very often in practice, especially in diffraction theory. Our aim is to provide a closed-form far-field asymptotic expansion of |$u$|⁠. In order to do so, we need to generalise the well-established complex analysis notion of contour indentation to integrals of functions of two complex variables. It is done by introducing the so-called bridge and arrow notation. Thanks to another integration surface deformation, we show that, to achieve our aim, we only need to study a finite number of real points in the Fourier space: the contributing points. This result is called the locality principle. We provide an extensive set of results allowing one to decide whether a point is contributing or not. Moreover, to each contributing point, we associate an explicit closed-form far-field asymptotic component of |$u$|⁠. We conclude the article by validating this theory against full numerical computations for two specific examples. [ABSTRACT FROM AUTHOR]

Earthquakes often occur in regions with complex material structure, such as sedimentary basins or mantle wedges. However, the majority of co-seismic modelling studies assume a simplified, often homogeneous elastic structure in order to expedite the process of model construction and speed up calculations. These co-seismic forward models are used to produce Green's functions for finite-fault inversions, so any assumptions made in the forward model may introduce bias into estimated slip models. In this study, we use a synthetic model of a sedimentary basin to investigate the impact of 3-D elastic structure on forward models of co-seismic surface deformation. We find that 3-D elastic structure can cause changes in the shape of surface deformation patterns. The magnitude of this effect appears to be primarily controlled by the magnitude of contrast in material properties, rather than the sharpness of contrast, the fault orientation, the location of the fault, or the slip orientation. As examples of real-world cases, we explore the impact of 3-D elastic structure with a model of the Taipei basin in Taiwan and a simulated earthquake on the Sanchaio fault, and with a 3-D geologic model of the San Francisco Bay Area and a slip model of the 1984 Morgan Hill earthquake on the Calaveras fault. Once again, we find that the presence of the basin leads to differences in the shape and amplitude of the surface deformation pattern, but we observe that the primary differences are in the magnitude of surface deformation and can be accounted for with a layered elastic structure. Our results imply that the use of homogeneous Green's functions may lead to bias in inferred slip models in regions with sedimentary basins, so, at a minimum, a layered velocity structure should be used. [ABSTRACT FROM AUTHOR]

Large thrust faults accommodate the convergence between India and Tibet along the southern margin of the Himalaya and have a history of producing great earthquakes that cause widespread damage. Along most parts of the Himalaya, there is geomorphological evidence that these thrusts can rupture to the surface in M w >8 earthquakes. However, in the Himalayan state of Jammu & Kashmir (NW India), the thrust faults are blind and large-scale folding is the only expression of active deformation at the surface, making it difficult to assess the seismic hazard in this region. In this paper, we use field, satellite, and seismological observations to determine the fault geometry in Jammu & Kashmir. We then estimate the ground motions from potential earthquakes in the region using models of the seismic wavefield that would be generated if the thrust fault beneath Jammu & Kashmir were to rupture. We find that earthquakes that rupture the buried, shallow part of the locked Main Himalayan Thrust could generate peak ground velocities that are >3 times larger than earthquakes of the same magnitude on its deeper portions. We also model the ground motions that would result from the thrust fault geometries representative of different parts of the Himalayan arc. These simulations show that even seemingly minor variations in the shallow fault geometry can lead to large differences in the expected ground motions, highlighting the importance of accurately determining the shallow geometry of thrust faults along the margins of mountain ranges for estimating seismic hazard. [ABSTRACT FROM AUTHOR]

The propagation of Tibetan Plateau high-elevation towards the Asian interior remains poorly understood. We conduct a series of 2-D viscoelastic plastic finite element models addressing the development of the Liupan Shan fold-and-thrust belt with multiple detachments in the NE Tibetan Plateau. We focus on the spatio–temporal relations of the Liupan Shan fold-and-thrust belt and dynamic mechanisms of growth in the NE margin of Tibetan Plateau during the late Cenozoic. Our models, consisting of a seismically-constrained two-layer-crust (of variable thicknesses) plus lithospheric mantle consider tectonic horizontal contraction and the influence of elastoplastic materials, including Model 1 (only the upper crust is elastoplastic), Model 2 (both the upper crust and the middle part of the lower crust are elastoplastic) and Model 3 (the whole lithosphere is elastoplastic). The results show that Model 1 can only reflect the properties of shallow detachment layer. Models 2 and 3 can better reflect the properties of multidetachment layers and approximately represent the process of fold-and-thrust belt formation and evolution; thus suggesting the importance of plastic rheology in the deep crust. In addition, Model 3 suggests that the pre-existing fault near the Moho below Liupan Shan is activated during the evolution process of ∼8 Ma. In summary, the results embody that the initiation and evolution of the fold-and-thrust belt depend on detachment structures: the deep and shallow detachments control the evolution of surface deformation and fold-and-thrust belt altogether. The simulated surface deformation is generally consistent with actual surface elevation distribution. Moreover, the growth of the Liupan Shan range may represent an expansion of the NE Tibetan Plateau toward the Ordos basin vicinity. Additionally, the tectonically driven uplift of Liupan Shan area supports the general crustal shortening tectonic mechanism of the Tibetan Plateau. [ABSTRACT FROM AUTHOR]

SUMMARY: Volcano seismology is an essential tool for monitoring volcanic processes in the advent and during eruptions. A variety of seismic signals can be recorded at volcanoes, of which some are thought to be related to the migration of fluids which is of primary importance for the anticipation of imminent eruptions. We investigate the volcanic crises at Villarrica volcano in 2015 and report on a newly discovered very-long-period (VLP) signal that accompanies phases of periodic long period (LP) signal burst. Despite their low amplitude emergent character, we can locate the source region of the 1 Hz LP signals to the close vicinity of the volcano using a network-based correlation method. The source of the VLP signal with a period of about 30–100 s appears to locate in the vicinity of two stations a few kilometres from the summit. Both stations record very similar VLP waveforms that are correlated with the envelope of the LP bursts. A shallow magma reservoir was inferred by Contreras from surface deformation as the source of inflation following the eruption in 2015. Cyclic volume changes of 6 m3 in this reservoir at 3 km depth can explain the observed amplitudes of the vertical VLP signal. We propose that the LP signal is generated by the migration of gas or gas-rich magma that is periodically released from the inflating reservoir through a non-linear valve structure which modulates the flux, and thereby causes bursts of flow-related LP signals and pressure changes observed as VLP deformation. Our model predicts that the correlated occurrence of LP bursts and VLP surface motion depends on the intensity of the fluid flux. A weaker flux of fluids may not exceed the opening pressure of valve structure, and higher rates might maintain pressure above the closing pressure. In both cases, the VLP signal vanishes. Our observation provides constrains for models of fluid transport inside volcanoes. At Villarrica the VLP signal, and its relation to the LP activity, reveal additional information about fluxes in the magmatic reservoir that might aide forecasting of volcanic activity. [ABSTRACT FROM AUTHOR]

On 29 December 2020, a shallow earthquake of magnitude M w 6.4 struck northern Croatia, near the town of Petrinja, more than 24 hr after a strong foreshock (M L 5). We formed a reconnaissance team of European geologists and engineers, from Croatia, Slovenia, France, Italy and Greece, rapidly deployed in the field to map the evidence of coseismic environmental effects. In the epicentral area, we recognized surface deformation, such as tectonic breaks along the earthquake source at the surface, liquefaction features (scattered in the fluvial plains of Kupa, Glina and Sava rivers), and slope failures, both caused by strong motion. Thanks to this concerted, collective and meticulous work, we were able to document and map a clear and unambiguous coseismic surface rupture associated with the main shock. The surface rupture appears discontinuous, consisting of multi-kilometre en échelon right stepping sections, along a NW–SE striking fault that we call the Petrinja-Pokupsko Fault. The observed deformation features, in terms of kinematics and trace alignments, are consistent with slip on a right lateral fault, in agreement with the focal solution of the main shock. We found mole tracks, displacement on faults affecting natural features (e.g. drainage channels), scarplets and more frequently breaks of anthropogenic markers (roads, fences). The surface rupture is observed over a length of ∼13 km from end-to-end, with a maximum displacement of 38 cm, and an average displacement of ∼10 cm. Moreover, the liquefaction extends over an area of nearly 600 km2 around the epicentre. Typology of liquefaction features include sand blows, lateral spreading phenomenon along the road and river embankments, as well as sand ejecta of different grain size and matrix. Development of large and long fissures along the fluvial landforms, current or ancient, with massive ejections of sediments is pervasive. These features are sometimes accompanied by small horizontal displacements. Finally, the environmental effects of the earthquake appear to be reasonably consistent with the usual scaling relationships, in particular the surface faulting. This rupture of the ground occurred on or near traces of a fault that shows clear evidence of Quaternary activity. Further and detailed studies will be carried out to characterize this source and related faults in terms of future large earthquakes potential, for their integration into seismic hazard models. [ABSTRACT FROM AUTHOR]

Large earthquakes can alter regional groundwater pressure, resulting in fluid flow, and the process of restoring hydrostatic equilibrium would in turn lead to observable surface deformation, termed poroelastic rebound, which is one of the most important post-seismic mechanisms for stress transfer and triggering. To constrain the poroelastic contributions to the early post-seismic deformation, we model the hydrologic response within 1.5 months following the 2015 M w 8.3 Illapel earthquake and remove its effects from the observed geodetic signals. Results demonstrate the post-seismic fluid-flow patterns from the co-seismic high-slip region to the north and south sides, and the northern poroelastic effects are remarkably stronger than those on the south side, verified by northern liquefaction phenomena. Therefore, previous pure afterslip models overestimate the asperities on both flanks of the co-seismic rupture zone and underestimate the middle region, with local errors of more than 50 per cent. It highlights the importance of considering the poroelastic effects, when modelling the transient post-seismic deformation. [ABSTRACT FROM AUTHOR]

Blade bearing is the key mechanism of the controllable pitch propellers (CPPs) to adjust the pitch to change the thrust. In the marine environment, wear and deformation of the contact surface of the blade bearing are quite large, which may result in leakage to reduce the service life of CPP. However, most studies on tolerance optimization of blade bearing disregard deformations and wear, and it could result in unreasonable tolerances and impair the CPP's service life. In this paper, to establish the statistical tolerance-cost-service life optimization model of the blade bearing of the CPP considering the marine environment effects in propellers, a prediction model about the relationship between the thrust and rotation speed of the CPP considering environment conditions of wind, wave, and ocean currents is constructed. Based on this, the wear and deformation are simulated according to the predicted rotation speed in the marine environment. Considering the deformation and wear in the marine environment, a tolerance analysis model of the assembly deviations of the blade bearing is introduced, and the assembly characteristics and functional requirements are established and converted into probability constraints to establish the service-life model based on statistical tolerance analysis and Monte Carlo simulation methods. Afterward, a statistical tolerance-cost-service life optimization model is established and transformed into a single-objective optimization problem through the normalization method and the weighted sum method, and solved by a meta-heuristic algorithm named cuckoo search. Finally, the application of the proposed model is demonstrated through a case involving the blade bearing of a CPP, and the optimized results show that the manufacturing cost of the blade bearing is reduced and its service life is extended. [ABSTRACT FROM AUTHOR]

We study compact hyperbolic surface laminations. These are a generalisation of closed hyperbolic surfaces, which appear to be more suited to the study of Teichmüller theory than arbitrary non-compact surfaces. We show that the Teichmüller space of any non-trivial hyperbolic surface lamination is infinite dimensional. In order to prove this result, we study the theory of deformations of hyperbolic surfaces, and we derive a new formula for the derivative of the length of a simple closed geodesic with respect to the action of grafting. This formula complements those derived by McMullen [ 24 ], in terms of the Weil–Petersson metric, and by Wolpert [ 34 ], for the case of earthquakes. [ABSTRACT FROM AUTHOR]

Estimates of line-of-sight (LOS) displacements from Interferometric Synthetic Aperture Radar (InSAR) observations serve as the basis of the long-term monitoring of an operating natural gas storage site at Honor Rancho in California. An inversion algorithm is used to estimate the portion of the signal that is attributable to deformation within the gas storage reservoir, located at a depth of around 3 km. Removing this contribution produces residuals that are used to characterize the background variation is surface deformation at the gas storage facility and to determine a threshold that can signify unusually large residuals. An application to almost 7 yr of InSAR data, from 2011 until 2018, indicates that there are intervals of heightened residuals as well as brief episodes of anomalously large misfits. An examination of the spatial distributions of the individual residual LOS displacements indicates larger displacements in an alluvial valley just south of the reservoir, with rapid spatial variations in sign, indicating a rather shallow origin. Furthermore, the two anomalous events also involve rapid spatial variations in the LOS displacement residuals directly above the storage facility. The results demonstrate that the technique of extracting residuals after removing the reservoir signal is a useful approach, even in the case of this deep reservoir, and is a promising method for long-term monitoring. [ABSTRACT FROM AUTHOR]

The Eastern Mediterranean is the most seismically active region in Europe due to the complex interactions of the Arabian, African, and Eurasian tectonic plates. Deformation is achieved by faulting in the brittle crust, distributed flow in the viscoelastic lower-crust and mantle, and Hellenic subduction, but the long-term partitioning of these mechanisms is still unknown. We exploit an extensive suite of geodetic observations to build a kinematic model connecting strike-slip deformation, extension, subduction, and shear localization across Anatolia and the Aegean Sea by mapping the distribution of slip and strain accumulation on major active geological structures. We find that tectonic escape is facilitated by a plate-boundary-like, trans-lithospheric shear zone extending from the Gulf of Evia to the Turkish-Iranian Plateau that underlies the surface trace of the North Anatolian Fault. Additional deformation in Anatolia is taken up by a series of smaller-scale conjugate shear zones that reach the upper mantle, the largest of which is located beneath the East Anatolian Fault. Rapid north–south extension in the western part of the system, driven primarily by Hellenic Trench retreat, is accommodated by rotation and broadening of the North Anatolian mantle shear zone from the Sea of Marmara across the north Aegean Sea, and by a system of distributed transform faults and rifts including the rapidly extending Gulf of Corinth in central Greece and the active grabens of western Turkey. Africa–Eurasia convergence along the Hellenic Arc occurs at a median rate of 49.8 mm yr–1 in a largely trench-normal direction except near eastern Crete where variably oriented slip on the megathrust coincides with mixed-mode and strike-slip deformation in the overlying accretionary wedge near the Ptolemy–Pliny–Strabo trenches. Our kinematic model illustrates the competing roles the North Anatolian mantle shear zone, Hellenic Trench, overlying mantle wedge, and active crustal faults play in accommodating tectonic indentation, slab rollback and associated Aegean extension. Viscoelastic flow in the lower crust and upper mantle dominate the surface velocity field across much of Anatolia and a clear transition to megathrust-related slab pull occurs in western Turkey, the Aegean Sea and Greece. Crustal scale faults and the Hellenic wedge contribute only a minor amount to the large-scale, regional pattern of Eastern Mediterranean interseismic surface deformation. [ABSTRACT FROM AUTHOR]

We present a theory of modern, thermally induced deformation in a realistic Earth. The heat conduction equation is coupled with standard elastic deformation theory to construct a boundary-value problem comprised of eighth-order differential equations. The accurate and stable dual variable and position propagating matrix technique is introduced to solve the boundary-value problem. The thermal load Love numbers are defined to describe the displacements and potential changes driven by thermally induced deformation. The proposed analytical method is validated by comparing the present results with exact solutions for a homogeneous sphere, which are also derived in this paper. The analytical method is then applied to a realistic Earth model to evaluate the effects of layering and self-gravitation of the Earth on displacement and changes of potential. Furthermore, the frequency dependence in the thermal load is illustrated by invoking different thermal periodicities in the computation. Thermal anisotropy is also considered by comparing the results obtained using isotropic and transversely isotropic Earth models. Results show that, when simulating thermally induced deformation, invoking a homogeneous spherical Earth leads to results that substantially differ from those obtained using a more realistic Earth model. [ABSTRACT FROM AUTHOR]

In this study, we devised a new set of analytical foundation solutions to compute the internal co-seismic displacement and strain changes caused by four independent point sources (strike-slip, dip-slip, horizontal tensile and vertical tensile) inside a homogeneous spherical earth model. Our model provides constraints on the deformation properties at depth and reveals that the internal co-seismic deformation is larger than that on the surface. The deformation near the source is convergent with our formulae. For the internal deformation at radial section plane, the patterns of horizontal displacements |${u_\theta },{u_\phi }$| and strain changes |${e_{{ rr}}},{e_{\theta \theta }},{e_{\phi \phi }},{e_{\theta \phi }}$| caused by strike-slip and tensile sources appear symmetric at the equidistance above and below the source. Their amplitudes are not identical but with a small discrepancy actually. Unlike these, the patterns of radial displacements |${u_r}$| for strike-slip and tensile sources exhibit point symmetry with the equidistance from the source. Also, the corresponding amplitudes are slightly different. The displacements |${u_\theta },{u_\phi }$| and strain changes |${e_{{ rr}}},{e_{\theta \theta }},{e_{\phi \phi }},{e_{\theta \phi }}$| caused by dip-slip also show the same properties as |${u_{ r}}$| of the strike-slip source. The magnitudes of the displacements and strain changes depend on the source types. The curvature effect on the near-field surface deformations is small, and it increases with the studied depth. However, for the far-field deformation caused by the strike-slip source (d s = 20 km), the curvature effect can be as large as 77 per cent when the epicentral distance approximates to 1778 km. [ABSTRACT FROM AUTHOR]

Geodetic observations from space continuously record surface deformation and global mass redistribution with an increasing accuracy. In parallel, surficial processes (oceanic, atmospheric and hydrological loading) are more and more precisely modelled. We propose a confrontation of the geodetic Global Positioning System (GPS) and gravity-field satellite laser ranging (SLR) observations at decadal and interannual timescales, in terms of resolution, correlation and comparison with surficial loading models. We focus on the largest global scale signals of degree 2. At interannual periods, surface deformations retrieved from GPS time-series do not exceed 0.8 mm. Our analysis does not reveal the presence of a dominant signal at a specific period, except perhaps for a signal of approximately 3 yr likely connected to the loading response to El Niño/Southern Oscillations. Contrary to the results of previous studies, we do not find in GPS time-series a clear 6-yr oscillation associated with a degree-2 order-2 pattern. Interannual variations in the degree-2 Stokes coefficients of the gravity field do not exceed 2 × 10−11. We do not detect a dominant gravity signal at one specific period but instead a broad spectrum of frequencies. The comparison between the degree-2 deformations built from GPS time-series with a prediction from SLR-derived gravity variations reveals some correlations, though their differences remain important. This highlights the present day limitations of these techniques in their ability to characterize global scale interannual variations. Hydrological loading models show some correlations with both GPS and SLR signals, but we cannot firmly establish that continental hydrology is dominantly responsible for the observed variations. Given the current limits in the resolution of both gravity and surface deformation and in the modelling of surface processes, we conclude that it will be a challenge to retrieve a geodetic signal of subdecadal period originating in the Earth's core. [ABSTRACT FROM AUTHOR]

While local earthquake tomography is typically used to image the crust, this technique has restricted depth penetration due to short receiver-source distances. Regional tomography however aims to image the upper mantle from teleseismic events but suffers from poor resolution from 0 down to 40 km depth. We present here a hybrid method that combines the two approaches taking advantage of the short-wavelength resolution within the crust to better constrain the ray path at depth, and thus to improve the lithospheric imaging. Using this new method enhances the continuity or disruption of mantle anomalies towards the surface. Such hybrid tomographic images of crust-to-upper mantle structures are then critical to understand the relation and interplay between the thermal and mechanical lithospheric processes and the role in the localization of the deformation at the surface. We apply our approach to the North Tanzanian Divergence (NTD), where those processes interact with a cold cratonic lithosphere. Our new tomographic images clearly demonstrate the impact of deep-seated processes on surface features. First, strong lateral velocity anomalies and clustered seismicity in the crust are consistent with the surface geology of the NTD (rifted basins, volcanoes and border faults). Then, at a lithospheric scale, the velocity distribution highlights the major role of inherited structures in guiding the rift opening. In particular, our study suggests a strong influence of the Masai cratonic block, south of the NTD, in the rift evolution. The transition from the north–south axial valley into three diverging rift arms (Eyasi, Natron-Manyara and Pangani) is likely due to the change in rheology and to the presence of magma along inherited sutures between the craton and the mobile belts. [ABSTRACT FROM AUTHOR]

Planetary scale interannual deformations of the Earth's surface, of millimetric amplitude, have recently been related to both geomagnetic field changes and motion within the fluid outer core. We calculate the temporal variations of the dynamical pressure at the surface of the core associated with core flow models inverted from geomagnetic observations. From these we compute predictions of the changes in Earth's topography in response to elastic deformations in the mantle. We show that at decadal periods, the predicted changes in Earth's topography are at most of the order of 0.3 mm. Focused at interannual periods between 4 and 9.5 yr, the predicted topography variations are smaller than 0.05 mm, at least an order of magnitude smaller than the reported observations. These amplitudes are only weakly sensitive to the choice of hypothesis used to reconstruct fluid motions at the core surface. We conclude that surface deformations induced by dynamical pressure changes in the core are below the detection level at present-day. Alternative geophysical sources must be sought to explain the observed millimetric interannual variations of the planetary scale topography, and its associated gravity variations. We currently see no justification for a physical relationship between interannual fluctuations of the geomagnetic field and of Earth's observed deformations. We conjecture that the largest gravity signal of core origin is potentially associated with decadal longitudinal oscillations of the inner core. It might be detectable as longer series will become available. [ABSTRACT FROM AUTHOR]

A noncommutative deformation of a quadric surface is usually described by a three-dimensional cubic Artin–Schelter regular algebra. In this paper we show that for such an algebra its bounded derived category embeds into the bounded derived category of a commutative deformation of the Hilbert scheme of two points on the quadric. This is the second example in support of a conjecture by Orlov. Based on this example we formulate an infinitesimal version of the conjecture and provide some evidence in the case of smooth projective surfaces. [ABSTRACT FROM AUTHOR]

Due to the steady moderate seismicity observed along the Briançon seismic arc, in the south-western French Alps, three temporary GNSS (Global Navigation Satellite System) surveys took place in 1996, 2006 and 2011, across a ∼50 × 60 km² wide area, to investigate the surface deformation field. The horizontal velocity field computed from these three surveys showed an east–west extension in the network. A fourth campaign was led in 2016, creating a 20 yr observation span, resulting in measurements which reach a sufficient accuracy to assess whether extension found within the Briançon network is localized onto any particular tectonic feature. Several faults in this area are known to be active normal faults. Assessing the localization of the deformation may lead to a better understanding of the active tectonics of the Alpine belt. To address this issue, a robust velocity field was computed from the combination of the different campaign and permanent GNSS data. Strain rate tensors were derived for the first time in this area on a 0.1 × 0.1 deg grid to assess the distribution of the deformation. The regional deformation appears localized in the Briançon area and reaches up to 20 ± 5 nanostrain yr−1 in the centre of the network. The observed velocities were projected on a profile across the network and compared with modelled interseismic deformation to characterize the behaviour of the major active faults known in the study zone. While a two-fault model provides the best fit to the data, a single fault model has only marginally higher residuals, with parameters which are more consistent with the seismotectonics of the region. The localization of the single modelled fault is consistent with the location of the High Durance Fault (HDF). Therefore, we used the known geological location of this structure as a priori information in a block model to compute a fault slip rate at the interface between the two blocks. The velocities on the interface indicate 0.4–0.5 mm yr−1 of extension, and therefore strain accumulates along the HDF throughout the seismic cycle. The geodetically derived fault slip rate is converted into an equivalent seismic moment release rate, which is consistent within its uncertainty bounds with the known historical and instrumental seismicity of the Briançon area. [ABSTRACT FROM AUTHOR]

The volcanism of the Eifel volcanic field (EVF), in west-central Germany, is often considered an example of hotspot volcanism given its geochemical signature and the putative mantle plume imaged underneath. EVF's setting in a stable continental area provides a rare natural laboratory to image surface deformation and test the hypothesis of there being a thermally buoyant plume. Here we use Global Positioning System (GPS) data to robustly image vertical land motion (VLM) and horizontal strain rates over most of intraplate Europe. We find a spatially coherent positive VLM anomaly over an area much larger than the EVF and with a maximum uplift of ∼1 mm yr−1 at the EVF (when corrected for glacial isostatic adjustment). This rate is considerably higher than averaged over the Late-Quaternary. Over the same area that uplifts, we find significant horizontal extension surrounded by a radial pattern of shortening, a superposition that strongly suggests a common dynamic cause. Besides the Eifel, no other area in NW Europe shows significant positive VLM coupled with extensional strain rates, except for the much broader region of glacial isostatic adjustment. We refer to this 3-D deformation anomaly as the Eifel Anomaly. We also find an extensional strain rate anomaly near the Massif Central volcanic field surrounded by radial shortening, but we do not detect a significant positive VLM signal there. The fact that the Eifel Anomaly is located above the Eifel plume suggests that the plume causes the anomaly. Indeed, we show that buoyancy forces induced by the plume at the bottom of the lithosphere can explain this remarkable surface deformation. Plume-induced deformation can also explain the relatively high rate of regional seismicity, particularly along the Lower Rhine Embayment. [ABSTRACT FROM AUTHOR]

An algorithm developed in Tanimoto & Wang for deriving shallow elasticity structure from colocated pressure and seismic instruments is applied to data at nine stations in the Piñon Flat Observatory (PFO) in Southern California. Depth resolution kernels indicate that this approach can recover near-surface shear-modulus structure in the upper 50–100 m of the Earth. Our estimate for Vs30 at the closest station (BPH01) agrees with an independent result by Yong et al. which was based on an active source survey with a geophone array. Vs30 derived from the borehole S-wave speed model by Fletcher et al. also agrees with estimates at nearby locations, considering the size of uncertainties. Derived structures are much smoother than their models, however. This is because our method analyses slow quasi-static deformation of the solid Earth in the frequency range 0.005–0.05 Hz that inherently limits depth resolution. Lack of detailed variations in structures does not seem to be the problem for estimating Vs30 as Vs30 is an averaged quantity from the upper 30 m of the Earth. Also our estimates are not affected by different choices of frequency range in the inversion. However, our method may have some difficulty at stiff, hard-rock sites because ground deformation caused by surface pressure can become small and the key observables in our method, the ratios between seismic amplitude and surface pressure change, become difficult to measure accurately. [ABSTRACT FROM AUTHOR]

A comprehensive GPS velocity field along the Dead Sea Fault System (DSFS) provides new constraints on along-strike variations of near-transform crustal deformation along this plate boundary, and internal deformation of the Sinai and Arabian plates. In general, geodetically derived slip rates decrease northwards along the transform (5.0 ± 0.2 to 2.2 ± 0.5 mm yr−1) and are consistent with geological slip rates averaged over longer time periods. Localized reductions in slip rate occur where the Sinai Plate is in ∼N–S extension. Extension is confined to the Sinai side of the fault and is associated with prominent changes in transform geometry, and with NW–SE striking, left-lateral splay faults, including the Carmel Fault in Israel and the Roum Fault in Lebanon. The asymmetry of the extensional velocity gradients about the transform reflects active fragmentation of the Sinai Plate along the continental margin. Additionally, elastic block modelling of GPS velocities requires an additional structure off-shore the northern DSF segment, which may correspond with a fault located along the continental margin, suggested by prior geophysical studies. [ABSTRACT FROM AUTHOR]

Using a numerical integral method, we deduced a set of formulae for the co-seismic internal deformation in a spherically symmetric earth model, simultaneously taking self-gravitation, compressibility and realistically stratified structure of the Earth into account. Using these formulae, we can calculate the internal deformation at an arbitrary depth caused by an arbitrary seismic source. To demonstrate the correctness of our formulae, we compared our numerical solutions of radial functions with analytical solutions reported by Dong & Sun based on a homogeneous earth model; we found that two sets of results agree well with each other. Our co-seismic internal Green's functions in the near field agree well with the results calculated by the formulae of Okada, which also verifies our Green's functions. Finally, we calculated the Coulomb stress changes on the Japanese Islands and Northeast China induced by the Tohoku-Oki M w 9.0 earthquake using the methods described above. We found that the effect of layered structure plays a leading role on the near field, while curvature occupies a dominant position on the deep region of the far field. Through a comparison of the Coulomb stress changes at a depth of 10 km on a layered earth model calculated by our method along with the corresponding results of Okada, we found that the discrepancy between them in near field was ∼31.5 per cent, and that of far field was >100 per cent of the signals. [ABSTRACT FROM AUTHOR]

A sequence of earthquake events consisting of three large shocks occurred in Central Italy from August to October in 2016 with the duration of almost 2 months. The preliminary study on the seismic mechanism suggests that the sequence of events is the result from the activity of the SW dipping Mt Bove–Mt Vettore–Mt Gorzano normal fault system. For investigation and understanding of the coseismic faulting of the seismogenic fault alignment, we collect a set of comprehensive satellite observations including the Sentinel-1A, ALOS-2/PALSAR-2 and GPS data to map the coseismic surface deformation and estimate the source models in this study. The derived faulting model for the first Amatrice event is characterized by two distinct slip asperities suggesting that it is a predominantly normal dip-slip motion with slight strike-slip component. The second event, Visso earthquake is almost a purely normal rupture. The third Norcia event is dominated by the normal dip-slip rupture of the seismogenic fault, and has propagated up to the ground with significant slip. The three faulting models are then utilized to quantify the Coulomb failure stress (CFS) change over the seismic zone. First, the CFS change on the subsequent two seismogenic faults of the earthquake sequence is estimated, and the derived positive CFS change induced by the preceding earthquakes suggests that the early events have positive effects on triggering the subsequent seismicity. We then explore the response relation of the aftershocks including 961 events with magnitudes larger than M 3.0 to the CFS change over the seismic zone. It suggests that the rupture pattern of the aftershocks is similar to the major shocks with predominantly normal dip-slip. To assess the risk of the future seismic hazard, we analyse quantitatively the spatial distribution of aftershock occurrence and CFS transfer at the seismogenic depth, indicating that the ruptures of the three major shocks do partly release the accumulated strain on the associated fault alignment as well as the dense aftershock, but the CFS increase zone with few aftershocks in the southwest of the eastern Quaternary fault alignment of Central Italy poses the potential of further rupture. In particular, the distribution of aftershock migration also suggests that the north extension of the Mt Bove fault is the potential zone with rupture risk. [ABSTRACT FROM AUTHOR]

Quasigeoid models can be determined from surface gravity anomalies, so are sensitive to changes in the shape of the topography as well as changes in gravity. Here we present results of forward modelling gravity/quasigeoid changes from synthetic aperture radar data following the 2016 M w 7.8 Kaikōura earthquake with land uplift of up to 10 m. We assess the impact of the topographic deformation on the reference surface of the New Zealand vertical datum in lieu of costly field gravity field measurements. The most significant modelled gravity and quasigeoid changes are—2.9 mGal and 5–7 mm, respectively. We compare our forward modelled gravity signal to terrestrial gravity observation data and show that differences between the data sets have a standard deviation of ±0.1 mGal. The largest modelled change in the quasigeoid is an order of magnitude smaller than the 57.7 mm estimated precision of the most recently computed NZGeoid model over the Kaikōura region. Modelled quasigeoid changes implied by this particular deformation event are not statistically significant with respect to estimated precision of the New Zealand quasigeoid model. [ABSTRACT FROM AUTHOR]

Teleseismic waveforms contain information on fault slip evolution during an earthquake, as well as on the fault geometry. A linear finite-fault inversion method is a tool for solving the slip-rate function distribution under an assumption of fault geometry as a single or multiple-fault-plane model. An inappropriate assumption of fault geometry would tend to distort the solution due to Green's function modelling errors. We developed a new inversion method to extract information on fault geometry along with the slip-rate function from observed teleseismic waveforms. In this method, as in most previous studies, we assumed a flat fault plane, but we allowed arbitrary directions of slip not necessarily parallel to the assumed fault plane. More precisely, the method represents fault slip on the assumed fault by the superposition of five basis components of potency-density tensor, which can express arbitrary fault slip that occurs underground. We tested the developed method by applying it to real teleseismic P waveforms of the M W 7.7 2013 Balochistan, Pakistan, earthquake, which is thought to have occurred along a curved fault system. The obtained spatiotemporal distribution of potency-density tensors showed that the focal mechanism at each source knot was dominated by a strike-slip component with successive strike angle rotation from 205° to 240° as the rupture propagated unilaterally towards the south-west from the epicentre. This result is consistent with Earth's surface deformation observed in optical satellite images. The success of the developed method is attributable to the fact that teleseismic body waves are not very sensitive to the spatial location of fault slip, whereas they are very sensitive to the direction of fault slip. The method may be a powerful tool to extract information on fault geometry along with the slip-rate function without requiring detailed assumptions about fault geometry. [ABSTRACT FROM AUTHOR]

The Democratic People's Republic of Korea (North Korea) conducted its sixth and largest affirmed underground nuclear test on 2017 September 3. Analysis of Interferometric Synthetic Aperture Radar (InSAR) data revealed detailed surface displacements associated with the nuclear explosion. The nuclear explosion produced large-scale surface deformation causing decorrelation of the InSAR data directly above the test site, Mt. Mantap, while the flanks of the Mountain experienced displacements up to 0.5 m along the Line-of-Sight of the Satellite. We determined source parameters of the explosion using the Bayesian inversion of the InSAR data. The explosive yield was estimated as 245–271 kiloton (kt) of TNT, while the previous yield estimations range from 70–400 kt. We determined the nuclear source at a depth of 542 ± 30 m below Mt. Mantap (129.0769°E, 41.0324°N). We demonstrated that the Bayesian modelling of the InSAR data reduces the uncertainties in the source parameters of the nuclear test, particularly the yield and source depth that are otherwise poorly resolved in seismic methods. [ABSTRACT FROM AUTHOR]

Most geodetic inversions of surface deformation rates consider the depth distribution of interseismic fault slip-rate to be time invariant. However, some numerical simulations show downdip penetration of dynamic rupture into regions with velocity-strengthening friction, with subsequent updip propagation of the locked-to-creeping transition. Recently, Bruhat and Segall developed a new method to characterize interseismic slip rates, that allows slip to penetrate up dip into the locked region. This simple model considered deep interseismic slip as a crack loaded at its downdip end, and provided analytical expressions for stress drop within the crack, slip and slip rate along the fault. This study extends this approach to strike-slip fault environments, and includes coupling of creep to viscoelastic flow in the lower crust and upper mantle. I use this model to investigate interseismic deformation rates along the Carrizo Plain section of the San Andreas fault. This study reviews possible models, elastic and viscoelastic, for fitting horizontal surface rates. Using this updated approach, I develop a physics-based solution for deep interseismic creep which accounts for possible slow vertical propagation, and investigate how it improves the fit of the horizontal deformation rates in the Carrizo Plain region. I found solutions for fitting the surface deformation rates that allow for reasonable estimates for earthquake rupture depth and coseismic displacement and improves the overall fit to the data. Best-fitting solutions present half-space relaxation time around 70 yr, and very low propagation speeds, less than a metre per year, suggesting a lack of creep propagation. [ABSTRACT FROM AUTHOR]

An organism's ability to control the timing and direction of energy flow both within its body and out to the surrounding environment is vital to maintaining proper function. When physically interacting with an external target, the mechanical energy applied by the organism can be transferred to the target as several types of output energy, such as target deformation, target fracture, or as a transfer of momentum. The particular function being performed will dictate which of these results is most adaptive to the organism. Chewing food favors fracture, whereas running favors the transfer of momentum from the appendages to the ground. Here, we explore the relationship between deformation, fracture, and momentum transfer in biological puncture systems. Puncture is a widespread behavior in biology requiring energy transfer into a target to allow fracture and subsequent insertion of the tool. Existing correlations between both tool shape and tool dynamics with puncture success do not account for what energy may be lost due to deformation and momentum transfer in biological systems. Using a combination of pendulum tests and particle tracking velocimetry (PTV), we explored the contributions of fracture, deformation and momentum to puncture events using a gaboon viper fang. Results on unrestrained targets illustrate that momentum transfer between tool and target, controlled by the relative masses of the two, can influence the extent of fracture achieved during high-speed puncture. PTV allowed us to quantify deformation throughout the target during puncture and tease apart how input energy is partitioned between deformation and fracture. The relationship between input energy, target deformation and target fracture is non-linear; increasing impact speed from 2.0 to 2.5 m/s created no further fracture, but did increase deformation while increasing speed to 3.0 m/s allowed an equivalent amount of fracture to be achieved for less overall deformation. These results point to a new framework for examining puncture systems, where the relative resistances to deformation, fracture and target movement dictate where energy flows during impact. Further developing these methods will allow researchers to quantify the energetics of puncture systems in a way that is comparable across a broad range of organisms and connect energy flow within an organism to how that energy is eventually transferred to the environment. [ABSTRACT FROM AUTHOR]

To study more fully the characteristic law of deformation and failure of tuff jointed rock mass of prefabricated parallel discontinuous joint test specimens, the uniaxial compression test was used. The stress–strain curve, peak intensity, deformation parameters, energy characteristics, etc. of the rock test specimens were systematically studied under different combinations of joint dip angle and joint spacing. The research found that: (1) during the failure process of tuff, the peak intensity and elastic modulus followed a U-shaped change pattern and the minimum value was reached when α = 60°; (2) the fracture modes of test specimens with different joint dip angles were different. When α = 30° and 45°, failure characteristics were mixed modes of tensile or tensile shear failure. When α = 60°, failure characteristics were shear. At α = 75°, the failure characteristic was tensile shear failure. (3) The absorbed and dissipated energy of the rock increased nonlinearly at each stage of deformation. (4) We quantified rock energy damage through a correlation between dissipated energy and absorbed energy of the rock in the process of energy evolution, and obtained an evolution of the relationship between the dissipated energy ratio, crack dip angle and crack spacing. Based on different fracture distribution methods and according to the strain equivalence principle, the constitutive equation of the pre-peak rock damage was obtained. [ABSTRACT FROM AUTHOR]

Motivation Cell shape provides both geometry for, and a reflection of, cell function. Numerous methods for describing and modeling cell shape have been described, but previous evaluation of these methods in terms of the accuracy of generative models has been limited. Results Here we compare traditional methods and deep autoencoders to build generative models for cell shapes in terms of the accuracy with which shapes can be reconstructed from models. We evaluated the methods on different collections of 2D and 3D cell images, and found that none of the methods gave accurate reconstructions using low dimensional encodings. As expected, much higher accuracies were observed using high dimensional encodings, with outline-based methods significantly outperforming image-based autoencoders. The latter tended to encode all cells as having smooth shapes, even for high dimensions. For complex 3D cell shapes, we developed a significant improvement of a method based on the spherical harmonic transform that performs significantly better than other methods. We obtained similar results for the joint modeling of cell and nuclear shape. Finally, we evaluated the modeling of shape dynamics by interpolation in the shape space. We found that our modified method provided lower deformation energies along linear interpolation paths than other methods. This allows practical shape evolution in high dimensional shape spaces. We conclude that our improved spherical harmonic based methods are preferable for cell and nuclear shape modeling, providing better representations, higher computational efficiency and requiring fewer training images than deep learning methods. Availability and implementation All software and data is available at http://murphylab.cbd.cmu.edu/software. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

In this paper, we derive analytical expressions for the dislocation Love numbers (DLNs) for a layered, spherical, transversely isotropic and self-gravitating Earth. This solution is based on the spherical system of vector functions (or the vector spherical harmonics) and a new propagating matrix method called the dual variable and position method. The DLNs can be obtained with high accuracy to an arbitrarily high degree, thereby allowing a wide range of applications based on high resolution earth models. Compared to the traditional numerical integration approach, the present analytical solution is at least 3 orders of magnitude faster. Our calculation also indicates that, compared to the isotropic PREM model, the effect of anisotropy is small, with a relative error about 5 per cent for DLN h 12 at degree n  = 69 with source located at depth 100 km. We provide MATLAB code for computing the analytical DLNs in the supplementary materials section, along with a user manual. [ABSTRACT FROM AUTHOR]

We study the surface deformation of a moderate size M 5+ earthquake swarm-type activity which burst at the tip of the Biga peninsula (western Turkey) in early 2017. No previous M 5+earthquakes have been recorded in the instrumental period on land, however, offshore normal faults, have ruptured to produce strong (M  > 6) earthquakes. We use the Interferometric Synthetic Aperture Radar (InSAR) technique and exploit a data set of Sentinel-1 and ALOS-2 images, to successfully detect the surface deformation caused by three M 5+ events of the sequence and constrain their source models, furtherly strengthened by seismic waveform analysis. The sequence occurred at shallow depths (<12 km) and is related with the activation of a normal fault. Our geodetic inversions constrained the rupture distribution of the main events and led us to conclude that they belong to a single fault plane, striking N110°E and dipping ∼40° to the SW, compatible with the regional tectonics. Furthermore, the InSAR analysis revealed that no aseismic transients occurred during the Biga swarm. InSAR applications in seismic swarms are limited mainly due to the low displacement signal they produce, which is not always easily detectable by radar satellites. Another obstacle is the time frequency of radar satellites acquisitions that often does not allow the temporal isolation of distinct seismic events. However, here we present a study that exploits InSAR data to isolate seismic sources of a swarm and proposes its use for the understanding of shallow seismic swarms in a systematic manner. [ABSTRACT FROM AUTHOR]

Western Slovenia is part of an actively deforming region accommodating anticlockwise rotation of Adria and its continuous collision with Eurasia. The geometry of the active faulting system in this plate boundary is not well defined. In this study, detailed analysis of earthquake activity was performed with relocation of earthquakes in the period between 2006 and 2017. With inspection of the waveform data, slight temporal clustering of activity was observed. To increase the detection rate of microearthquakes we used a matched filter detection algorithm method. Templates of earthquakes were created and a database of continuous waveform data within the period 2006–2017 was investigated. As a result, high temporal correlation allowed us to identify swarms and earthquake sequences that affected the active fault system in the study region. Relocated seismicity allowed us to constrain the geometry of 5 nearly parallel faults, namely: Ravne, Idrija, Predjama, Selce and Raša faults. All these faults do have an expression in the geomorphology and reach a seismogenic depth of up to 20 km. Vertical and along strike extents of these active faults can favour earthquakes of moment magnitude equal to 7 or larger. The most recent large earthquake that occurred in this region is the 1511 earthquake with a magnitude 6.8. The leading fault in the system being the Idrija right-lateral strike-slip fault, experiences earthquake activity from 5 to 20 km on its northern segment, while on its southern segment no earthquake activity is detected over the decade of observations. We show that the interseismic loading on the southern segment of Idrija fault is likely unclamping the locked adjacent faults promoting the observed bursts of seismicity. Moreover, in 2009 the Predjama fault accommodated a sudden increase of the surface deformation at the extensometer accompanied by a simultaneous swarm activity at its seismogenic depth. This behaviour might correspond to velocity strengthening and weakening processes taking place at both the surface and depth terminations of a locked vertical fault. These processes can be driven by a slow-slip event on the deeper part of Idrija fault that would generate a temporary acceleration of the interseismic loading rate along with a change within the fluid circulation. [ABSTRACT FROM AUTHOR]

Knowledge of deformation at plate boundaries has been improved greatly by the development of observational techniques in space geodesy. However, most theoretical and numerical models of coseismic deformation have remained very simple and do not include realistic Earth structure. 3-D material heterogeneity and topography are often neglected because simple models are assumed to be sufficient and available tools cannot easily accommodate complex heterogeneity. In this study, we demonstrate the importance of 3-D heterogeneity using a spectral element method that incorporates topography and 3-D material properties. Using a parabolic hill model and a topographic model of the 2010 Maule earthquake, we show that topographic features can alter the shape of observed surface deformation patterns. We also estimate the coseismic surface deformation due to a model of the slip distribution of the 2015 Gorkha earthquake using realistic topography and 3-D elastic structure, and find that the presence of topography causes changes in the shape of observed surface displacement patterns, while material heterogeneities primarily affect the magnitude of observed displacements. Our results show that the inclusion of topography in particular can affect predictions of coseismic deformation modelling. [ABSTRACT FROM AUTHOR]

Magnitudes of differential stress in the lithosphere, especially in the crust, are still disputed. Earthquake-based stress drop estimates indicate median values <10 MPa whereas the lateral variation of gravitational potential energy per unit area (GPE) across significant relief indicates stress magnitudes of ca. 100 MPa in average across a 100 km thick lithosphere between the Indian lowland and the Tibetan plateau. These standard GPE -based stress estimates correspond to membrane stresses because they are associated with a deformation that is uniform with depth. We show here with new analytical results that lateral variations in GPE can also cause bending moments and related bending stresses of several hundreds of MPa. Furthermore, we perform 2-D thermomechanical numerical simulations (1) to evaluate estimates for membrane and bending stresses based on GPE variations, (2) to quantify minimum crustal stress magnitudes that are required to maintain the topographic relief between Indian lowland and Tibetan plateau for ca. 10 Ma and (3) to quantify the corresponding relative contribution of crustal strength to the total lithospheric strength. The numerical model includes viscoelastoplastic deformation, gravity and heat transfer. The model configuration is based on density fields from the CRUST1.0 data set and from a geophysically and petrologically constrained density model based on in situ field campaigns. The numerical results indicate that values of differential stress in the upper crust must be > ca. 180 MPa, corresponding to a friction angle of ca. 10° to maintain the topographic relief between lowland and plateau for >10 Ma. The relative contribution of crustal strength to total lithospheric strength varies considerably laterally. In the region between lowland and plateau and inside the plateau the depth-integrated crustal strength is approximately equal to the depth-integrated strength of the mantle lithosphere. Simple analytical formulae predicting the lateral variation of depth-integrated stresses agree with numerically calculated stress fields, which show both the accuracy of the numerical results and the applicability of simple, rheology-independent, analytical predictions to highly variable, rheology-dependent stress fields. Our results indicate that (1) crustal strength can be locally equal to mantle lithosphere strength and that (2) crustal stresses must be at least one order of magnitude larger than median stress drops in order to support the plateau relief over a duration of ca. 10 Ma. [ABSTRACT FROM AUTHOR]

Progressive deformation of upper mantle rocks via dislocation creep causes their constituent crystals to take on a non-random orientation distribution (crystal preferred orientation or CPO). The resulting anisotropy of the rock's elastic properties can be observed by seismic techniques, and provides a means to constrain mantle convective flow patterns. Existing methods for calculating the evolution of CPO in deforming olivine-dominated rocks rely on unwieldy discrete representations of the crystal orientation distribution in terms of a large number (103–104) of individual grains. Here we propose a more efficient method in which CPO is represented using three continuous analytical functions (structured basis functions or SBFs), each of which represents a virtual CPO produced by the action of just one of the three dominant slip systems of olivine. The SBFs are then combined using an appropriate weighting scheme to represent a realistic CPO that results from the simultaneous activity of all three slip systems. We assume that olivine CPO is a unique function of the finite strain experienced by the aggregate, which implies that the weights of the SBFs depend only on the two ratios of the lengths of the axes of the finite strain ellipsoid (FSE) and the two ratios of the strengths (critical resolved shear stresses) of the slip systems. Our preferred set of weighting coefficients is obtained by least-squares fitting of the SBF expansion to the predictions of a kinematic model (solved by the method of characteristics) in which the amplitudes of the crystallographic spins do not increase with strain. Calculation of CPO using this model is ≈107 times faster than full homogenization approaches such as the second-order self-consistent model, and the result fits the characteristics-based solution with a variance reduction |$\ge 88.6\hbox{ per cent}$| for equivalent strains up to 0.9. Finally, we propose a simple modification of the FSE that prevents the CPO from becoming singular at large strains. [ABSTRACT FROM AUTHOR]

The geodynamics of the Mediterranean comprises a transitional setting in which slab rollback and plate convergence compete to shape the region. In the central Mediterranean, where the balance of driving and resisting forces changes continuously and rapidly since the Miocene, both kinematic and seismo-tectonic observations display a strong variation in deformation style and, therefore possibly, lithospheric forces. We aim to understand the current kinematics in southern Italy and Sicily in terms of lithospheric forces that cause them. The strong regional variation of geodetic velocities appears to prohibit such simple explanation. We use mechanical models to quantify the deformation resulting from large-scale Africa–Eurasia convergence, ESE retreat of the Calabrian subduction zone, pull by the Aegean slab, and regional variations in gravitational potential energy (topography). A key model element is the resistance to slip on major regional fault zones. We show that geodetic velocities, seismicity and sense of slip on regional faults can be understood to result from lithospheric forces. Our most important new finding is that regional variations in resistive tractions are required to fit the observations, with notably very low tractions on the Calabrian subduction contact, and a buildup towards a significant earthquake in the Calabrian forearc. We also find that the Calabrian net slab pull force is strongly reduced (compared to the value possible in view of the slab's dimensions) and that trench suction tractions are negligible. Such very small contributions to the present-day force balance in the south-central Mediterranean suggest that the Calabrian arc is now further transitioning towards a setting dominated by Africa–Eurasia plate convergence, whereas during the past 30 Myr slab retreat continually was the dominant factor. [ABSTRACT FROM AUTHOR]

Monitoring surface deformation in real-time helps at planning and protecting infrastructures and populations, manages sensitive production (i.e. SEVESO-type) and mitigates long-term consequences of modifications implemented. We present RT-SHAKE, an algorithm developed to detect ground motions associated with landslides, subsurface collapses, subsidences, earthquakes or rock falls. RT-SHAKE detects first transient changes in individual GPS time-series before investigating for spatial correlation(s) of observations made at neighbouring GPS sites and eventually issues a motion warning. In order to assess our algorithm on fast (seconds to minute), large (from 1 cm to metres) and spatially consistent surface motions, we use the 1 Hz GEONET GNSS network data of the Tohoku-Oki M w9.0 2011 as a test scenario. We show that the delay of detection of seismic wave arrival by GPS records is of ∼10 s with respect to an identical analysis based on strong-motion data and this time delay depends on the level of the time-variable noise. Nevertheless, based on the analysis of the GPS network noise level and ground motion stochastic model, we show that RT-SHAKE can narrow the range of earthquake magnitude, by setting a lower threshold of detected earthquakes to M w6.5–7, if associated with a real-time automatic earthquake location system. [ABSTRACT FROM AUTHOR]