Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26Linked to original sources

Analysis of geophysical well logs obtained in the State 2-14 borehole, Salton Sea geothermal area, California

A complete suite of conventional geophysical well logs was obtained in the upper part of a 3220-m- deep borehole drilled into geothermally altered alluvial sediments on the southeastern edge of the Salton Sea. In situ temperatures greater than 300°C and an inability to cool parts of the borehole by circulation limited the suite of logs run below 2000 m in depth to deep induction, spontaneous potential, un-calibrated neutron, natural gamma, and temperature. Bottom-hole temperature trends given by repeat temperature logs were extrapolated to undisturbed temperatures approaching 355°C at a depth of 3220 m. Geophysical logs obtained in the State 2–14 borehole indicate that neutron porosity, gamma-gamma, and deep-induction logs provide useful information on lithologic trends with depth. The natural gamma log contains almost continuous, high-frequency fluctuations that obscure lithologic trends and that may be related to recent radioisotope redistribution and departure from radiometric equilibrium. Acoustic transit time logs give unrealistically low in situ compressional velocities ranging from 1.8 to 3.0 km/s, whereas acoustic waveform logs indicate that sediment compressional velocities range from less than 3.0 km/s shallower than 1000 m in depth to almost 5.0 km/s at depths greater than 2000 m. Analyses indicate that most log values lie between two lithologic end points: an electrically conductive claystone with moderate neutron porosity, but no effective porosity, and an electrically nonconductive, fully cemented siltstone that has small but finite porosity. A limited number of clean sandstones depart from this trend; geophysical logs from these sandstones indicate an effective porosity ranging from 5 to 14%, and saturation with brines having equivalent NaCl concentrations greater than 100,000 mg/L. Depth- averaged trends in neutron porosity and deep-induction logs, along with trends in acoustic velocity determined from acoustic waveform logs, demonstrate that major changes in the properties of alluvial sediments occur within the depth range from 1200 to 1800 m. Although caliper logs were not obtained deeper than 2000 m, resistivity values less than 10 ohm m at those depths probably correspond to borehole enlargements in production zones rather than local increases in effective porosity. The transition in sediment properties indicated by the geophysical logs in the depth interval from 1200 to 2000 m apparently represents a detailed vertical profile of the transition from relatively unaltered clay minerals in alluvial sediments to electrically nonconductive alteration products such as epidote and feldspar.

Journal of Geophysical Research Solid Earth↗

Geologic mapping of Europa

Galileo data enable the major geological units, structures, and surface features to be identified on Europa. These include five primary units (plains, chaos, band, ridge, and crater materials) and their subunits, along with various tectonic structures such as faults. Plains units are the most widespread. Ridged plains material spans a wide range of geological ages, including the oldest recognizable features on Europa, and appears to represent a style of tectonic resurfacing, rather than cryovolcanism. Smooth plains material typically embays other terrains and units, possibly as a type of fluid emplacement, and is among the youngest material units observed. At global scales, plains are typically mapped as undifferentiated plains material, although in some areas differences can be discerned in the near infrared which might be related to differences in ice grain size. Chaos material is composed of plains and other preexisting materials that have been severely disrupted by inferred internal activity; chaos is characterized by blocks of icy material set in a hummocky matrix. Band material is arrayed in linear, curvilinear, wedge-shaped, or cuspate zones with contrasting albedo and surface textures with respect to the surrounding terrain. Bilateral symmetry observed in some bands and the relationships with the surrounding units suggest that band material forms by the lithosphere fracturing, spreading apart, and infilling with material derived from the subsurface. Ridge material is mapped as a unit on local and some regional maps but shown with symbols at global scales. Ridge material includes single ridges, doublet ridges, and ridge complexes. Ridge materials are considered to represent tectonic processes, possibly accompanied by the extrusion or intrusion of subsurface materials, such as diapirs. The tectonic processes might be related to tidal flexing of the icy lithosphere on diurnal or longer timescales. Crater materials include various interior (smooth central, rough inner, and annular massif) and exterior (continuous ejecta) subunits. Structural features and landforms are shown with conventional symbols. Type localities for the units are identified, along with suggestions for portraying the features on geological maps, including colors and letter abbreviations for material units. Implementing these suggestions by the planetary mapping community would facilitate comparisons of maps for different parts of Europa and contribute to an eventual global synthesis of its complex geology. On the basis of initial mapping results, a stratigraphic sequence is suggested in which ridged plains form the oldest unit on Europa, followed by development of band material and individual ridges. Band materials tend to be somewhat older than ridges, but in many areas the two units formed simultaneously. Similarly, the formation of most chaos follows the development of ridged plains; although chaos is among the youngest materials on Europa, some chaos units might have formed contemporaneously with ridged plains. Smooth plains generally embay all other units and are late-stage in the evolution of the surface. C1 craters are superposed on ridged plains but are crosscut by other materials, including bands and ridges. Most c2 craters postdate all other units, but a few c2 craters are cut by ridge material. C3 craters constitute the youngest recognizable material on Europa. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Structure, stratigraphy, and origin of Husband Hill, Columbia Hills, Gusev Crater, Mars

