Tectonic Displacement Research Articles

Structure and Composition of Basement and Sedimentary Cover in the Southwestern Part of the Siljan Ring, Central Sweden: New Data from the C-C-1 Drill Core

The results from the geological and geophysical investigations of the Siljan Ring impact structure (central Sweden) have shown that the Paleozoic sedimentary succession and the Precambrian basement were strongly affected by complex deformational processes. Studies of a new drill core from the C-C-1 well provide valuable additional information necessary for the reconstruction of the geological setting in the southwestern part of the Siljan Ring. It was found that the contact between the basement and the sedimentary cover is tectonic, not normal sedimentary, in origin. The basement interval comprises Precambrian metavolcanic and metasedimentary rocks with a single mafic intrusion (gabbro-dolerite) in the upper part. The rocks have only been partially metamorphosed. The intercalation of calcareous mudstones, skeletal wackstones, and black shales in the sedimentary cover interval is not consistent with the regional lithostratigraphy scheme. Thus, more likely that the sedimentary sequence is not complete as a result of tectonic displacements, and a significant part of the Lower and Middle Ordovician succession is missing. The Post-Proterozoic tectonic reactivation and impact event also caused the formation of four types of fracture. The third type of fracture is accompanied by cataclastic zones and probably have an impact-related nature. In the highly fractured basement rocks, a dissolution along the second type of fracture has resulted in the development of open vugs. Open vugs and microporosity in cataclastic zones have been considered to be an effective storage space for hydrocarbons.

Seismic stratigraphy of the north-western Sea of Marmara shelf along the North Anatolian Fault system

The Ganos Fault, a part of the Northern strand of the North Anatolian Fault system, is an active- strike slip fault and divides the narrow NW shelf of the Sea of Marmara into two parts near the town of Gaziköy. This paper presents recently collected shallow high-resolution seismic data to discriminate the sedimentary successions, each characterized by distinctive stratigraphic patterns on both sides of the Ganos Fault. Two main units, namely U1 and U2, and three para-sequences (U1a, U1b and U1c) were identified, depending on their internal reflection patterns, accommodation depths as well as the presence of conformity and the unconformity surfaces. The thickness of Unit U1 reaches its maximum at the northern side of the Ganos Fault; it is much thinner to the south. The para-sequences of U1b and U1c have “progradation” and “aggradation to progradation” depositional characters, respectively. This probably implies fluvial deposition controlled by sea- level fluctuations. Unit U1b can only be observed at the northern side of the Ganos Fault, while Unit U1c at the southern side. Units U1a and U1b were deposited during the transgressive system tract, while Unit U1c was deposited during a sea-level fall and/or a lowstand phase marked by an erosional surface. The marine terraces in the study area are shallower than those along the northern shelf of the Sea of Marmara, possibly due to successive tectonic displacements along the Ganos Fault, which also controls the distribution and thickness of the parasequences identified in this study.

Preliminary numerical study on seismic response of ordinary long-span suspension bridges crossing active faults

The objective of this paper is to expediently expose the seismic performance pertinent to demand and capacity of general long-span suspension bridges crossing active faults. Firstly three dimensional finite element model of the ordinary long-span suspension bridge is established based on the powerful and attractive finite element software ANSYS. Secondly a series of appropriate fault ground motions with different target final permanent displacements (Tectonic displacements or ground offset) in the direction perpendicular to the fault plane are assumed and applied to the employed long-span suspension bridge. And then the Newmark method is utilized to solve the equation of motion of the long-span suspension bridge structure subjected to fault ground motions in the elastic range. Finally some important conclusions are drawn that the final permanent displacements in the direction perpendicular to the fault plane has significant influence on the seismic responses and demands of general long-span suspension bridges crossing active faults. And the resultant conclusions deliver explicitly and directly specifications and guidelines for seismic design of ordinary long-span suspension bridges across fault-rupture zones.

Overturned Eocene – Lower Pliocene alluvial stratum on the southern coast of Lake Baikal and its neotectonic significance

