Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Applied Geophysics”

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 181 records · Page 10Linked to original sources

Reoccupation of floodplains by rivers and its relation to the age structure of floodplain vegetation

River channel dynamics over many decades provide a physical control on the age structure of floodplain vegetation as a river occupies and abandons locations. Floodplain reoccupation by a river, in particular, determines the interval of time during which vegetation can establish and mature. A general framework for analyzing floodplain reoccupation and a time series model are developed and applied to five alluvial rivers in the United States. Channel dynamics in these rivers demonstrate time-scale dependence with short-term oscillation in active channel area in response to floods and subsequent vegetation growth and progressive lateral movement that accounts for much of the cumulative area occupied by the rivers over decades. Rivers preferentially reoccupy locations recently abandoned causing a decreasing probability of reoccupation with time since abandonment. For a typical case, a river is 10 times more likely to reoccupy an area it abandoned in the past decade than it is to reoccupy an area it abandoned 30 yrs ago. The decreasing probability of reoccupation over time is consistent with observations of persistent stands of late seral stage floodplain forest. A power function provides a robust approach for estimating the cumulative area occupied by a river and the age structure of riparian forests resulting from a specific historical sequence of streamflow in comparison to either linear or exponential alternatives.

Journal of Geophysical Research

Reconstructing the velocity and deformation of a rapid landslide using multiview video

Noncontact measurements of spatially varied ground surface deformation during landslide motion can provide important constraints on landslide mechanics. Here, we present and test a new method for extracting measurements of rapid landslide surface displacement and velocity (accelerations of approximately 1 m/s 2 ) using sequences of stereo images obtained from a pair of inexpensive, stationary 4K video cameras with nominal frame rates of 29.97 Hz. The method combines elements of Structure from Motion with those of optical flow to extract data on 3‐D evolution of the ground surface during slope failure. We apply the method to an experiment at the U.S. Geological Survey debris‐flow flume in which a high‐speed, liquefying landslide was triggered by gradually adding water to a 6‐m 3 prism of loosely packed sediment on a 31° slope. Strip‐scanning lidar measurements made during the experiment corroborate our video‐based measurements, but the latter encompassed the entire landslide surface and were much lower in cost. Our video‐based measurements enabled computation of depth‐integrated landslide dilation/contraction rates. The range of computed rates was within the ranges inferred from independent measurements of evolving pore water pressures and reasonable estimates of the hydraulic permeability of the sediment. Dilation and contraction rates play a crucial role in landslide mechanics. The dilation and contraction we observe contradict the incompressible flow assumption used in many studies that have employed noncontact methods to infer landslide properties.

Oregon

Distributed energy balance modeling of South Cascade Glacier, Washington and assessment of model uncertainty

We have developed a physically based, distributed surface energy balance model to simulate glacier mass balance under meteorological and climatological forcing. Here we apply the model to estimate summer ablation on South Cascade Glacier, Washington, for the 2004 and 2005 mass balance seasons. To arrive at optimal mass balance simulations, we investigate and quantify model uncertainty associated with selecting from a range of physical parameter values that are not commonly measured in glaciological mass balance field studies. We optimize the performance of the model by varying values for atmospheric transmissivity, the albedo of surrounding topography, precipitation-elevation lapse rate, surface roughness for turbulent exchange of momentum, and snow albedo aging coefficient. Of these the snow aging parameter and precipitation lapse rates have the greatest influence on the modeled ablation. We examined model sensitivity to varying parameters by performing an additional 103 realizations with parameters randomly chosen over a ??5% range centered about the optimum values. The best fit suite of model parameters yielded a net balance of -1.69??0.38 m water equivalent (WE) for the 2004 water year and -2.10??0.30 m WE up to 11 September 2005. The 2004 result is within 3% of the measured value. These simulations account for 91% and 93% of the variance in measured ablation for the respective years. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface

Spatial-temporal variation of low-frequency earthquake bursts near Parkfield, California