The strike and dip of lithologic units imaged in stereo by the Spirit rover in the Columbia Hills using three-dimensional imaging software shows that measured dips (15-32??) for bedding on the main edifice of the Columbia Hill are steeper than local topography (???8-10??). Outcrops measured on West Spur are conformable in strike with shallower dips (7-15??) than observed on Husband Hill. Dips are consistent with observed strata draping the Columbia Hills. Initial uplift was likely related either to the formation of the Gusev Crater central peak or ring or through mutual interference of overlapping crater rims. Uplift was followed by subsequent draping by a series of impact and volcaniclastic materials that experienced temporally and spatially variable aqueous infiltration, cementation, and alteration episodically during or after deposition. West Spur likely represents a spatially isolated depositional event. Erosion by a variety of processes, including mass wasting, removed tens of meters of materials and formed the Tennessee Valley primarily after deposition. This was followed by eruption of the Adirondack-class plains basalt lava flows which embayed the Columbia Hills. Minor erosion, impact, and aeolian processes have subsequently modified the Columbia Hills. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Radar and photoclinometric studies of wrinkle ridges on Mars

Earth-based radar altimetry and image derived photoclinometric profiles were analyzed to examine both the long- and short-wavelength topography associated with wrinkle ridges on Mars. Photoclinometrically derived elevation data across wrinkle ridges were evaluated to determine the sensitivity of profiles to two empirical photoclinometric parameters, the horizontal digital number (HDN) and the scattered light value (SLV). The photoclinometric profiles are extremely sensitive to small variations in HDN. The sense of slope of a profile can be completely reversed over a range in HDN of as little as ±1. Comparably small variations in the SLV have relatively minor effects on the photoclinometrically derived elevations. The existence of elevation offsets from one side of the ridge to the other, reported in previous photoclinometric studies of martian wrinkles, were not confirmed through photoclinometry. In addition, no evidence of elevation offsets were found in Earth-based radar altimetry profiles across wrinkle ridges. In order to more accurately model wrinkle ridge topography, we controlled photoclinometrically derived elevations with long-wavelength topography obtained from the radar altimetry. The results of this study do not support kinematic models for the origin of planetary wrinkle ridges that involve deeply rooted thrust faults which separate crustal blocks at different elevations. A kinematic model involving buckling of shallow crustal layers into concentric folds that close, leading to the development of thrust faults, is consistent with wrinkle ridge morphology and terrestrial analogs. Recent geophysical studies of terrestrial analogs and the influence of shallow subsurface structures, particularly buried craters, on the localization of many wrinkle ridges on Mars suggest that thrust faults associated with the ridges are confined to the ridged plains material and do not extend into the lithosphere.

Journal of Geophysical Research E: Planets↗

Flow of variably fluidized granular masses across three-dimensional terrain 2. Numerical predictions and experimental tests

