Rupture Model Research Articles

Revealing the slip distribution and seismogenic structure of the magnitude 6.2 Jishishan (Gansu) earthquake in 2023 through multi-source InSAR observations

ABSTRACT On 18 December 2023, Jishishan County in Gansu Province, China, was struck by an earthquake with magnitude 6.2, resulting in considerable loss of life and property. The seismogenic fault and rupture model of the event has sparked significant controversy, leading to diverse opinions among researchers. Some attribute the earthquake to the northern segment of the west-dipping Lajishan fault, whereas others consider an east-dipping blind fault on the Jishishan fault to be the seismogenic fault. This study integrates InSAR observations, aftershock relocation, and active tectonic analysis to determine the structure and rupture characteristics of seismogenic faults. Considering the active fault distribution, a listric fault model has emerged as the most plausible source of this earthquake event. The model of this study suggests a hidden branch of the southern segment of the Lajishan fault, whose reverse movement caused an uplift in the shallow crust. Supported by InSAR-derived coseismic deformation data and relocated aftershocks, this seismogenic model provides valuable insights into the regional tectonic framework. This study emphasizes the complexity of fault systems involved in earthquake genesis and contributes to a deeper understanding of seismic risk in the area.

TEAD1-Mediated Trans-Differentiation of Vascular Smooth Muscle Cells into Fibroblast-Like Cells Contributes to the Stabilization and Repair of Disrupted Atherosclerotic Plaques.

Atherosclerotic plaque rupture mainly contributes to acute coronary syndrome (ACS). Insufficient repair of these plaques leads to thrombosis and subsequent ACS. Central to this process is the modulation of vascular smooth muscle cells (VSMCs) phenotypes, emphasizing their pivotal role in atherosclerotic plaque stability and healing post-disruption. Here, an expansion of FSP1+ cells in a tandem stenosis (TS) model of atherosclerotic mice is unveiled, predominantly originating from VSMCs through single-cell RNA sequencing (scRNA-seq) analyses and VSMC lineage tracing studies. Further investigation identified TEA domain transcription factor 1 (TEAD1) as the key transcription factor driving the trans-differentiation of VSMCs into fibroblast-like cells. In vivo experiments using a TS model of plaque rupture demonstrated that TEAD1 played a crucial role in promoting plaque stability and healing post-rupture through pharmacological or TEAD1-AAV treatments. Mechanistically, it is found that TEAD1 promoted the expression of fibroblast markers through the Wnt4/β-Catenin pathway, facilitating the trans-differentiation. Thus, this study illustrated that TEAD1 played a critical role in promoting the trans-differentiation of VSMCs into fibroblast-like cells and subsequent extracellular matrix production through the Wnt4/β-Catenin pathway. Consequently, this process enhanced the healing mechanisms following plaque rupture, elucidating potential therapeutic avenues for managing atherosclerotic instability.

Earthquake Fault Rupture Modeling and Ground-Motion Simulations for the Southwest Iceland Transform Zone Using CyberShake

ABSTRACT CyberShake is a high-performance computing workflow for kinematic fault-rupture and earthquake ground-motion simulation developed by the Statewide California Earthquake Center to facilitate physics-based probabilistic seismic hazard assessment (PSHA). CyberShake exploits seismic reciprocity for wave propagation by computing strain green tensors along fault planes, which in turn are convolved with rupture models to generate surface seismograms. Combined with a faultwide hypocentral variation of each simulated rupture, this procedure allows for generating ground-motion synthetics that account for realistic source variability. This study validates the platform’s kinematic modeling of physics-based seismic wave propagation simulations in Southwest Iceland as the first step toward migrating CyberShake from its original study region in California. Specifically, we have implemented CyberShake workflows to model 2103 fault ruptures and simulate the corresponding two horizontal components of ground-motion velocity on a 5 km grid of 625 stations in Southwest Iceland. A 500-yr-long earthquake rupture forecast consisting of 223 hypothetical finite-fault sources of Mw 5–7 was generated using a physics-based model of the bookshelf fault system of the Southwest Iceland transform zone. For each station, every reciprocal simulation uses 0–1 Hz Gaussian point sources polarized along two horizontal grid directions. Comparison of the results in the form of rotation-invariant synthetic pseudoacceleration spectral response values at 3, 4, and 5 s periods are in good agreement with the Icelandic strong motion data set and a suite of empirical Bayesian ground-motion prediction equations (GMPEs). The vast majority of the physics-based simulations fall within one standard deviation of the mean GMPE predictions, previously estimated for the area. At large magnitudes for which no data exist in Iceland, the synthetic data set may play an important role in constraining GMPEs for future applications. Our results comprise the first step toward comprehensive and physics-based PSHA for Southwest Iceland.