Tectonic tremor (TT) and low-frequency earthquakes (LFEs) have been found in the deeper crust of various tectonic environments globally in the last decade. The spatial-temporal behaviour of LFEs provides insight into deep fault zone processes. In this study, we examine recurrence times from a 12-yr catalogue of 88 LFE families with ∼730 000 LFEs in the vicinity of the Parkfield section of the San Andreas Fault (SAF) in central California. We apply an automatic burst detection algorithm to the LFE recurrence times to identify the clustering behaviour of LFEs (LFE bursts) in each family. We find that the burst behaviours in the northern and southern LFE groups differ. Generally, the northern group has longer burst duration but fewer LFEs per burst, while the southern group has shorter burst duration but more LFEs per burst. The southern group LFE bursts are generally more correlated than the northern group, suggesting more coherent deep fault slip and relatively simpler deep fault structure beneath the locked section of SAF. We also found that the 2004 Parkfield earthquake clearly increased the number of LFEs per burst and average burst duration for both the northern and the southern groups, with a relatively larger effect on the northern group. This could be due to the weakness of northern part of the fault, or the northwesterly rupture direction of the Parkfield earthquake.

California

Uncovering the nonadiabatic response of geosynchronous electrons to geomagnetic disturbance

We describe an energy spectrum method for scaling electron integral flux, which is measured at a constant energy, to phase space density at a constant value of the first adiabatic invariant which removes much of the variation due to reversible adiabatic effects. Applying this method to nearly a solar cycle (1995 - 2006) of geosynchronous electron integral flux (E>2.0MeV) from the GOES satellites, we see that much of the diurnal variation in electron phase space density at constant energy can be removed by the transformation to phase space density at constant μ (4000 MeV/G). This allows us a clearer picture of underlying non-adiabatic electron population changes due to geomagnetic activity. Using scaled phase space density, we calculate the percentage of geomagnetic storms resulting in an increase, decrease or no change in geosynchronous electrons as 38%, 7%, and 55%, respectively. We also show examples of changes in the electron population that may be different than the unscaled fluxes alone suggest. These examples include sudden electron enhancements during storms which appear during the peak of negative Dst for μ-scaled phase space density, contrary to the slow increase seen during the recovery phase for unscaled phase space density for the same event.

Journal of Geophysical Research

Coupled semivariogram uncertainty of hydrogeological and geophysical data on capture zone uncertainty analysis

This study investigates capture zone uncertainty that relates to the coupled semivariogram uncertainty of hydrogeological and geophysical data. Semivariogram uncertainty is represented by the uncertainty in structural parameters (range, sill, and nugget). We used the beta distribution function to derive the prior distributions of structural parameters. The probability distributions of structural parameters were further updated through the Bayesian approach with the Gaussian likelihood functions. Cokriging of noncollocated pumping test data and electrical resistivity data was conducted to better estimate hydraulic conductivity through autosemivariograms and pseudo-cross-semivariogram. Sensitivities of capture zone variability with respect to the spatial variability of hydraulic conductivity, porosity and aquifer thickness were analyzed using ANOVA. The proposed methodology was applied to the analysis of capture zone uncertainty at the Chicot aquifer in Southwestern Louisiana, where a regional groundwater flow model was developed. MODFLOW-MODPATH was adopted to delineate the capture zone. The ANOVA results showed that both capture zone area and compactness were sensitive to hydraulic conductivity variation. We concluded that the capture zone uncertainty due to the semivariogram uncertainty is much higher than that due to the kriging uncertainty for given semivariograms. In other words, the sole use of conditional variances of kriging may greatly underestimate the flow response uncertainty. Semivariogram uncertainty should also be taken into account in the uncertainty analysis. ?? 2008 ASCE.

Journal of Hydrologic Engineering

Using corrected and imputed polarity measurements to improve focal mechanisms in a regional earthquake catalog near the Mt. Lewis Fault Zone, California

