A Fluid-Diffusion-Hybrid Limiting Approximation for Priority Systems with Fast and Slow Customers

The paper “Fluid-Diffusion-Hybrid (FDH) Approximation” proposes a new heavy-traffic asymptotic regime for a two-class priority system in which the high-priority customers require substantially larger service times than the low-priority customers. In the FDH limit, the high-priority queue is a diffusion, whereas the low-priority queue operates as a (random) fluid limit, whose dynamics are driven by the former diffusion. A characterizing property of our limit process is that, unlike other asymptotic regimes, a non-negligible proportion of the customers from both classes must wait for service. This property allows us to study the costs and benefits of de-pooling, and prove that a two-pool system is often the asymptotically optimal design of the system.

Fast Evaluation of Caprock Strength Sensitivity to Different CO2 Solutions Using Small Sample Techniques

The Draupne shale is a rock formation functioning as overburden for gas reservoirs in the Norwegian Sea and potentially as caprock for future CO2 storage locations as well. In this paper, the Draupne shale was exposed to several fluids: CO2 gas, supercritical CO2, CO2 gas dissolved in brine, supercritical CO2 dissolved in brine, as well as brine and dry air. The motivation for the exposure tests was to investigate whether injected CO2 in a reservoir coming into contact with the caprock could change the caprock’s mechanical properties and increase the risk for leakage of the stored CO2. In addition, a systematic exposure study will provide more insight into the various processes susceptible of altering the shale’s shear strength and acoustic velocity, such as clay hydration, mineral dissolution, and capillary forces. Due to the low permeability of the shale, experiments were conducted on mm-sized disk samples, reducing fluid diffusion into the shale, and allowing for many repeated tests on disks close by in the original core. The punch method, where a small circle is punched out of the shale, was used to assess shear strength, while continuous wave technique was used to assess ultrasonic velocity. Results show that the shale is not noticeably sensitive to CO2, in the sense that no additional weakening is observed in the presence of CO2 as compared to brine exposure. This last weakening effect is probably due to poor matching between pore fluid salinity and exposure brine strength.

Near-Fault Monitoring Reveals Combined Seismic and Slow Activation of a Fault Branch within the Istanbul–Marmara Seismic Gap in Northwest Turkey

Various geophysical observations show that seismic and aseismic slip on a fault may occur concurrently. We analyze microseismicity recordings from a temporary near-fault seismic network and borehole strainmeter data from the eastern Marmara region in northwest Turkey to track seismic and aseismic deformation around the hypocentral region of an Mw 4.5 earthquake in 2018. A slow transient is observed that lasted about 30 days starting at the time of the Mw 4.5 event. We study about 1200 microseismic events that occurred during 417 days after the Mw 4.5 event around the mainshock fault rupture. The seismicity reveals a strong temporal clustering, including four episodic seismic sequences, each containing more than 30 events per day. Seismicity from the first two sequences displayed typical characteristics driven by aseismic slip and/or fluids, such as the activation of a broader region around the mainshock and swarm-like topology. The third and fourth sequences correspond to typical mainshock–aftershock sequences. These observations suggest that slow slip and potentially fluid diffusion along the fault plane could have controlled the seismicity during the initial 150 days following the Mw 4.5 event. In contrast, stress redistribution and breaking of remaining asperities may have caused the activity after the initial 150 days. Our observation from a newly installed combined dense seismic and borehole strainmeter network follows an earlier observation of a slow transient occurring in conjunction with enhanced local seismic moment release in the same region. This suggests a frequent interaction of seismic and aseismic slip in the Istanbul–Marmara seismic gap.

Effect of Modification on the Fluid Diffusion Coefficient in Silica Nanochannels

