Search USGS⌕ Search

SEARCH · Search USGS

Results for “Geophysical Journal”

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

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

At least 289 records · Page 16Linked to original sources
USGS2006

Effect of tree line advance on carbon storage in NW Alaska

We investigated the size, distribution, and temporal dynamics of ecosystem carbon (C) pools in an area of recent tree line advance, northwest Alaska. Repeat aerial photographs show forest cover increased ???10% in our study area since 1949. We sampled C pools of four principal ecosystem types, tussock tundra, shrub tundra, woodland, and forest, all located on a 600-800 year old river terrace. Significant differences between ecosystem C pools, both above ground and below ground existed. Tundra sites store >22.2 kg C/m2, shrub tundra sites and woodland sites store 9.7 kg C/m2 and 14.3 kg C/m2, respectively, and forest sites store 14.4 kg C/m2. Landscape variation of total ecosystem C was primarily due to organic soil C and was secondarily due to C stored in trees. Soil C/N profiles of shrub tundra sites and woodland sites showed similarities with forest site soils at surface and tundra site soils at depth. We hypothesize that tundra systems transformed to forest systems in this area under a progression of permafrost degradation and enhanced drainage. On the basis of C pool estimates for the different ecosystem types, conversion of tundra sites to forest may have resulted in a net loss of > 7.8 kg C/m2, since aboveground C gains were more than offset by belowground C losses to decomposition in the tundra sites. Tree line advance therefore might not increase C storage in high-latitude ecosystems and thus might not, as previously suggested, act as a negative feedback to warming. Key to this hypothesis and to its projection to future climate response is the fate of soil carbon upon warming and permafrost drainage. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2008

Constraining fault constitutive behavior with slip and stress heterogeneity

We study how enforcing self-consistency in the statistical properties of the preshear and postshear stress on a fault can be used to constrain fault constitutive behavior beyond that required to produce a desired spatial and temporal evolution of slip in a single event. We explore features of rupture dynamics that (1) lead to slip heterogeneity in earthquake ruptures and (2) maintain these conditions following rupture, so that the stress field is compatible with the generation of aftershocks and facilitates heterogeneous slip in subsequent events. Our three-dimensional fmite element simulations of magnitude 7 events on a vertical, planar strike-slip fault show that the conditions that lead to slip heterogeneity remain in place after large events when the dynamic stress drop (initial shear stress) and breakdown work (fracture energy) are spatially heterogeneous. In these models the breakdown work is on the order of MJ/m2, which is comparable to the radiated energy. These conditions producing slip heterogeneity also tend to produce narrower slip pulses independent of a slip rate dependence in the fault constitutive model. An alternative mechanism for generating these confined slip pulses appears to be fault constitutive models that have a stronger rate dependence, which also makes them difficult to implement in numerical models. We hypothesize that self-consistent ruptures could also be produced by very narrow slip pulses propagating in a self-sustaining heterogeneous stress field with breakdown work comparable to fracture energy estimates of kJ/M2. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Constitutive relationships and physical basis of fault strength due to flash heating

We develop a model of fault strength loss resulting from phase change at asperity contacts due to flash heating that considers a distribution of contact sizes and nonsteady state evolution of fault strength with displacement. Laboratory faulting experiments conducted at high sliding velocities, which show dramatic strength reduction below the threshold for bulk melting, are well fit by the model. The predicted slip speed for the onset of weakening is in the range of 0.05 to 2 m/s, qualitatively consistent with the limited published observations. For this model, earthquake stress drops and effective shear fracture energy should be linearly pressure-dependent, whereas the onset speed may be pressure-independent or weakly pressure-dependent. On the basis of the theory, flash weakening is expected to produce large dynamic stress drops, small effective shear fracture energy, and undershoot. Estimates of the threshold slip speed, stress drop, and fracture energy are uncertain due to poor knowledge of the average ontact dimension, shear zone thickness and gouge particle size at seismogenic depths. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Constraining earthquake source inversions with GPS data: 2. A two-step approach to combine seismic and geodetic data sets