We utilized relative polarity measurements and machine learning techniques to better resolve focal mechanisms and stress orientations considering a catalog of ∼29,000 relocated earthquakes that occurred during 1984–2021 in the southeastern San Francisco Bay Area. Earthquake focal mechanisms are commonly produced using P wave first motion polarities, which traditionally requires events to be well-recorded across a seismic network with good focal sphere coverage. We adapted recently developed approaches that are less dependent on high signal-to-noise records and exploit similar waveforms to produce relative polarity and amplitude measurements between earthquake pairs. These techniques were previously only applied on localized earthquake sequences, and we further developed these approaches so that they can be utilized for regional catalogs. We validated or corrected manually identified polarities by performing polarity consensuses using earthquake pairs. Missing and unreliable polarity measurements were imputed using iterative random forests, an unsupervised ensemble machine learning method. Relative P and S wave amplitude measurements were made between earthquakes, constraining S / P ratios for low signal-to-noise waveforms. Using these techniques, we were able to reduce focal mechanism uncertainties by an average of ∼13° and produced well-constrained focal mechanisms for ∼6 times as many earthquakes than those produced using only the traditionally derived polarities. We performed stress inversions using the focal mechanisms by grouping the focal mechanism results into a quadtree structure. Our stress results are consistent with previous work, albeit at a higher spatial resolution, and demonstrate these techniques can aid our understanding of fault structures and kinematics in more detail than was previously possible.

California

Noise in two-color electronic distance meter measurements revisited

Frequent, high-precision geodetic data have temporally correlated errors. Temporal correlations directly affect both the estimate of rate and its standard error; the rate of deformation is a key product from geodetic measurements made in tectonically active areas. Various models of temporally correlated errors are developed and these provide relations between the power spectral density and the data covariance matrix. These relations are applied to two-color electronic distance meter (EDM) measurements made frequently in California over the past 15-20 years. Previous analysis indicated that these data have significant random walk error. Analysis using the noise models developed here indicates that the random walk model is valid for about 30% of the data. A second 30% of the data can be better modeled with power law noise with a spectral index between 1 and 2, while another 30% of the data can be modeled with a combination of band-pass-filtered plus random walk noise. The remaining 10% of the data can be best modeled as a combination of band-pass-filtered plus power law noise. This band-pass-filtered noise is a product of an annual cycle that leaks into adjacent frequency bands. For time spans of more than 1 year these more complex noise models indicate that the precision in rate estimates is better than that inferred by just the simpler, random walk model of noise.

Journal of Geophysical Research B: Solid Earth

Effect of wave-current interaction on wind-driven circulation in narrow, shallow embayments

The effect of wind waves on the steady wind-driven circulation in a narrow, shallow bay is investigated with a two-dimensional ( y , z ) circulation model and the Grant and Madsen [1979] bottom-boundary layer model, which includes wave-current interaction. A constant wind stress is applied in the along-channel x direction to a channel with a constant cross-sectional profile h ( y ). The wind-induced flushing of shallow bays is shown to be sensitive to both the shape of the cross section and the effects of surface waves. The flushing increases with increasing , where h ′ is the standard deviation of cross-channel depth and is the mean depth. This is consistent with the findings of Hearn et al. [1987]. The flushing decreases, however, with the inclusion of surface wave effects which act to increase the bottom drag felt by the currents. Increasing effective bottom friction reduces the strength of the circulation, while the along-bay surface slope, bottom stress and the structure of current profiles remain nearly unchanged. An implication of the circulation dependence on wave-current interaction is that low-frequency oscillatory winds may drive a mean circulation when the wave field changes with wind direction. x

Journal of Geophysical Research C: Oceans

Scarp degraded by linear diffusion: Inverse solution for age