Clinical-Radiomics Nomogram Model Based on CT Angiography for Prediction of Intracranial Aneurysm Rupture: A Multicenter Study.

Risk estimation of intracranial aneurysm rupture is critical in determining treatment strategy. There is a scarcity of multicenter studies on the predictive power of clinical-radiomics models for aneurysm rupture. This study aims to develop a clinical-radiomics model and explore its additional value in the discrimination of aneurysm rupture. A total of 516 aneurysms, including 273 (52.9%) with ruptured aneurysms, were retrospectively enrolled from four hospitals between January 2019 and August 2020. Relevant clinical features were collected, and radiomic characteristics associated with aneurysm were extracted. Subsequently, three models, including a clinical model, a radiomics model, and a clinical-radiomics model were constructed using multivariate logistic regression analysis to effectively classify aneurysm rupture. The performance of models was analyzed through operating characteristic curves, decision curve, and calibration curves analysis. Different models' comparison used DeLong tests. To offer an understandable and intuitive scoring system for assessing rupture risk, we developed a comprehensive nomogram based on the developed model. Three clinical risk factors and fourteen radiomics features were explored to establish three models. The area under the receiver operating curve (AUC) for the radiomics model was 0.775 (95% CI,0.719-0.830), 0.752 (95% CI,0.663-0.841), 0.747 (95% CI,0.658-0.835) in the training, internal and external test datasets, respectively. The AUC for clinical model was 0.802 (95% CI, 0.749-0.854), 0.736 (95% CI, 0.644-0.828), 0.789 (95% CI, 0.709-0.870) in these three sets, respectively. The clinical-radiomics model showed an AUC of 0.880 (95% CI,0.840-0.920), 0.807 (95% CI,0.728-0.887), 0.815 (95% CI,0.740-0.891) in three datasets respectively. Compared with the radiomics and clinical models, the clinical-radiomics model demonstrated better diagnostic performance (DeLong' test P < 0.05). The clinical-radiomics model represents a promising approach for predicting rupture of intracranial aneurysms.

Spatiotemporal Clustering of Large Earthquakes Along the Central‐Eastern Sections of the Altyn Tagh Fault, China

AbstractThe understanding of the spatial‐temporal distribution of past earthquakes is essential to assess the event recurrence behavior and to estimate the size of potential earthquakes along active strike‐slip fault systems. However, the scarcity of paleoseismic data remains a major hurdle in this endeavor. This is the case of the longest strike‐slip fault in Asia, the Altyn Tagh Fault (ATF). We documented six very likely large earthquakes that potentially ruptured the Aksay section of the ATF. Employing a Bayesian approach, we present modeled date ranges of 6339–5220 BC, 5296–4563 BC, 3026–2677 BC, 1324–808 BC, 314–632 AD, and 915–1300 AD. The mean recurrence time is 1,329 ± 588 years with a coefficient of variation (COV) of ∼0.44. In the same fault section, 90 horizontal offsets record an average coseismic slip of 5.1 ± 1.4 m for the last event and suggest four older earthquakes plausibly with a similar slip distribution. Although at the local‐scale the COV indicates quasi‐periodic rupture behavior, the individual interevent times exhibit significant irregularity, a pattern also observed in adjacent fault sections (Xorxoli, Annanba and Tashi sections). We found that such irregularities are a natural consequence of long‐term fault interactions, which allow for synchronized ruptures along the northern and southern strands of the central‐eastern ATF. Our rupture model highlights bursty periods of seismic activity with mean interevent times of 475 ± 108 years separated by long‐lull periods of 1.1–1.6 kyr. Based on this temporal organization and considering the 401‐year elapsed time since the most recent event on the Xorxoli section, there exists a possibility of a forthcoming large earthquake occurring within the next century.