We present a new method to combine static and wavefield data to image earthquake ruptures. Our combined inversion is a two-step procedure, following the work of Hernandez et al. (1999), and takes into account the differences between the resolutions of the two data sets. The first step consists of an inversion of the static field, which yields a map of slip amplitude. This inversion exploits a special irregular grid that takes into account the resolution of the static data. The second step is an inversion of the radiated wavefield; it results in the determination of the time evolution of slip on the fault. In the second step, the slip amplitude is constrained to resemble the static slip amplitude map inferred from the GPS inversion. Using this combined inversion, we study the source process of the 2004 M6 Parkfield, California, earthquake. We conclude that slip occurred in two main regions of the fault, each of which displayed distinct rupture behaviors. Slip initiated at the hypocenter with a very strong bilateral burst of energy. Here, slip was localized in a narrow area approximately 10 km long, the rupture velocity was very fast (???3.5 km/s), and slip only lasted a short period of time (<1 s). Then the rupture proceeded to a wider region 12-20 km northwest of the hypocenter. Here, the earthquake developed in a more moderated way: the rupture velocity slowed to ???3.0 km/s and slip lasted longer (1-2 s). The maximum slip amplitude was 0.45 m. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Generalized analytical model for benthic water flux forced by surface gravity waves

A generalized analytical model for benthic water flux forced by linear surface gravity waves over a series of layered hydrogeologic units is developed by adapting a previous solution for a hydrogeologic unit with an infinite thickness (Case I) to a unit with a finite thickness (Case II) and to a dual-unit system (Case III). The model compares favorably with laboratory observations. The amplitude of wave-forced benthic water flux is shown to be directly proportional to the amplitude of the wave, the permeability of the hydrogeologic unit, and the wave number and inversely proportional to the kinematic viscosity of water. A dimensionless amplitude parameter is introduced and shown to reach a maximum where the product of water depth and the wave number is 1.2. Submarine groundwater discharge (SGD) is a benthic water discharge flux to a marine water body. The Case I model estimates an 11.5-cm/d SGD forced by a wave with a 1 s period and 5-cm amplitude in water that is 0.5-m deep. As this wave propagates into a region with a 0.3-m-thick hydrogeologic unit, with a no-flow bottom boundary, the Case II model estimates a 9.7-cm/d wave-forced SGD. As this wave propagates into a region with a 0.2-m-thick hydrogeologic unit over an infinitely thick, more permeable unit, the Case III quasi-confined model estimates a 15.7-cm/d wave-forced SGD. The quasi-confined model has benthic constituent flux implications in coral reef, karst, and clastic regions. Waves may undermine tracer and seepage meter estimates of SGD at some locations. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2009

Constraints on deep moonquake focal mechanisms through analyses of tidal stress

[1] A relationship between deep moonquake occurrence and tidal forcing is suggested by the monthly periodicities observed in the occurrence times of events recorded by the Apollo Passive Seismic Experiment. In addition, the typically large S wave to P wave arrival amplitude ratios observed on deep moonquake seismograms are indicative of shear failure. Tidal stress, induced in the lunar interior by the gravitational influence of the Earth, may influence moonquake activity. We investigate the relationship between tidal stress and deep moonquake occurrence by searching for a linear combination of the normal and shear components of tidal stress that best approximates a constant value when evaluated at the times of moonquakes from 39 different moonquake clusters. We perform a grid search at each cluster location, computing the stresses resolved onto a suite of possible failure planes, to obtain the best fitting fault orientation at each location. We find that while linear combinations of stresses (and in some cases stress rates) can fit moonquake occurrence at many clusters quite well; for other clusters, the fit is not strongly dependent on plane orientation. This suggests that deep moonquakes may occur in response to factors other than, or in addition to, tidal stress. Several of our inferences support the hypothesis that deep moonquakes might be related to transformational faulting, in which shear failure is induced by mineral phase changes at depth. The occurrence of this process would have important implications for the lunar interior. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2009