Numerical solutions of the equations describing flow of variably fluidized Coulomb mixtures predict key features of dry granular avalanches and water-saturated debris flows measured in physical experiments. These features include time-dependent speeds, depths, and widths of flows as well as the geometry of resulting deposits. Threedimensional (3-D) boundary surfaces strongly influence flow dynamics because transverse shearing and cross-stream momentum transport occur where topography obstructs or redirects motion. Consequent energy dissipation can cause local deceleration and deposition, even on steep slopes. Velocities of surge fronts and other discontinuities that develop as flows cross 3-D terrain are predicted accurately by using a Riemann solution algorithm. The algorithm employs a gravity wave speed that accounts for different intensities of lateral stress transfer in regions of extending and compressing flow and in regions with different degrees of fluidization. Field observations and experiments indicate that flows in which fluid plays a significant role typically have high-friction margins with weaker interiors partly fluidized by pore pressure. Interaction of the strong perimeter and weak interior produces relatively steep-sided, flat-topped deposits. To simulate these effects, we compute pore pressure distributions using an advection-diffusion model with enhanced diffusivity near flow margins. Although challenges remain in evaluating pore pressure distributions in diverse geophysical flows, Riemann solutions of the depthaveraged 3-D Coulomb mixture equations provide a powerful tool for interpreting and predicting flow behavior. They provide a means of modeling debris flows, rock avalanches, pyroclastic flows, and related phenomena without invoking and calibrating Theological parameters that have questionable physical significance.

Journal of Geophysical Research B: Solid Earth↗

Compact Reconnaissance Imaging Spectrometer for Mars investigation and data set from the Mars Reconnaissance Orbiter's primary science phase

The part of the Compact Reconnaissance Imaging Spectrometer (CRISM) for Mars investigation conducted during the Mars Reconnaissance Orbiter's (MRO's) primary science phase was a comprehensive investigation of past aqueous environments, structure of the planet's crust, past climate, and current meteorology. The measurements to implement this investigation include over 9500 targeted observations of surface features taken at spatial resolutions of better than 40 m/pixel, monitoring of seasonal variations in atmospheric aerosols and trace gases, and acquisition of a 200 m/pixel map covering over 55% of Mars in 72 selected wavelengths under conditions of relatively low atmospheric opacity. Key results from these data include recognition of a diversity of aqueous mineral-containing deposits, discovery of a widespread distribution of phyllosilicates in early to middle Noachian units, the first definitive detection of carbonates in bedrock, new constraints on the sequence of events that formed Hesperian-aged, sulfate-rich layered deposits, characterization of seasonal polar processes, and monitoring of the 2007 global dust event. Here we describe CRISM's science investigations during the Primary Science Phase, the data sets that were collected and their calibration and uncertainties, and how they have been processed and made available to the scientific community. We also describe the ongoing investigation during MRO's extended science phase. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Residence time controls on the fate of nitrogen in flow‐through lakebed sediments

For many glacial lakes with highly permeable sediments, water exchange rates control hydrologic residence times within the sediment‐water interface (SWI) and the removal of reactive compounds such as nitrate, a common pollutant in lakes and groundwater. Here we conducted a series of focused tracer injection experiments in the upper 20 cm of the naturally downwelling SWI in a flow‐through lake on Cape Cod, MA. We systematically varied residence time and reactant controls on nitrate processing, using isotopically labeled 15 N nitrate to monitor the effect of these changes on nitrate removal via denitrification. The addition of acetate, a labile carbon compound, triggered the lake SWI to switch from net production to net removal of nitrate. When acetate was combined with increased residence time created by controlled reductions in water flux, we observed a fivefold increase in nitrate removal, a 26‐fold increase in N 2 production, and a 42‐fold increase in N 2 O production. We demonstrate that water residence time is an important control on the fate of nitrate in these lake SWIs and illustrate that seasonal conditions that alter lake exchange rates and variability in lake carbon may predict dynamic nitrate removal across the SWI. Additionally, observed N 2 O production during the oxic pore water experiments paired with geophysical characterization of the sediment porosity revealed that the lake SWI has less mobile pores occupying upward of 50% of the total porosity volume, which function as reactive microzones for nitrate processing.

Journal of Geophysical Research: Biogeosciences↗

Temporal stress changes caused by earthquakes: A review

Earthquakes can change the stress field in the Earth’s lithosphere as they relieve and redistribute stress. Earthquake-induced stress changes have been observed as temporal rotations of the principal stress axes following major earthquakes in a variety of tectonic settings. The stress changes due to the 2011 Mw9.0 Tohoku-Oki, Japan, earthquake were particularly well documented. Earthquake stress rotations can inform our understanding of earthquake physics, most notably addressing the long-standing problem of whether the Earth’s crust at plate boundaries is “strong” or “weak.” Many of the observed stress rotations, including that due to the Tohoku-Oki earthquake, indicate near-complete stress drop in the mainshock. This implies low background differential stress, on the order of earthquake stress drop, supporting the weak crust model. Earthquake stress rotations can also be used to address other important geophysical questions, such as the level of crustal stress heterogeneity and the mechanisms of postseismic stress reloading. The quantitative interpretation of stress rotations is evolving from those based on simple analytical methods to those based on more sophisticated numerical modeling that can capture the spatial-temporal complexity of the earthquake stress changes.