Exploring Fission Dynamics through Fission Fragment Spectroscopy

Abstract Fission Fragment Spectroscopy (FFS) measurement technique has been utilized to measure the relative isotopic yield distributions of the even-even correlated fission fragment nuclei produced from 236U* at two different values of excitation energies (E ex ). The fissioning compound nucleus, 236U* was populated at E ex = 6.5 and 21.5 MeV through two different sets of experiments by using the reactions, 235U(n th ,f) and 232Th(α,f), respectively. The experimentally measured relative isotopic yield distributions were further used to obtain the corresponding mass yield distributions. The extracted mass yield distributions of 236U* indicates the transition from asymmetric distribution to bimodal distribution with increase in the values of E ex from 6.5 to 21.5 MeV. The experimental mass yield distribution pattern for the reaction, 232Th(α,f) at E ex = 21.5 MeV could be fitted well through three Gaussian fits as per the prescription based on the Multimodal Random-Neck Rupture Model (MM-RNRM).

CD1530, selective RARγ agonist, facilitates Achilles tendon healing by modulating the healing environment including less chondrification in a mouse model.

Heterotopic ossification (HO) in Achilles tendon often arises due to endochondral ossification during the healing process following trauma. Retinoic acid receptor γ (RARγ) plays a critical role in this phenomenon. This study aims to elucidate the therapeutic effects of CD1530, an RARγ selective agonist, along with the contributing cells, in Achilles tendon healing, utilizing a cell lineage tracing system. Local injection of CD1530 facilitated histological tendon healing by inhibiting chondrification in a mouse Achilles rupture model. Resident Scleraxis (Scx)+ cells in Achilles tendon were not found to be actively involved in HO or tendon healing following injury. Instead, these processes were primarily driven by tendon stem/progenitor cells (TSPC)-like cells. Furthermore, an in vitro assay revealed that CD1530 attenuated inflammation in injured Achilles tendon-derived tendon fibroblasts (iATF) and inhibited the chondrogenesis of iATF. This dual effect suggests the potential of CD1530 in effectively modulating the healing environment during tendon healing. Together, the present study demonstrated that the local administration of CD1530 accelerated tendon healing by modulating the healing environment, including reducing chondrification via targeting TSPC-like cells in a mouse Achilles tendon rupture model. These results suggest that CD1530 may have the potential to be a novel tendon therapy that offers benefits via the inhibition of chondrogenesis.

Characterizing the spatial correlation of coseismic slip distributions: a data driven Bayesian approach

SUMMARY The spatial correlation of coseismic slip is a necessary input for generating stochastic seismic rupture models, which are commonly used in seismic and tsunami hazard assessments. To date, the spatial correlation of individual earthquakes is characterized using finite fault models by finding the combination of parameters of a von Kármán autocorrelation function that best fits the observed autocorrelation function of the finite fault model. However, because a priori spatial correlation conditions (i.e. not in the data) are generally applied in finite fault model generation, the results obtained using this method may be biased. Additionally, robust uncertainty estimates for spatial correlations of coseismic slip are generally not performed. Considering these limitations in the classic method, here, a method is developed based on a Bayesian formulation of Finite Fault Inversion (FFI) with positivity constraints. This method allows for characterizing the spatial correlation of coseismic slip and its uncertainties for an earthquake by using samples of coseismic slip from a posterior probability density function (PDF). Furthermore, a Bayesian model selection criterion called Akaike Bayesian Information Criterion (ABIC) is applied to objectively choose between different prior spatial correlation schemes before computing the posterior, to reduce subjectivity due to this prior condition. The ABIC is calculated using an approximate analytical expression of Bayesian evidence. The method is applied to simulated P waves, demonstrating that model selection allows for objectively estimating the most suitable prior spatial correlation scheme in FFI. Additionally, the target spatial correlation of coseismic slip is accurately recovered using samples from the posterior PDF, as well as their uncertainties. Moreover, in the simulated experiment, it is shown that a non-robust choice of the prior spatial correlation scheme can significantly bias the estimated spatial correlations of coseismic slip. We apply our method to observed P waves from the 2015, Illapel earthquake ($M_{\rm w} = 8.3$), finding that the spatial correlation of coseismic slip of this earthquake is better described by a von Kármán ACF, with mean correlation lengths of around 47 km and Hurst parameter of 0.58. We conclude that using our method reduces biases associated with prior spatial correlation conditions and allows for robust estimation of spatial correlations of coseismic slip and their uncertainties.

