Search NASA⌕ Search
Home/Archive/DOE OSTI

DOE OSTI · 2423205

Multiphysics modelling in PyLith: poroelasticity

Walker, Robert L.·Knepley, Matthew G.·Aagaard, Brad T.·Williams, Charles A.
View original source ↗Export metadata ↓

Abstract

SUMMARY PyLith, a community, open-source code for modelling quasi-static and dynamic crustal deformation with an emphasis on earthquake faulting, has recently been updated with a flexible multiphysics implementation. We demonstrate the versatility of the multiphysics implementation by extending the code to model fully coupled continuum poromechanics. We verify the newly incorporated physics using standard benchmarks for a porous medium saturated with a slightly compressible fluid. The benchmarks include the 1-D consolidation problem as outlined by Terzaghi, Mandel’s problem for the 2-D case, and Cryer’s problem for the 3-D case. All three benchmarks have been added to the PyLith continuous integration test suite. We compare the closed form analytical solution for each benchmark against solutions generated by our updated code, and lastly, demonstrate that the poroelastic material formulation may be used alongside the existing fault implementation in PyLith.

Explore related subjects

Geochemistry & Geophysics

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Walker, Robert L., Knepley, Matthew G., Aagaard, Brad T., Williams, Charles A.. 2023-09-18. Multiphysics modelling in PyLith: poroelasticity. https://doi.org/10.1093/gji%2Fggad370

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

DOE OSTI2024

Formation of a Large Cold Groundwater Mantle Helium Anomaly and High Temperature Geothermal Resources in Response to Bimodal Magmatism Near Roosevelt Hot Springs and Utah FORGE, Milford Valley, Southwest Utah

Abstract A large mantle helium anomaly and separate domains of high heat flow are the predominant manifestations of bimodal magmatic activity in the Milford valley. The mantle helium anomaly (1.9–2.6 R/Ra) covers 270 km 2 and is subdivided into two separated domains: a cold shallow groundwater regime and high temperature hydrothermal activity. The zone of anomalous heat flow covers >100 km 2 and is also subdivided into two adjacent domains, comprising hydrothermal activity at Roosevelt Hot Springs (RHS) (3–7 W/m 2 ) and conductive heat flow (100–180 mW/m 2 ). While the transfer of heat and mantle helium at RHS are coupled, heat and helium transfer are decoupled in the adjacent cold groundwater regime to the west. Both the mantle helium and geothermal anomalies are attributed to recent mafic‐felsic magmatic intrusions of >400 km 3 , however, the absence of volcanic eruptions <500,000 years indicates magmas stall before rising to shallow crustal level <10 km depth. Deep level magmatism produces a felsic composition melt, which is inferred to be responsible for the widespread and near uniform range of diluted mantle helium values. A thick and impermeable mass of crystalline granitic basement rock at the mid‐crustal level divides the ascent of mantle helium into separate flow paths. It may also impede the rise of buoyant magma trapping thermal energy that facilitates partial melting, slow cooling, and development of a thick thermal aureole. Partitioning of convective and conductive thermal regimes and independent flow paths supplying deeply derived helium characterize the development of a large long‐lived magma‐related geothermal system.

Geochemistry & Geophysics↗
DOE OSTI2024

Geodetic Evidence for Distributed Shear Below the Brittle Crust of the Walker Lane, Western United States

Abstract Models of active deformation of the Earth's crust are predominantly represented with dislocations having a downdip continuation into the lower crust, where the fault slips continuously. This model predicts surface strain accumulation concentrated near the fault during the interseismic period. In an alternative model, faults do not extend beneath the elastic portion of the crust and are accompanied by a wide zone of distributed shear underneath, predicting a more constant strain rate lacking concentrations at the faults. We use high‐precision GPS data collected across the northern and central Walker Lane, USA— a region of complex faulting near the western edge of the Basin and Range Province to evaluate which model is appropriate. Despite the existence of dense continuous and semi‐continuous geodetic networks that have been surveyed for ∼20 years, the horizontal velocities reveal no evidence of localized strain accumulation across the fault surface expressions. Instead, deformation within the Walker Lane is uniformly linear, suggesting that the surface deformation reflects distributed shear within the ductile crust rather than focused deformation at faults. This suggests no downdip extension of the faults below the seismogenic layer. The shear zone is 172 ± 6 km wide in the northernmost Walker Lane narrowing to 116 ± 4 km in the central Walker Lane. The total velocity budget across the shear zone is 7.2 ± 0.1 mm/yr in the north, increasing to 10.1 ± 0.1 mm/yr in the central Walker Lane. We conclude that assuming the presence of lower crustal dislocations when estimating geodetic faults slip rates may be inappropriate.

Geochemistry & Geophysics↗