USGS Science⌕ Search

SEARCH · USGS Science

Results for “The R 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 1,009 records · Page 56Linked to original sources
USGS2021

Carbon fluxes and microbial activities from boreal peatlands experiencing permafrost thaw

Permafrost thaw in northern ecosystems may cause large quantities of carbon (C) to move from soil to atmospheric pools. Because soil microbial communities play a critical role in regulating C fluxes from soils, we examined microbial activity and greenhouse gas production soon after permafrost thaw and ground collapse (into collapse-scar bogs), relative to the permafrost plateau or older thaw features. Using multiple field and laboratory-based assays at a field site in interior Alaska, we show that the youngest collapse-scar bog had the highest CH 4 production potential from soil incubations, and, based upon temporal changes in porewater concentrations and 13 C-CH 4 and 13 C-CO 2 , had greater summer in situ rates of respiration, methanogenesis, and surface CH 4 oxidation. These patterns could be explained by greater C and N availability in the young bog, while alternative terminal electron accepting processes did not play a significant role. Field diffusive CH 4 fluxes from the young bog were 4.1 times greater in the shoulder season and 1.7–7.2 times greater in winter relative to older bogs, but not during summer. Greater relative CH 4 flux rates in the shoulder season and winter could be due to reduced CH 4 oxidation relative to summer, magnifying the importance of differences in production. Both the permafrost plateau and collapse-scar bogs were sources of C to the atmosphere due in large part to winter C fluxes. In collapse scar bogs, winter is a critical period when differences in thermokarst age translates to differences in surface fluxes.

Alaska↗
USGS2022

Insights on multistage rock avalanche behavior from runout modeling constrained by seismic inversions

Inversion of low-frequency regional seismic records to solve for a time series of bulk forces exerted on the earth by a landslide (a force-time function) is increasingly being used to infer volumes and dynamics of large, highly energetic landslides, such as rock avalanches and flowslides, and to provide calibration information on event dynamics and volumes for numerical landslide runout models. Much of the work to date using landslide runout modeling constrained by seismic data has focused on using single-phase models with frictional or velocity-weakening rheologies. Awareness of multistage landslide initiations is increasing, with discrete failures separated in time contributing to the final impact of an event. Our work utilizes a method for incorporating seismic data as a calibration constraint for landslide runout models, considering variable rheologies and different initiation conditions. This study presents a systematic examination of multiple rheologies and initiation conditions, and shows how these factors affect the force-time function derived from the landslide runout model. Our work confirms that, while rheology and fragmenting or initially coherent initiations affect the force-time function, multiple collapses separated by tens of seconds have the greatest impact on the shape and amplitude. We apply this method to the analysis of three real rock avalanches to better constrain plausible initiation conditions and rheology parameters using both seismic and field data. This study provides insights on how assumptions about the initiation dynamics of the source zone and the runout model definition can aid in the interpretation of seismic inversions for multistage rock avalanches.

Journal of Geophysical Research: Solid Earth↗
USGS1996

Gravity field over the Sea of Galilee: Evidence for a composite basin along a transform fault

The Sea of Galilee (Lake Kinneret) is located at the northern portion of the Kinneret-Bet Shean basin, in the northern Dead Sea transform. Three hundred kilometers of continuous marine gravity data were collected in the lake and integrated with land gravity data to a distance of more than 20 km around the lake. Analyses of the gravity data resulted in a free-air anomaly map, a variable density Bouguer anomaly map, and a horizontal first derivative map of the Bouguer anomaly. These maps, together with gravity models of profiles across the lake and the area south of it, were used to infer the geometry of the basins in this region and the main faults of the transform system. The Sea of Galilee can be divided into two units. The southern half is a pull-apart that extends to the Kinarot Valley, south of the lake, whereas the northern half was formed by rotational opening and transverse normal faults. The deepest part of the basinal area is located well south of the deepest bathymetric depression. This implies that the northeastern part of the lake, where the bathymetry is the deepest, is a young feature that is actively subsiding now. The pull-apart basin is almost symmetrical in the southern part of the lake and in the Kinarot Valley south of the lake. This suggests that the basin here is bounded by strike-slip faults on both sides. The eastern boundary fault extends to the northern part of the lake, while the western fault does not cross the northern part. The main factor controlling the structural complexity of this area is the interaction of the Dead Sea transform with a subperpendicular fault system and rotated blocks.

Sea of Galilee↗
USGS1997

Continuous monitoring of surface deformation at Long Valley Caldera, California, with GPS

