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,153 records · Page 64Linked to original sources

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↗

Chemical, multispectral, and textural constraints on the composition and origin of rocks at the Mars Pathfinder landing site

Rocks at the Mars Pathfinder site are probably locally derived. Textures on rock surfaces may indicate volcanic, sedimentary, or impact-generated rocks, but aeolian abration and dust coatings prevent unambiguous interpretation. Multispectral imaging has resolved four spectral classes of rocks: gray and red, which occur on different surfaces of the same rocks; pink, which is probably soil crusts; and maroon, which occurs as large boulders, mostly in the far field. Rocks are assigned to two spectral trends based on the position of peak reflectance: the primary spectral trend contains gray, red, and pink rocks; maroon rocks constitute the secondary spectral trend. The spatial pattern of spectral variations observed is oriented along the prevailing wind direction. The primary spectral trend arises from thin ferric coatings of aeolian dust on darker rocks. The secondary spectral trend is apparently due to coating by a different mineral, probably maghemite or ferrihydrite. A chronology based on rock spectra suggests that rounded maroon boulders constitute the oldest petrologic unit (a flood deposit), succeeded by smaller cobbles possibly deposited by impact, and followed by aeolian erosion and deposition. Nearly linear chemical trends in alpha proton X-ray spectrometer rock compositions are interpreted as mixing lines between rock and adhering dust, a conclusion supported by a correlation between sulfur abundance and red/blue spectral ratio. Extrapolations of regression lines to zero sulfur give the composition of a presumed igneous rock. The chemistry and normative mineralogy of the sulfur-free rock resemble common terrestrial volcanic rocks, and its classification corresponds to andesite. Igneous rocks of this composition may occur with clastic sedimentary rocks or impact melts and breccias. However, the spectral mottling expected on conglomerates or breccias is not observed in any APXS-analyzed rocks. Interpretation of the rocks as andesites is complicated by absence of a "1 μm" pyroxene absorption band. Plausible explanations include impact glass, band masking by magnetite, or presence of calcium- and iron-rich pyroxenes and olivine which push the absorption band minimum past the imager's spectral range. The inferred andesitic composition is most similar to terrestrial anorogenic icelandites, formed by fractionation of tholeiitic basaltic magmas. Early melting of a relatively primitive Martian mantle could produce an appropriate parent magma, supporting the ancient age of Pathfinder rocks inferred from their incorporation in Hesperian flood deposits. Although rocks of andesitic composition at the Pathfinder site may represent samples of ancient Martian crust, inferences drawn about a necessary role for water or plate tectonics in their petrogenesis are probably unwarranted.

Journal of Geophysical Research E: Planets↗

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↗

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↗

Origin of lead from green glass of Apollo 15426: A search for primitive lunar lead.

A major obstacle in lunar U-Pb chronology is the elusiveness of the primary Pb isotopic composition and U/Pb ratio and therefore the Pb evolution for the early history of the moon. In an attempt to seek the primitive lunar Pb isotopic composition, green glass from lunar clod 15426,49 was studied for U-Th-Pb systematics because it is extremely Mg-rich and known to be the most primitive among sampled lunar volcanic rocks. Green glass and yellow impact glass were handpicked, washed with water, then washed in an acid mixture (0.3 N HBr + 0.2 N HF), and the washes and residues were analyzed. In contrast to the Pb from Apollo 17 orange glass, Pb from the surface of the green glass spherules is totally different from that of the spherule interior, but is similar to lunar surface lead. Data for the interior Pb yield a 207 Pb/ 206 Pb age of 3.41 Ga. Because of the low Pb concentration and high U/Pb ratio observed for the interior, the initial Pb was poorly defined. Nevertheless, the data indicate that lunar Pb evolved in an environment with 238 U/ 204 Pb = 19–55, which is considerably lower than those for mare basalts (around 300) but higher than values for the earth (6–8).

Journal of Geophysical Research Solid Earth↗

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↗

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↗

One-dimensional thermal modelling of Acadian metamorphism in southern Vermont, USA

One-dimensional thermal (1DT) modelling of an Acadian (Devonian) tectonothermal regime in southern Vermont, USA, used measured metamorphic pressures and temperatures and estimated metamorphic cooling ages based on published thermobarometric and geochronological studies to constrain thermal and tectonic input parameters. The area modelled lies within the Vermont Sequence of the Acadian orogen and includes: (i) a western domain containing garnet-grade pre-Silurian metasedimentary and metavolcanic rocks from the eastern flank of an Acadian composite dome structure (Rayponda-Sadawga Dome); and (ii) an eastern domain containing similar, but staurolite- or kyanite-grade, rocks from the western flank of a second dome structure (Athens Dome), approximately 10 km farther east. Using reasonable input parameters based on regional geological, petrological and geochronological constraints, the thermal modelling produced plausible P-T paths, and temperature-time (T-t) and pressure-time (P-t) curves. Information extracted from P-T-t modelling includes values of maximum temperature and pressure on the P-T paths, pressure at maximum temperature, predicted Ar closure ages for hornblende, muscovite and K-feldspar, and integrated exhumation and cooling rates for segments of the cooling history. The results from thermal modelling are consistent with independently obtained pressure, temperature and Ar cooling age data on regional metamorphism in southern Vermont. Modelling results provide some important bounding limits on the physical conditions during regional metamorphism, and indicate that the pressure contemporaneous with the attainment of peak temperature was probably as much as 2.5 kbar lower than the actual maximum pressure experienced by rocks along various particle paths. In addition, differences in peak metamorphic grade (garnet-grade versus staurolite-grade or kyanite-grade) and peak temperature for rocks initially loaded to similar crustal depths, differences in calculated exhumation rates, and differences in 40Ar/39Ar closure ages are likely to have been consequences of variations in the duration of isobaric heating (or 'crustal residence periods') and tectonic unroofing rates. Modelling results are consistent with a regional structural model that suggests west to east younging of specific Acadian deformational events, and therefore diachroneity of attainment of peak metamorphic conditions and subsequent 40Ar/39Ar closure during cooling. Modelling is consistent with the proposition that regional variations in timing and peak conditions of metamorphism are the result of the variable depths to which rocks were loaded by an eastward-thickening thrust-nappe pile rooted to the east (New Hampshire Sequence), as well as by diachronous structural processes within the lower plate rocks of the Vermont Sequence.

Journal of Metamorphic Geology↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Effects of uranium mining discharges on water quality in the Puerco River basin, Arizona and New Mexico

From 1967 until 1986, uranium mine dewatering increased dissolved gross alpha, gross beta, uranium and radium activities and dissolved selenium and molybdenum concentrations in the Puerco River as indicated by time trends, areal patterns involving distance from the mines and stream discharge. Additionally, increased dissolved uranium concentrations were identified in groundwater under the Puerco River from where mine discharges entered the river to approximately the Arizona-New Mexico State line about 65 km downstream. Total mass of uranium and gross alpha activity released to the Puerco River by mine dewatering were estimated as 560 Mg (560 × 10 6 g) and 260 Ci, respectively. In comparison, a uranium mill tailings pond spill on 16 July 1979, released an estimated 1.5 Mg of uranium and 46 Ci of gross alpha activity. Mass balance calculations for alluvial ground water indicate that most of the uranium released did not remain in solution. Sorption of uranium on sediments and uptake of uranium by plants probably removed the uranium from solution.

Arizona, New Mexico↗