The study is focused on a section of sediments exposed on the right bank of Mishikha River, Russia. These sediments have a wide range of ages, from the Eocene to the Lower Pliocene. The stratigraphic subdivision of the section is based on the lithogeochemical data and X-ray phase analysis of the mineral compositions. The particle-size analysis shows the alluvial origin of the deposits. Their ages are constrained by spore-pollen spectra in three palynozones: I – Eocene – Oligocene, II – Early – Middle Miocene (subzone a – Tsuga, Picea in the lower part, and Quercus, Taxodiaceae, Momipites, Carya in the upper part; subzone b – Fagus, Quercus, Tsuga), and III – the Late Miocene – beginning of the Pliocene (subzone ν – Ulmus, Juglans, Carya; subzone g – Carya, Alnus). The section shows a combination of normal and overturned sedimentary layers. The tectonic displacement of the block with its flip was accompanied by the entry into contact of the unlithified Pliocene sediments with a rigid bed and the development of a landslide. The lower age limit of deformations is constrained from the youngest (beginning of the Pliocene) spore and pollen spectrum extracted from deformed layers. It is suggested that the overturned layers result from strike-slip deformations of the sediments at the beginning of the late orogenic stage of the Baikal rift development. The regional correlations of the sedimentary strata give grounds to conclude that the Mishikha section is characteristic of alluvial sedimentation that dominated at the eastern end of the Tankhoi tectonic step (Mishikha-Klyuevka paleovalley), in contrast to the Tankhoi block in the central part of the step, wherein a thick Lower Miocene stratum of swampy-oxbow sediments accumulated. The stratons of the Mishikha section correlate with sedimentary units detected by drilling in the Selenga delta at the central part of the South Baikal basin.

Role of landslides on the volume balance of the Nepal 2015 earthquake sequence

The 7.8 Mw earthquake that struck Nepal on April 25th, 2015 triggered over 21,000 landslides over an area of more than 25,000 km2. These landslides contributed to mass wasting, partially compensating the tectonic uplift by the earthquake. In this paper we quantify the volume balance resulting from the 2015 earthquake uplift (or subsidence) and landslide erosion. Starting from a new complete earthquake-induced landslide inventory, we calculated landslide volume by adopting different strategies for low-mobility and high-mobility landslides, considering also the potential supply of sediments to the drainage network. The results show that the contribution of earthquake-induced landslides to erosion is about one order of magnitude smaller than the vertical coseismic displacement. We found landslide volume values, due to the 2015 Nepal earthquake, ranging between 251 (− 15/ + 16) Mm3 up to 1503 (− 183/ + 210) Mm3 based on the adopted method, and a volume due to coseismic vertical displacement of 2134 (± 1269) Mm3 for the whole area. The volume balance of the 2015 Nepal earthquake is strongly dominated by tectonic displacement. We show that these estimates depend on several uncertainties. We identified and quantified uncertainties related to: (1) the choice of empirical volume-area scaling relationships and their parameters; (2) the completeness and quality of landslide inventory through comparison with available inventories; (3) the approach adopted for the assessment of elongated landslide volume; (4) the InSAR displacement data.

THE CENOZOIC CRUSTAL STRESS STATE OF MONGOLIA ACCORDING TO GEOLOGICAL AND STRUCTURAL DATA (REVIEW)

This article gives a chronological review of the main published research results concerning the Cenozoic crustal stress-strain state in Mongolia and adjacent territories. The studies commenced in the southern Baikal rift zone in the 1970s and were extended further southwards to cover mobile regions neighbouring the Siberian platform. Geological, structural and morphostructural data were collected and analysed to define the crustal stress types and spatial characteristics. The authors have consolidated their reconstructions of the crustal stress-strain state of Mongolia, which were based on tectonic fracturing data and displacements along fractures in fault zones active in the Cenozoic. We consolidated a database of reconstructed stress tensors, which now contains more than 750+ solutions. The Late Cenozoic stress field was mapped. The map shows domains differing in types of the paleostress state of the crust. The reconstructions were compared to our calculations of the present-day crustal stress state, which were based on earthquake focal mechanisms, and to calculations by other authors. At the Late Cenozoic and current stages, the maximum horizontal compression axis (S Hmax ) has varying orientations, from submeridional (Western Mongolia) to NE and ENE (Eastern Mongolia). The role of compression increases from the northern domains, where the reconstructions show shear and transtension, to the southern domains with prevailing transpression and compression. Regular changes occur in the stress state and rupture parageneses along the largest latitudinal faults, North Khangai and Dolinoozersky; such changes are related to left-lateral strike-slip faulting. We analysed the sequence of the occurrence of stress fields differing in types and spatial characteristics, and revealed the main regularities in the evolution of the crustal stress-strain state in time. In the Cenozoic history of crust deformation in Mongolia, we can distinguish several episodes that differ in the dominant impacts of various tectonic force sources or combinations of such impacts. At the beginning of the Cenozoic, tectonic structures developed mainly under the influence of the interaction of East Asia and the Pacific Plate, which was manifested in the southeastern domains of the study area. The long-term SE-trending asthenospheric flow caused crustal stretching, which initiated the formation of tectonic structures comprising the Baikal rift system. Starting from the Pliocene, crustal stretching took place in combination with NNE compression caused by the India–Eurasia convergence. As a result, shearing occurred along the large faults. At this background, the Khangai and Khentei uplifts (including crust extension zones at their crests) are large structures that developed due to the dynamic effect of local mantle anomalies.