Journal of Geophysical Research B: Solid Earth↗

Mg-spinel lithology: A new rock type on the lunar farside

High-resolution compositional data from Moon Mineralogy Mapper (M 3) for the Moscoviense region on the lunar farside reveal three unusual, but distinctive, rock types along the inner basin ring. These are designated "OOS" since they are dominated by high concentrations of orthopyroxene, olivine, and Mg-rich spinel, respectively. The OOS occur as small areas, each a few kilometers in size, that are widely separated within the highly feldspathic setting of the basin rim. Although the abundance of plagioclase is not well constrained within the OOS, the mafic mineral content is exceptionally high, and two of the rock types could approach pyroxenite and harzburgite in composition. The third is a new rock type identified on the Moon that is dominated by Mg-rich spinel with no other mafic minerals detectable (<5% pyroxene, olivine). All OOS surfaces are old and undisturbed since basin formation. They are effectively invisible in image data and are only recognized by their distinctive composition identified spectroscopically. The origin of these unusual lithologies appears to be linked to one or more magmatic intrusions into the lower crust, perhaps near the crust-mantle interface. Processes such as fractional crystallization and gravity settling within such intrusions may provide a mechanism for concentrating the mafic components within zones several kilometers in dimension. The OOS are embedded within highly anorthositic material from the lunar crust; they may thus be near contemporaneous with crustal products from the cooling magma ocean. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Late Cenozoic volcanism, geochronology, and structure of the Coso Range, Inyo County, California

The Coso Range lies at the west edge of the Great Basin, adjacent to the southern part of the Sierra Nevada. A basement complex of pre‐Cenozoic plutonic and metamorphic rocks is partly buried by ∼35 km 3 of late Cenozoic volcanic rocks that were erupted during two periods, as defined by K‐Ar dating: (1) 4.0–2.5 m.y., ∼31 km 3 of basalt, rhyodacite, dacite, andesite, and rhyolite, in descending order of abundance, and (2) ≤1.1 m.y., nearly equal amounts of basalt and rhyolite, most of the rhyolite being ≤0.3 m.y. old. Vents for the volcanic rocks of the younger period are localized on and near a horst of basement rocks within a concavity defined by the distribution of vents of the older period. The alignment of many vents and the presence of a considerable number of roughly north‐trending normal faults of late Cenozoic age reflect basin and range tectonics dominated by roughly east‐west lithospheric extension. Fumaroles, intermittently active thermal springs, and associated altered rocks occur within and immediately east of the central part of the field of Quaternary rhyolite, in an area characterized by various geophysical anomalies that are evidently related to an active hot‐water geothermal system. This system apparently is heated by a reservoir of silicic magma at ≥8‐km depth, itself produced and sustained through partial melting of crustal rocks by thermal energy contained in mantle‐derived basaltic magma that intrudes the crust in response to lithospheric extension.

California↗

Evidence for a Battle Mountain-Eureka crustal fault zone, north-central Nevada, and its relation to Neoproterozoic-Early Paleozoic continental breakup

Combined evidence from gravity, radiogenic isotope, and magnetotelluric (MT) data indicates a crustal fault zone that coincides with the northwest-trending Battle Mountain-Eureka (BME) mineral trend in north-central Nevada, USA. The BME crustal fault zone likely originated during Neoproterozoic-Early Paleozoic rifting of the continent and had a large influence on subsequent tectonic events, such as emplacement of allochthons and episodic deformation, magmatism, and mineralization throughout the Phanerozoic. MT models show the fault zone is about 10 km wide, 130-km long, and extends from 1 to 5 km below the surface to deep crustal levels. Isotope data and gravity models imply the fault zone separates crust of fundamentally different character. Geophysical evidence for such a long-lived structure, likely inherited from continental breakup, defies conventional wisdom that structures this old have been destroyed by Cenozoic extensional processes. Moreover, the coincidence with the alignment of mineral deposits supports the assertion by many economic geologists that these alignments are indicators of buried regional structures.

