Towards a poroelastodynamics framework for induced earthquakes

Earthquakes can be triggered after pore pressure perturbations activate critically stressed seismogenic faults, where the perturbations can originate from natural causes like earth tides, rainfall, snowfall or anthropogenic causes like wastewater disposal, CO$_2$ injection, oil production, or groundwater extraction. As the faults slip under the action of the induced stress field, seismic waves are spawned from the hypocenter location. The waves propagate through the domain with a velocity that evolves with the evolving pressure and stress fields. The effect of these waves on the surrounding rock and the seismic velocity recorded on the seismograph can be modeled accurately only by incorporating elastodynamics in the deformation model coupled with flow-induced pressure perturbations. Hitherto, most of the literature in the realm has been limited to elastostatics coupled with flow within a prescribed/kinematic or quasi-dynamic fault slip framework. In this work, we provide a framework for coupling of wave propagation with pore pressure perturbations using one-way coupled poroelastodynamics in the presence of faults in which the pore pressure is specified apriori as a spatiotemporal function.We present results from analysis of displacement and velocity fields in the domain and tractions and slip evolution on the fault. The rendition of two-way coupled poroelastodynamics in which the flow problem is also solved is proposed as future work.

Search for Infrasound Signals in InSight Data Using Coupled Pressure/Ground Deformation Methods

ABSTRACT The unprecedented quality and sampling rate of seismometer and pressure sensors of the InSight Mars mission allow us to investigate infrasound through its pressure and ground deformation signals. This study focuses on compliance effects induced by acoustic waves propagating almost horizontally close to the surface. The compliance of acoustic waves is first estimated using the compliance estimates from pressure perturbations moving at wind speed. Then, a marker of compliance events is used to detect events of ground deformation induced by pressure variations, in three frequency bands from 0.4 to 3.2 Hz, from InSight sol 180 to 690. Additional selection criteria are imposed on the detected events to focus on acoustic waves and to remove various noise sources (e.g., wind effects or seismometer artifacts). After an automated selection, the visual inspection of the records allows us to validate two infrasound candidates that cannot be related to pressure perturbations moving at wind speed nor to known noise sources. For our highest quality infrasound candidate, the relation between this event and a convective vortex occuring 10 s later is tested. The azimuth of the vortex position at the time of infrasound detection is not consistent with the arrival azimuth of the suspected infrasound inferred from the polarization of seismometer records, thus the link between these two phenomena cannot be demonstrated. Further investigations would require a better understanding of wind-related noise impacting InSight sensors and of the effects of lateral variations of subsurface mechanical properties on the ground deformations induced by atmospheric pressure variations.

Detecting Wave Features in Doppler Radial Velocity Radar Observations

Mesoscale, wave-like perturbations in horizontal air motions in the troposphere (velocity waves) are associated with vertical velocity, temperature, and pressure perturbations that can initiate or enhance precipitation within clouds. The ability to detect velocity waves from horizontal wind information is an important tool for atmospheric research and weather forecasting. This paper presents a method to routinely detect velocity waves using Doppler radial velocity data from a scanning weather radar. The method utilizes the difference field between consecutive PPI scans at a given elevation angle. Using the difference between fields a few minutes apart highlights small scale perturbations associated with waves because the larger scale wind field changes more slowly. Image filtering retains larger contiguous velocity bands and discards noise. Wave detection scales are limited by the size of the temporal difference relative to the wave motion and the radar resolution volume size.

Impact of Foreshock Transients on the Flank Magnetopause and Magnetosphere and the Ionosphere

Mesoscale (on the scales of a few minutes and a few RE) magnetosheath and magnetopause perturbations driven by foreshock transients have been observed in the flank magnetotail. In this paper, we present the 3D global hybrid simulation results to show qualitatively the 3D structure of the flank magnetopause distortion caused by foreshock transients and its impacts on the tail magnetosphere and the ionosphere. Foreshock transient perturbations consist of a low-density core and high-density edge(s), thus, after they propagate into the magnetosheath, they result in magnetosheath pressure perturbations that distort magnetopause. The magnetopause is distorted locally outward (inward) in response to the dip (peak) of the magnetosheath pressure perturbations. As the magnetosheath perturbations propagate tailward, they continue to distort the flank magnetopause. This qualitative explains the transient appearance of the magnetosphere observed in the flank magnetosheath associated with foreshock transients. The 3D structure of the magnetosheath perturbations and the shape of the distorted magnetopause keep evolving as they propagate tailward. The transient distortion of the magnetopause generates compressional magnetic field perturbations within the magnetosphere. The magnetopause distortion also alters currents around the magnetopause, generating field-aligned currents (FACs) flowing in and out of the ionosphere. As the magnetopause distortion propagates tailward, it results in localized enhancements of FACs in the ionosphere that propagate anti-sunward. This qualitatively explains the observed anti-sunward propagation of the ground magnetic field perturbations associated with foreshock transients.

