SALSA3D: Updated Tomographic Velocity Models for Improved Travel-Time Prediction and Uncertainty .
Abstract not provided.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Abstract not provided.
This poster will be shown at the 2023 Seismological Society of America conference in San Juan, Puerto Rico. The audience will consist of international members of academia and U.S. lab and government members. Due to time limits, full review is being undertaken in parallel with HQ review. The full review will not be accepted until HQ approval comes through.
Abstract Location algorithms have historically relied on simple, one-dimensional (1D) velocity models for fast seismic event locations. 1D models are generally used as travel-time lookup tables, one for each seismic phase, with travel-times pre-calculated for event distance and depth. These travel-time lookup tables are extremely fast to use and this fast computational speed makes them the preferred type of velocity model for operational needs. Higher-dimensional (i.e., three-dimensional—3D) seismic velocity models are becoming readily available and provide more accurate event locations over 1D models. The computational requirements of these 3D models tend to make their operational use prohibitive. Additionally, comparing location accuracy for 3D seismic velocity models tends to be problematic, as each model is determined using different ray-tracing algorithms. Attempting to use a different algorithm than the one used to develop a model usually results in poor travel-time prediction. We demonstrate and test a framework to create first-P and first-S 3D travel-time correction surfaces using an open-source framework ( PCalc + GeoTess , https://www.sandia.gov/salsa3d/software/geotess ) that easily stores 3D travel-time and uncertainty data. This framework produces fast travel-time and uncertainty predictions and overcomes the ray-tracing algorithm hurdle because the lookup tables can be generated using the exact ray-tracing algorithm that is preferred for a model.
Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.
Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.
Leaching of excess nitrogen (N) to groundwater in fertilized landscapes can overwhelm natural biogeochemical processes and cause long-term eutrophication of aquatic systems. We investigated N fate and transport from an intensively managed short-rotation woody crop (Pinus taeda) plantation through the riparian zone of an intermittent, low-gradient blackwater stream. Fertilization of the P. taeda plantation on the uplands resulted in contamination of groundwater with nitrate concentrations between 0.9 and 1.9 mg N L -1 . No corresponding increase in nitrate was observed in stream water or shallow groundwater in the riparian zone. Groundwater travel-time modeling predicted that N from near-stream, upland plantation areas should have reached streams during the monitoring period. Two years of measuring N species in well water in contrasting landscape positions (within the plantation, swale, riparian edge, forested hillslope, and valley), indicated rapid nitrate transformation and denitrification within the forested wetland valleys. Denitrification in the shallow groundwater system within the toeslopes and the riparian zone was estimated to have removed > 90% of nitrate. Finally, these results highlight the importance of riparian zones as pathways for the removal of N and for controlling downstream N loads.
In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.
Recently, Duvall and Hanasoge (Solar Phys. 287, 71, 2013) found that large distance separation [delta] travel-time differences from a center to an annulus [deltat(sub oi)] implied a model of the average super granular cell that has a peak upflow of 240 ms(exp −1) at a depth of 2.3 Mm and a corresponding peak outward horizontal flow of 700 ms(exp −1) at a depth of 1.6 Mm. In the present work, this effect is further studied by measuring and modeling center-to-quadrant travel-time differences [deltat(sub qu)], which roughly agree with this model. Simulations are analyzed that show that such a model flow would lead to the expected travel-time differences. As a check for possible systematic errors, the center-to-annulus travel-time differences [deltat(sub oi)] are found not to vary with heliocentric angle. A consistency check finds an increase of deltat(sub oi) with the temporal frequency [ν] by a factor of two, which is not predicted by the ray theory.
We report the results of simultaneous Chandra and RXTE observations of the Seyfert 1 galaxy Mkn 509. We deconvolve the broad and narrow Fe-K emission-line components for which we measure rest-frame equivalent widths of 119+/-18 eV and 57+/-13 eV respectively. The broad line has a FWHM of 57,600((sup 14,400)(sub -21,000)) km/s and the narrow line is unresolved, with an upper limit on the FWHM of 4,940 km/s. Both components must originate in cool matter since we measure rest-frame center energies of 6.36((sup +0.13)(sub -0.12)) keV and 6.42+/-0.01 keV for the broad and narrow line respectively. This rules out He-like and H-like Fe for the origin of both the broad and narrow lines. If, as is widely accepted, the broad Fe-K line originates in Thomson-thick matter (such as an accretion disk), then one expects to observe spectral curvature above approximately 10 keV, (commensurate with the observed broad line), characteristic of the Compton-reflection continuum. However our data sets very stringent limits on deviations of the observed continuum from a power law. Light travel-time delays cannot be invoked to explain anomalies in the relative strengths of the broad Ferry line and Compton-reflection continuum since they are supposed to originate in the same physical location. We are forced to conclude that both the broad and narrow Fe-K lines had to originate in Thomson-thin matter during our observation. This result, for a single observation of just one source, means that our understanding of Fe K line emission and Compton reflection from accreting X-ray sources in general needs to be re-examined. For example, if an irradiated accretion disk existed in Mkn 509 at the time of the observations, the lack of spectral curvature above approximately 10 keV suggests two possibilities. Either the disk was Thomson-thick and highly ionized, having negligible Fe-K line emission and photoelectric absorption or the disk was Thomson-thin producing some or all of the broad Fe-K line emission. In the former case, the broad Fe-K line had to have produced in a Thomson-thin region elsewhere. In both cases the predicted spectral curvature above approximately 10 keV is negligible. An additional implication of our results is that any putative obscuring torus in the system, required by unification models of active galaxies, must also be Thomson-thin. The same applies to the optical broad line region (BLR) if it has a substantial covering factor.