Under the assumption that landforms unaffected by drainage channels are degraded according to the linear diffusion equation, a procedure is developed to invert a scarp profile to find its “diffusion age.” Diffusion age, having dimension [length] 2 , is the product of diffusivity times chronological age. The second moment of scarp slope grows linearly with age. This fact, together with an assumption about initial scarp shape, allows the inverse determination of diffusion age. Age found assuming a vertical initial scarp is termed “apparent age”; any nonvertical initial scarp profile has a nonzero initial apparent age. True diffusion age differs from apparent age by a fraction of scarp offset squared. The inverse procedure applied to synthetic data yields the following rules of thumb. Evidence of initial scarp shape has been lost when apparent age reaches twice its initial value. If a scarp is formed by two events, the inversion gives their offset-weighted-mean age with an error that is a fraction of offset squared. A scarp that appears to have been formed by one event may have been formed by two with an interval between them as large as apparent age. After scarps of two fault traces have diffused to appear as one, the error in inferred age may be as large as half the apparent age. Variation of apparent age along strike would indicate multiple fault traces. The simplicity of scarp profile measurement and this inversion makes profile analysis attractive. If linearity of the flow law, time for a free face to be reduced to the angle of repose, and variation of diffusivity with climate and material could be established, profile analysis would become a reliable dating technique.

Journal of Geophysical Research Solid Earth

3-D joint geodetic and strong-motion finite fault inversion of the 2008 May 12, Wenchuan, China Earthquake

We present a source inversion of the 2008 Wenchuan, China earthquake, using strong-motion waveforms and geodetic offsets together with three-dimensional synthetic ground motions. We applied the linear multiple time window technique considering geodetic and dynamic Green's functions computed with the finite element method and the reciprocity and Strain Green’s Tensor formalism. All ground motion estimates, valid up to 1 Hz, accounted for three-dimensional effects, including the topography and the geometry of the Beichuan and Pengguan faults. Our joint inversion has a higher moment (M0) than a purely geodetic inversion and the slip distribution presents differences when compared to one-dimensional model source inversions. The moment is estimated to be M0=1.2x1021 Nm, slightly larger than other works. Our results show that considering a complex 3D structure reduces the size of large areas of 10 m slip or greater by distributing it in wider zones, with reduced slips, in the central portion of the Beichuan and the Pengguan faults. Finally, we compare our source with a relocated aftershock catalog and conclude that the 4-5 m slip contours approximately bound the absence or presence of aftershocks.

Wenchuan

Coseismic slip resolution along a plate boundary megathrust: the Nankai Trough, southwest Japan

Geodetic survey measurements are used to estimate the coseismic slip distribution in the 1944 Tonankai (M w =8.1) and 1946 Nankaido (M w =8.3) earthquakes and to assess quantitatively the degree to which this slip is resolved on the plate boundary megathrust. Data used include 798 angle changes from triangulation surveys, 328 leveling section differences, and 5 coseismic tidal gage offsets. Many of the nominally coseismic triangulation data span ∼50 years, nearly half the earthquake cycle, and correction for interseismic deformation using post-1950 observations is applied. Microseismicity is used to define the configuration of the plate boundary interface and approximate it with a continuous, multisegment fault model. Because the onshore geodetic data have very limited resolving power for offshore fault segments, offshore coseismic slip was constrained by Satctke's [1993] estimation based on tsunami data. The majority of the coseismic slip occurs between 15 and 25 km depth. Although resolution declines toward the trench axis, it is sufficiently good to define two distinct high-slip regions, one off southeastern Shikoku Island (11 m maximum) and the other offshore of Kii Peninsula (3 m maximum). The slip magnitude off southeastern Shikoku, coupled with the plate convergence rate, would imply an recurrence interval of about 270 years, much-longer than the average repeat time of ∼120 years for historical great earthquakes on the Nankai Trough. However, the maximum coseismic slip is sensitive to the assumed fault geometry, and slippage on trough-parallel splay faults could significantly decrease the maximum slip to about 6 m.

Journal of Geophysical Research B: Solid Earth

Formation of fine sediment deposit from a flash flood river in the Mediterranean Sea