Source mechanism of the 2023 Ms 5.5 earthquake in Subei, Gansu Province revealed by relocated aftershocks and InSAR: complement to the ‘shallow slip deficit’ of the eastern boundary of the Altyn Tagh fault

The Ms 5.5 earthquake struck on 24 October 2023, in Subei County, Gansu Province, China, occurring along the eastern segment of the Altyn Tagh fault. It raises the question of whether this earthquake is linked to the ongoing shortening slip rate along this segment or triggered by other seismic events. Analyzing the fault geometry of the Subei earthquake and understanding the significance of the weakening activity rate for seismic hazards in neighboring regions is crucial. The surface deformation from small- and medium-sized earthquakes (magnitudes less than Mw5.5) is often subtle, and the coseismic deformation detected by interferometric synthetic aperture radar (InSAR) is vulnerable to atmospheric disturbances, leading to significant measurement errors. Moreover, inaccuracies in the regional crustal velocity structure can cause errors in earthquake localization based on seismic data. These challenges complicate the establishment of a rupture model for seismogenic faults and hinder the inversion of fault slip models. To overcome these limitations, we employed the time-series InSAR stacking method and aftershock relocation to determine the fault geometry of the Subei earthquake. A two-step inversion method was utilized to ascertain both the fault geometry and slip distribution. Our modeling indicates that the 2023 Subei earthquake had a thrust mechanism with a component of strike-slip. The rupture did not reach the surface, with the maximum fault slip measuring 0.45 m at a depth of 2.5–3.5 km. The fault dips westward, and the moment magnitude is calculated at 5.4. This earthquake is associated with the ongoing weakening of the left-lateral strike-slip rupture along the Altyn Tagh fault in the Subei region. Furthermore, retrograde thrust tectonics significantly contribute to the absorption of accumulated stress during this process.Our findings highlight the potential of utilizing time-series InSAR images to enhance earthquake catalogs with geodetic observations, offering valuable data for further studies of the earthquake cycle and active tectonics. This approach is also applicable in other tectonically active regions, enhancing understanding of seismic hazards and risk assessment.

Exploring the Dynamic Interactions Between the Southern San Andreas Fault and a Normal Fault Under the Salton Sea

AbstractWe investigate the dynamic interactions between the Southern San Andreas Fault (SSAF) and a proximal normal fault (NF) beneath the Salton Sea in southern California. The NF, positioned near the SSAF terminus at Bombay Beach, exhibits 11–15 displacement events across 14 stratigraphic sequences, with a range of 0.2–1.4 m of vertical offset since ∼2–3 ka. Notably, four of these events may align temporally with SSAF earthquakes, raising questions about the possible interplay between the two faults. Utilizing dynamic rupture models, we analyze the coseismic interactions between the SSAF and NF, addressing under what conditions the SSAF induces slip on the NF. Our findings reveal that a suite of SSAF ruptures, particularly those propagating from north to south, can trigger slip on the normal fault and replicate observed vertical offsets. If the SSAF extends beneath the Salton Sea, earthquakes originating south of the NF intersection are less likely to trigger normal fault slip, although we cannot exclude this possibility. Some SSAF ruptures do not trigger discernible slip on the NF, rendering such events undetectable in the stratigraphic record. Our research contributes toward discussions regarding the seismic hazard in southern California, shedding light on the interplay between the SSAF and NF.

Dynamic Rupture Modeling of Large Earthquake Scenarios at the Hellenic Arc Toward Physics‐Based Seismic and Tsunami Hazard Assessment