Simulating and understanding sand wave variation: A case study of the Golden Gate sand waves

In this paper we present a detailed comparison between measured features of the Golden Gate sand wave field and the results of a nonlinear sand wave model. Because the Golden Gate sand waves exhibit large variation in their characteristics and in their environmental physics, this area gives us the opportunity to study sand wave variation between locations, within one well-measured, large area. The nonlinear model used in this paper is presently the only tool that provides information on the nonlinear evolution of large-amplitude sand waves. The model is used to increase our understanding of the coupling between the variability in environmental conditions and the sand wave characteristics. Results show that the model is able to describe the variation in the Golden Gate sand waves well when both the local oscillating tidal current and the residual current are taken into account. Current and water depth seem to be the most important factors influencing sand wave characteristics. The simulation results give further confidence in the underlying model hypothesis and assumptions. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2010

A rain splash transport equation assimilating field and laboratory measurements

Process-based models of hillslope evolution require transport equations relating sediment flux to its major controls. An equation for rain splash transport in the absence of overland flow was constructed by modifying an approach developed by Reeve (1982) and parameterizing it with measurements from single-drop laboratory experiments and simulated rainfall on a grassland in East Africa. The equation relates rain splash to hillslope gradient, the median raindrop diameter of a storm, and ground cover density; the effect of soil texture on detachability can be incorporated from other published results. The spatial and temporal applicability of such an equation for rain splash transport in the absence of overland flow on uncultivated hillslopes can be estimated from hydrological calculations. The predicted transport is lower than landscape-averaged geologic erosion rates from Kenya but is large enough to modify short, slowly eroding natural hillslopes as well as microtopographic interrill surfaces between which overland flow transports the mobilized sediment. Copyright 2010 by the American Geophysical Union. Copyright 2010 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2011

The magnitude distribution of earthquakes near Southern California faults

We investigate seismicity near faults in the Southern California Earthquake Center Community Fault Model. We search for anomalously large events that might be signs of a characteristic earthquake distribution. We find that seismicity near major fault zones in Southern California is well modeled by a Gutenberg-Richter distribution, with no evidence of characteristic earthquakes within the resolution limits of the modern instrumental catalog. However, the b value of the locally observed magnitude distribution is found to depend on distance to the nearest mapped fault segment, which suggests that earthquakes nucleating near major faults are likely to have larger magnitudes relative to earthquakes nucleating far from major faults. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Process-based modeling of tsunami inundation and sediment transport

The infrequent and unpredictable nature of tsunamis precludes the use of field experiments to measure the hydrodynamic and sediment transport processes that occur. Instead, these processes are often approximated from laboratory, numerical, and theoretical studies or inferred from observations of the resultant sediment deposits. Here Delft3D, a three-dimensional numerical model, is used to simulate the inundation and sediment transport of a tsunami similar in magnitude to the 26 December 2004 Indian Ocean tsunami over one measured and three idealized morphologies. The model is first shown to match well the observations taken at Kuala Meurisi, Sumatra, and then used to examine in detail the processes that occur during the tsunami. The model predicts that at a given cross-shore location the onshore flow accelerates rapidly to a maximum as the wavefront passes, and then gradually decelerates before reversing direction and flowing offshore. The onshore flow does not tend to zero everywhere at maximum inundation, but instead flow reversal occurs near the shoreline even as the wavefront continues to inundate landward. While some sediment is eroded by the passing wavefront, the suspension of sandy sediment is dominated by the long-duration, high-velocity backwash that occurs along the beach face and offshore of the shoreline. Some of the sediment suspended during backwash is advected shoreward by the subsequent wave, creating large spatial gradients in the suspended sediment concentrations, which may not be in equilibrium with the local hydrodynamics. The inundation and transport of sediment during a tsunami can be affected by complexities in the morphological profile and interactions between multiple waves, and many of the hydrodynamic and sediment transport processes predicted here are similar to analogous processes previously observed in the swash zone. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2011