Journal of Geophysical Research B: Solid Earth↗

GCIP water and energy budget synthesis (WEBS)

As part of the World Climate Research Program's (WCRPs) Global Energy and Water-Cycle Experiment (GEWEX) Continental-scale International Project (GCIP), a preliminary water and energy budget synthesis (WEBS) was developed for the period 1996-1999 fromthe "best available" observations and models. Besides this summary paper, a companion CD-ROM with more extensive discussion, figures, tables, and raw data is available to the interested researcher from the GEWEX project office, the GAPP project office, or the first author. An updated online version of the CD-ROM is also available at http://ecpc.ucsd.edu/gcip/webs.htm/. Observations cannot adequately characterize or "close" budgets since too many fundamental processes are missing. Models that properly represent the many complicated atmospheric and near-surface interactions are also required. This preliminary synthesis therefore included a representative global general circulation model, regional climate model, and a macroscale hydrologic model as well as a global reanalysis and a regional analysis. By the qualitative agreement among the models and available observations, it did appear that we now qualitatively understand water and energy budgets of the Mississippi River Basin. However, there is still much quantitative uncertainty. In that regard, there did appear to be a clear advantage to using a regional analysis over a global analysis or a regional simulation over a global simulation to describe the Mississippi River Basin water and energy budgets. There also appeared to be some advantage to using a macroscale hydrologic model for at least the surface water budgets. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Evidence for nonuniform permafrost degradation after fire in boreal landscapes

Fire can be a significant driver of permafrost change in boreal landscapes, altering the availability of soil carbon and nutrients that have important implications for future climate and ecological succession. However, not all landscapes are equally susceptible to fire-induced change. As fire frequency is expected to increase in the high latitudes, methods to understand the vulnerability and resilience of different landscapes to permafrost degradation are needed. We present a combination of multiscale remote sensing, geophysical, and field observations that reveal details of both near-surface (<1 m) and deeper (>1 m) impacts of fire on permafrost. Along 11 transects that span burned-unburned boundaries in different landscape settings within interior Alaska, subsurface electrical resistivity and nuclear magnetic resonance data indicate locations where permafrost appears to be resilient to disturbance from fire, areas where warm permafrost conditions exist that may be most vulnerable to future change, and also areas where permafrost has thawed. High-resolution geophysical data corroborate remote sensing interpretations of near-surface permafrost and also add new high-fidelity details of spatial heterogeneity that extend from the shallow subsurface to depths of about 10 m. Results show that postfire impacts on permafrost can be variable and depend on multiple factors such as fire severity, soil texture, soil moisture, and time since fire.

Alaska↗

Externally triggered renewed bubble nucleation in basaltic magma: the 12 October 2008 eruption at Halema‘uma‘u Overlook vent, Kīlauea, Hawai‘i, USA

From October 2008 until present, dozens of small impulsive explosive eruptions occurred from the Overlook vent on the southeast side of Halema&lsquo;uma&lsquo;u Crater, at Kīlauea volcano, USA. These eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column. Here we use microtextural observations and data from clasts erupted during the well-characterized 12 October 2008 explosive eruption at Halema&lsquo;uma&lsquo;u to extend existing models of eruption triggering. We present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets. We measure the size and number density of bubbles preserved in juvenile ejecta using 2D images and X-ray microtomography. Our data suggest that accumulations of large bubbles with diameters of >50 &mu; m to at least millimeters existed at shallow levels within the conduit prior to the 12 October 2008 explosion. Furthermore, a high number density of small bubbles <50 &mu; m is measured in the clasts, implying very rapid nucleation of bubbles. Visual observations, combined with preexisting geophysical data, suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that (1) nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and (2) expanded the preexisting bubbles driving upward acceleration. The trigger of eruption and bubble nucleation is thus external to the degassing system.

Hawaii↗

Microgravity measurements spanning 2019–2022 indicate no change in magma supply or storage during Kilauea's reawakening