Traces of Seismic Activity in Architectural Monuments in Derbent, Dagestan

As a result of studies conducted in the city of Derbent, it is established that its monumental fortification structures erected at the end of the 560s AD underwent numerous strong seismic shocks during their existence. The number and approximate ages of the strong seismic events that affected Derbent throughout its history are determined by the age of the deformed stone masonry. Five episodes of seismogenic destructions are identified. Clockwise rotations around the vertical axis along the interblock joints at the segment of the Sasanian masonry of the curtain wall and the tower of the southern wall of the citadel are noted to occur in late Sasanian and early Arab times (from the 6th century to the beginning of the 8th century). They were large in value and caused by seismic action. Numerous cases of rebuilding and relaying in the defensive walls of the city are also revealed. The building works during the Arab period, especially in the 8th century, described in written sources, were probably caused not only by the military confrontation of the Arab Caliphate and the Khazar Khaganate (the Arab–Khazar wars), which led to the damage of some segments in the defensive walls, but also by seismic activity. At the boundary of the Arab period and the rise of the Derbent Emirate (8th–11th centuries), a significant seismogenic rupture occurred in the upper part of the northern wall of the city between the citadel and the Dzharchi-Kapi Gate. The character of deformations corresponds to seismic shocks from the northwest. It is established that the repair and building works on reinforcing the defensive walls in the Sasanian citadel by later masonry (including rustication) were performed during the existence of the Derbent Emirate and in the Seljuk period. These works could be not only a consequence of destruction during the confrontation between the Derbent emirs with the city elite and the Shirvanshah state, which is described in the Tarikh al-Bab chronicle, but also as a result of the seismic destructions in the fortifications. The destruction and repair of the Juma Mosque occurred in the middle of 14th century (the 1360s). The seismic shocks, resulting in the destruction, arrived from the north. This destruction apparently caused the subsequent building activity in the Derbendi Shirvanshah period in the 15th century, especially during the reign of Khalilullah I (1417–1462), which was described in the building inscriptions and the typical features of Shirvan architecture (the Orta-Kapi Gate, Minaret Mosque). The women’s bathhouse, which was damaged by seismic oscillations propagating in submeridional direction, was repaired at the end of 18th century. The tectonic displacement of the basement blocks in the dugout of Peter the Great dates to this period as well.

Geomorphological and paleoseismological evidence of capable faulting in the Northern Apennines (Italy): Insights into active tectonics and seismic hazard of the Lunigiana basin

