USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Physics of the Earth”

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 55 records · Page 3Linked to original sources

In situ Raman spectroscopic investigation of the structure of subduction-zone fluids

In situ Raman spectra of synthetic subduction-zone fluids (KAlSi3O8-H2O system) were measured to 900?? and 2.3 GPa using a hydrothermal diamond-anvil cell. The structures of aqueous fluid and hydrous melt become closer when conditions approach the second critical endpoint. Almost no three-dimensional network was observed in the supercritical fluid above 2 GPa although a large amount of silicate component is dissolved, suggesting that the physical and chemical properties of these phases change drastically at around the second critical endpoint. Our experimental results indicate that the fluids released from a subducting slab change from aqueous fluid to supercritical fluid with increasing depth under the volcanic arcs. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

The birth of a Hawaiian fissure eruption

Most basaltic explosive eruptions intensify abruptly, allowing little time to document processes at the start of eruption. One opportunity came with the initiation of activity from fissure 8 (F8) during the 2018 eruption on the lower East Rift Zone of Kīlauea, Hawaii. F8 erupted in four episodes. We recorded 28 min of high‐definition video during a 51‐min period, capturing the onset of the second episode on 5 May. From the videos, we were able to analyze the following in‐flight parameters: frequency and duration of explosions; ejecta heights; pyroclast exit velocities; in‐flight total mass and estimated mass eruption rates; and the in‐flight total grain size distributions. The videos record a transition from initial pulsating outgassing, via spaced, but increasingly rapid, discrete explosions, to quasisustained, unsteady fountaining. This transition accompanied waxing intensity (mass flux) of the F8 eruption. We infer that all activity was driven by a combination of the ascent of a coupled mixture of small bubbles and melt, and the buoyant rise of decoupled gas slugs and/or pockets. The balance between these two types of concurrent flow determined the exact form of the eruptive activity at any point in time, and changes to their relative contributions drove the transition we observed at early F8. Qualitative observations of other Hawaiian fountains at Kīlauea suggest that this physical model may apply more generally. This study demonstrates the value of in‐flight parameters derived from high‐resolution videos, which offer a rapid and highly time‐sensitive alternative to measurements based on sampling of deposits posteruption.

Island of Hawai'i↗

Observations and a model of undertow over the inner continental shelf

Onshore volume transport (Stokes drift) due to surface gravity waves propagating toward the beach can result in a compensating Eulerian offshore flow in the surf zone referred to as undertow. Observed offshore flows indicate that wave-driven undertow extends well offshore of the surf zone, over the inner shelves of Martha’s Vineyard, Massachusetts, and North Carolina. Theoretical estimates of the wave-driven offshore transport from linear wave theory and observed wave characteristics account for 50% or more of the observed offshore transport variance in water depths between 5 and 12 m, and reproduce the observed dependence on wave height and water depth. During weak winds, wave-driven cross-shelf velocity profiles over the inner shelf have maximum offshore flow (1–6 cm s −1 ) and vertical shear near the surface and weak flow and shear in the lower half of the water column. The observed offshore flow profiles do not resemble the parabolic profiles with maximum flow at middepth observed within the surf zone. Instead, the vertical structure is similar to the Stokes drift velocity profile but with the opposite direction. This vertical structure is consistent with a dynamical balance between the Coriolis force associated with the offshore flow and an along-shelf “Hasselmann wave stress” due to the influence of the earth’s rotation on surface gravity waves. The close agreement between the observed and modeled profiles provides compelling evidence for the importance of the Hasselmann wave stress in forcing oceanic flows. Summer profiles are more vertically sheared than either winter profiles or model profiles, for reasons that remain unclear.

Massachusetts, North Carolina↗

Biocrusts: The living skin of the Earth

