Alexander S. Yates, Corentin Caudron, Aurélien Mordret, Philippe Lesage, Virginie Pinel, Thomas Lecocq, Craig A. Miller, Oliver D. Lamb, Nicolas Fournier
{"title":"Seasonal Snow Cycles and Their Possible Influence on Seismic Velocity Changes and Eruptive Activity at Ruapehu Volcano, New Zealand","authors":"Alexander S. Yates, Corentin Caudron, Aurélien Mordret, Philippe Lesage, Virginie Pinel, Thomas Lecocq, Craig A. Miller, Oliver D. Lamb, Nicolas Fournier","doi":"10.1029/2024JB029568","DOIUrl":"10.1029/2024JB029568","url":null,"abstract":"<p>Understanding volcanic eruption triggers is critical toward anticipating future activity. While internal magma dynamics typically receive more attention, the influence of external processes remains less understood. In this context, we explore the relationship between seasonal snow cycles and eruptive activity at Ruapehu, New Zealand. This is motivated by apparent seasonality in the eruptive record, where a higher than expected proportion of eruptions (post-1960) occur in spring (including the two previous eruptions of 2006 and 2007). Employing recent advancements in passive seismic interferometry, we compute sub-surface seismic velocity changes between 2005 and 2009 using the cross-wavelet transform approach. Stations on the volcano record a higher velocity in winter, closely correlated with the presence of snow. Inverting for depth suggests these changes occur within the upper 300 m. Notably, we observe that the timing of the previous two eruptions coincides with a period associated with an earlier velocity decrease at approximately 200–300 m depth relative to the surface. Reduced water infiltration (as precipitation falls as snow) is considered a likely control of seasonal velocities, while modeling also points to a contribution from snow-loading. We hypothesize that this latter process may play a role toward explaining seasonality in the eruptive record. Our findings shed light on the complex interactions between volcanoes and external environmental processes, highlighting the need for more focused research in this area. Pursuing this line of inquiry has significant implications toward improved risk and hazard assessments at not just Ruapehu, but also other volcanoes globally that experience seasonal snow cover.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2024JB029568","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142820929","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Melting of B1-Phase MgO From Simultaneous True Radiative Shock Temperature and Sound Speed Measurements to 250 GPa on Samples Preheated to 2300 K","authors":"O. V. Fat’yanov, P. D. Asimow","doi":"10.1029/2024JB029137","DOIUrl":"10.1029/2024JB029137","url":null,"abstract":"<p>To refine the melting curve, equation of state, and physical properties of MgO we performed plate impact experiments spanning <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>∼</mo>\u0000 </mrow>\u0000 <annotation> ${sim} $</annotation>\u0000 </semantics></math>170–250 GPa on <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo><</mo>\u0000 <mn>100</mn>\u0000 <mo>></mo>\u0000 </mrow>\u0000 <annotation> ${< } 100 > $</annotation>\u0000 </semantics></math> MgO single crystals, preheated to 2300 K. A controlled thermal gradient in <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>∼</mo>\u0000 </mrow>\u0000 <annotation> ${sim} $</annotation>\u0000 </semantics></math>20 mm long samples enabled radiative temperature (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>±</mo>\u0000 </mrow>\u0000 <annotation> $pm $</annotation>\u0000 </semantics></math>3%–4%) and rarefaction overtake observations (yielding sound speed <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>±</mo>\u0000 </mrow>\u0000 <annotation> $pm $</annotation>\u0000 </semantics></math>2%) close to the hot Mo driver with a free surface below <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>∼</mo>\u0000 </mrow>\u0000 <annotation> ${sim} $</annotation>\u0000 </semantics></math>2000 K that minimized evaporation. Ta flyers were launched by two-stage light-gas gun up to 7.6 km/s and sample radiance was recorded with a 6-channel (500–850 nm) pyrometer. Shock front reflectivity was measured at 198 and 243 GPa using <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>∼</mo>\u0000 </mrow>\u0000 <annotation> ${sim} $</annotation>\u0000 </semantics></math>50/50 sapphire beam-splitters. Most experiments show monotonic increases of shock temperature with pressure, from (168 GPa, 7100 K) to (243 GPa, 9400 K), in good agreement with predictions of our MgO B1 phase equation of state. Measured sound speeds are parallel to but <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>∼</mo>\u0000 </mrow>\u0000 <annotation> ${sim} $</annotation>\u0000 </semantics></math>10% higher than model predictions for bulk sound speed of solid B1 MgO, confirming ductile behavior of preheated MgO. Two experiments, at 238 and 246 GPa, showed anomalously low shock temperature and sound speed, suggesting melting. Using reported MgO melting data up to 120 GPa and our data at 232–246 GPa, we constructed a maximum-likelihood Simon-Glatzel fit. At Earth's core-mantle boundary pressure (135 GPa), our best-fit interpo","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142816028","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Harison S. Wiesman, Thomas Breithaupt, David Wallis, Lars N. Hansen
{"title":"Microstructural and Micromechanical Evolution of Olivine Aggregates During Transient Creep","authors":"Harison S. Wiesman, Thomas Breithaupt, David Wallis, Lars N. Hansen","doi":"10.1029/2024JB029812","DOIUrl":"10.1029/2024JB029812","url":null,"abstract":"<p>To examine the microstructural evolution that occurs during transient creep, we deformed samples of polycrystalline olivine to different strains that spanned the initial transient deformation. Two sets of samples with different initial grain sizes of 5 μm and 20 μm were deformed in torsion at <i>T</i> = 1,523 K, <i>P</i> = 300 MPa, and a constant shear strain rate of 1.5 × 10<sup>−4</sup> s<sup>−1</sup>, during which both sets of samples experienced strain hardening. We characterized the microstructures at the end of each experiment using high-angular resolution electron backscatter diffraction (HR-EBSD) and dislocation decoration. In the coarse-grained samples, dislocation density increased from 1.5 × 10<sup>11</sup> m<sup>−2</sup> to 3.6 × 10<sup>12</sup> m<sup>−2</sup> with strain. Although the same final dislocation density was reached in the fine-grained samples, it did not vary significantly at small strains, potentially due to concurrent grain growth during deformation. In both sets of samples, HR-EBSD analysis revealed that intragranular stress heterogeneity increased in magnitude with strain and that elevated stresses are associated with regions of high geometrically necessary dislocation density. Further analysis of the stresses and their probability distributions indicate that the stresses are imparted by dislocations and cause long-range elastic interactions among them. These characteristics indicate that dislocation interactions were the primary cause of strain hardening during transient creep in our samples. A comparison of the results to the predictions of three recent models reveals that the models do not correctly predict the evolution in stress and dislocation density with strain in our experiments due to a lack of previous such data in their calibrations.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2024JB029812","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142816029","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Permeability Development During Fault Growth and Slip in Granite","authors":"Franciscus M. Aben, Ado Farsi, Nicolas Brantut","doi":"10.1029/2024JB029057","DOIUrl":"10.1029/2024JB029057","url":null,"abstract":"<p>In tight crystalline rocks faults are known to be substantially more hydraulically conductive than the rock matrix. However, most of our knowledge relies on static measurements, or before/after failure data sets. The spatio-temporal evolution of the permeability field during faulting remains unknown. Here, we determine at which stage of faulting permeability changes most, and the degree of permeability heterogeneity along shear faults. We conducted triaxial deformation experiments in intact Westerly granite, where faulting was stabilized by monitoring acoustic emission rate. At repeated stages during deformation and faulting we paused deformation and imposed macroscopic fluid flow to characterize the overall permeability of the material. The pore pressure distribution was measured along the prospective fault to estimate apparent hydraulic transmissivity, and propagation of the macroscopic shear fault was monitored by locating acoustic emissions. We find that average permeability increases dramatically (by two orders of magnitude) with increasing deformation up to peak stress, where the fault is not yet through-going. Post-peak stress, overall permeability increases by a factor of three. However, at this stage we observed local heterogeneities in permeability by up to factors of eight, ascribed to a partially connected fracture network. This heterogeneity decreases with fault completion at residual shear stress. With further slip on the newly formed fault, the average hydraulic transmissivity remains mostly stable. Our results show that permeability enhancement during shear rupture mostly occurs ahead of the rupture tip, and that strongly heterogeneous permeability patterns are generated in the fault cohesive zone due to partial fracture connectivity.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2024JB029057","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142816465","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Teleseismic Full-Waveform Tomography Based on a 3D FK-LTSOS Hybrid Method for Local Imaging","authors":"Weijuan Meng, Dinghui Yang, Ling Chen, Xingpeng Dong, Xiaobing Xu, Jiandong Huang","doi":"10.1029/2024JB029544","DOIUrl":"10.1029/2024JB029544","url":null,"abstract":"<p>Teleseismic waveforms contain abundant interpretable information about Earth's properties. They can be used to explore the refined structure of Earth's interior, especially in the regions with imbalanced spatial distribution of seismic activity. However, it's technically infeasible to numerically simulate high-frequency (>1 Hz) teleseismic wave propagation within a whole domain iteratively in full-waveform inversion due to its vast computational costs. We develop a 3D FK-LTSOS (Frequency-Wavenumber, Layered Time-Space Optimized Symplectic) hybrid method and then apply it to teleseismic full-waveform tomography to tackle this computational challenge. The 3D FK-LTSOS hybrid method combines the semi-analytical solution computed by the FK method rapidly in a 1D background model and the numerical solution calculated by the 3D LTSOS method accurately in 3D localized heterogeneous media with topography to simulate teleseismic wave propagation efficiently and accurately. The comparison of synthetic seismograms shows its accuracy and stability when simulating wave propagation in the topographic model. Based on this hybrid method, the teleseismic full-waveform tomographic method is developed to efficiently resolve the detailed structure of the local research domain utilizing high-frequency teleseismic data. The essential contents of full-waveform tomography are presented, including misfit function, Fréchet kernels, smoothing strategy, and nonlinear conjugate gradient method. Synthetic data application for spherical anomaly models with planer surface and Gaussian topography, and observed data application for the crust-upper mantle structure beneath eastern Tibet confirm the validity of the teleseismic full-waveform tomography and demonstrate that our tomographic method can image the localized structures speedily from full-waveform information.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142797831","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Weiwei Zhu, Zhiqiang Chen, Xupeng He, Jingyao Liu, Songfeng Guo, Bowen Zheng, Ali Yousef, Shengwen Qi, Moran Wang
{"title":"Numerical Analysis of the Dynamic Mechanisms in Hydraulic Fracturing With a Focus on Natural Fractures","authors":"Weiwei Zhu, Zhiqiang Chen, Xupeng He, Jingyao Liu, Songfeng Guo, Bowen Zheng, Ali Yousef, Shengwen Qi, Moran Wang","doi":"10.1029/2024JB029487","DOIUrl":"10.1029/2024JB029487","url":null,"abstract":"<p>The hydraulic fracturing process is a prominent example of fracture network evolution under stress. However, the interactions between hydraulic fractures and natural fracture networks, along with the connectivity evolution of the resultant fracture networks, require more research. This research incorporates discrete fracture networks to characterize subsurface structures and employs the Discrete Element - Lattice Boltzmann Method to simulate the hydraulic fracturing process. The dynamic evolution of subsurface structures, including the initiation of hydraulic fractures and their interaction with natural fractures, is systematically investigated. Results indicate that natural fractures significantly impact fracture initiation, propagation, and connectivity evolution, which in turn affects fluid production. Fracture strength is key for the interaction, and hydraulic fractures tend to propagate along weak natural fractures with low resistance. Fracture strength variability determines the final fracture networks, with low-strength fractures breaking due to the altered in-situ stress and forming local clusters. High injection rates and fluid viscosity result in a large pressure buildup and exaggerate the influential region. A multi-cluster system is thus formed during the hydraulic fracturing process, and its connectivity can be well quantified with a novel connectivity metric. In low-permeable reservoirs, fracture clusters connected to production wells can contribute instantly, while local clusters may contribute to production from a long-term perspective. Injection rate, fluid viscosity, fracture orientation, and clustering effects have consistent positive correlations with the total connectivity and production. Heterogeneity has a weak positive correlation with fluid production, while a moderate negative correlation with total connectivity.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142797833","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Xi Liu, Fanfan Tian, Juan Li, Yang Li, Weidong Sun
{"title":"The Evolution of Dense ULVZs Originating Outside LLSVPs and Implications for Dynamics at LLSVP Margins","authors":"Xi Liu, Fanfan Tian, Juan Li, Yang Li, Weidong Sun","doi":"10.1029/2024JB028972","DOIUrl":"10.1029/2024JB028972","url":null,"abstract":"<p>Interactions between multiple-scale thermochemical heterogeneities in the lowermost mantle, specifically ultralow velocity zones (ULVZs) and large low shear velocity provinces (LLSVPs), are critical in lower mantle dynamics. However, the evolution of ULVZs formed outside LLSVPs has not been thoroughly explored. Here we perform two-dimensional numerical experiments to examine the evolution of highly dense ULVZs originating beneath cold downwellings and their interactions with the LLSVP. We find that ULVZs with an intrinsic density anomaly more than 500 kg/m<sup>3</sup> compared with the ambient lowermost mantle cannot fully enter the LLSVP and would dwell at LLSVP margins for an indefinitely long time. This suggests that dense ULVZ within LLSVPs might have different sources from those outside LLSVPs. The buoyancy number and compositional viscosity of ULVZs are controlling factors on their dynamics and imprints on the core-mantle boundary (CMB), such as how much the ULVZ protrudes into the LLSVP and the CMB topography beneath the ULVZ. The excess density of ULVZs dictates their width but not their thickness. The oscillations of ULVZ morphology suggest that various types of plumes occur at the LLSVP margin. The mobility of ULVZ implies that the bottom margin of the LLSVP moves much faster than its center. Hot zones exist within the LLSVP near its margins, which may affect the evolution of ULVZs and subducted material nearby. The CMB topography under dense ULVZs are positive unless the buoyancy number of ULVZs exceeds 6.0. These results have intriguing implications for the distribution of ULVZs as well as thermochemical evolution in the lowermost mantle.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142797828","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Tao Xia, Lingling Ye, Yefei Bai, Thorne Lay, Shiqing Xu, Hiroo Kanamori, Luis Rivera, Sesar Prabu Dwi Sriyanto
{"title":"The 2022 MW 7.3 Southern Sumatra Tsunami Earthquake: Rupture Up-Dip of the 2007 MW 8.4 Bengkulu Event","authors":"Tao Xia, Lingling Ye, Yefei Bai, Thorne Lay, Shiqing Xu, Hiroo Kanamori, Luis Rivera, Sesar Prabu Dwi Sriyanto","doi":"10.1029/2024JB030284","DOIUrl":"10.1029/2024JB030284","url":null,"abstract":"<p>On 18 November 2022, a large earthquake struck offshore southern Sumatra, generating a tsunami with 25 cm peak amplitude recorded at tide gauge station SBLT. Our <i>W</i>-phase solution indicates a shallow dip of 6.2°, compatible with long-period surface wave radiation patterns. Inversion of teleseismic body waves indicates a shallow slip distribution extending from about 10 km deep to near the trench with maximum slip of ∼4.1 m and seismic moment of <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>1.05</mn>\u0000 <mo>×</mo>\u0000 <msup>\u0000 <mn>10</mn>\u0000 <mn>20</mn>\u0000 </msup>\u0000 </mrow>\u0000 <annotation> $1.05times {10}^{20}$</annotation>\u0000 </semantics></math> Nm (<i>M</i><sub><i>W</i></sub> 7.3). Joint modeling of seismic and tsunami data indicates a shallow rigidity of ∼23 GPa. We find a low moment-scaled radiated energy of <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>4.15</mn>\u0000 <mo>×</mo>\u0000 <msup>\u0000 <mn>10</mn>\u0000 <mrow>\u0000 <mo>−</mo>\u0000 <mn>6</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation> $4.15times {10}^{-6}$</annotation>\u0000 </semantics></math>, similar to that of the 2010 <i>M</i><sub><i>W</i></sub> 7.8 Mentawai event (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>3.1</mn>\u0000 <mo>×</mo>\u0000 <msup>\u0000 <mn>10</mn>\u0000 <mrow>\u0000 <mo>−</mo>\u0000 <mn>6</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation> $3.1times {10}^{-6}$</annotation>\u0000 </semantics></math>) and other tsunami earthquakes. These characteristics indicate that the 2022 event should be designated as a smaller moment magnitude tsunami earthquake compared to the other 12 well-documented global occurrences since 1896. The 2022 event ruptured up-dip of the 2007 <i>M</i><sub><i>W</i></sub> 8.4 Bengkulu earthquake, demonstrating shallow seismogenic capability of a megathrust that had experienced both a deeper seismic event and adjacent shallow aseismic afterslip. We consider seismogenic behavior of shallow megathrusts and concern for future tsunami earthquakes in subduction zones globally, noting a correlation between tsunami earthquake occurrence and subducting seafloor covered with siliceous pelagic sediments. We suggest that the combination of pelagic clay and siliceous sediments and rough seafloor topography near the trench play important roles in controlling the genesis of tsunami earthquakes along Sumatra and other regions, rather than the subduction tectonic framework of accretionary or erosive mar","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142797832","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nicola Tisato, Kyle T. Spikes, Nishank Saxena, Ronny Hofmann
{"title":"Scattering and Frequency Effects on Ultrasonic Velocities of Carbonates","authors":"Nicola Tisato, Kyle T. Spikes, Nishank Saxena, Ronny Hofmann","doi":"10.1029/2024JB029491","DOIUrl":"10.1029/2024JB029491","url":null,"abstract":"<p>Scattering of elastic waves causes velocity dispersion, which increases uncertainty in seismic analysis. Understanding the sources of scattering and the degree of velocity dispersion are critical to improve subsurface imaging in efforts to locate resources and study subsurface processes. In addition to scattering, other mechanisms, such as the wave-induced fluid flow in saturated rocks cause velocity dispersion. To study the effect of scattering on velocity dispersion, we conducted laboratory measurements of ultrasonic velocities on dry rock samples and performed wave-propagation simulations on CT-scanned 3D volumes of those samples. The set of samples consists of homogeneous and heterogeneous carbonate rocks with porosities between 3% and 26%. Ultrasonic velocities were measured at frequencies between 0.3 and 1 MHz, and numerical wave propagation simulations on the digital volumes were performed using an elastic approximation and a finite-difference method. The homogeneous sample and the corresponding numerical simulations exhibit negligible velocity dispersion. On the other hand, heterogeneous samples exhibit significant dispersion, and the corresponding numerical simulations accurately reproduce the observed dispersion in terms of magnitude and frequency shift. We conclude that scattering has a first-order effect on the velocities of the elastic waves in heterogeneous samples. This effect should be considered in conjunction with laboratory measurements in heterogeneous carbonates similar to those studied here. Furthermore, we illustrate a method to characterize frequency-dependent ultrasonic velocities (i.e., dispersion) and show that finite-difference modeling can reproduce the laboratory-observed dispersion.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142797829","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sebastian Noe, Dirk-Philip van Herwaarden, Solvi Thrastarson, Marta Pienkowska, Neda Masouminia, Jincheng Ma, Hans-Peter Bunge, Deborah Wehner, Nicholas Rawlinson, Yajian Gao, Frederik Tilmann, Artie Rodgers, Andreas Fichtner
{"title":"The Collaborative Seismic Earth Model: Generation 2","authors":"Sebastian Noe, Dirk-Philip van Herwaarden, Solvi Thrastarson, Marta Pienkowska, Neda Masouminia, Jincheng Ma, Hans-Peter Bunge, Deborah Wehner, Nicholas Rawlinson, Yajian Gao, Frederik Tilmann, Artie Rodgers, Andreas Fichtner","doi":"10.1029/2024JB029656","DOIUrl":"10.1029/2024JB029656","url":null,"abstract":"<p>Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.</p>","PeriodicalId":15864,"journal":{"name":"Journal of Geophysical Research: Solid Earth","volume":"129 12","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2024JB029656","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142793504","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}