The Market Weighton High in the 21st century – new understandings of a long standing problem

The sequence of Jurassic and Cretaceous strata laid down in north Lincolnshire is traced northwards onto the Market Weighton High, and compared with the equivalent Cleveland Basin sequence as it is traced south onto the High. Understanding of the manner of operation of the High has long been hindered by the amount of erosion of strata over the High, particularly in the Early Cretaceous. Facies analysis of the remaining strata present on either side of the High is used to estimate the thickness and the facies of strata originally deposited over the High. Although there were periodic uplifts of the High, leading to increasing omissions of strata approaching the High, there is little evidence that during the Jurassic the High was emergent for significant periods. Most strata were either originally present at Market Weighton in marginal marine facies, or continued across the High without interruption. The operation of the High in the light of modern understanding of crustal stresses during the Mesozoic is reviewed.

Episodic fluid pressure cycling controls earthquake sequences on subduction megathrusts

<p>A major goal in earthquake physics is to derive a constitutive framework for fault slip that captures the dependence of friction on lithology, sliding velocity, temperature, and pore fluid pressure. Here, we present a newly-developed two-phase flow numerical model — which couples solid rock deformation and pervasive fluid flow — to show how crustal stresses and fluid pressures within subducting megathrust evolve before and during slow slip and fast events. This unified 2D numerical framework couples inertial mechanical deformation and fluid flow by using finite difference methods, marker-in-cell technique, and poro-visco-elasto-plastic rheology. An adaptive time stepping allows the correct resolution of both long- and short-time scales, ranging from years to milliseconds during the dynamic propagation of dynamic rupture.</p><p>We investigate how permeability and its spatial distribution control the interseismic coupling along the megathrust interface, the interplay between seismic and aseismic slip, and the nucleation of large earthquakes. While a constant permeability leads to more regular seismic cycles, a depth dependent permeability contributes substantially to the development of two distinct megathrust zones: a shallow, locked seismogenic zone and a deep, narrow aseismic segment characterized by slow-slip events. Furthermore, we show that without requiring any specific friction law, our models reveal that permeability, episodic stress transfer and fluid pressure cycling control the predominant slip mode along the subduction megathrust. Furthermore, we analyze how rate dependent strength and dilatation affect rupture propagation and arrest. Our preliminary results show that fluid-solid poro-visco-elasto-plastic coupling behaves similarly to rate- and state-dependent friction. In this context, fluid pressure plays the role of state parameter whose time evolution is governed by: (i) the short-term elasto-plastic collapse of pores inside faults during the rupture (coseismic self-pressurization of faults) and (ii) the long-term pore-pressure diffusion from the faults into surrounding rocks (post- and interseismic relaxation of fluid pressure). This newly-developed numerical framework contributes to improve our understanding of the physical mechanisms underlying large megathrust earthquakes, and demonstrate that fluid play a key role in controlling the interplay between seismic and aseismic slip.</p>

Propagation velocity of magma intrusions, a new 2D numerical approach

<p>The tortuous travel of magma through the crust may sometimes result in volcanic eruptions at the surface. In the brittle crust, magma propagation usually occurs by fracturing the rock and opening its own way through them. This process of diking is controlled by the interaction of many complex physical processes including rock fracture, flow of compressible fluids, phase transitions, heat exchange. Current models of dikes consider either a fracturing-dominated approach, that neglects the viscous flow and allow to estimate the trajectory of dike propagation, or a viscous-dominated approach that neglects the fracturing at the dike tip allowing to infer the propagation velocity of the dike. Here we propose a new numerical approach aiming at modeling both the magma path and velocity. We start from a two-dimensional Boundary Element model solving for the trajectory of a quasi-static crack in an elastic medium subjected to external stress (Maccaferri et al, 2011), and implement the effects of the viscous fluid flow assuming a Poiseuille flow. We build on the previous work by Dahm (2000) but relaxing the assumption of stationarity, and thus allowing to take into account heterogeneous crustal stresses, complex dike paths, and dike velocity variations. The fluid flow results in a viscous pressure drop applied to the crack wall, which modifies the crack shape and contributes to the energy balance of the propagating dike. In fact, the energy dissipated by viscous flow is linearly dependant on the viscosity of the fluid and the crack velocity. It follows that the velocity can be inferred from the total energy budget by imposing that the viscous energy dissipation and the energy spent to fracture the rocks equals the strain-plus-gravitational energy release. However, the viscous dissipation strongly depends on the opening of each dislocation element, causing numerical instabilities in the calculation of the dike velocity due to the fracture closure at the dike tail. We will present first results of velocities derived with this approach considering only a static crack shape (that is to say neglecting the modification of the crack shape induced by the flow). We will discuss the influence of various parameters (crack size, Young’s modulus value...), and will compare the numerical velocities obtained with observations, first considering velocities measured in analogue experiments when injecting fluids of various viscosity (air and oils) in gelatin tanks, and secondly using diking events documented at basaltic volcanoes (such as Piton de la Fournaise (Réunion) and Mount Etna (Sicily)).</p>

Delineating active faults across the Narmada-Son Lineament (NSL), India using the technique of Fracture Induced Electromagnetic Radiation (FEMR)