Biological soil crusts (biocrusts) form a “living skin” at the soil surface in many low productivity ecosystems around the world including water- and cold-limited environments, and early successional seres (Belnap et al. 2003). They may be composed of any configuration of soil surface-dwelling cyanobacteria, eukaryotic algae, lichens, mosses or liverworts, and support assemblages of decomposers and a faunal food web (Belnap et al. 2003). These soil surface communities have global relevance, as it has been recently estimated that they cover about 12% of the terrestrial surface currently (Rodriguez-Caballero et al. 2018). Biocrust communities are perhaps an ideal subject for the journal Plant and Soil, because they are simultaneously plant-like, due to their dominance by autotrophs, yet biocrusts are also clearly a physical feature of the soil given that component organisms are enmeshed in, adherent to, or otherwise in direct contact with the soil surface. The activity of the organisms is what engineers the well-aggregated thin layer at the soil surface that we recognize as a biocrust (Belnap et al. 2003). The contributions of biocrusts to ecosystem function has fueled much research interest, initially in the observation of biocrusts’ soil aggregating and erosion-resisting nature, and later as a multifunctional, globally-relevant ecosystem element instrumental in: 1. building or otherwise altering soil nutrient stocks through N-fixation (Elbert et al. 2012), dust trapping (Reynolds et al. 2001) and nutrient cycling (Strauss et al. 2012), 2. influencing hydrological properties of soil such as the water balance (Chamizo et al. 2016), and 3. The thermal energy balance of the ecosystem (Coradeau et al. 2016, Rutherford et al. 2017).

Plant and Soil↗

Morphometric variability within the axial zone of the southern Juan de Fuca Ridge: Interpretation from Sea MARC II, Sea MARC I, and deep-sea photography

The morphometric characteristics of the axial regions of oceanic spreading centers are determined by (1) the type of volcanic flows, (2) the relation between primary volcanic relief (on a scale of a few meters to tens of meters) and degree of sediment cover, and (3) the extent of surficial expression and timing of tectonic disruption of the young oceanic crust. Even within a single, continuous, linear spreading-ridge segment with relatively uniform axial valley dimensions over a distance of 50 or more kilometers, such as along the southern Juan de Fuca Ridge, the changes in morphometric characteristics along axis within the youngest crust indicate distinct variation in tectonic and volcanic activity over short distances within short time periods. An integrated analysis of Sea MARC I, Sea MARC II, and photographic data for the southernmost continuous segment of the Juan de Fuca Ridge shows that generalizations about tectonic and volcanic processes at spreading ridges must consider both the temporal scale of processes as well as the physical scales of observations if predictive models are to be successful. Comparison of the morphometric expression within the major hydrothermal vent area and the rest of the southernmost ridge segment suggests that the mapped distribution of hydrothermal vents may reflect the extent of survey effort rather than uniqueness of geologic setting.

Journal of Geophysical Research B: Solid Earth↗

Physics-based forecasts of eruptive vent locations at calderas

Constraining stresses in the Earth's crust in volcanic regions is critical for understanding many mechanical processes related to eruptive activity. Dike pathways, in particular, are shaped by the orientation of principal stress axes. Therefore, accurate models of dike trajectories and future vent locations rely on accurate estimates of stresses in the subsurface. This work presents a framework for probabilistic constraint of the stress state of calderas by combining three-dimensional physics-based dike pathway models with observed past vent locations using a Monte Carlo approach. The retrieved stress state is then used to produce probability maps of future vent opening across a caldera. We test our stress inversion and vent forecast approach on synthetic scenarios, and find it successful depending on the distribution of the available vents and the complexity of the volcano's structural history. We explore the potential and limitations of the approach, show how its performance is sensitive to the assumptions in the models and available prior information, and discuss how it may be applied to real calderas.

Journal of Geophysical Research: Solid Earth↗

Computation of free oscillations of the earth

Although free oscillations of the Earth may be computed by many different methods, numerous practical considerations have led us to use a Rayleigh-Ritz formulation with piecewise cubic Hermite spline basis functions. By treating the resulting banded matrix equation as a generalized algebraic eigenvalue problem, we are able to achieve great accuracy and generality and a high degree of automation at a reasonable cost.