The diffusion behavior of fluid water in nanochannels with hydroxylation of silica gel and silanization of different modified chain lengths was simulated by the equilibrium molecular dynamics method. The diffusion coefficient of fluid water was calculated by the Einstein method and the Green–Kubo method, so as to analyze the change rule between the modification degree of nanochannels and the diffusion coefficient of fluid water. The results showed that the diffusion coefficient of fluid water increased with the length of the modified chain. The average diffusion coefficient of fluid water in the hydroxylated nanochannels was 8.01% of the bulk water diffusion coefficient, and the diffusion coefficients of fluid water in the –(CH2)3CH3, –(CH2)7CH3, and –(CH2)11CH3 nanochannels were 44.10%, 49.72%, and 53.80% of the diffusion coefficients of bulk water, respectively. In the above four wall characteristic models, the diffusion coefficients in the z direction were smaller than those in the other directions. However, with an increase in the silylation degree, the increased self-diffusion coefficient due to the surface effect could basically offset the decreased self-diffusion coefficient owing to the scale effect. In the four nanochannels, when the local diffusion coefficient of fluid water was in the range of 8 Å close to the wall, Dz was greater than Dxy, and beyond the range of 8 Å of the wall, the Dz was smaller than Dxy.

Recurrence of Fault Valve Behavior in a Continental Collision Area: Evidence From Tilt/Strain Measurements in Northern Adria

We analyzed the data recorded by the NE-Italy subsurface tilt and strainmeter network evidencing a coherent transient signal in the recordings of four tiltmeter sites in the 1984–1990 period that produced a tilt along the main fractures. Borrowing from classical seismology techniques, we used the uprise times to locate the transient signal source. The propagation velocity is compatible with a fluid diffusion process that starts from a source located close to the hypocenter of the February 10, 1983 Uccea earthquake, MD = 4.2 at the Italy-Slovenia border, at an estimated depth of 10.8 km. Our results add to the previous interpretation of a transient signal recorded by several global navigation satellite system (GNSS) stations in the 2006–2009 period in terms of fluid diffusion below the Bovec basin (Slovenia). That source was located upon continuation to the northwest of the Ravne fault, few kilometers to the northeast from the present one, and about 6 km from the July 12, 2004 Bovec–Krn earthquake, Mw = 5.1, depth ~6.1 km. These observations suggest that the area is subject to fault valve behavior episodes that released fluids trapped at depth to the surrounding region as pore-pressure bulges. The convergence between Alpine and External Dinarides structures in this area puts highly permeable dolomitic limestones in contact with low-permeable fine-grained limestones and flysch formations. Therefore, the conditions for overpressure generation can be created, whereas fault movements, from time to time, in close relation with seismic events, can enable fluid diffusion in the surroundings. We also estimated the possible fluid influx needed to maintain overpressure and possible discharge across both the faults. The study provides insights on pore–fluid pressure variations related to slow slip events from a context different from subduction or transform margins, i.e., in a continental collision area.

The tilt and strainmeter network of NE Italy: multi-decadal observations of crustal deformation as ground truth for DinSAR.