A Bayesian network to predict coastal vulnerability to sea level rise

Sea level rise during the 21st century will have a wide range of effects on coastal environments, human development, and infrastructure in coastal areas. The broad range of complex factors influencing coastal systems contributes to large uncertainties in predicting long-term sea level rise impacts. Here we explore and demonstrate the capabilities of a Bayesian network (BN) to predict long-term shoreline change associated with sea level rise and make quantitative assessments of prediction uncertainty. A BN is used to define relationships between driving forces, geologic constraints, and coastal response for the U.S. Atlantic coast that include observations of local rates of relative sea level rise, wave height, tide range, geomorphic classification, coastal slope, and shoreline change rate. The BN is used to make probabilistic predictions of shoreline retreat in response to different future sea level rise rates. Results demonstrate that the probability of shoreline retreat increases with higher rates of sea level rise. Where more specific information is included, the probability of shoreline change increases in a number of cases, indicating more confident predictions. A hindcast evaluation of the BN indicates that the network correctly predicts 71% of the cases. Evaluation of the results using Brier skill and log likelihood ratio scores indicates that the network provides shoreline change predictions that are better than the prior probability. Shoreline change outcomes indicating stability (-1 < rate < 1 m/yr) or erosion (rate < -1 m/yr) tend to occur for two sets of input scenarios. Stable shoreline change rates occur mainly for low rates of relative sea level rise and occur in low-vulnerability geomorphic settings. Rates indicating erosion result for cases where the rate of relative sea level rise is high and moderate-to-high vulnerability geomorphic settings occur. In contrast, accretion (rate > 1 m/yr) was not well predicted. We find that BNs can assimilate important factors contributing to coastal change in response to sea level rise and can make quantitative, probabilistic predictions that can be applied to coastal management decisions. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2012

A review of the regional geophysics of the Arizona Transition Zone

A review of existing geophysical information and new data presented in this special section indicate that major changes in crustal properties between the Basin and Range and Colorado Plateau occur in, or directly adjacent to, the region defined as the Arizona Transition Zone. Although this region was designated on a physiographic basis, studies indicate that it is also the geophysical transition between adjoining provinces. The Transition Zone displays anomalous crustal and upper mantle seismic properties, shallow Curie isotherms, high heat flow, and steep down-to-the-plateau Bouguer gravity gradients. Seismic and gravity studies suggest that the change in crustal thickness, from thin crust in the Basin and Range to thick crust in the Colorado Plateau, may occur as a series of steps rather than a planar surface. Anomalous P wave velocities, high heat flow, shallow Curie isotherms, and results of gravity modeling suggest that the upper mantle is heterogeneous in this region. A relatively shallow asthenosphere beneath the Basin and Range and Transition Zone contrasted with a thick lithosphere beneath the Colorado Plateau would be one explanation that would satisfy these geophysical observations.

Journal of Geophysical Research↗
USGS2003

Fault slip and seismic moment of the 1700 Cascadia earthquake inferred from Japanese tsunami descriptions