<p>The Fracture Induced Electromagnetic Radiation (FEMR) technique has gradually progressed in the past decade as a useful geophysical tool to determine the direction and magnitude of recent crustal stresses, visualize the modification and realignment of stresses inside tunnels thus proving to be an important precursor for geohazards, earthquake forecasting, as well as delineate landslide-prone slip planes in unstable regions. Its working principle is based on the generation of geogenic electromagnetic radiation emanating from the brittle rock bodies that are fractured being subjected to an incremental increase of the differential stress in the near-surface of the Earth’s crust. The “Process zone” at the fractured crack tip contains numerous microcracks which subsequently creates dipoles due to the polarization of charges on such microcrack tips which rapidly oscillates emitting FEMR waves of frequencies between KHz to MHz range. The coalescence of the microcracks eventually leads to a macro failure dampening the amplitude of the FEMR pulses. The attenuation of FEMR pulses is comparatively lesser than seismic waves making it a more efficient precursor to potential tectonic activities indicating an upcoming earthquake a few hours/days before the actual event. In the current study, we have attempted to exploit this technique to identify the locations of the potential active faults across the tectonically active Narmada-Son Lineament (NSL), Central India. Although the first tectonic stage involved rifting and formation of the NSL during the Precambrian time, the rifting continued at least till the time of Gondwana deposition. Later, tectonic inversion took place as a result of the collision between the Indian and the Eurasian plate resulting in reverse reactivation of the faults. Episodic reverse movement along NSL caused recurrent earthquakes and linear disposition of the sediments that were deposited at the foothills of the Satpura Horst. Although the origin of East-West trending NSL dates back to the Precambrian time, it is very much tectonically active as manifested by recent earthquakes. The study has been conducted by taking linear FEMR readings across 3 traverses along the NSL which on analysis provides an idea about the potential active faults, their locations, and frequency of occurrence. The accumulation of strain in the brittle rocks that can eventually lead to a macro failure is demarcated as an anomalous increase in the amplitude of the FEMR pulses indicative of an upcoming tectonic episode in the region. To further corroborate the analysis, we have attempted to determine the neo-tectonic activity in the region by calculating the morphometric parameters across the Khandwa-Itarsi-Jabalpur region, Central India. Finally, we attempt to comment on the tectonic evolution of Central India in the recent past. We also encourage researchers to adapt the novel technique of FEMR which is swift, affordable, and feasible compared to conventional techniques deployed to survey the active tectonics of a region.</p>

Fracture Mechanism on Supercritical Carbon Dioxide Blasting in Low-Permeability Coal Under Various Crustal Stresses

The low permeability of coalbed inhabits exploring coalbed methane (CBM) in China. In order to enhance CBM production, supercritical carbon dioxide (SC-CO2) blasting is used as an efficient and promising method. In the paper, the experiments of SC-CO2 blasting were performed for coal and concrete specimens in the large size (1 m × 1 m × 0.5 m) under different crustal stresses, and then the fracture mechanism was analyzed. The results show that Crushing zone, rupture zone, and vibration zone were formed in sequence around the blasting hole, which provide the experimental and theoretical foundation for the field blasting.

Episodic fluid pressure cycling controls the interplay between Slow Slip Events and Large Megathrust Earthquakes

<p>Slow slip events (SSEs) are part of a spectrum of aseismic processes that relieve tectonic stress on faults. Their occurrence in subduction zones have been suggested to trigger megathrust earthquakes due to perturbations in fluid pressure. However, examples to date have been poorly recorded and physical observations of temporal fluid pressure fluctuations through slow slip cycles remain elusive. Here, we use a newly developed two-phase flow numerical model — which couples solid rock deformation and pervasive fluid flow — to show how crustal stresses and fluid pressures within subducting megathrust evolve before and during slow slip and regular events. This unified 2D numerical framework couples inertial mechanical deformation and fluid flow by using finite difference methods, marker-in-cell technique, and poro-visco-elasto-plastic rheologies. Furthermore, an adaptive time stepping allows the correct resolution of both long- and short-time scales, ranging from years to milliseconds during the dynamic propagation of earthquake rupture.</p><p>Here we show how permeability and its spatial distribution control the degree of locking along the megathrust interface and the interplay between seismic and aseismic slip. While a constant permeability leads to more regular seismic cycles, a depth dependent permeability contributes substantially to the development of two distinct megathrust zones: a shallow, locked seismogenic zone and a deep, narrow aseismic segment characterized by SSEs. Furthermore, we show that without requiring any specific friction law, our model shows that permeability, episodic stress transfer and fluid pressure cycling control the predominant slip mode along the subduction megathrust. Specifically, we find that the up-dip propagation of episodic SSEs systematically decreases the fault strength due to a continuous accumulation and release of fluid pressure within overpressured subducting interface, thus affecting the timing of large megathrust earthquakes. These results contribute to improve our understanding of the physical driving forces underlying the interplay between seismic and aseismic slip, and demonstrate that slow slip events may prove useful for short-term earthquake forecasts.</p>

Sedimentary loading-unloading cycles and faulting in intermontane basins: insights from numerical modeling and field observations from the broken foreland basin of NW Argentine Andes

<p>The removal, redistribution, and transient storage of sediments in tectonically active mountain belts is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but unambiguous evidence for this potential link has been elusive so far.</p><p>Here we combine 2D numerical experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks) with field data from intermontane basins in the NW Argentine Andes. Our modelling results suggest that thicker sedimentary fills (> 0.7-1.0 km) may suppress basinal faulting processes, while thinner fills (< 0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and basin excavation promotes renewed intrabasinal faulting.</p><p>These results help to better understand the tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean Plateau in northwestern Argentina. For example, the Santa María and the Humahuaca basins record intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal deformation patterns in intermontane basins of tectonically active broken forelands.</p>