AbstractThe Mediterranean Hellenic Arc subduction zone (HASZ) has generated several 8 earthquakes and tsunamis. Seismic‐probabilistic tsunami hazard assessment typically utilizes uniform or stochastic earthquake models, which may not represent dynamic rupture and tsunami generation complexity. We present an ensemble of ten 3D dynamic rupture earthquake scenarios for the HASZ, utilizing a realistic slab geometry. Our simplest models use uniform along‐arc pre‐stresses or a single circular initial stress asperity. We then introduce progressively more complex models varying initial shear stress along‐arc, multiple asperities based on scale‐dependent critical slip weakening distance, and a most complex model blending all aforementioned heterogeneities. Thereby, regional initial conditions are constrained without relying on detailed geodetic locking models. Varying epicentral locations in the simplest, homogeneous model leads to different rupture speeds and moment magnitudes. We observe dynamic fault slip penetrating the shallow slip‐strengthening region and affecting seafloor uplift. Off‐fault plastic deformation can double vertical seafloor uplift. A single‐asperity model generates a 8 scenario resembling the 1303 Crete earthquake. Using along‐strike varying initial stresses results in 8.0–8.5 dynamic rupture scenarios with diverse slip rates and uplift patterns. The model with the most heterogeneous initial conditions yields a 7.5 scenario. Dynamic rupture complexity in prestress and fracture energy tends to lower earthquake magnitude but enhances tsunamigenic displacements. Our results offer insights into the dynamics of potential large Hellenic Arc megathrust earthquakes and may inform future physics‐based joint seismic and tsunami hazard assessments.

The Postseismic Deformation of the 6 July 2019 Mw 7.1 Ridgecrest Earthquake from Burst Overlap Interferometry, InSAR, and GNSS

Abstract The July 2019 Ridgecrest earthquake sequence consists of an Mw 6.4 foreshock and an Mw 7.1 mainshock, which ruptured a complex set of orthogonal faults in the eastern California shear zone. We measure the co- and postseismic deformation associated with this sequence using the Burst Overlap Interferometry (BOI) technique in addition to the commonly used Interferometric Synthetic Aperture Radar (InSAR). The BOI technique, which provides displacement in the satellite’s along-track direction, offers important information on the postseismic deformation that cannot be measured by traditional InSAR and is only sparsely measured by the Global Navigation Satellite System networks. The BOI data reveal up to 4 cm displacement in the along-track direction, 10 km north of the northern tip of the seismic rupture, and up to 3 cm displacement along the coseismically active faults. These results rule out the possibility of significant shallow afterslip near the mainshock hypocenter, suggesting that the previously suggested poroelastic rebound is likely to be the cause for the significant postseismic line of sight deformation near the mainshock rupture. Based on the aftershocks’ moment tensor distribution, surface rupture, and simple forward modeling, we propose that the postseismic deformation north of the Ridgecrest rupture is caused by an aseismic slip along a north-trending normal fault of the Ridgecrest rupture that was induced by the Ridgecrest earthquake. Furthermore, we observed that both deformation and seismic activity decay slower over time as the distance from the Coso geothermal area increases. This decay is influenced by the mechanical properties of the crust, which are affected by the increased heat flow at Coso and thus suppress deformation and seismicity, ultimately controlling their temporal evolution.

Effects of Dip Angle on Rupture Propagation Along Branch Fault Systems

ABSTRACT An important consideration in assessing seismic hazards is determining what is likely to happen when an earthquake rupture encounters a geometric complexity such as a branch fault. Previous studies showed parameters such as branch angle, stress orientation, and stress heterogeneity as key factors in the self-determined rupture path on branch faults. However, most of these studies were conducted in 2D or 3D with perfectly vertical faults. Therefore, in this study, we investigate the effects of dipping angle on rupture propagation along a branch fault system. We construct 3D finite-element meshes where we vary the dip angles (nine geometries in total) of the main and secondary faults, the stressing angle (Ψ=20°, 40°, and 65°), and the hypocenter location with nucleation on both the main and secondary segments. We find that for Ψ=40°, a rupture on the main fault is most likely to propagate across the branch intersection when the secondary fault is dipping. In addition, for Ψ=65°, a rupture on the secondary fault is most likely to propagate to the main fault when the secondary fault is shallowly dipping. This is caused by a fast rupture speed on the secondary fault and the dynamic stress effect that develops with the interaction of the free surface and the dipping secondary fault. These results indicate that dip angle is an important parameter in the determination of rupture path on branch fault systems, with potentially significant impact for seismic hazard, and should be considered in future dynamic rupture modeling studies.