The 1700 Cascadia earthquake attained moment magnitude 9 according to new estimates based on effects of its tsunami in Japan, computed coseismic seafloor deformation for hypothetical ruptures in Cascadia, and tsunami modeling in the Pacific Ocean. Reports of damage and flooding show that the 1700 Casscadia tsunami reached 1-5 m heights at seven shoreline sites in Japan. Three sets of estimated heights express uncertainty about location and depth of reported flooding, landward decline in tsunami heights from shorelines, and post-1700 land-level changes. We compare each set with tsunami heights computed from six Cascadia sources. Each source is vertical seafloor displacement calculated with a three-dimensional elastic dislocation model, for three sources the rupture extends the 1100 km length of the subduction zone and differs in width and shallow dip; for the other sources, ruptures of ordinary width extend 360-670 km. To compute tsunami waveforms, we use a linear long-wave approximation with a finite difference method, and we employ modern bathymetry with nearshore grid spacing as small as 0.4 km. The various combinations of Japanese tsunami heights and Cascadia sources give seismic moment of 1-9 ?? 1022 N m, equivalent to moment magnitude 8.7-9.2. This range excludes several unquantified uncertainties. The most likely earthquake, of moment magnitude 9.0, has 19 m of coseismic slip on an offshore, full-slip zone 1100 km long with linearly decreasing slip on a downdip partial-slip zone. The shorter rupture models require up to 40 m offshore slip and predict land-level changes inconsistent with coastal paleoseismological evidence. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Constraints on the viscosity of the continental crust and mantle from GPS measurements and postseismic deformation models in western Mongolia

We use GPS measurements and models of postseismic deformation caused by seven M6.8 to 8.4 earthquakes that occurred in the past 100 years in Mongolia to assess the viscosity of the lower crust and upper mantle. We find an upper mantle viscosity between 1 ?? 1018 and 4 ?? 1018 Pa s. The presence of such a weak mantle is consistent with results from independent seismological and petrological studies that show an abnormally hot upper mantle beneath Mongolia. The viscosity of the lower crust is less well constrained, but a weak lower crust (3 ?? 1016 to 2 ?? 1017 Pa s) is preferred by the data. Using our best fit upper mantle and lower crust viscosities, we find that the postseismic effects of viscoelastic relaxation on present-day horizontal GPS velocities are small (<2 mm yr-1) but still persist 100 years after the 1905, M8.4, Bolnay earthquake. This study shows that the GPS velocity field in the Baikal-Mongolia area can be modeled as the sum of (1) a rigid translation and rotation of the whole network, (2) a 3-5 mm yr-1 simple shear velocity gradient between the Siberian platform to the north and northern China to the south, and (3) the contribution of postseismic deformation, mostly caused by the 1905 Bolnay-Tsetserleg sequence and by the smaller, but more recent, 1957 Bogd earthquake. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Fault interaction and stress triggering of twentieth century earthquakes in Mongolia

A cluster of exceptionally large earthquakes in the interior of Asia occurred from 1905 to 1967: the 1905 M7.9 Tsetserleg and M8.4 Bolnai earthquakes, the 1931 M8.0 Fu Yun earthquake, the 1957 M8.1 Gobi-Altai earthquake, and the 1967 M7.1 Mogod earthquake (sequence). Each of the larger (M ??? 8) earthquakes involved strike-slip faulting averaging more than 5 m and rupture lengths of several hundred kilometers. Available geologic data indicate that recurrence intervals on the major source faults are several thousands of years and distances of about 400 km separate the respective rupture areas. We propose that the occurrences of these and many smaller earthquakes are related and controlled to a large extent by stress changes generated by the compounded static deformation of the preceding earthquakes and subsequent viscoelastic relaxation of the lower crust and upper mantle beneath Mongolia. We employ a spherically layered viscoelastic model constrained by the 1994-2002 GPS velocity field in western Mongolia [Vergnolle et al., 2003]. Using the succession of twentieth century earthquakes as sources of deformation, we then analyze the time-dependent change in Coulomb failure stress (????f). At remote interaction distances, static ????f values are small. However, modeled postseismic stress changes typically accumulate to several tenths of a bar over time intervals of decades. Almost all significant twentieth century regional earthquakes (M ??? 6) with well-constrained fault geometry lie in positive ????f lobes of magnitude about +0.5 bar. Our results suggest that significant stress transfer is possible among continental faults separated by hundreds of kilometers and on timescales of decades. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Geodetic and seismic constraints on some seismogenic zone processes in Costa Rica