The Lunigiana basin is a NW trending tectonic depression located in the Northern Apennines of Italy. The basin is bounded by active normal faults, characterized by moderate to strong seismicity, with historical earthquakes reaching up to magnitude 6 (February 14, 1834) and instrumental events with magnitude up to 5.4 (June 21, 2013).Using the classical methodologies of active tectonic studies (i.e. remote sensing and digital topography analysis, field mapping, structural geology, and tectonic geomorphology), the main active faults have been mapped for their geometries, kinematics and evidence of activity. Both fault systems bounding the basin to the northeast and southwest have revealed clear evidence of fault capability, i.e., tectonic displacement of the latest Pleistocene-Holocene sediments and geomorphic markers. Then, the study focused on the southwestern Arzelato-Mulazzo-Tresana fault system (AMT), depicting its Late Quaternary tectono-sedimentary evolution and linkage relationship between important fault segments. The throw measurement of key geomorphic markers on the central fault of the AMT allowed to derive medium- (last 350–250 ka) and short-term (18–13 ka) vertical slip rates: 0.44–0.62 and 0.45–0.68 mm/yr, respectively.Lastly, an exploratory trench for paleoseismological analysis was dug across the fault scarp where geomorphic indicators of active surface faulting were most pronounced, just south of the village of Mulazzo.Its paleoseismological analysis revealed at least two events of surface faulting, with a cumulative throw in excess of 70 cm, occurred in historical times. Radiocarbon dating of the exposed alluvial and slope deposits allowed to bracket the surface rupture events to 14th to 19th century CE and after the 18th century CE. Typically, in the active extensional belt of the Apennines, the observed offsets correspond to earthquakes of magnitude >6.3, also confirmed by the empirical fault-scaling laws. This magnitude is, therefore, higher than that currently inferred for this seismic source (Mw 5.7÷6.0 in the DISS catalog) and is comparable to that experienced in the adjoining Garfagnana basin (Mw 6.5 on September 7, 1920). The findings of this research suggests a need to revise the seismic hazard of the Lunigiana basin, including greater consideration of the surface faulting hazard in the seismic microzonation studies.More generally, our findings provide new insight into the comprehension of a) the still unclear deep geometry (listric or not?) and the hierarchy of bounding faults, and b) provide new elements to decipher the puzzling seismotectonic relationships of Lunigiana with the adjoining basins on both its ends, characterized by quite different structural settings.

Shallow Subsurface Structure in the Hualien Basin and Relevance to the Damage Pattern and Fault Rupture during the 2018 Hualien Earthquake

ABSTRACTThe 2018 Mw 6.4 Hualien earthquake generated a large peak-to-peak velocity of over 2 m/s, with a period of 3 s at the south end of the Milun fault, which resulted in the collapse of five buildings. To investigate the shallow subsurface soil structure and evaluate possible effects on the ground motion and building damage, we performed microtremor measurements in the Hualien basin. Based on the velocity structure jointly inverted from both Rayleigh-wave dispersion curves and microtremor horizontal-to-vertical spectral ratio data, we found that the shallow subsurface structure generally deepens from west to east. Close to the Milun fault, the structure becomes shallower, which is consistent with faulting during the 2018 earthquake and the long-term tectonic displacement. There is no significant variation for the site conditions in the north–south direction that can explain the large peak ground velocity in the south. As a result of the dense measurements in the heavily damaged area, where three high-rise buildings totally collapsed, these locations have the average S-wave velocity of the upper 30 m (AVS30) values and are relatively high compared to the more distant area from the Meilun River. This is somewhat unusual, because lower AVS30 values indicating softer ground conditions are expected close to the river. We did not find any characteristic subsurface soil structure that may contribute to the building collapses. The large 3 s pulse was probably generated by source effects, rather than subsurface soil amplification.

New Insights for Understanding the Structural Deformation Style of the Strike-Slip Regime along the Wadi Shueib and Amman-Hallabat Structures in Jordan Based on Remote Sensing Data Analysis

This paper presents new findings that contribute to the understanding of the deformational style of the Wadi Shueib Structure (WSS) and the Amman-Halabat Structure (AHS) and their relationship with the regional tectonic regime of the Dead Sea Transform Fault (DSTF). Our research utilized Landsat-8 OLI imagery for the automatic extraction of lineaments, and our lineament mapping was facilitated by processing and digital image enhancement using principal component analysis (PCA). Our data revealed a relatively higher density of lineaments along the extension of the major faults of the WSS and AHS. However, a relatively lower density of lineaments was shown in areas covered by recent deposits. Two major lineament trends were observed (NNE-SSW and NW-SE) in addition to a minor one (NE-SW), and most of these lineaments are parallel to the orientation of the WSS and AHS. We offer the supposition that the DSTF has merged into the major faults of the WSS and AHS. We further suppose that these faults were reactivated as a restraining bend composed of active strike-slip fault branches that developed due to the NNW-SSE-trending Dead Sea transpressional stress field. Depending on the relationship between the direction of the WSF and AHF strands and the regional tectonic displacement along the DSTF, thrust components are present on faults with horsetail geometry, and these movements are accompanied by folding and uplifting. Thus, the major faults of the WSS and AHS represent a contractional horsetail geometry with associated folding and thrusting deformation.