After more than 2 years of quiescence following its major 2018 flank lava effusion and summit caldera collapse, Kīlauea, Hawaiʻi, returned to frequent eruptive activity in late 2020. Microgravity measurements from the summit region in 2019 and 2022 span this reawakening and indicate a gravity increase of up to hundreds of microgals near the center of the summit caldera. Much of this increase is caused by lava that partially filled the 2018 collapse pit. We accounted for this mass by differencing digital elevation models of the summit area to determine the geometry and amount of accumulated lava and its gravity contribution at each station in the network. The remaining gravity change is best approximated by a mass increase of 6.9 × 10 10 kg (5.1–23.0 × 10 10 kg at 95% credible interval) at a depth of 480 m (121–795 m at 95% credible interval) beneath the locus of the 2018 collapse. The modeled source depth reflects the combined impacts of magma accumulation in a reservoir known from previous geophysical studies to exist at about 1–2 km depth, intrusions between that reservoir and the surface, and mass accumulation in the heavily fractured caldera floor. The spatial and temporal patterns of subsurface mass accumulation during 2019–2022 imply that there were no significant changes in magma reservoir geometry or magma supply to the volcano in the years following the 2018 collapse.

Hawaii↗

Does Europa have a subsurface ocean? Evaluation of the geological evidence

It has been proposed that Jupiter's satellite Europa currently possesses a global subsurface ocean of liquid water. Galileo gravity data verify that the satellite is differentiated into an outer H 2 O layer about 100 km thick but cannot determine the current physical state of this layer (liquid or solid). Here we summarize the geological evidence regarding an extant subsurface ocean, concentrating on Galileo imaging data. We describe and assess nine pertinent lines of geological evidence: impact morphologies, lenticulae, cryovolcanic features, pull‐apart bands, chaos, ridges, surface frosts, topography, and global tectonics. An internal ocean would be a simple and comprehensive explanation for a broad range of observations; however, we cannot rule out the possibility that all of the surface morphologies could be due to processes in warm, soft ice with only localized or partial melting. Two different models of impact flux imply very different surface ages for Europa; the model favored here indicates an average age of ∼50 Myr. Searches for evidence of current geological activity on Europa, such as plumes or surface changes, have yielded negative results to date. The current existence of a global subsurface ocean, while attractive in explaining the observations, remains inconclusive. Future geophysical measurements are essential to determine conclusively whether or not there is a liquid water ocean within Europa today.

Journal of Geophysical Research: Planets↗

Introduction and digest to the Special Issue on Chemical Effects of Water on the Deformation and Strengths of Rocks

The important role of pore pressure in promoting such brittle processes as cataclasis, hydraulic fracturing, large‐scale faulting, and earthquakes within the crust is widely accepted in geology and geophysics [, 1957; , 1959; , 1958; ., 1963; , 1968; ., 1968; ., 1976; , 1973, 1980; , 1981]. Provided that fluid pressure is fully communicated with rock pore space, the effective normal stresses that control crack growth, macroscopic fracture, and friction are reduced by the magnitude of the fluid pressui'e. Beyond this physical effect of pore fluids, there are chemical effects of water on the strength of rocks that are also important in governing differential stresseg and flow in the continental crust. Some of these chemical effects of water on rock deformation have long been recognized.

Journal of Geophysical Research B: Solid Earth↗

Damped regional-scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence

We present a new focal mechanism stress inversion technique to produce regional-scale models of stress orientation containing the minimum complexity necessary to fit the data. Current practice is to divide a region into small subareas and to independently fit a stress tensor to the focal mechanisms of each subarea. This procedure may lead to apparent spatial variability that is actually an artifact of overfitting noisy data or nonuniquely fitting data that does not completely constrain the stress tensor. To remove these artifacts while retaining any stress variations that are strongly required by the data, we devise a damped inversion method to simultaneously invert for stress in all subareas while minimizing the difference in stress between adjacent subareas. This method is conceptually similar to other geophysical inverse techniques that incorporate damping, such as seismic tomography. In checkerboard tests, the damped inversion removes the stress rotation artifacts exhibited by an undamped inversion, while resolving sharper true stress rotations than a simple smoothed model or a moving-window inversion. We show an example of a spatially damped stress field for southern California. The methodology can also be used to study temporal stress changes, and an example for the Coalinga, California, aftershock sequence is shown. We recommend use of the damped inversion technique for any study examining spatial or temporal variations in the stress field.

Journal of Geophysical Research B: Solid Earth↗