Journal of Computational Physics↗

Introduction to special issue on gas hydrate in porous media: Linking laboratory and field‐scale phenomena

The proliferation of drilling expeditions focused on characterizing natural gas hydrate as a potential energy resource has spawned widespread interest in gas hydrate reservoir properties and associated porous media phenomena. Between 2017 and 2019, a Special Section of this journal compiled contributed papers elucidating interactions between gas hydrate and sediment based on laboratory, numerical modeling, and field studies. Motivated mostly by field observations in the northern Gulf of Mexico and offshore Japan, several papers focus on the mechanisms for gas hydrate formation and accumulation, particularly with vapor phase gas, not dissolved gas, as the precursor to hydrate. These studies rely on numerical modeling or laboratory experiments using sediment packs or benchtop micromodels. A second focus of the Special Section is the role of fines in inhibiting production of gas from methane hydrate, controlling the distribution of hydrate at a pore scale, and influencing the bulk behavior of seafloor sediments. Other papers fill knowledge gaps related to the physical properties of hydrate-bearing sediments and advance new approaches in coupled thermal-mechanical modeling of these sediments during hydrate dissociation. Finally, one study addresses the long-standing question about the fate of methane hydrate at the molecular level when CO2 is injected into natural reservoirs under hydrate-forming conditions.

Journal of Geophysical Research B: Solid Earth↗

Venus cartography

The entire surface of the planet Venus is being mapped at global and regional scales (1:50,000,000 through 1:1,500,000) with synthetic aperture radar (SAR), radar altimeter, and radiometer measurements of physical properties from the Magellan spacecraft. The mapping includes SAR image mosaics, shaded relief maps, and topographic contour overlays made from altimetry data and by radargrammetric methods. Methods used include new techniques of radar image processing that became operational as a result of the Magellan mission. The area being mapped (at a resolution of 75 m/pixel) is roughly equivalent to that of Earth, including seafloors. The mapping is designed to support geologic and geophysical investigations.

Journal of Geophysical Research E: Planets↗

Putting it all together: Exhumation histories from a formal combination of heat flow and a suite of thermochronometers

A suite of new techniques in thermochronometry allow analysis of the thermal history of a sample over a broad range of temperature sensitivities. New analysis tools must be developed that fully and formally integrate these techniques, allowing a single geologic interpretation of the rate and timing of exhumation and burial events consistent with all data. We integrate a thermal model of burial and exhumation, (U-Th)/He age modeling, and fission track age and length modeling. We then use a genetic algorithm to efficiently explore possible time-exhumation histories of a vertical sample profile (such as a borehole), simultaneously solving for exhumation and burial rates as well as changes in background heat flow. We formally combine all data in a rigorous statistical fashion. By parameterizing the model in terms of exhumation rather than time-temperature paths (as traditionally done in fission track modeling), we can ensure that exhumation histories result in a sedimentary basin whose thickness is consistent with the observed basin, a physically based constraint that eliminates otherwise acceptable thermal histories. We apply the technique to heat flow and thermochronometry data from the 2.1 -km-deep San Andreas Fault Observatory at Depth pilot hole near the San Andreas fault, California. We find that the site experienced <1 km of exhumation or burial since the onset of San Andreas fault activity

Journal of Geophysical Research B: Solid Earth↗

Time-dependent earthquake probabilities

