USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research A: Space Physics”

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.

15 recordsLinked to original sources
USGS2024

Challenging ring-current models of the Carrington storm

A detailed analysis is made of horizontal-component geomagnetic-disturbance data acquired at the Colaba observatory in India recording the Carrington magnetic storm of September 1859. Prior to attaining its maximum absolute value, disturbance at Colaba increased with an e -folding timescale of 0.46 hr (28 min). Following its maximum, absolute disturbance at Colaba decreased as a trend having an e -folding timescale of 0.31 hr (19 min). Both of these timescales are much shorter than those characterizing the drift period of ring-current ions. Furthermore, over one 28-min interval when absolute disturbance was increasing, the data indicate an absolute rate of change of ≥2,436 nT/hr. If this is representative of disturbance generated by a symmetric magnetospheric ring current, then, assuming a standard and widely used parameterization, an interplanetary electric field of ≥451 mV/m is indicated. An idealized and extreme solar-wind dynamic pressure could, conceivably, reduce this bound on the interplanetary electric field to ≥202 mV/m. If the parameterization for electric-field extrapolation is accurate, but the field strengths obtained are deemed implausible, then it can be concluded that the Colaba disturbance data were significantly affected by partial-ring, field-aligned, or ionospheric currents. The same conclusion is supported by the shortness of the e -folding timescales characterizing the Colaba data. Several prominent studies of the Carrington event need to be reconsidered.

Journal of Geophysical Research Space Physics↗
USGS2008

Lightning‐driven electric fields measured in the lower ionosphere: Implications for transient luminous events

Transient luminous events above thunderstorms such as sprites, halos, and elves require large electric fields in the lower ionosphere. Yet very few in situ measurements in this region have been successfully accomplished, since it is typically too low in altitude for rockets and satellites and too high for balloons. In this article, we present some rare examples of lightning‐driven electric field changes obtained at 75–130 km altitude during a sounding rocket flight from Wallops Island, Virginia, in 1995. We summarize these electric field changes and present a few detailed case studies. Our measurements are compared directly to a 2D numerical model of lightning‐driven electromagnetic fields in the middle and upper atmosphere. We find that the in situ electric field changes are smaller than predicted by the model, and the amplitudes of these fields are insufficient for elve production when extrapolated to a 100 kA peak current stroke. This disagreement could be due to lightning‐induced ionospheric conductivity enhancement, or it might be evidence of flaws in the electromagnetic pulse mechanism for elves.

Journal of Geophysical Research A: Space Physics↗
USGS2011

An empirical model of the quiet daily geomagnetic field variation

An empirical model of the quiet daily geomagnetic field variation has been constructed based on geomagnetic data obtained from 21 stations along the 210 Magnetic Meridian of the Circum‐pan Pacific Magnetometer Network (CPMN) from 1996 to 2007. Using the least squares fitting method for geomagnetically quiet days ( Kp ≤ 2+), the quiet daily geomagnetic field variation at each station was described as a function of solar activity SA , day of year DOY , lunar age LA , and local time LT . After interpolation in latitude, the model can describe solar‐activity dependence and seasonal dependence of solar quiet daily variations (S) and lunar quiet daily variations (L). We performed a spherical harmonic analysis (SHA) on these S and L variations to examine average characteristics of the equivalent external current systems. We found three particularly noteworthy results. First, the total current intensity of the S current system is largely controlled by solar activity while its focus position is not significantly affected by solar activity. Second, we found that seasonal variations of the S current intensity exhibit north‐south asymmetry; the current intensity of the northern vortex shows a prominent annual variation while the southern vortex shows a clear semi‐annual variation as well as annual variation. Thirdly, we found that the total intensity of the L current system changes depending on solar activity and season; seasonal variations of the L current intensity show an enhancement during the December solstice, independent of the level of solar activity.

Journal of Geophysical Research A: Space Physics↗
USGS2012

Problem of the Love‐Gannon relation between the asymmetric disturbance field and Dst

Love and Gannon (2009) discovered that statistically, over a fifty year period the difference in the dawn and dusk disturbance‐field H component at low latitudes (hourly averaged) is linearly proportional to Dst. If the difference is designated by δ DD in units of nT/R E , then the Love‐Gannon (L‐G) relation is δ DD = − 0.2 Dst. At any time departures from the relation can be large. Nonetheless, the relation is evident for all values of Dst and persists throughout magnetic storms, both the main phase and the recovery phase. The Love‐Gannon discovery presents a problem to current understanding of the relation between the causes of δ DD and Dst because the dawn dusk asymmetry in the disturbance field is presumably governed by a long‐established magnetosphere‐ionosphere coupling theory which predicts a characteristic time scale (the shielding time) of less than an hour whereas the characteristic time scale for Dst (the ring current decay time) is more like ten hours. Thus, without forcing both time scales toward each other to the limits of their ranges, a linear proportionality between δ DD and Dst cannot be derived from the current understanding of the causes of the asymmetry and the ring current. This conclusion is the paper's main contribution. In addition, we attempt to get around the conflict of time scales by looking at other possibilities for generating δ DD that depend directly on the ring current. The most promising of these is the possibility that the ring current decay mechanism creates a quasi‐permanent, local‐time modification of the ring current compared to what it would be in the absence of the decay mechanism and that this modification causes a field‐aligned current that closes through the ionosphere and generates the asymmetry δ DD . This idea has the virtue of coupling the asymmetry directly to the ring current and of accounting for the persistence of the L‐G proportionality through the recovery phase of magnetic storms.