Examples of Simultaneous Appearance of Deformations of the Earth’s Crust with Hour-long Period and Earthquakes

Records from broadband seismic station AAK located on Kirgiz Ridge, Tien Shan, are presented. The episodes of discrepancy between tectonic displacements and tidal oscillations are registered and compared with the timing of local earthquakes of energy class of K > 11 which occurred during 2000–2018. Several times the tectonic anomalies and the earthquakes occurred on the same day. Comparison with the meteorological data is carried out. The following physical mechanism of this phenomenon is suggested. The records of the broadband seismic station reflect gravity variations and vertical ground deformations caused by the displacements of air masses. This process initiates slow movements and earthquakes near unstable faults.

Deformation detection through the realization of reference frames

Abstract Hellenic Military Geographical Service (HMGS) has established and measured various networks in Greece which constitute the geodetic infrastructure of the country. One of them is the triangulation network consisting of about 26.000 pillars all over Greece. Classical geodetic measurements that held by the Hellenic Military Geographic Service (HMGS) through the years have been used after adjustment for the state reference frame which materializes the current Hellenic Geodetic Reference System of 1987 (HGRS87). The aforementioned Reference System (RS) is a static one and is in use since 1990. Through the years especially in the era of satellite navigation systems many Global Navigation Satellite System (GNSS) networks have been established. The latest such network materialized by HMGS is ongoing and covers until now more than the 2/3 of the country. It is referenced by International GNSS Service (IGS) permanent stations and consists a local densification IGS08 Reference Frame. Firstly, this gives the opportunity to calculate transformation parameters between the two systems and a statistical analysis of the residuals leads to intermediate conclusions. After that and in conjunction with existing past transformations, tectonic deformations and their directions are concluded. Moreover past GPS observations on the same pillars in compare to the newer ones give also a sense of tectonic displacements. Greece is one of the most tectonically active countries in Europe and the adoption of a modern kinematic or semi-kinematic geodetic datum is a necessity as it should incorporate a deformation model like 3d velocities on the reference frame realization. The detection of geodynamic changes is a continuous need and should be taken into consideration at each epoch.

Using the method of boundary integral equations in solving of the geological mapping problems with faults

Methods for creating digital grid models of geological surfaces based on approximation by bicubic B-splines are widely used in solving problems of mathematical geology. The variational-gridding method of geological mapping is a flexible and powerful tool that allows to use a large amount of source data for mapping, as well as a priori information about the spatial distribution of the mapped parameter. The smoothness of the basis functions does not allow the direct use of this effective method for mapping geological surfaces complicated by faults. This fact requires adaptation of the variational-gridding method to mapping with faults. The article discusses the technology of separate construction of the fault and plicative (smooth) components of the structural map. The fault component is represented in the form of an antiplane shear field of an elastic membrane described in a stationary two-dimensional Laplace equation. The fault network is modeled by narrow contours, at the boundaries of which the values of tectonic displacements are set. The Laplace equation is solved by the method of boundary integral equations, which allows one to calculate the displacement field at an arbitrary point of the mapped region, as well as to most accurately approximate the complex geometry of faults. Modeling of the plicative component of the structural surface takes place on the basis of a variational-gridding approach with adjustment of tectonic displacements. The approach combines the advantages of the spline approximation method and the accuracy of the semi-analytical solution for the fault component. It does not impose restrictions on the geometry of faults, and also allows for more efficient mathematical operations with structural surfaces, and their differential and integral characteristics.

Time-dependent modelling of a mountain front retreat due to a fold-to-fault controlled lateral spreading