We have attempted to provide a careful examination of a class of approaches for estimating the conditional probability of failure of a single large earthquake, particularly approaches that account for static stress perturbations to tectonic loading as in the approaches of Stein et al. (1997) and Hardebeck (2004). We have loading as in the framework based on a simple, generalized rate change formulation and applied it to these two approaches to show how they relate to one another. We also have attempted to show the connection between models of seismicity rate changes applied to (1) populations of independent faults as in background and aftershock seismicity and (2) changes in estimates of the conditional probability of failures of different members of a the notion of failure rate corresponds to successive failures of different members of a population of faults. The latter application requires specification of some probability distribution (density function of PDF) that describes some population of potential recurrence times. This PDF may reflect our imperfect knowledge of when past earthquakes have occurred on a fault (epistemic uncertainty), the true natural variability in failure times, or some combination of both. We suggest two end-member conceptual single-fault models that may explain natural variability in recurrence times and suggest how they might be distinguished observationally. When viewed deterministically, these single-fault patch models differ significantly in their physical attributes, and when faults are immature, they differ in their responses to stress perturbations. Estimates of conditional failure probabilities effectively integrate over a range of possible deterministic fault models, usually with ranges that correspond to mature faults. Thus conditional failure probability estimates usually should not differ significantly for these models. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Ionospheric current source modeling and global geomagnetic induction using ground geomagnetic observatory data

Long-period global-scale electromagnetic induction studies of deep Earth conductivity are based almost exclusively on magnetovariational methods and require accurate models of external source spatial structure. We describe approaches to inverting for both the external sources and three-dimensional (3-D) conductivity variations and apply these methods to long-period ( T ≥1.2 days) geomagnetic observatory data. Our scheme involves three steps: (1) Observatory data from 60 years (only partly overlapping and with many large gaps) are reduced and merged into dominant spatial modes using a scheme based on frequency domain principal components. (2) Resulting modes are inverted for corresponding external source spatial structure, using a simplified conductivity model with radial variations overlain by a two-dimensional thin sheet. The source inversion is regularized using a physically based source covariance, generated through superposition of correlated tilted zonal (quasi-dipole) current loops, representing ionospheric source complexity smoothed by Earth rotation. Free parameters in the source covariance model are tuned by a leave-one-out cross-validation scheme. (3) The estimated data modes are inverted for 3-D Earth conductivity, assuming the source excitation estimated in step 2. Together, these developments constitute key components in a practical scheme for simultaneous inversion of the catalogue of historical and modern observatory data for external source spatial structure and 3-D Earth conductivity.

Journal of Geophysical Research↗

Rheology of the lithosphere inferred from postseismic uplift following the 1959 Hebgen Lake earthquake

We have modeled the broad postseismic uplift measured by geodetic leveling in the epicentral area of the 1959 Mw = 7.3 Hebgen Lake, Montana earthquake, a normal faulting event in the northern Basin and Range province. To fit the observed uplift we calculate synthetic postseismic deformation using the relaxation response of a gravitational viscoelastic Earth to the earthquake. For a model with an elastic plate overlying a viscoelastic half-space, we find that the elastic thickness is 38 ?? 8 km, which isclose to the local crustal thickness. The half-space viscosity is estimated at 4 ?? 1018??0.5 Pa s. The leveling data do not require a viscous lower crust but permit a lower bound viscosity of 1020 Pa s. The observed broad uplift cannot be explained by physically plausible afterslip on and below the coseismic fault. However, local deformation across the coseismic surface rupture requires shallow afterslip reaching the surface. The postseismic deformation induced by the estimated viscoelastic structure decays exponentially with a time constant of ???15 years. Because of coupling between the elastic layer and the viscoelastic substrate, this relaxation time is significantly longer than the 2 year Maxwell relaxation time of the viscous half-space itself. Our result suggests the importance of postseismic relaxation in interpreting high-precision global positioning system velocities. For example, our model results suggest that postseismic transient velocities from both the 1959 Hebgen Lake and the 1983 Mw = 6.9 Borah Peak earthquakes are currently as large as 1-2 mm/yr.

Journal of Geophysical Research B: Solid Earth↗

Olivine friction at the base of oceanic seismogenic zones