<p>Decade-long series of tilt- and strain-meter observations in NE Italy allow monitoring the crustal deformation from short transient to long-term phenomena. These recordings, some of them started in 1960, are generated by sources spanning a wide spectrum of spatial scales, such as sudden underground flooding due to extreme rainfall [1, 2], years-long fluid diffusion transients due to fault behavior [3], the free oscillation arising from megathrust earthquakes (e.g. Chile 1960, Sumatra 2004, Tohoku 2011).<br>The instrumental sites lie on karst formations, in an area of continental collision and active seismicity, the northeastern portion of the Adria microplate, where the south-directed thrusts of the Alpine system merge with the NW-SE transpressive regime of the External Dinarides. Measurements include the ongoing interseismic strain accumulation processes, including the peculiar observation of episodic disturbances and southward tilting in the three years preceding the 1976 Mw6.4 Friuli earthquake [4].<br><br>The channel systems of Karst hydrology, which undergo complete flooding and overpressure buildup in extremely short time spans (e.g. near-simultaneous flooding over a distance of 30 km) result in observable surface deformation and a change in the gravity field. Tilt time series allow to extract and model this type of hydrology-forced uplift and associated deformation [2,5].<br><br>Tilt- and strain-meters allow for accuracy and precision in measuring crustal deformation, to a level which space-borne geodesy cannot provide. The main drawback, however, is that only point measurements are provided, in locations where stations could be set up.<br>On the other hand, the thousands of points on the surface that DInSAR can provide are affected by coarser accuracy and influenced by atmospheric effects - resulting in LoS displacements uncorrelated to the actual surface deformations. We aim at enabling the transfer of knowledge from tilt- and strain-meters observations to DInSAR-derived data, thus allowing a first assessment of ground-truth constrained displacement models.<br><br>[1] Braitenberg C. (2018). The deforming and rotating Earth - A review of the 18th International Symposium on Geodynamics and Earth Tide, Trieste 2016 , Geodesy and Geodynamics, 187-196, doi::10.1016/j.geog.2018.03.003</p><p>[2] Braitenberg C., Pivetta T., Barbolla D. F., Gabrovsek F., Devoti R., Nagy I. (2019). Terrain uplift due to natural hydrologic overpressure in karstic conduits. Scientific Reports, 9:3934, 1-10, doi.:10.1038/s41598-019-38814-1</p><p>[3] Rossi, G., Fabris, P. & Zuliani, D. Overpressure and Fluid Diffusion Causing Non-hydrological Transient GNSS Displacements. Pure Appl. Geophys. 175, 1869–1888 (2018). https://doi.org/10.1007/s00024-017-1712-x</p><p>[4] Dragoni M., Bonafede M., and Boschi E. (1985). On the interpretation of slow ground deformation precursory to the 1976 Friuli earthquake. Pure and Applied Geophysics 122, 781–792. doi:10.1007/978-3-0348-6245-5_3</p><p>[5] Grillo B., Braitenberg C., Nagy I., Devoti R., Zuliani D., Fabris P. (2018). Cansiglio Karst-Plateau: 10 years of geodetic-hydrological observations in seismically active northeast Italy. Pure and Applied Geophysics, 175, 5, 1765-1781, doi:10.1007/s00024-018-1860-7.</p><p> </p>

Dual seismic migration velocities reveal imbricated fluid diffusion and aseismic slip In a Corinth Gulf swarm (Greece)

<p>Fluid induced earthquake sequences generally appear as expanding swarms activating a particular fault. Such swarms are generally interpreted as fluid diffusion, which ignores the possibility of static, dynamic or aseismic triggering, and the existence of rapid migration. Here, we study the temporal evolution of a seismic swarm that occurred over a 10-day period in October 2015 in the extensional rift of the Corinth Gulf (Greece) using high-resolution earthquakes relocations. The seismicity radially migrates on a normal fault at a fluid diffusion velocity (~125 m/day). However, this migration occurs intermittently, with periods of fast expansion (2-to-10 km/day) during short seismic bursts alternating with quiescent periods. Moreover, the growing phases of the swarm illuminate a high number of repeaters. Therefore, we propose a new model to explain the combination of multiple driving processes for such swarms.  Fluid up flow in the fault may induce aseismic slip episodes, separated by phases of fluid pressure build-up. The stress perturbation due to aseismic slip may activate small asperities in the fault that produce bursts of seismicity during the most intense phase of the swarm. We then validated this model through hydro-mechanical modeling, where earthquakes consist in the failure of asperities on a creeping fault infiltrated by fluid. For that, we couple rate‐and‐state friction, non‐linear diffusivity and elasticity along a 1D interface. This model reproduces the dual migration speeds observed in real swarms. We show that migration speeds increase linearly with the mean pressurization, and are not dependent on the hydraulic diffusivity, as traditionally suggested.</p>

The 2020 Haenam earthquake sequence: The first observation of a seismic front on the Korean Peninsula migrating in a manner similar to fluid diffusion