Continuous Global Positioning System (GPS) measurements at Long Valley Caldera, an active volcanic region in east central California, have been made on the south side of the resurgent dome since early 1993. A site on the north side of the dome was added in late 1994. Special adaptations for autonomous operation in remote regions and enhanced vertical precision were made. The data record ongoing volcanic deformation consistent with uplift and expansion of the surface above a shallow magma chamber. Measurement precisions (1 standard error) for “absolute” position coordinates, i.e., relative to a global reference frame, are 3–4 mm (north), 5–6 mm (east), and 10–12 mm (vertical) using 24 hour solutions. Corresponding velocity uncertainties for a 12 month period are about 2 mm/yr in the horizontal components and 3–4 mm/yr in the vertical component. High precision can also be achieved for relative position coordinates on short (less than 10 km) baselines using broadcast ephemerides and observing times as short as 3 hours, even when data are processed rapidly on site. Comparison of baseline length changes across the resurgent dome between the two GPS sites and corresponding two-color electronic distance measurements indicates similar extension rates within error (∼2 mm/yr) once we account for a random walk noise component in both systems that may reflect spurious monument motion. Both data sets suggest a pause in deformation for a 3.5 month period in mid-1995, when the extension rate across the dome decreased essentially to zero. Three dimensional positioning data from the two GPS stations suggest a depth (5.8±1.6 km) and location (west side of the resurgent dome) of a major inflation center, in agreement with other geodetic techniques, near the top of a magma chamber inferred from seismic data. GPS systems similar to those installed at Long Valley can provide a practical method for near real-time monitoring and hazard assessment on many active volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Radar and photoclinometric studies of wrinkle ridges on Mars

Earth-based radar altimetry and image derived photoclinometric profiles were analyzed to examine both the long- and short-wavelength topography associated with wrinkle ridges on Mars. Photoclinometrically derived elevation data across wrinkle ridges were evaluated to determine the sensitivity of profiles to two empirical photoclinometric parameters, the horizontal digital number (HDN) and the scattered light value (SLV). The photoclinometric profiles are extremely sensitive to small variations in HDN. The sense of slope of a profile can be completely reversed over a range in HDN of as little as ±1. Comparably small variations in the SLV have relatively minor effects on the photoclinometrically derived elevations. The existence of elevation offsets from one side of the ridge to the other, reported in previous photoclinometric studies of martian wrinkles, were not confirmed through photoclinometry. In addition, no evidence of elevation offsets were found in Earth-based radar altimetry profiles across wrinkle ridges. In order to more accurately model wrinkle ridge topography, we controlled photoclinometrically derived elevations with long-wavelength topography obtained from the radar altimetry. The results of this study do not support kinematic models for the origin of planetary wrinkle ridges that involve deeply rooted thrust faults which separate crustal blocks at different elevations. A kinematic model involving buckling of shallow crustal layers into concentric folds that close, leading to the development of thrust faults, is consistent with wrinkle ridge morphology and terrestrial analogs. Recent geophysical studies of terrestrial analogs and the influence of shallow subsurface structures, particularly buried craters, on the localization of many wrinkle ridges on Mars suggest that thrust faults associated with the ridges are confined to the ridged plains material and do not extend into the lithosphere.

Journal of Geophysical Research E: Planets↗
USGS2012

Two-mica granites of northeastern Nevada

The field settings are described and analytical data are presented for six two-mica granites from north-eastern Nevada. High δ 18 O and 87 Sr/ 86 Sr values indicate that all are S-type granites, derived from continental crust. The major element chemistry and accessory mineral contents of these rocks also are characteristic of S-type granites. Chemical, X ray, and other data are presented for the micas recovered from these granites. The muscovites are notably high in Fe 2 O 3 , FeO, and MgO. Except for one hydrobiotite, each of the biotites has an MgO content near 6.0 weight percent. Two different types of two-mica granites are recognized in the area of this study. One type is distinguished by the presence of many biotite euhedra within muscovite phenocrysts and by an unusual suite of accessory minerals completely devoid of opaque oxides. This type probably resulted from anatexis of late Precambrian argillites under conditions of relatively low oxygen fugacity, along a line that roughly coincides with the westward disappearance of continental basement. In the other textural type of two-mica granite the micas are equigranular and there is a greater variety of accessory minerals. The magmatic evolution of this type also appears to reflect the influence of late Precambrian argillites; there may be age differences between the two types of two-mica granites.

Journal of Geophysical Research Solid Earth↗
USGS2012

Paleomagnetic evidence for the timing of collapse and resurgence of the Lake City Caldera, San Juan Mountains, Colorado