Source Modeling of Deep Plate Boundary 2021 Miyagi Earthquake (Mw7.0) Employing Modified Semi-Empirical Technique with Site Effects: A Step Forward Towards Hazard Mitigation

ABSTRACT The source modeling of 2021 Miyagi earthquake (Mw7.0) offers a significant tool to understand the initial evaluation of M9 earthquake cycle in this region. This article seeks to simulate the 2021 Miyagi earthquake (Mw7.0) using the modified semi-empirical technique (MSET) after incorporating site effects determined using the Horizontal to Vertical spectral ratio technique. We propose the best-fitting source model of this earthquake from a spectrum of rupture model parameters using MSET. We believe that this effort is the first to use MSET to model this earthquake and will provide significant contribution for seismic hazard assessment of the Miyagi region.

Kinematic rupture modeling of broadband ground motion from the 2022 MS6.9 Menyuan earthquake

We propose a novel kinematic rupture modeling procedure for synthesizing broadband ground motions derived from the frequency-wavenumber integration algorithm. This procedure addresses two key issues in characterizing the rupture processes relevant to broadband seismic radiation: an accurate Green's function and a well-constrained kinematic source model. For the first issue, we derive the theoretical Green's function based on an improved dynamic stiffness matrix approach that effectively handles wave propagation in a 1D crustal velocity structure across a broad frequency band. For the second issue, we generate the hybrid source model that combines asperity slip and random slip over the fault plane to effectively implement constraints on the radiated energy during the whole rupture process. The accuracy and effectiveness of the proposed methodology are verified by comparing with the surface acceleration traces and Fourier spectra calculated by spectral element method. With the hybrid source model and crustal velocity structure applicable to the target area, the broadband (0–10 Hz) ground motion of the 2022 MS6.9 Menyuan earthquake is synthesized. The amplitude, duration, and frequency content of the synthetic motions are systematically compared with those of the available observed records and ground motion attenuation relationships, as well as the spatial distribution characteristics of the near-field ground motions from the earthquake scenarios are presented. In conclusion, the case study of the Menyuan MS6.9 earthquake demonstrates that the presented modeling procedure can estimate broadband ground motions rapidly and reliably from a physics-based kinematic rupture perspective.

Non‐Typical Supershear Rupture: Fault Heterogeneity and Segmentation Govern Unilateral Supershear and Cascading Multi‐Fault Rupture in the 2021 Mw ${M}_{w}$7.4 Maduo Earthquake

AbstractPrevious geodetic and teleseismic observations of the 2021 7.4 Maduo earthquake imply surprising but difficult‐to‐constrain complexity, including rupture across multiple fault segments and supershear rupture. Here, we present an integrated analysis of multi‐fault 3D dynamic rupture models, high‐resolution optical correlation analysis, and joint optical‐InSAR slip inversion. Our preferred model, validated by the teleseismic multi‐peak moment rate release, includes unilateral eastward double‐onset supershear speeds and cascading rupture dynamically triggering two adjacent fault branches. We propose that pronounced along‐strike variation in fracture energy, complex fault geometries, and multi‐scale variable prestress drives this event's complex rupture dynamics. We illustrate how supershear transition has signatures in modeled and observed off‐fault deformation. Our study opens new avenues to combine observations and models to better understand complex earthquake dynamics, including local and potentially repeating supershear episodes across immature faults or under heterogeneous stress and strength conditions, which are potentially not unusual.

Shallow Reverse Moderate Earthquakes in the Weiyuan Shale Gas Field, Sichuan Basin, China, Related to Hydraulic Fracturing

Abstract The Weiyuan shale gas field in the stable southern Sichuan basin, China, has experienced increasing seismicity since systematic hydraulic fracturing (HF) operations in 2015. Three Mw≥4.4 shallow earthquakes occurred in the Weiyuan area between September 2019 and February 2020, yet their seismogenic faults, rupture models, and relationship with HF are unknown. In this study, we first obtain the high-resolution coseismic deformation fields of these three events and then invert their slip distribution. The result shows that all three events are shallow high-dip reverse events under the contractional Weiyuan anticline environment with peak slips of 158, 68, and 34 cm and at depths of 4, 3, and 1.6 km, respectively. The spatial relationship between seismogenic faults, horizontal wells, as well as geological data reveals that pore-pressure diffusion due to the HF may be the main mechanism of the 8 September 2019 and the 18 December 2019 events, whereas the 16 February 2020 event may be attributed to the poroelastic stress perturbation caused by the HF. Our study highlights that HF activities and regional geological characteristics jointly influence the properties of earthquakes in the Sichuan basin.