<p>Low-magnitude earthquakes (maximum Mw: 3.2) were recorded from late April 2020 onward in the county of Haenam, southwestern South Korea. Moderate to strong earthquakes had not previously been documented in instrumental, historical, or geological records. We identified 226 hypocentres in this earthquake sequence from April 25 to May 11, 2020. The seismic front of this sequence migrated in a manner similar to a diffusing fluid, with a hydraulic diffusivity of 0.012 m<sup>2</sup>/s. This is the first observation of natural seismicity on the Korean Peninsula imitating fluid diffusion. We applied a cross-correlation approach to detect unrecorded events, and relocated the hypocentres of the 71 previously recorded and 155 newly detected events using data collected at permanent seismic stations; clear linearity was observed at the metre scale. Spatially, the hypocentres were distributed within a 0.3 km × 0.3 km fault plane at a depth of ~20 km, trending west-northwest–east-southeast with a dip of ~70° in the south-southwestern direction. The moment tensor solution of the largest event had a strike of 98°, dip of 65°, and rake of 7°, which correspond to the fault geometry of the relocated hypocentres. The hypocentres progressed toward the upper eastern edge of the lineament. The largest event occurred at a shallow region of the fault plane, in the direction of hypocentre migration. Together, these results showed that the migration sequence of the 2020 Haenam earthquake mimicked the flow of a diffusing fluid. The structural data indicate that a fault–fracture mesh geometry channelled fluid flow, supporting the concept of a “fluid-driven earthquake swarm” for the 2020 Haenam earthquake sequence. Regarding the final parts of the sequence, there appeared to be a second intrusion at the western end, and a permeability barrier at the eastern end, of the fault plane. The well-constrained hypocentre locations in our study provide essential data for future research, and our interpretations of hypocentre migration during this earthquake sequence may help to elucidate the mechanisms driving earthquake swarms under conditions of intraplate stress.</p>

Bridging scales in disordered porous media by mapping molecular dynamics onto intermittent Brownian motion

Owing to their complex morphology and surface, disordered nanoporous media possess a rich diffusion landscape leading to specific transport phenomena. The unique diffusion mechanisms in such solids stem from restricted pore relocation and ill-defined surface boundaries. While diffusion fundamentals in simple geometries are well-established, fluids in complex materials challenge existing frameworks. Here, we invoke the intermittent surface/pore diffusion formalism to map molecular dynamics onto random walk in disordered media. Our hierarchical strategy allows bridging microscopic/mesoscopic dynamics with parameters obtained from simple laws. The residence and relocation times – tA, tB – are shown to derive from pore size d and temperature-rescaled surface interaction ε/kBT. tA obeys a transition state theory with a barrier ~ε/kBT and a prefactor ~10−12 s corrected for pore diameter d. tB scales with d which is rationalized through a cutoff in the relocation first passage distribution. This approach provides a formalism to predict any fluid diffusion in complex media using parameters available to simple experiments.

Conjugate faulting and structural complexity on the young fault system associated with the 2000 Tottori earthquake

Young faults display unique complexity associated with their evolution, but how this relates to earthquake occurrence is unclear. Unravelling the fine-scale complexity in these systems could lead to a greater understanding of ongoing strain localization in young fault zones. Here we present high-spatial-resolution images of seismic sources and structural properties along a young fault zone that hosted the Tottori earthquake (Mw 6.8) in southwest Japan in 2000, based on data from a hyperdense network of ~1,000 seismic stations. Our precise micro-earthquake catalog reveals conjugate faulting over multiple length scales. These conjugate faults are well developed in zones of low seismic velocity. A vertically dipping seismic cluster of about 200 m length occurs within a width of about 10 m. Earthquake migrations in this cluster have a speed of about 30 m per day, which suggests that fluid diffusion plays a role. We suggest that fine structural complexities influence the pattern of seismicity in a developing fault system.