Rocks of the 23.1‐m.y.‐old Lake City caldera consist of the compositionally zoned Sunshine Peak Tuff, postcollapse intracaldera lava flows, and resurgent quartz syenite intrusions. Declinations of reversely magnetized ( I = −45° to −75°) Sunshine Peak Tuff change from easterly ( D = 93°–130°) throughout most of the tuff to southerly ( D = 195°–207°) within the late eruptive phases. The postcollapse lava flows are also reversely magnetized, and their emplacement appears to have spanned resurgence of the caldera. In contrast, the uppermost levels of the resurgent intrusion record two distinct components of magnetization: a reversed component isolated by thermal demagnetization typically at temperatures greater than 400°C and a normal component defined by thermal demagnetization at lower temperatures. A baked‐contact test implies that emplacement and initial cooling of the intrusion occurred during a time of reversed polarity and that continued cooling at lower blocking temperatures took place during normal polarity. Assuming that the reversed directions of the extrusive and intrusive rocks record a single period of reversed polarity as suggested by the paleomagnetic and radiometric age data, the time for caldera development (from eruption of the ash flow tuffs to emplacement of the resurgent intrusion) was less than about 300,000 years on the basis of best estimates of durations of reversed polarities in the early Miocene. The estimated 300,000‐year time span for the development of the Lake City caldera is very close in duration to the radiometrically determined time spans of much younger (Pleistocene) calderas having similar dimensions and similar volumes of erupted material.

Colorado↗
USGS2012

Pillow basalts of the Angayucham terrane: Oceanic plateau and island crust accreted to the Brooks Range

The Angayucham Mountains (north margin of the Yukon-Koyukuk province) are made up of an imbricate stack of four to eight east-west trending, steeply dipping, fault slabs composed of Paleozoic (Devonian to Mississippean), Middle to Late Triassic, and Early Jurassic oceanic upper crustal rocks (pillow basalt, subordinate diabase, basaltic tuff, and radiolarian chert). Field relations and geochemical characteristics of the basaltic rocks suggest that the fault slabs were derived from an oceanic plateau or island setting and were emplaced onto the Brooks Range continental margin. The basalts are variably metamorphosed to prehnite-pumpellyite and low-greenschist facies. Major element analyses suggest that many are hypersthene-normative olivine tholeiites. Classification based on immobile trace elements confirms the tholeiitic character of most of the basalts but suggests that some had primary compositions transitional to alkali basalt. Although field and petrographic features of the basalts are similar, trace element characteristics allow definition of geographically distinct suites. A central outcrop belt along the crest of the mountains is made up of basalt with relatively flat rare earth element (REE) patterns. This belt is flanked to the north and south by LREE (light rare earth element)-enriched basalts. Radiolarian and conodont ages from interpillow and interlayered chert and limestone indicate that the central belt of basalts is Triassic in age, the southern belt is Jurassic in age, and the northern belt contains a mixture of Paleozoic and Mesozoic ages. Data for most of the basalts cluster in the “within-plate basalt” fields of trace element discriminant diagrams; none have trace-element characteristics of island arc basalt. The Triassic and Jurassic basalts are geochemically most akin to modern oceanic plateau and island basalts. Field evidence also favors an oceanic plateau or island setting. The great composite thickness of pillow basalt probably resulted from obduction faulting, but the lack of fault slabs of gabbro or peridotite suggests that obduction faults did not penetrate below oceanic layer 2, a likely occurrence if layer 2 were anomalously thick, as in the vicinity of an oceanic island. The presence of basaltic tuff interbeds indicates proximity to an explosive basaltic eruptive center. The juxtaposition of submarine basalts of differing chemical affinity and age, adjacent to higher-grade Paleozoic metamorphic rocks of the Brooks Range to the north, may be explained by obduction of internally complex (thickened) oceanic crust formed in an ocean plateau setting. Emplacement and rotation of thrust plates to steep attitudes occurred during accretion of the Brooks Range passive margin, probably beginning in the Late to Middle Jurassic.

Journal of Geophysical Research Solid Earth↗
USGS2002

Migration and habitats of diadromous Danube River sturgeons in Romania: 1998-2000

