Editorial: Exploring new frontiers with JGR–Solid Earth
No abstract available.
Journal of Geophysical Research B: Solid Earth↗
SEARCH · USGS Science
Results for “JGR - Atmospheres”
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.
Introduction to special section on the California-Arizona crustal transect: CACTIS, Part 3
The CACTIS (California‐Arizona Crustal Transect Interim Synthesis) workshop in May 1988 brought researchers together in Flagstaff, Arizona, to discuss the geologic evolution and crustal structure of the southern Cordillera between the San Andreas fault in southeastern California and the Colorado Plateau in Arizona [ Sass et al ., 1988]. The first set of papers resulting from the workshop appeared in the Journal of Geophysical Research (JGR) as the special CACTIS 1 section in January 1990 [ Howard et al ., 1990], and a second as CACTIS 2 in November 1990 [ Haxel et al ., 1990]. Twelve papers in this issue form part 3 of the CACTIS series. A final collection of papers will appear in the future as a joint special section with the California Consortium for Crustal Studies (CALCRUST).
Journal of Geophysical Research B: Solid Earth↗
Inundation of a barrier island (Chandeleur Islands, Louisiana, USA) during a hurricane: Observed water-level gradients and modeled seaward sand transport
Large geomorphic changes to barrier islands may occur during inundation, when storm surge exceeds island elevation. Inundation occurs episodically and under energetic conditions that make quantitative observations difficult. We measured water levels on both sides of a barrier island in the northern Chandeleur Islands during inundation by Hurricane Isaac. Wind patterns caused the water levels to slope from the bay side to the ocean side for much of the storm. Modeled geomorphic changes during the storm were very sensitive to the cross-island slopes imposed by water-level boundary conditions. Simulations with equal or landward sloping water levels produced the characteristic barrier island storm response of overwash deposits or displaced berms with smoother final topography. Simulations using the observed seaward sloping water levels produced cross-barrier channels and deposits of sand on the ocean side, consistent with poststorm observations. This sensitivity indicates that accurate water-level boundary conditions must be applied on both sides of a barrier to correctly represent the geomorphic response to inundation events. More broadly, the consequence of seaward transport is that it alters the relationship between storm intensity and volume of landward transport. Sand transported to the ocean side may move downdrift, or aid poststorm recovery by moving onto the beach face or closing recent breaches, but it does not contribute to island transgression or appear as an overwash deposit in the back-barrier stratigraphic record. The high vulnerability of the Chandeleur Islands allowed us to observe processes that are infrequent but may be important at other barrier islands.
Louisiana↗
Euler-vector clustering of GPS velocities defines microplate geometry in southwest Japan
I have used Euler-vector clustering to assign 469 GEONET stations in southwest Japan to k clusters ( k = 2, 3,..., 9) so that, for any k , the velocities of stations within each cluster are most consistent with rigid-block motion on a sphere. That is, I attempt to explain the raw (i.e., uncorrected for strain accumulation), 1996–2006 velocities of those 469 Global Positioning System stations by rigid motion of k clusters on the surface of a spherical Earth. Because block geometry is maintained as strain accumulates, Euler-vector clustering may better approximate the block geometry than the values of the associated Euler vectors. The microplate solution for each k is constructed by merging contiguous clusters that have closely similar Euler vectors. The best solution consists of three microplates arranged along the Nankaido Trough-Ryukyu Trench between the Amurian and Philippine Sea Plates. One of these microplates, the South Kyushu Microplate (an extension of the Ryukyu forearc into the southeast corner of Kyushu), had previously been identified from paleomagnetic rotations. Relative to ITRF2000 the three microplates rotate at different rates about neighboring poles located close to the northwest corner of Shikoku. The microplate model is identical to that proposed in the block model of Wallace et al. (2009, https://doi.org/10.1130/G2522A.1) except in southernmost Kyushu. On Shikoku and Honshu, but not Kyushu, the microplate model is consistent with that proposed in the block models of Nishimura and Hashimoto (2006, https://doi.org/10.1016/j.tecto.2006.04.017) and Loveless and Meade (2010, https://doi.org/10.1029/2008JB006248) without the low-slip-rate boundaries proposed in the latter.
southwest Japan↗
A newly identified role of the deciduous forest floor in the timing of green‐up
Plant phenology studies rarely consider controlling factors other than air temperature. We evaluate here the potential significance of physical and chemical properties of soil (edaphic factors) as additional important controls on phenology. More specifically, we investigate causal connections between satellite‐observed green‐up dates of small forest watersheds and soil properties in the Adirondack Mountains of New York, USA. Contrary to the findings of previous studies, where edaphic controls of spring phenology were found to be marginal, our analyses show that at least three factors manifest themselves as significant controls of seasonal patterns of variation in vegetated land surfaces observed from remote sensing: (1) thickness of the forest floor, (2) concentration of exchangeable soil potassium, and (3) soil acidity. For example, a thick forest floor appears to delay the onset of green‐up. Watersheds with elevated concentrations of potassium are associated with early surface greening. We also found that trees growing in strongly acidified watersheds demonstrate delayed green‐up dates. Overall, our work demonstrates that, at the scale of small forest watersheds, edaphic factors can explain a significant percentage of the observed spatial variation in land surface phenology that is comparable to the percentage that can be explained by climatic and landscape factors. We conclude that physical and chemical properties of forest soil play important roles in forest ecosystems as modulators of climatic drivers controlling the rate of spring soil warming and the transition of trees out of winter dormancy.
New York↗
Revisiting earthquakes in the Los Angeles, California, basin during the early instrumental period: Evidence for an association with oil production
A total of seven independent M L ≥ 4.0 earthquakes occurred in the Los Angeles, California, basin, during the early instrumental period between 1932 and 1952, the largest of which was the 1933 Long Beach earthquake. Revising available macroseismic and instrumental data for a total of 6 4.0 ≤ M L ≤ 5.1 events between 1938 and 1944, we conclude that early instrumental locations can be grossly inconsistent with detailed macroseismic data. We use available macroseismic data to revisit event locations. We further present evidence that most if not all of these moderate earthquakes may have been induced by oil production. We quantify the predicted stress change associated with production from eight oil fields in the southwestern Los Angeles basin and show that frictional failure would have been encouraged beneath and at the periphery of high-volume fields, with stress changes upward of 0.1 MPa at 5-km depth. The results suggest that if earthquakes are induced by stress changes associated with production, the magnitudes of events might tend to be limited by the limited spatial extent of lobes of increased Coulomb failure stress. It further appears that the advent of fluid injection recovery methods ( water-flooding ) around 1960 mitigated induced earthquake risk considerably.
California↗
Less fine particle retention in a restored versus unrestored urban stream: Balance between hyporheic exchange, resuspension, and immobilization
Stream restoration goals include reducing erosion and increasing hyporheic exchange to promote biogeochemical processing and improve water quality. Little is known, however, about fine particle dynamics in response to stream restoration. Fine particles (<100 μm) are exchanged with transient storage areas near and within streambeds and banks. Fine particle retention directly impacts carbon and nutrient cycling by supporting benthic and hyporheic primary production, but overaccumulation of fine particle deposits can impair metabolism by burying benthic biofilms and reducing streambed permeability. We analyzed the transport and retention of water and fine particles at both the reach and local scales in a restored urban stream, 9 years postrestoration. We injected conservative solute and fine particle tracers under summer baseflow conditions and monitored their distribution between surface water, porewaters, and storage areas (i.e., biofilms, hyporheic zones, and slow surface waters). Comparison of the results to a nearby unrestored stream demonstrate that the restored reach had 10–45 times greater exchange of fine particles with transient storage zones, but 5 times lower rate of net particle immobilization. Local-scale results showed that restoration increased fine particle exchange with short-term storage areas but did not increase long-term particle retention. Thus, the restored stream rapidly exchanged fine sediments with transient storage areas, but did not store fine sediments as efficiently as the unrestored stream. The decreased retention of particulate organic matter in the restored stream may reduce biogeochemical processes, such as denitrification, by not providing sufficient organic carbon or the surface area required for microbial colonization.
JGR Biogeosciences↗
Morphodynamic resilience of intertidal mudflats on a seasonal time scale
Intertidal mudflats are morphodynamic features present in many estuaries worldwide. Often located between vegetated shores and deep channels they comprise valuable ecosystems and serve to protect the hinterland by attenuating waves. Although mudflats are persistently present on yearly to decadal time scales, little is known on their morphodynamic adaptation to short‐term variations in forcing such as storms, spring‐neap tidal cycles, and sediment supply. This study aims to explore the morphodynamic resilience of mudflats to seasonal variations in forcing. First, we compare transects observed in South Bay, California, at 3‐ to 6‐monthly intervals. Second, we present the results of a process‐based, morphodynamic profile model (Mflat). Mflat is an open source, Matlab code that describes both cross‐shore and alongshore tidal hydrodynamics as well as a stationary wave model. An advection‐diffusion equation solves sediment transport while bed level changes occur by the divergence of the sediment transport field. Mflat reproduces the observed South San Francisco Bay profile in equilibrium with significant skill. Short‐term variations in hydrodynamic forcing and sediment characteristics disturb the profile mainly at the channel‐shoal edge. The modeled profile disturbance is consistent with observations. The modeled profile is remarkably resilient since it recovers to the equilibrium profile within weeks to months. The model results suggest that 3‐monthly observation intervals are probably too long to discriminate processes responsible for the profile disturbance. These processes may include variations in sediment supply, mudflat erodibility, and wave action as well as the spring‐neap tidal cycle.
California↗
Diagenesis revealed by fine-scale features at Vera Rubin ridge, Gale crater, Mars
Fine-scale (submillimeter to centimeter) depositional and diagenetic features encountered during the Curiosity rover's traverse in Gale crater provide a means to understand the geologic history of Vera Rubin ridge (VRR). VRR is a topographically high feature on the lower north slope of Aeolis Mons, a 5-km high stratified mound within Gale crater. We use high-spatial resolution images from the Mars Hand Lens Imager (MAHLI) as well as grain sizes estimated with the Gini index mean score technique that uses ChemCam Laser-Induced Breakdown Spectroscopy (LIBS) chemical data to constrain the postdepositional history of the strata exposed on this ridge. MAHLI images were used to examine the color, grain size, and style of lamination of the host rocks, as well as to explore the occurrence of nodules, diagenetic crystals, pits, and a variety of dark-gray iron-rich features. This survey revealed abundant and widespread diagenetic features within the rocks exposed on VRR and demonstrated that rock targets estimated to be coarser generally contain more diagenetic features than those estimated to have finer grains, which indicate that grain size may have influenced the degree and type of diagenesis. A subset of rocks within VRR are gray in color and exhibit the highest proportion of diagenetic features. We suggest that these targets experienced a different diagenetic history than the other rocks on VRR and hypothesize that redistribution and recrystallization of iron within specific intervals may have resulted in both the gray color and the abundance of dark-gray iron-rich diagenetic features.
JGR Planets↗
Diagenesis of Vera Rubin ridge, Gale crater, Mars from Mastcam multispectral images
Images from the Mars Science Laboratory (MSL) mission of lacustrine sedimentary rocks of Vera Rubin ridge on “Mt. Sharp” in Gale crater, Mars, have shown stark color variations from red to purple to gray. These color differences crosscut stratigraphy and are likely due to diagenetic alteration of the sediments after deposition. However, the chemistry and timing of these fluid interactions is unclear. Determining how diagenetic processes may have modified chemical and mineralogical signatures of ancient Martian environments is critical for understanding the past habitability of Mars and achieving the goals of the MSL mission. Here we use visible/near-infrared spectra from Mastcam and ChemCam to determine the mineralogical origins of color variations in the ridge. Color variations are consistent with changes in spectral properties related to the crystallinity, grain size, and texture of hematite. Coarse-grained gray hematite spectrally dominates in the gray patches and is present in the purple areas, while nanophase and fine-grained red crystalline hematite are present and spectrally dominate in the red and purple areas. We hypothesize that these differences were caused by grain-size coarsening of hematite by diagenetic fluids, as observed in terrestrial analogs. In this model, early primary reddening by oxidizing fluids near the surface was followed during or after burial by bleaching to form the gray patches, possibly with limited secondary reddening after exhumation. Diagenetic alteration may have diminished the preservation of biosignatures and changed the composition of the sediments, making it more difficult to interpret how conditions evolved in the paleolake over time.
JGR Planets↗
The application of ensemble wave forcing to quantify uncertainty of shoreline change predictions
Reliable predictions and accompanying uncertainty estimates of coastal evolution on decadal to centennial time scales are increasingly sought. So far, most coastal change projections rely on a single, deterministic realization of the unknown future wave climate, often derived from a global climate model. Yet, deterministic projections do not account for the stochastic nature of future wave conditions across a variety of temporal scales (e.g., daily, weekly, seasonally, and interannually). Here, we present an ensemble Kalman filter shoreline change model to predict coastal erosion and uncertainty due to waves at a variety of time scales. We compare shoreline change projections, simulated with and without ensemble wave forcing conditions by applying ensemble wave time series produced by a computationally efficient statistical downscaling method. We demonstrate a sizable (site-dependent) increase in model uncertainty compared with the unrealistic case of model projections based on a single, deterministic realization (e.g., a single time series) of the wave forcing. We support model-derived uncertainty estimates with a novel mathematical analysis of ensembles of idealized process models. Here, the developed ensemble modeling approach is applied to a well-monitored beach in Tairua, New Zealand. However, the model and uncertainty quantification techniques derived here are generally applicable to a variety of coastal settings around the world.
JGR Earth Surface↗
Recovery of soils from acidic deposition may exacerbate nitrogen export from forested watersheds
Effects of ambient decreases in N deposition on forest N cycling remain unclear as soils recover from acidic deposition. To investigate, repeated soil sampling data were related to deposition, vegetation, and stream data, for 2000–2015 in North and South Buck Creek watersheds, in the Adirondack region of New York, USA. In 63 other Adirondack streams, NO 3 − concentrations were also compared between 2004–2005 and 2014–2015, and a link between soil calcium and stream NO 3 − was investigated using data from 387 Adirondack streams that were sampled in either 2003–2005 or 2010–2011. No trends in N export or NO 3 − concentrations were observed in either Buck watershed despite a 45% decrease in N deposition, although South Buck N export was 2 to 3 times higher than in North Buck, where 48% of deposited N was accounted for by accumulation in the upper soil. In marked contrast, the upper profile in South Buck showed a net loss of N. Increased decomposition appeared likely in South Buck as those soils are adjusted to lower levels of acidifying S deposition, whereas decomposition increases in North Buck were likely suppressed by high levels of natural organic acidity. Stream NO 3 − concentrations in Buck watersheds bracketed regional results and were consistent with the regional streams that showed no overall change in NO 3 − concentrations between 2004 and 2014. A negative correlation observed between NO 3 − concentration and watershed buffering capacity expressed as the ratio of Ca 2+ to SO 4 2− also suggested that stream NO 3 − concentrations were elevated where soil Ca depletion had occurred.
JGR: Biogeosciences↗
Constraints on the geometry of the subducted Gorda Plate with converted phases generated by local earthquakes
The largest slip in great megathrust earthquakes often occurs in the 10–30 km depth range, yet seismic imaging of the material properties in this region has proven difficult. We utilize a dense onshore-offshore passive seismic dataset from the southernmost Cascadia subduction zone where seismicity in the mantle of the subducted Gorda Plate produces S -to- P and P -to- S conversions generated within a few km of the plate interface. These conversions typically occur in the 10–20 km depth range at either the top or bottom of a ∼5 km thick layer with a high Vp/Vs that we infer to be primarily the subducted crust. We use their arrival times and amplitudes to infer the location of the top and bottom of the subducted crust as well as the velocity contrasts across these discontinuities. Comparing with both the Slab1.0 and the updated Slab2 interface models, the Slab2 model is generally consistent with the converted phases, while the Slab1.0 model is 1–2 km deeper in the 2–20 km depth range and ∼6–8 km too deep in the 10–20 km depth range between 40.25°N and 40.4°N. Comparing the amplitudes of the converted phases to synthetics for simplified velocity structures, the amplitude of the converted phases requires models containing a ∼5 km thick zone with at least a ∼10%–20% reduction in S wave velocity. Thus, the plate boundary is likely contained within or at the top of this low velocity zone, which potentially indicates a significant porosity and fluid content within the seismogenic zone.
California↗
Seismic reflection imaging of the low-angle Panamint normal fault system, eastern California
Shallowly dipping (<30°) low‐angle normal faults (LANFs) have been documented globally; however, examples of active LANFs in continental settings are limited. The western margin of the Panamint Range in eastern California is defined by a LANF that dips west beneath Panamint Valley and has evidence of Quaternary motion. In addition, high‐angle dextral‐oblique normal faults displace middle to late Quaternary alluvial fans near the range front. To image shallow (<1 km depth), crosscutting relationships between the low‐ and high‐angle faults along the range front, we acquired two high‐resolution P wave seismic reflection profiles. The northern, 4.6‐km‐long profile crosses the 2‐km‐wide Wildrose graben and the southern, 0.8‐km‐long profile extends onto the Panamint Valley playa, ~7.5 km S of Ballarat, CA. The profile across the Wildrose graben reveals a robust, low‐angle reflector interpreted to represent the LANF separating Plio‐Pleistocene alluvial fanglomerate and Proterozoic metasedimentary deposits. High‐angle faults interpreted in the seismic profile correspond to fault scarps on Quaternary alluvial fan surfaces. Interpretation of the reflection data suggests that the high‐angle faults vertically displace the LANF up to 80 m within the Wildrose graben. Similarly, the profile south of Ballarat reveals a low‐angle reflector, which appears both rotated and displaced up to 260 m by high‐angle faults. These results suggest that near the Panamint range front, the high‐angle faults are the dominant active structures. We conclude that at least at shallow (<1 km) depths, the LANF we imaged is not active today.
California↗
Fluid-earthquake and earthquake-earthquake interactions in southern Kansas, USA
An increase in injection activity associated with energy production in southern Kansas starting in 2013 has been linked to the occurrence of more than 130,000 earthquakes (M −1.5 to 4.9) between 2014 and 2017. Studies suggest that the dramatic increase in seismicity rate is related to wastewater injection into the highly permeable Arbuckle formation. Most of the seismicity is located in the underlying crystalline basement, for which hydrological properties and specific fault geometries are unknown. Additionally, some earthquake clusters occurred relatively far (tens of kilometers) from the main injection wells. Therefore, the effect of pore pressure diffusion may be insufficient to explain the relation between the volume of injected fluids and the spatiotemporal evolution of seismicity. Combining physical models (static stress and poroelasticity) and a statistical cluster analysis applied to a high-resolution relocated catalog, we analyze the evolution of seismicity in southern Kansas. We find that pore pressure changes (Δ p ) and Coulomb stress changes (ΔCFS) due to fluid diffusion smaller than 0.1 MPa are enough to initiate seismic sequences, which then evolve depending on their distance from the major injection wells. However, we find that earthquake sequences have different seismogenic responses to Δ p and ΔCFS in terms of triggering threshold. In regions located close to disposal wells (Harper area) our cluster analysis suggests that both earthquake interactions and fluid diffusion control the evolution of seismicity. On the other hand, at greater distances (Milan area), where clustering behavior suggests greater earthquake interactions, we find that coseismic ΔCFS are larger than Δ p .
Kansas, Oklahoma↗
A geology and geodesy based model of dynamic earthquake rupture on the Rodgers Creek‐Hayward‐Calaveras Fault System, California
The Hayward fault in California's San Francisco Bay area produces large earthquakes, with the last occurring in 1868. We examine how physics‐based dynamic rupture modeling can be used to numerically simulate large earthquakes on not only the Hayward fault, but also its connected companions to the north and south, the Rodgers Creek and Calaveras faults. Equipped with a wealth of images of this fault system, including those of its 3D geology and 3D geometry, in addition to inferences about its interseismic creep‐rate pattern and rock‐friction behavior, we use a finite‐element computer code to perform 3D dynamic earthquake rupture simulations. We find that the rock properties affect the locations and amount of slip produced in our simulated large earthquakes. Crucial factors that control rupture behavior in our modeling are the earthquake nucleation locations, the fault geometry, and the data that reveal where the fault system is creeping or locked. Our findings suggest that large Rodgers Creek‐Hayward‐Calaveras‐Northern Calaveras (RC‐H‐C‐NC) fault‐system earthquakes may result from dynamic rupture that starts in a locked part of the fault system, but is then stopped by the creeping parts, leading to high‐magnitude‐6 earthquakes; or, from dynamic rupture that starts in a locked part of the fault system, then cascades through some of the creeping parts, leading to magnitude‐7 earthquakes.
California↗
Geophysical constraints on the crustal architecture of the transtensional Warm Springs Valley fault zone, northern Walker Lane, western Nevada, USA
The Walker Lane is a zone of distributed transtension where normal faults are overprinted by strike-slip motion. We use two newly-acquired high-resolution seismic reflection profiles and a reprocessed Consortium for Continental Reflection Profiling (COCORP) deep crustal reflection profile to assess the subsurface geometry of the Holocene-active, transtensional Warm Springs Valley fault zone (WSVFZ) near Reno, Nevada, USA. Our multi-scale observations extend to 12 km depth and suggest that the WSVFZ is more complex in the subsurface than implied by late Pleistocene surface fault traces. Two ~4-km-long high-resolution profiles image to a depth of ~2 km and reveal moderately dipping reflections and truncations, some of which project to mapped scarps formed in late Pleistocene surfaces. The shallow lines are co-located with COCORP profile NV 08 along ~40° N latitude. Re-analysis of the COCORP data reveals previously unidentified coherent reflections to a depth of ~12 km and a previously mapped ~30 west-dipping fault at 8-12 km. From these seismic profiles, the WSVFZ is not a simple, sub-vertical fault zone extending through the entire seismogenic crust. Instead, the reflections are consistent with a zone of steep- and moderately-dipping faults that simplify and steepen with depth before intersecting a mid-crustal, low angle (~25-30°) fault. The complex fault geometry of the WSVFZ implies that crustal shear is accommodated by a mix of dipping and subvertical faults in the transtensional northern Walker Lane. If so, transtensional fault zones may present challenges to paleoseismic and geodetic studies and require careful treatment when included in seismic hazard analyses.
Nevada↗
A unified model of crustal stress heterogeneity from borehole breakouts and earthquake focal mechanisms
Observations of crustal stress orientation from the regional inversion of earthquake focal mechanisms often conflict with those from borehole breakouts, possibly indicating local stress heterogeneity, either laterally or with depth. To investigate this heterogeneity, we compiled SHmax estimates from previous studies for 57 near‐vertical boreholes with measured breakout azimuths across the Los Angeles region. We identified subsets of earthquake focal mechanisms from established earthquake catalogs centered around each borehole with various criteria for maximum depth and maximum lateral distance from the borehole. Each subset was independently inverted for 3‐D stress orientation and corresponding SHmax probability distributions, then compared with the corresponding borehole breakout‐derived estimate. We find good agreement when both methods sample the basement stress (breakouts are close to the sediment‐basement interface), or when both methods sample the mid‐basin stress (sufficient earthquakes are present within a sedimentary basin). Along sedimentary basin margins, in contrast, we find acceptable agreement only when focal mechanisms are limited to shallow and close earthquakes, implying short‐length‐scale heterogeneity of <20 km. While the region as a whole shows evidence of both lateral and vertical stress orientation heterogeneity, we find a more homogeneous stress state within basement rock, over length scales of 1–35 km. These results reconcile the apparently conflicting observations of short‐length‐scale heterogeneity observed in boreholes, which sample primarily the basins, with the relative homogeneity of stress inferred from focal mechanisms, which sample primarily the basement, and imply distinct regimes for the appropriate use of each type of stress indicator.
JGR Solid Earth↗