The Eastern Betic Cordillera (Spain) exemplifies at the Sierra de Aitana an impressive case of lateral spreading which involves a mountain anticline ridge dislodged by a post-orogenic extensional tectonic. The deep trenches testify the ongoing gravitational slope deformation featured by a rock mass lateral spreading. The Sierra de Aitana anticline ridge was dislodged along a normal fault zone which generated a scarp with a kilometric lateral continuity. The lateral spreading originated because Eocene limestone overlays Eocene marls, inducing a continuous evolution over time due to visco-plastic deformations and isolating huge prismatic blocks. To quantify the creep-driven stress-strain effects at Sierra de Aitana ridge, the mechanical and rheological parameter values of the involved rocks were quantified through field and laboratory tests which allowed to constrain an engineering-geological model based on equivalent-continuum approach. A stress-strain numerical modelling aimed at back-analysing the sequential morpho-structural evolution of the Sierra de Aitana anticline-ridge was performed by reproducing the retreat of the ridge front and the time succession of faults (re)activations. A time-dependent solution was approached by assuming a time-dependent creep configuration of the modelling and a parametric solution adopted to calibrate rock mass rheology, also taking into account: tectonic displacement due to the Aitana normal fault system by simulating these element as interface in the discretized domain; influence of pre-existing joint sets by assuming a main anisotropy in the elasto-plastic failure solution; admissible variations of the regional stress-field related to compressive and extensional phases and a generalised visco-plastic behaviour representing the ductile marls. The numerical modelling outputs that the lateral spreading mainly evolved as a creep-driven process while stress-release induced by tectonic activity is not sufficient to justify the observed landforms. The obtained results highlight a not negligible role of inherited structural elements for regulating the time evolution of the ongoing gravitational process at the Sierra the Aitana anticline.

On spreading modes and magma supply at slow and ultraslow mid-ocean ridges

The ultraslow eastern Southwest Indian Ridge (SWIR) offers an opportunity to study the effect of magma supply on an ultraslow mid-ocean ridge starting from quasi-melt-free detachment-dominated spreading, and transitioning to volcanic spreading as one nears prominent axial volcanos. Detachments in the quasi-melt-free mode extend along-axis 60 to 95 km and have a lifetime of 0.6 to 1.5 myrs. They cut into their predecessor's footwall with an opposite polarity, causing part of the footwall lithosphere to experience further deformation, hydrothermal alteration, sparse magmatism and possibly thermal rejuvenation, in a hanging wall position. The accretion of the oceanic lithosphere in this context therefore occurs in two distinct stages over the lifetime of two successive detachment faults. We examine the transition from this nearly amagmatic detachment-dominated mode to the more common volcanic mode of spreading, showing that it occurs along-axis over distances ≤30 km. It involves a significant thinning of the axial lithosphere and a gradual decrease of the amount of tectonic displacement on faults, as the magmatic contribution to the divergence of the two plates increases. We develop a conceptual model of this transition, in which magma plays a double role: it fills the space between the diverging plates, thus reducing the need for displacement along faults, and it modifies the thermal state and the rheology of the plate boundary, affecting its thickness and its tectonic response to plate divergence. Based on a comparison of the ultraslow eastern SWIR, with the faster spreading Mid-Atlantic Ridge, we show that the activation of the volcanic, or of the detachment-dominated modes of spreading is connected with the volume of magma supplied per increment of plate separation, over a range of axial lithosphere thickness, and therefore over a range of the M ratio defined by (Buck et al., 2005) as the relative contribution of magma and faults to plate divergence (M is smaller, for a given volume of melt per increment of plate separation, if the plate is thicker). We therefore propose that M does not fully explain the variability in faulting styles observed at slow and ultraslow ridges and propose that rheological changes induced by magma also play a key role (melt itself is weak, hydrothermally altered gabbro-peridotite mixtures are weak, and melt heat sustains more vigorous hydrothermal circulation), resulting in contrasted potentials for strain localization, footwall flexure on faults and the development of detachment faults.

A Macroscale Hydrogeological Numerical Model of the Suio Hydrothermal System (Central Italy)

The complex behaviour of the Suio hydrothermal system (central Italy) and its potential exploitation as a renewable energy source are still unclear. To quantitatively evaluate the geothermal resource, the Suio hydrothermal system has been investigated with a hydrogeological numerical model that couples fluid flow, thermal convection, and transport of diluted species inside a hybrid continuum-discrete medium. The numerical model, calibrated and validated with available and new experimental data, unveiled the complex behaviour of the hydrothermal system. The normal tectonic displacements, the fracturing of the karst hydrostructure, and the aquitard distribution strongly influence the hydrothermal basin. In particular, a dual fluid circulation, sustained by steady-state thermal and pressure gradients, modulates the hydrothermalism at the several springs and wells. The presence of a medium to a low-temperature reservoir allows for potential exploitation of the geothermal resource.

A journey on plate tectonics sheds light on European crayfish phylogeography.