Journal of Geophysical Research A: Space Physics↗
USGS2012

Interseismic uplift at the Nankai subduction zone, southwest Japan, 1951–1990

Uplift as a function of time from 1951 through 1990 has been deduced from annual mean sea level measured at 15 tide gages along the Nankai subduction zone, southwest Japan. The recurrence interval for rupture of the Nankai subduction zone is about 100 years, and the most recent rupture was in late 1946. Thus the 1951–1990 uplift record covers most of the first half of the earthquake cycle. The precision of the uplift record is better than could be obtained currently by annual Global Positioning System (GPS) measurements. The pre‐1960 uplift record shows rapid deformation that appears to terminate in 1959. After 1959 the uplift record is remarkably linear in time: Significant curvature in the uplift‐versus‐time plots is detected at only three of the 15 tide gages. The inferred uplift rates are not quantitatively consistent with the predictions of either the viscoelastic coupling or elastic half‐space models of subduction, but the agreement is sufficient to suggest that the causative physical mechanisms have been identified. The immediate postseismic response is controlled by the propagation of slip downward along the downdip extension of the coseismic rupture, and the interseismic response is caused by the accumulation of a slip deficit on the main thrust zone. Asthenosphere relaxation is not required to explain the observations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

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↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 2. Experimental tests

Scaling considerations indicate that miniature experiments can be used to test models of granular avalanches in which the effects of intergranular fluid and cohesion are negligible. To test predictions of a granular avalanche model described in a companion paper, we performed bench top experiments involving avalanches of dry sand across irregular basal topography that mimicked the complexity of natural terrain. The experiments employed a novel method of laser-assisted cartography to map the three-dimensional morphology of rapidly moving avalanches, thereby providing high-resolution data for comparison with model output. Model input consisted of two material properties (angles of internal and basal Coulomb friction of the sand), which were measured in independent tests, and of initial and boundary conditions that characterized the geometry of the experimental apparatus. Experimental results demonstrate that the model accurately predicts not only the gross behavior but also many details of avalanche motion, from initiation to deposition. We attribute this accuracy to a mathematical and computational formulation that conserves mass and momentum in three-dimensional physical space and satisfies the Coulomb equation in three-dimensional stress space. Our results support the hypothesis that a Coulomb proportionality between shear and normal stresses applies in moderately rapid granular flows and that complicated constitutive postulates are unnecessary if momentum conservation is strictly enforced in continuum avalanche models. Furthermore, predictions of our Coulomb continuum model contrast with those of a Coulomb point mass model, illustrating the importance of multidimensional modeling and model testing.

Journal of Geophysical Research Earth Surface↗
USGS2007

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.

Journal of Geophysical Research F: Earth Surface↗
USGS2013

Evaluation of a new model of aeolian transport in the presence of vegetation

Aeolian transport is an important characteristic of many arid and semiarid regions worldwide that affects dust emission and ecosystem processes. The purpose of this paper is to evaluate a recent model of aeolian transport in the presence of vegetation. This approach differs from previous models by accounting for how vegetation affects the distribution of shear velocity on the surface rather than merely calculating the average effect of vegetation on surface shear velocity or simply using empirical relationships. Vegetation, soil, and meteorological data at 65 field sites with measurements of horizontal aeolian flux were collected from the Western United States. Measured fluxes were tested against modeled values to evaluate model performance, to obtain a set of optimum model parameters, and to estimate the uncertainty in these parameters. The same field data were used to model horizontal aeolian flux using three other schemes. Our results show that the model can predict horizontal aeolian flux with an approximate relative error of 2.1 and that further empirical corrections can reduce the approximate relative error to 1.0. The level of error is within what would be expected given uncertainties in threshold shear velocity and wind speed at our sites. The model outperforms the alternative schemes both in terms of approximate relative error and the number of sites at which threshold shear velocity was exceeded. These results lend support to an understanding of the physics of aeolian transport in which (1) vegetation's impact on transport is dependent upon the distribution of vegetation rather than merely its average lateral cover and (2) vegetation impacts surface shear stress locally by depressing it in the immediate lee of plants rather than by changing the bulk surface's threshold shear velocity. Our results also suggest that threshold shear velocity is exceeded more than might be estimated by single measurements of threshold shear stress and roughness length commonly associated with vegetated surfaces, highlighting the variation of threshold shear velocity with space and time in real landscapes.