Rupture direction of paleoearthquakes on the Alpine Fault, New Zealand, as recorded by curved slickenlines

Abstract We observed and further exhumed curved slickenlines on fault planes associated with paleo-surface rupture of the Alpine Fault, New Zealand's ∼30 mm/yr continental transform plate boundary. Dynamic rupture modeling indicates that the geometry of such curvature provides a record of past earthquake rupture directions. We focused our efforts on three sites that span a region known to variably halt or allow passage of past earthquakes (an “earthquake gate”) to contribute rupture direction constraints to the fault's spatiotemporally rich paleoseismic record. At Hokuri Creek and Martyr River, we observed both convex-up and convex-down curved slickenlines on and adjacent to principal slip surfaces, indicating past ruptures from both the northeast and southwest of these locations. At Martyr River, relationships suggest that the most recent event (inferred to correlate to 1717 CE) ruptured from the southwest. Our results demonstrate the utility of curved slickenlines as a valuable new paleoseismological tool for determining past rupture directions, applicable to surface-rupturing faults globally.

Panda Rope Bridge Technique promoted Achilles tendon regeneration in a novel rat tendon defect model.

This study aimed to determine whether the Achilles tendon tissue can undergo the pathological process of Achilles tendon regeneration after the Panda Rope Bridge Technique (PRBT). Rats (n = 120) that operated with Achilles tendon rupture were divided into three treatment groups: Defect group (D group), PRBT group and Defect + Fix group (DF group). The D group represented natural healing with no treatment, the PRBT group represented healing receiving PRBT treatment and the DF group represented healing through conservative treatment by ankle fixation. The morphological, histological and biomechanical properties of the defective Achilles tendon were assessed at 7, 10, 12, 14, 28 and 56 days postoperatively. Compared to that observed in the other two groups, defected rat Achilles tendons that underwent PRBT recruited more cells earlier, eventually forming mature tendons, as revealed by histological analysis. PRBT also enabled defected tendons to regain stronger mechanical properties, thereby improving the prognosis. This improvement may be related to the earlier polarization of macrophages. By establishing and using a novel surgical model of Achilles tendon rupture in rats, most injured Achilles tendons can regenerate and regain normal histological properties, whereas tendons with other interventions formed fibrotic scar tissue. The strong regenerative capacity of tendon tissue enabled us to describe the pathological process of tendon regeneration after PRBT surgery in detail, which would aid in the treatment of tendon injuries. PRBT promotes Achilles tendon regeneration and has the potential to become a standard treatment. Not applicable.

Comparative transcriptomic analysis of articular cartilage of post-traumatic osteoarthritis models.

Animal models of post-traumatic osteoarthritis (PTOA) recapitulate the pathological changes observed in human PTOA. Here, skeletally mature C57Bl6 mice were subjected to either rapid-onset non-surgical mechanical rupture of the anterior cruciate ligament (ACL) or to surgical destabilisation of the medial meniscus (DMM). Transcriptome profiling of micro-dissected cartilage at day7 or day42 following ACL or DMM procedure, respectively, showed that the two models were comparable and highly correlative. Gene ontology (GO) enrichment analysis identified similarly enriched pathways that were overrepresented by anabolic terms. To address the transcriptome changes more completely in the ACL model, we also performed small RNA sequencing, describing the first microRNA profile of this model. miR-199-5p was amongst the most abundant, yet differentially expressed, microRNAs, and its inhibition in primary human chondrocytes led to a transcriptome response that was comparable to that observed in both human 'OA damaged vs intact cartilage' and murine DMM cartilage datasets. We also experimentally verified CELSR1, GIT1, ECE1 and SOS2 as novel miR-199-5p targets. Together, these data support the use of the ACL rupture model as a non-invasive companion to the DMM model.