Upstream migrant adults of stellate sturgeon, Acipenser stellatus (10 in 1998, 43 in 1999) and Russian sturgeon, A. gueldenstaedtii (three in 1999) were captured at river km (rkm) 58-137, mostly in the spring, and tagged with acoustic tags offering a reward for return. The overharvest was revealed by tag returns (38% in 1998, 28% in 1999) and by harvest within 26 days (and before reaching spawning grounds) of the six stellate sturgeon tracked upstream. A drop-back of > 50% of the tagged sturgeon, some to the Black Sea, shows a high sensitivity to interruption of migration and capture/handling/holding. Harvesting and dropback prevented tracking of sturgeon to spawning sites. Gillnetting and tracking of stellate sturgeon showed that the autumn migration ended in early October (river temperature 16??C) and identified a likely wintering area at river km (rkm) 75-76 (St George Branch). Thus, fishery harvesting after early October captures wintering fish, not migrants. Rare shoreline cliffs in the lower river likely create the only rocky habitat for sturgeon spawning. A survey for potential spawning habitats found five sites with rocky substrate and moderate water velocity, all ???rkm 258. Drift netting caught early life-stages of 17 fish species and one sturgeon, a beluga, Huso huso, larva likely spawned at ???rkm 258. All diadromous Danube sturgeons likely spawn at ???rkm 258.

Journal of Applied Ichthyology↗
USGS1996

Numerical simulation of widening and bed deformation of straight sand-bed rivers. II: Model evaluation

In this paper the numerical model presented in the companion paper is tested and applied. Assessment of model accuracy was based on two approaches. First, predictions of evolution of a 13.5 km reach of the South Fork of the Forked Deer River, in west Tennessee, were compared to observations over a 24-yr period. Results suggest that although the model was able to qualitatively predict trends of widening and deepening, quantitative predictions were not reliable. Simulated widths and depths were within 15% of the corresponding observed values, but observed change in these parameters at the study sites were also close to these values. Simulated rates of depth adjustment were within 15% of observed rates, but observed rates of channel widening at the study sites were approximately three times those simulated by the model. In the second approach, the model was used to generate relationships between stable channel width and bank-full discharge. The model was able to successfully replicate the form of empirically derived regime-width equations. Simulations were used to demonstrate the model's ability to obtain more realistic predictions of bed evolution in widening channels.

Journal of Hydraulic Engineering↗
USGS2003

Bedload transport in alluvial channels

Hydraulic, sediment, land-use, and rock-erosivity data of 22 alluvial streams were used to evaluate conditions of bedload transport and the performance of selected bedload-transport equations. Transport categories of transport-limited (TL), partially transport-limited (PTL), and supply-limited (SL) were identified by a semiquantitative approach that considers hydraulic constraints on sediment movement and the processes that control sediment availability at the basin scale. Equations by Parker et al. in 1982, Schoklitsch in 1962, and Meyer-Peter and Muller in 1948 adequately predicted sediment transport in channels with TL condition, whereas the equations of Bagnold in 1980, and Schoklitsch, in 1962, performed well for PTL and SL conditions. Overall, the equation of Schoklitsch predicted well the measured bedload data for eight of 22 streams, and the Bagnold equation predicted the measured data in seven streams.

Journal of Hydraulic Engineering↗
USGS2008

Computing nonhydrostatic shallow-water flow over steep terrain

Flood and dambreak hazards are not limited to moderate terrain, yet most shallow-water models assume that flow occurs over gentle slopes. Shallow-water flow over rugged or steep terrain often generates significant nonhydrostatic pressures, violating the assumption of hydrostatic pressure made in most shallow-water codes. In this paper, we adapt a previously published nonhydrostatic granular flow model to simulate shallow-water flow, and we solve conservation equations using a finite volume approach and an Harten, Lax, Van Leer, and Einfeldt approximate Riemann solver that is modified for a sloping bed and transient wetting and drying conditions. To simulate bed friction, we use the law of the wall. We test the model by comparison with an analytical solution and with results of experiments in flumes that have steep (31°) or shallow (0.3°) slopes. The law of the wall provides an accurate prediction of the effect of bed roughness on mean flow velocity over two orders of magnitude of bed roughness. Our nonhydrostatic, law-of-the-wall flow simulation accurately reproduces flume measurements of front propagation speed, flow depth, and bed-shear stress for conditions of large bed roughness.

Journal of Hydraulic Engineering↗
USGS2002

Dynamic fuzzy modeling of storm water infiltration in urban fractured aquifers