Journal of Geophysical Research F: Earth Surface↗
USGS2011

Detecting aseismic strain transients from seismicity data

Aseismic deformation transients such as fluid flow, magma migration, and slow slip can trigger changes in seismicity rate. We present a method that can detect these seismicity rate variations and utilize these anomalies to constrain the underlying variations in stressing rate. Because ordinary aftershock sequences often obscure changes in the background seismicity caused by aseismic processes, we combine the stochastic Epidemic Type Aftershock Sequence model that describes aftershock sequences well and the physically based rate‐ and state‐dependent friction seismicity model into a single seismicity rate model that models both aftershock activity and changes in background seismicity rate. We implement this model into a data assimilation algorithm that inverts seismicity catalogs to estimate space‐time variations in stressing rate. We evaluate the method using a synthetic catalog, and then apply it to a catalog of M ≥ 1.5 events that occurred in the Salton Trough from 1990 to 2009. We validate our stressing rate estimates by comparing them to estimates from a geodetically derived slip model for a large creep event on the Obsidian Buttes fault. The results demonstrate that our approach can identify large aseismic deformation transients in a multidecade long earthquake catalog and roughly constrain the absolute magnitude of the stressing rate transients. Our method can therefore provide a way to detect aseismic transients in regions where geodetic resolution in space or time is poor.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

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↗
USGS2022

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗
USGS2017

Decaying lava extrusion rate at El Reventador Volcano, Ecuador measured using high-resolution satellite radar

Lava extrusion at erupting volcanoes causes rapid changes in topography and morphology on the order of tens or even hundreds of meters. Satellite radar provides a method for measuring changes in topographic height over a given time period to an accuracy of meters, either by measuring the width of radar shadow cast by steep sided features, or by measuring the difference in radar phase between two sensors separated in space. We measure height changes, and hence estimate extruded lava volume flux, at El Reventador, Ecuador, between 2011 and 2016, using data from the RADARSAT-2 and TanDEM-X satellite missions. We find that 39 new lava flows were extruded between 9 February 2012 and 24 August 2016, with a cumulative volume of 44.8M m 3 dense rock equivalent and a gradually decreasing eruption rate. The average dense rock rate of lava extrusion during this time is 0.31 ± 0.02 m 3 s −1 , which is similar to the long-term average from 1972 to 2016. Apart from a volumetrically small dyke opening event between 9 March and 10 June 2012, lava extrusion at El Reventador is not accompanied by any significant magmatic ground deformation. We use a simple physics-based model to estimate that the volume of the magma reservoir under El Reventador is greater than 3 km 3 . Our lava extrusion data can be equally well fit by models representing a closed reservoir depressurising during the eruption with no magma recharge, or an open reservoir with a time-constant magma recharge rate of up to 0.35 ± 0.01 m 3 s −1 .

El Reventador Volcano↗
USGS2001

Mars Color Imager (MARCI) on the Mars Climate Orbiter

The Mars Color Imager, or MARCI, experiment on the Mars Climate Orbiter (MCO) consists of two cameras with unique optics and identical focal plane assemblies (FPAs), Data Acquisition System (DAS) electronics, and power supplies. Each camera is characterized by small physical size and mass (∼6×6×12 cm, including baffle; <500 g), low power requirements (<2.5 W, including power supply losses), and high science performance (1000×1000 pixel, low noise). The Wide Angle (WA) camera will have the capability to map Mars in five visible and two ultraviolet spectral bands at a resolution of better than 8 km/pixel under the worst case downlink data rate. Under better downlink conditions the WA will provide kilometer-scale global maps of atmospheric phenomena such as clouds, hazes, dust storms, and the polar hood. Limb observations will provide additional detail on atmospheric structure at 1/3 scale-height resolution. The Medium Angle (MA) camera is designed to study selected areas of Mars at regional scale. From 400 km altitude its 6° FOV, which covers ∼40 km at 40 m/pixel, will permit all locations on the planet except the poles to be accessible for image acquisitions every two mapping cycles (roughly 52 sols). Eight spectral channels between 425 and 1000 nm provide the ability to discriminate both atmospheric and surface features on the basis of composition. The primary science objectives of MARCI are to (1) observe Martian atmospheric processes at synoptic scales and mesoscales, (2) study details of the interaction of the atmosphere with the surface at a variety of scales in both space and time, and (3) examine surface features characteristic of the evolution of the Martian climate over time. MARCI will directly address two of the three high-level goals of the Mars Surveyor Program: Climate and Resources. Life, the third goal, will be addressed indirectly through the environmental factors associated with the other two goals.

Journal of Geophysical Research E: Planets↗