Search

Dynamic modeling of sequences of earthquakes and aseismic slip (SEAS) provides a self-consistent, physics-based framework to connect, interpret, and predict diverse geophysical observations across spatial and temporal scales. Amid growing applications of SEAS models, numerical code verification is essential to ensure reliable simulation results but is often infeasible due to the lack of analytical solutions. Here, we develop two benchmarks for three-dimensional (3D) SEAS problems to compare and verify numerical codes based on boundary-element, finite-element, and finite-difference methods, in a community initiative. Our benchmarks consider a planar vertical strike-slip fault obeying a rate- and state-dependent friction law, in a 3D homogeneous, linear elastic whole-space or half-space, where spontaneous earthquakes and slow slip arise due to tectonic-like loading. We use a suite of quasi-dynamic simulations from 10 modeling groups to assess the agreement during all phases of multiple seismic cycles. We find excellent quantitative agreement among simulated outputs for sufficiently large model domains and small grid spacings. However, discrepancies in rupture fronts of the initial event are influenced by the free surface and various computational factors. The recurrence intervals and nucleation phase of later earthquakes are particularly sensitive to numerical resolution and domain-size-dependent loading. Despite such variability, key properties of individual earthquakes, including rupture style, duration, total slip, peak slip rate, and stress drop, are comparable among even marginally resolved simulations. Our benchmark efforts offer a community-based example to improve numerical simulations and reveal sensitivities of model observables, which are important for advancing SEAS models to better understand earthquake system dynamics.

To what extent mechanical anisotropy is required to explain the dynamics of the lithosphere is an important yet unresolved question. If anisotropy affects stress and deformation, and hence processes such as fault loading, how can we quantify its role from observations? Here, we derive analytical solutions and build a theoretical framework to explore how a shear zone with linear anisotropic viscosity can lead to deviatoric stress heterogeneity, strain-rate enhancement, as well as non-coaxial principal stress and strain rate. We develop an open-source finite-element software based on FEniCS for more complicated scenarios in both 2-D and 3-D. Mechanics of shear zones with transversely isotropic and orthorhombic anisotropy subjected to misoriented shortening and simple shearing are explored. A simple regional example for potential non-coaxiality for the Leech River Schist above the Cascadia subduction zone is presented. Our findings and these tools may help to better understand, detect and evaluate mechanical anisotropy in natural settings, with potential implications including the transfer of lithospheric stress and deformation through fault loading. [ABSTRACT FROM AUTHOR]

This study investigates the use of calcium-based slurry for simultaneous removal NO and SO2 from oxy-fuel combustion flue gas, and recovery of the sulfur and nitrogen species in resulting solutions. The experiments were performed in a bubbling reactor in a transient mode under the pressure of 20 bar. The various influencing factors including the CaO amount, carrier gas (N2/CO2), and absorption time on the simultaneous NO and SO2 removal process, and the solution products were studied comprehensively. The results show that the NO2 removal efficiency can be improved by the presence of CO2, and the gas phase HNO2 produces in this process. The addition of CaO has positive effects not only on the NO2 removal efficiency but also on the formation of stable HNO3. With the presence of CO2, CaCO3 is formed in a solution initially. With the decrease of pH, CaCO3 is gradually converted to CaSO4, and in particular CaCO3 can be fully avoided through decreasing the pH of an absorption solution to 1.14. At the same time, the formation of unstable S(IV) and NO2- can be prevented when the solution pH is lower than 1.37. The nitrogen and sulfur compounds in the absorption solution (at pH 1.14) were further separated by the addition of different amounts of CaO. In particular, 95% of SO42- finally can be recovered in the form of CaSO4.2H2O with nitrogen in solution existing as NO3- by controlling the Ca/S ratio at 4.70. The effectiveness of calcium-based slurry on the removal and recovery of SO2 and NO is confirmed.

A bubbling reactor is an important type of gas scrubber to reduce SO2 emissions in maritime shipping. Both experiments and simulations were conducted to study the relationship between the periodic gas bubbling process and SO2 concentration at the outlet of the reactor, and the entrainment of liquid droplets on SO2 absorption. The accuracy of the model was verified by comparing the bubble size, the depth of bubbles injected into the water, and the SO2 concentration obtained in both experiments and simulations. The gas bubbling process is accompanied by bubble formation, rise, and collapse. The gas bubbling period is affected by the disturbance of the liquid level. The period of the SO2 concentration at the outlet of the gas bubbling reactor is smaller than that at the gas jar outlet which acts as the gas buffering region. The amounts of water carried out of the bubbling reactor by the gas bubbling process increase with the gas flow rates. The droplets and liquid film in the gas jar and the connecting tube play an important role in the absorption of SO2. This study encourages more research to reduce the fluctuation of SO2 concentration and consider droplet entrainment in the design of bubbling reactors. [ABSTRACT FROM AUTHOR]