We identify the mechanisms controlling fine deposits on the inner-shelf in front of the Besòs River, in the northwestern Mediterranean Sea. This river is characterized by a flash flood regime discharging large amounts of water (more than 20 times the mean water discharge) and sediment in very short periods lasting from hours to few days. Numerical model output was compared with bottom sediment observations and used to characterize the multiple spatial and temporal scales involved in offshore sediment deposit formation. A high-resolution (50 m grid size) coupled hydrodynamic-wave-sediment transport model was applied to the initial stages of the sediment dispersal after a storm-related flood event. After the flood, sediment accumulation was predominantly confined to an area near the coastline as a result of preferential deposition during the final stage of the storm. Subsequent reworking occurred due to wave-induced bottom shear stress that resuspended fine materials, with seaward flow exporting them toward the midshelf. Wave characteristics, sediment availability, and shelf circulation determined the transport after the reworking and the final sediment deposition location. One year simulations of the regional area revealed a prevalent southwestward average flow with increased intensity downstream. The circulation pattern was consistent with the observed fine deposit depocenter being shifted southward from the river mouth. At the southern edge, bathymetry controlled the fine deposition by inducing near-bottom flow convergence enhancing bottom shear stress. According to the short-term and long-term analyses, a seasonal pattern in the fine deposit formation is expected.

Journal of Geophysical Research: Oceans

Dislocation pileup as a representation of strain accumulation on a strike-slip fault

The conventional model of strain accumulation on a vertical transform fault is a discrete screw dislocation in an elastic half-space with the Burgers vector of the dislocation increasing at the rate of relative plate motion. It would be more realistic to replace that discrete dislocation by a dislocation distribution, presumably a pileup in which the individual dislocations are in equilibrium. The length of the pileup depends upon the applied stress and the amount of slip that has occurred at depth. I argue here that the dislocation pileup (the transition on the fault from no slip to slip at the full plate rate) occupies a substantial portion of the lithosphere thickness. A discrete dislocation at an adjustable depth can reproduce the surface deformation profile predicted by a pileup so closely that it will be difficult to distinguish between the two models. The locking depth (dislocation depth) of that discrete dislocation approximation is substantially (???30%) larger than that (depth to top of the pileup) in the pileup model. Thus, in inverting surface deformation data using the discrete dislocation model, the locking depth in the model should not be interpreted as the true locking depth. Although dislocation pileup models should provide a good explanation of the surface deformation near the fault trace, that explanation may not be adequate at greater distances from the fault trace because approximating the expected horizontally distributed deformation at subcrustal depths by uniform slip concentrated on the fault is not justified.

Journal of Geophysical Research B: Solid Earth

Earthquake-derived seismic velocity changes during the 2018 caldera collapse of Kīlauea volcano

The 2018 Kīlauea caldera collapse produced extraordinary sequences of seismicity and deformation, with 62 episodic collapse events which significantly altered the landscape of the summit region. Despite decades of focused scientific studies at Kīlauea, detailed information about the internal structure of the volcano is limited. Recently developed techniques in seismic interferometry can be used to monitor the internal structure of an active volcano more directly by detecting subtle spatiotemporal changes in seismic wave velocity, but their utility relies on accurate interpretations of the underlying phenomena causing those velocity changes. Here, we retrospectively apply repeating-earthquake-based seismic interferometry to the 2018 Kīlauea eruption sequence. We find that seismic velocities changed over two distinct time scales: a sudden increase followed by a slower decrease in velocity in the hours following each collapse event, and a gradual, long-term decrease in velocity over several weeks that ceased approximately 1 month prior to the end of the eruption. Modeling suggests that short-term changes can be explained by magma reservoir pressurization which specifically closed vertical ring fractures. Long-term changes are related to subsidence of the caldera and likely include the influence of inelastic strain from the formation of new fractures. These observations provide new insights into the evolution of Kīlauea during its progressive collapse and will inform future interpretations for near-real-time monitoring at hazardous volcanoes around the world using similar techniques, especially where a dominant fracture orientation is present.

Hawaii

Forecast model for moderate earthquakes near Parkfield, California