In an urban fractured-rock aquifer in the Mt. Eden area of Auckland, New Zealand, disposal of storm water is via "soakholes" drilled directly into the top of the fractured basalt rock. The dynamic response of the groundwater level due to the storm water infiltration shows characteristics of a strongly time-varying system. A dynamic fuzzy modeling approach, which is based on multiple local models that are weighted using fuzzy membership functions, has been developed to identify and predict groundwater level fluctuations caused by storm water infiltration. The dynamic fuzzy model is initialized by the fuzzy clustering algorithm and optimized by the gradient-descent algorithm in order to effectively derive the multiple local models-each of which is associated with a locally valid model that represents the groundwater level state as a response to different intensities of rainfall events. The results have shown that even if the number of fuzzy local models derived is small, the fuzzy modeling approach developed provides good prediction results despite the highly time-varying nature of this urban fractured-rock aquifer system. Further, it allows interpretable representations of the dynamic behavior of the groundwater system due to storm water infiltration.

Journal of Hydrologic Engineering↗
USGS1966

Variational method of determining effective moduli of polycrystals with tetragonal symmetry

Variational principles have been applied to aggregates of randomly oriented pure‐phase polycrystals having tetragonal symmetry. The bounds of the effective elastic moduli obtained in this way show a substantial improvement over the bounds obtained by means of the Voigt and Reuss assumptions. The Hill average is found to be a good approximation in most cases when compared to the bounds found from the variational method. The new bounds reduce in their limits to the Voigt and Reuss values.

Journal of Applied Physics↗
USGS2009

Forecasting distributions of large federal-lands fires utilizing satellite and gridded weather information

The current study presents a statistical model for assessing the skill of fire danger indices and for forecasting the distribution of the expected numbers of large fires over a given region and for the upcoming week. The procedure permits development of daily maps that forecast, for the forthcoming week and within federal lands, percentiles of the distributions of (i) number of ignitions; (ii) number of fires above a given size; (iii) conditional probabilities of fires greater than a specified size, given ignition. As an illustration, we used the methods to study the skill of the Fire Potential Index an index that incorporates satellite and surface observations to map fire potential at a national scale in forecasting distributions of large fires. ?? 2009 IAWF.

International Journal of Wildland Fire↗
USGS2011

An evaluation of the Bayesian approach to fitting the N-mixture model for use with pseudo-replicated count data

The N-mixture model proposed by Royle in 2004 may be used to approximate the abundance and detection probability of animal species in a given region. In 2006, Royle and Dorazio discussed the advantages of using a Bayesian approach in modelling animal abundance and occurrence using a hierarchical N-mixture model. N-mixture models assume replication on sampling sites, an assumption that may be violated when the site is not closed to changes in abundance during the survey period or when nominal replicates are defined spatially. In this paper, we studied the robustness of a Bayesian approach to fitting the N-mixture model for pseudo-replicated count data. Our simulation results showed that the Bayesian estimates for abundance and detection probability are slightly biased when the actual detection probability is small and are sensitive to the presence of extra variability within local sites.

Journal of Statistical Computation and Simulation↗
USGS2020

Influence of hydropower outflow characteristics affecting riverbank stability: The lower Osage River case (Missouri, USA)

This research examined the influences of outflow characteristics affecting riverbank stability. The 130 km stretch of the lower Osage River downstream from Bagnell Dam (Missouri, USA) provided an excellent case study for this purpose. The integrated BSTEM model with the HEC-RAS model was accurately calibrated and validated with data from the U.S. Geological Survey (USGS). Then, the outflow characteristics (peak flow duration, flow drawdown rate, and low flow duration) were investigated individually. The results of this study showed that: 1) Riverbank stability is little affected by the duration time of the peak flow, especially on the reaches far from the dam. 2) Sudden flow drawdown significantly reduces riverbank stability. However, the impact of the drawdown rate decreases with distance from the dam. 3) The duration of the low flow after peak flow influences the riverbank stability value proportional to the distance from the dam. The time of low flow before failure increases as the distance from the dam increases.

Missouri↗
USGS2024

New insights on the origin of the Richardson-Richards equation

The Richardson-Richards equation (RRE), despite known shortcomings especially in regard to preferential flow, provides the basis of the vast majority of unsaturated flow models in use today. L.F. Richardson published this equation in 1922, nine years before L.A. Richards. Whereas Richards approached this problem directly from the groundbreaking developments of Edgar Buckingham, Richardson, surprisingly, cited as his starting point only the earlier work of L.J. Briggs. Collectively, these four scientists’ published and unpublished work reveals that: (1) Briggs’ work, though qualitative, captured the essential physical principles needed for quantifying unsaturated flow; (2) Buckingham came very close to deriving the RRE and explained why he stopped short of doing so; (3) derivation of the RRE from the work of either Briggs or Buckingham required only modest developmental work; and (4) besides deriving the RRE, Richards carried through much of the experimental agenda that Buckingham considered a necessary precursor to mathematical treatment.

Hydrological Sciences Journal↗