Differential responsiveness of glabrous and non-glabrous skin to local transmural pressure elevations; the impact of 5 weeks of iterative local pressure loading

We examined the in vivo pressure-flow relationship in human cutaneous vessels during acute and repeated elevations of local transmural pressure. In 10 healthy men, red blood cell flux was monitored simultaneously on the non-glabrous skin of the forearm and the glabrous skin of a finger during a vascular pressure provocation, wherein the blood vessels of an arm were exposed to a wide range of stepwise increasing distending pressures. Forearm skin blood flux was relatively stable at slight and moderate elevations of distending pressure, whereas it increased ~3-4-fold at the highest levels (P = 0.004). Finger blood flux on the contrary, dropped promptly and consistently throughout the provocation (P < 0.001). Eight of the subjects repeated the provocation trial after a 5-week pressure-training regimen, during which the vasculature in one arm was exposed intermittently (40 min, 3 times･week-1) to increased transmural pressure (from +65 mmHg week-1 to +105 mmHg week-5). The training regimen diminished the pressure-induced increase in forearm blood flux by ~34% (P = 0.02), whereas it inhibited the reduction in finger blood flux (P < 0.001) in response to slight and moderate distending pressure elevations. The present findings demonstrate that, during local pressure perturbations, the cutaneous autoregulatory function is accentuated in glabrous compared to in the non-glabrous skin regions. Prolonged intermittent regional exposures to augmented intravascular pressure blunt the responsiveness of the glabrous skin, but enhance arteriolar pressure resistance in the non-glabrous skin.

Ocean glider flight in the presence of surface waves

Autonomous, buoyancy-driven ocean gliders are increasingly used as a platform for the measurement of turbulence microstructure. In the processing of such measurements, there is a sensitive (quartic) dependence of the turbulence dissipation rate, ϵ, on the speed of flow past the sensors, or alternatively, the speed of the glider through the ocean water column. The mechanics of glider flight is therefore examined by extending previous flight models to account for the effects of ocean surface waves. It is found that due to the relatively small buoyancy changes used to drive gliders, the surface wave-induced motion, superimposed onto the steady-state motion, follows to a good approximation the motion of the wave orbitals. Errors expected in measuring ϵ at the ocean near-surface due to wave-induced relative velocities are generally less than 10%. However, pressure perturbations associated with the wave motion can be significant when using the glider-measured pressure signal to infer the glider vertical velocity. This effect of surface waves is only present in the shallow water regime, and can also affect glider depth measurements. It arises from an incomplete cancellation of the wave-induced pressure perturbation with the hydrostatic component due to vertical glider displacements, whereas for deep-water waves this cancellation is complete.

Hydrodynamic instability at impact interfaces and planetary implications

Impact-induced mixing between bolide and target is fundamental to the geochemical evolution of a growing planet, yet aside from local mixing due to jetting – associated with large angles of incidence between impacting surfaces – mixing during planetary impacts is poorly understood. Here we describe a dynamic instability of the surface between impacting materials, showing that a region of mixing grows between two media having even minimal initial topography. This additional cause of impact-induced mixing is related to Richtmyer-Meshkov instability (RMI), and results from pressure perturbations amplified by shock-wave refraction through the corrugated interface between impactor and target. However, unlike RMI, this new impact-induced instability appears even if the bodies are made of the same material. Hydrocode simulations illustrate the growth of this mixing zone for planetary impacts, and predict results suitable for experimental validation in the laboratory. This form of impact mixing may be relevant to the formation of stony-iron and other meteorites.

A fault slip model to study earthquakes due to pore pressure perturbations

The burgeoning need to sequester anthropogenic CO_2 for climate mitigation and the need for energy sustenance leading upto enhanced geothermal energy production has made it incredibly critical to study potential earthquakes due to fluid activity in the subsurface. These earthquakes result from reactivation of faults in the subsurface due to pore pressure perturbations. In this work, we provide a framework to model fault slip due to pore pressure change leading upto quantifying the earthquake magnitude.

A kinetic formation model of foreshock transients

&lt;p&gt;Foreshock transients are ion kinetic structures in the ion foreshock. Due to their dynamic pressure perturbations, they can disturb the bow shock, magnetosheath, magnetopause, and magnetosphere-ionosphere system. Recent studies found that they can also accelerate particles through shock drift acceleration, Fermi acceleration, betatron acceleration, and magnetic reconnection. Although foreshock transients are important, how they form is still not fully understood. Using particle-in-cell simulations and MMS observations, we propose a physical formation process that the positive feedback of demagnetized foreshock ions on the varying magnetic field caused by the foreshock ion Hall current enables an &amp;#8220;instability&amp;#8221; and the growth of the structure.&amp;#160;&amp;#160;&amp;#160;&amp;#160;&amp;#160;&amp;#160;&lt;/p&gt;