New seismic and geodetic data from Costa Rica provide insight into seismogenic zone processes in Central America, where the Cocos and Caribbean plates converge. Seismic data are from combined land and ocean bottom deployments in the Nicoya peninsula in northern Costa Rica and near the Osa peninsula in southern Costa Rica. In Nicoya, inversion of GPS data suggests two locked patches centered at 14 ± 2 and 39 ± 6 km depth. Interplate microseismicity is concentrated in the more freely slipping intermediate zone, suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone, the rupture zone for future large earthquakes, at least over the short (∼1 year) observation period. We also estimate northwest motion of a coastal “sliver block” at 8 ± 3 mm/yr, probably related to oblique convergence. In the Osa region to the south, convergence is orthogonal to the trench. Cocos-Caribbean relative motion is partitioned here, with ∼8 cm/yr on the Cocos-Panama block boundary (including a component of permanent shortening across the Fila Costeña fold and thrust belt) and ∼1 cm/yr on the Panama block–Caribbean boundary. The GPS data suggest that the Cocos plate–Panama block boundary is completely locked from ∼10–50 km depth. This large locked zone, as well as associated forearc and back-arc deformation, may be related to subduction of the shallow Cocos Ridge and/or younger lithosphere compared to Nicoya, with consequent higher coupling and compressive stress in the direction of plate convergence.

Nicoya Peninsula, Osa Peninsula↗
USGS2008

First in situ investigation of a dark wind streak on Mars

Prominent low‐albedo wind streaks issue from embayments at the north end of Victoria crater in Meridiani Planum, the site of surface investigations by the Mars Exploration Rover Opportunity. Opportunity was sent to examine the darkest of these wind streaks and compare it to the adjacent soil in the brighter corridor between the streaks. Two nearby sites inside and outside of the streak were selected for visible and infrared remote sensing and for elemental abundance measurements and images of the microscopic soil morphology. The rover was next sent to study the interaction of the wind streak with an obstacle, a small rock near the source of the streak. Opportunity's observations suggest that the Victoria wind streaks are deposits of basaltic sand blown out of the crater from the dark dunes nestled below the crater rim, particularly at the base of the alcove leading up to the darkest streak. No local sources of sand have been identified within the Victoria crater wall rock, suggesting that the sand is supplied from outside the crater and is presently escaping its temporary topographic trap. This process presents a possible explanation for the serrated margin of Victoria crater through abrasion of the soft rock as trapped sands are blown out of the crater and carve alcoves under various seasonal winds. Some indication of the rate of activity of these wind streaks will be provided by continued monitoring of the erasure of the rover's tracks.

Journal of Geophysical Research E: Planets↗
USGS2001

Moss and lichen cover mapping at local and regional scales in the boreal forest ecosystem of central Canada

Mosses and lichens are important components of boreal landscapes [Vitt et al., 1994; Bubier et al., 1997]. They affect plant productivity and belowground carbon sequestration and alter the surface runoff and energy balance. We report the use of multiresolution satellite data to map moss and lichens over the BOREAS region at a 10 m, 30 m, and 1 km scales. Our moss and lichen classification at the 10 m scale is based on ground observations of associations among soil drainage classes, overstory composition, and cover type among four broad classes of ground cover (feather, sphagnum, and brown mosses and lichens). For our 30 m map, we used field observations of ground cover-overstory associations to map mosses and lichens in the BOREAS southern study area (SSA). To scale up to a 1 km (AVHRR) moss map of the BOREAS region, we used the TM SSA mosaics plus regional field data to identify AVHRR overstory-ground cover associations. We found that: 1) ground cover, overstory composition and density are highly correlated, permitting inference of moss and lichen cover from satellite-based land cover classifications; 2) our 1 km moss map reveals that mosses dominate the boreal landscape of central Canada, thereby a significant factor for water, energy, and carbon modeling; 3) TM and AVHRR moss cover maps are comparable; 4) satellite data resolution is important; particularly in detecting the smaller wetland features, lakes, and upland jack pine sites; and 5) distinct regional patterns of moss and lichen cover correspond to latitudinal and elevational gradients. Copyright 2001 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