We investigate the strength and frictional behavior of olivine aggregates at temperatures and effective confining pressures similar to those at the base of the seismogenic zone on a typical ridge transform fault. Triaxial compression tests were conducted on dry olivine powder (grain size ???60 ??m) at effective confining pressures between 50 and 300 MPa (using Argon as a pore fluid), temperatures between 600??C and 1000??C, and axial displacement rates from 0.06 to 60 ??m/s (axial strain rates from 3 ?? 10-6 to 3 ?? 10-3 s-1). Yielding shows a negative pressure dependence, consistent with predictions for shear enhanced compaction and with the observation that samples exhibit compaction during the initial stages of the experiments. A combination of mechanical data and microstructural observations demonstrate that deformation was accommodated by frictional processes. Sample strengths were pressure-dependent and nearly independent of temperature. Localized shear zones formed in initially homogeneous aggregates early in the experiments. The frictional response to changes in loading rate is well described by rate and state constitutive laws, with a transition from velocity-weakening to velocity-strengthening at 1000??C. Microstructural observations and physical models indicate that plastic yielding of asperities at high temperatures and low axial strain rates stabilizes frictional sliding. Extrapolation of our experimental data to geologic strain rates indicates that a transition from velocity weakening to velocity strengthening occurs at approximately 600??C, consistent with the focal depths of earthquakes in the oceanic lithosphere. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Close temporal correspondence between geomagnetic anomalies and earthquakes during the 2002-2003 eruption of Etna volcano

The early stages of the 2002-2003 lateral eruption at Mount Etna were accompanied by slow changes (over some hours) and some rapid step offsets in the local magnetic field. At five monitoring locations, the total magnetic field intensity has been measured using continuously operating Overhauser magnetometers at a sampling rate of 10 s. The very unique aspect of these observations is the close temporal correspondence between magnetic field offsets and earthquakes that occurred in the upper northern flank of the volcano on 27 October 2002 prior to a primary eruption. Rapid coseismic changes of the magnetic field were clearly identified for three of the most energetic earthquakes, which were concentrated along the Northeast Rift at a depth of about 1 km below sea level. Coseismic magnetic signals, with amplitudes from 0.5 to 2.5 nT, have been detected for three of the largest seismic events located roughly midway between the magnetic stations. We quantitatively examine possible geophysical mechanisms, which could cause the magnetic anomalies. The comparison between magnetic data, seismicity and surface phenomena implies that piezomagnetic effects are the primary physical mechanism responsible for the observed magnetic anomalies although the detailed cause of the rapid high stress change required is not clear. The modeling of the observed coseismic magnetic changes in terms of piezomagnetic mechanism provides further evidence of the complex interaction between volcanic and tectonic processes during dike propagation along the Northeast Rift. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

When and where the aftershock activity was depressed: Contrasting decay patterns of the proximate large earthquakes in southern California

Seismic quiescence has attracted attention as a possible precursor to a large earthquake. However, sensitive detection of quiescence requires accurate modeling of normal aftershock activity. We apply the epidemic-type aftershock sequence (ETAS) model that is a natural extension of the modified Omori formula for aftershock decay, allowing further clusters (secondary aftershocks) within an aftershock sequence. The Hector Mine aftershock activity has been normal, relative to the decay predicted by the ETAS model during the 14 months of available data. In contrast, although the aftershock sequence of the 1992 Landers earthquake (M = 7.3), including the 1992 Big Bear earthquake (M = 6.4) and its aftershocks, fits very well to the ETAS up until about 6 months after the main shock, the activity showed clear lowering relative to the modeled rate (relative quiescence) and lasted nearly 7 years, leading up to the Hector Mine earthquake (M = 7.1) in 1999. Specifically, the relative quiescence occurred only in the shallow aftershock activity, down to depths of 5-6 km. The sequence of deeper events showed clear, normal aftershock activity well fitted to the ETAS throughout the whole period. We argue several physical explanations for these results. Among them, we strongly suspect aseismic slips within the Hector Mine rupture source that could inhibit the crustal relaxation process within "shadow zones" of the Coulomb's failure stress change. Furthermore, the aftershock activity of the 1992 Joshua Tree earthquake (M = 6.1) sharply lowered in the same day of the main shock, which can be explained by a similar scenario.