Crayfish can be used as model organisms in phylogeographic and divergence time studies if reliable calibrations are available. This study presents a comprehensive investigation into the phylogeography of the European stone crayfish (Austropotamobius torrentium) and includes samples from previously unstudied sites. Two mitochondrial markers were used to reveal evolutionary relationships among haplogroups throughout the species’ distributional range and to estimate the divergence time by employing both substitution rates and geological calibration methods. Our haplotype network reconstruction and phylogenetic analyses revealed the existence of a previously unknown haplogroup distributed in Romania's Apuseni Mountains. This haplogroup is closely related to others that are endemic in the Dinarides, despite their vast geographical separation (~600 km). The separation is best explained by the well‐dated tectonic displacement of the Tisza–Dacia microplate, which started in the Miocene (~16 Ma) and possibly carried part of the A. torrentium population to the current location of the Apuseni Mountains. This population may thus have been isolated from the Dinarides for a period of ca. 11 m.y. by marine and lacustrine phases of the Pannonian Basin. The inclusion of this geological event as a calibration point in divergence time analyses challenges currently accepted crayfish evolutionary time frames for the region, constraining the evolution of this area's crayfish to a much earlier date. We discuss why molecular clock calibrations previously employed to date European crayfish species divergences should therefore be reconsidered.

Strength of Strained Two‐Phase Mixtures: Application to Rapid Creep and Stress Amplification in Subduction Zone Mélange

AbstractAseismic creep may occur by distributed deformation in mélange shear zones comprising weak matrix and stronger clast materials. Slow slip events or steady tectonic displacement can be distributed over <100‐m thick shear zones if weak matrix controls bulk shear zone deformation. We use 2‐D numerical models to quantify the rheology of moderately strained (shear strain <1.75) mélange for various volumetric proportions of competent clasts. Mélange deformation with <50% clasts is matrix dominated and can accommodate steady creep. At higher clast proportions mélange viscosity increases more than tenfold after small strains, because strong clasts interact and form force chains. Clast shear stress is amplified above the imposed shear stress, by a factor of <14 where force chains develop. Slow slip events may occur due to a temporary absence of force chains, while localized regions of high shear stress generate coincident fracturing and potentially tremor events.

A landslide induced by the 2016 Kumamoto Earthquake adjacent to tectonic displacement - Generation mechanism and long-term monitoring

The 2016 Kumamoto earthquake, which occurred in the central Kyushu region of Southwest Japan, caused numerous landslides. Mass movements from the events damaged many structures in the Minami-aso village residential areas. We conducted geological and geophysical surveys and monitoring of landslides in one of these residential areas. It was found that the landslide occurred directly adjacent to the earthquake's main surface fault-ruptures. While the landslide slip surface formed within a layer of hydrothermally altered rhyolite at about 15 m depth for the downslope side, the slip surface for the upslope side formed within intact tephra layers due to the shallow dip of the rhyolite layer. Geophysical survey results show that fault ruptures in the tephra layers along with the tectonic faulting induce a decrease in soil strength. It was inferred that the landslide head scarp and upper slip surface formed inside this weak soil. Long-term monitoring of the landslide demonstrated that deformation continued for at least four months after the main slide. Although the initial and the major slide was activated by the earthquake, the landslide did not immediately stabilize and continued to pose risks after the main event. Toward the mitigation for earthquake-induced landslides, post seismic monitoring as well as the relationship between surface fault traces and landslide blocks should be taken into consideration.

GPU-based Implementation of Discrete Element Method for Simulation of the Geological Fault Geometry and Position

We present an algorithm for numerical simulation of the geological fault formation. The approach is based on the discrete elements method, which allows modeling of the deformations and structural discontinuity of the Upper part of the Earth crust. In the discrete elements method, the medium is represented as an combination of discrete particles which interact as elastic or viscoelastic bodies. Additionally, external potential forces, for example gravitational forces, may be introduced. At each time step the full set of forces acting at each particle is computed, after that the position of the particle is evaluated on the base of Newtonian mechanics. We implement the algorithm using CUDA technology to simulate single statistical realization of the model, whereas MPI is used to parallelize with respect to different statistical realizations. Obtained numerical results show that for low dip angles of the tectonic displacements relatively narrow faults form, whereas high dip angles of the tectonic displacements lead to a wide V-shaped deformation zones.