Earthquake instability models have possible application to earthquake forecasting because the models simulate both preseismic and coseismic changes of fault slip and ground deformation. In the forecast procedure proposed here, repeated measurements of preseismic fault slip and ground deformation constrain the values of model parameters. The early part of the model simulation corresponds to the available field data, and the subsequent part constitutes an estimate of future faulting and ground deformation. In particular, the time, location, and size of unstable faulting are estimates of the pending earthquake parameters. The forecast accuracy depends on the model realism and parameter resolution. The forecast procedure is applied to fault creep and trilateration data measured near Parkfield, California, where at least five magnitude 5.5 to 6 earthquakes have occurred regularly since 1881, the last in 1966. The quasi-static model consists of a flat vertical plane embedded in an elastic half space. Spacially variable fault slip of strike-slip sense is driven by an increasing regional shear stress but is impeded by a relatively strong patch of brittle, strain-softening fault. The field data are consistent with these approximate values of patch parameters: radius of 3 km, patch center 5 km deep and 8 km southeast of the 1966 epicenter, and maximum brittle strength of 26 bars. Fluctuations in the available field data prevent estimating the earthquake time with any more precision than use of the 21±8 year recurrence interval. However, the model may later give a more precise estimate of the earthquake time if the fault slip rate near the inferred patch increases before the earthquake, as predicted by the model.

Journal of Geophysical Research Solid Earth

Structure of the California Coast Ranges and San Andreas Fault at SAFOD from seismic waveform inversion and reflection imaging

[1] A seismic reflection and refraction survey across the San Andreas Fault (SAF) near Parkfield provides a detailed characterization of crustal structure across the location of the San Andreas Fault Observatory at Depth (SAFOD). Steep-dip prestack migration and frequency domain acoustic waveform tomography were applied to obtain highly resolved images of the upper 5 km of the crust for 15 km on either side of the SAF. The resulting velocity model constrains the top of the Salinian granite with great detail. Steep-dip reflection seismic images show several strong-amplitude vertical reflectors in the uppermost crust near SAFOD that define an ∼2-km-wide zone comprising the main SAF and two or more local faults. Another prominent subvertical reflector at 2–4 km depth ∼9 km to the northeast of the SAF marks the boundary between the Franciscan terrane and the Great Valley Sequence. A deep seismic section of low resolution shows several reflectors in the Salinian crust west of the SAF. Two horizontal reflectors around 10 km depth correlate with strains of seismicity observed along-strike of the SAF. They represent midcrustal shear zones partially decoupling the ductile lower crust from the brittle upper crust. The deepest reflections from ∼25 km depth are interpreted as crust-mantle boundary.

Journal of Geophysical Research B: Solid Earth

Submarine permafrost map in the arctic modelled using 1D transient heat flux (SuPerMAP)

Offshore permafrost plays a role in the global climate system, but observations of permafrost thickness, state, and composition are limited to specific regions. The current global permafrost map shows potential offshore permafrost distribution based on bathymetry and global sea level rise. As a first‐order estimate, we employ a heat transfer model to calculate the subsurface temperature field. Our model uses dynamic upper boundary conditions that synthesize Earth System Model air temperature, ice mass distribution and thickness, and global sea level reconstruction and applies globally distributed geothermal heat flux as a lower boundary condition. Sea level reconstruction accounts for differences between marine and terrestrial sedimentation history. Sediment composition and pore water salinity are integrated in the model. Model runs for 450 ka for cross‐shelf transects were used to initialize the model for circumarctic modeling for the past 50 ka. Preindustrial submarine permafrost (i.e., cryotic sediment), modeled at 12.5‐km spatial resolution, lies beneath almost 2.5 ×10 6 km 2 of the Arctic shelf. Our simple modeling approach results in estimates of distribution of cryotic sediment that are similar to the current global map and recent seismically delineated permafrost distributions for the Beaufort and Kara seas, suggesting that sea level is a first‐order determinant for submarine permafrost distribution. Ice content and sediment thermal conductivity are also important for determining rates of permafrost thickness change. The model provides a consistent circumarctic approach to map submarine permafrost and to estimate the dynamics of permafrost in the past.

Journal of Geophysical Research C: Oceans