Journal of Geophysical Research B: Solid Earth↗

Radiometric cross-calibration of EO-1 ALI with L7 ETM+ and Terra MODIS sensors using near-simultaneous desert observations

The Earth Observing-1 (EO-1) satellite was launched on November 21, 2000, as part of a one-year technology demonstration mission. The mission was extended because of the value it continued to add to the scientific community. EO-1 has now been operational for more than a decade, providing both multispectral and hyperspectral measurements. As part of the EO-1 mission, the Advanced Land Imager (ALI) sensor demonstrates a potential technological direction for the next generation of Landsat sensors. To evaluate the ALI sensor capabilities as a precursor to the Operational Land Imager (OLI) onboard the Landsat Data Continuity Mission (LDCM, or Landsat 8 after launch), its measured top-of-atmosphere (TOA) reflectances were compared to the well-calibrated Landsat 7 (L7) Enhanced Thematic Mapper Plus (ETM+) and the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) sensors in the reflective solar bands (RSB). These three satellites operate in a near-polar, sun-synchronous orbit 705 km above the Earth's surface. EO-1 was designed to fly one minute behind L7 and approximately 30 minutes in front of Terra. In this configuration, all the three sensors can view near-identical ground targets with similar atmospheric, solar, and viewing conditions. However, because of the differences in the relative spectral response (RSR), the measured physical quantities can be significantly different while observing the same target. The cross-calibration of ALI with ETM+ and MODIS was performed using near-simultaneous surface observations based on image statistics from areas observed by these sensors over four desert sites (Libya 4, Mauritania 2, Arabia 1, and Sudan 1). The differences in the measured TOA reflectances due to RSR mismatches were compensated by using a spectral band adjustment factor (SBAF), which takes into account the spectral profile of the target and the RSR of each sensor. For this study, the spectral profile of the target comes from the near-simultaneous EO-1 Hyperion data over these sites. The results indicate that the TOA reflectance measurements for ALI agree with those of ETM+ and MODIS to within 5% after the application of SBAF.

Arabia↗

Ultraviolet and visible remote sensing of volcanic gases

As magma rises in volcanic systems, volatile species exsolve from the silicate melt and are emitted as gases into the atmosphere. Measuring the magnitude and composition of gas emissions from volcanoes provides insights into processes occurring deep within the Earth and helps constrain the impact of volcanic degassing on atmospheric chemistry. Optical remote sensing techniques allow volcanic gas emissions to be characterized without the need to access hazardous areas near active volcanic vents. This paper reviews the state of the art in ultraviolet and visible volcanic gas remote sensing from the ground, air, and space. Special attention is given to discussing the physics of atmospheric radiative transfer on which these techniques are based. The functionality and limitations of different remote sensing instruments are examined, making clear that the ideal choice of instrumentation will depend on the volcanic system to which it is applied and the sought measurement parameters. Common algorithms for determining trace gas column densities, gas burdens, and volcanic emission rates from measurements of spectral radiance are outlined and compared, showing how some algorithms attempt to model the physics of the measurement while others maximize sensitivity. Several examples demonstrate how remote sensing measurements continue to advance our understanding of volcanic systems and their impact on the atmosphere. Finally, a few promising directions of inquiry are suggested that could lead to improvements in remote sensing instrumentation and analysis techniques. By combining spectroscopic and imaging techniques, improving our understanding of atmospheric radiative transfer, expanding the suite of target gases, and increasing the coverage and frequency of observations, we stand to significantly improve our ability to detect and quantify volcanic gas emissions and gain new insights into important Earth-system processes.

Journal of Volcanology and Geothermal Research↗