USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geological Research”

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

Geodetic measurement of crustal deformation on the San Andreas, Hayward, and Calaveras faults near San Francisco, California

Analysis of a geodetic network of 115 lines crossing the San Andreas, Hayward, and Calaveras faults in the vicinity of San Francisco Bay and measured repeatedly between 1970 and 1980 has revealed details about the accommodation of relative plate motion in this area. The most striking result is that the deformation is not uniformly distributed across the area. In the east bay, along the Hayward and Calaveras faults, all motion appears to take place as slip directly on the fault, with no accumulation of strain in the adjacent crust. On both the Calaveras and the Hayward faults the rate obtained for the 1970 to 1980 period agrees with geologic rates spanning a few million years and with creep rates spanning a few decades. The Hayward fault slip rate is 7 ± 1 mm/yr. The Calaveras fault slip rate is 7 ± 1 mm/yr, with perhaps half of this slip distributed across a zone a few kilometers wide, probably as inelastic deformation of weak near‐surface material. The absence of strain accumulation in the east bay is surprising since the Hayward and Calaveras faults have been the site of large earthquakes in the past. A block located east of the Calaveras fault and south of the Las Positas fault has been rotating clockwise at a rate of 0.3 ± 0.1 μrad/yr with very little internal deformation. Along the San Francisco peninsula no detectable slip occurs (less than 1.5 mm/yr) at the surface, but appreciable strain is accumulating. Near fault shear strain rates are 0.6 ± 0.1 μstrain/yr (engineering) with direction N47°W ± 9. The slip rate near the San Andreas fault is 12.2 ± 3.9 mm/yr distributed across a broad zone. The relative motion across the whole region during the period 1970–1980 is 32.1 ± 7.4 mm/yr.

California↗

U.S. Geological Survey Trace Elements and related reports through 1953

This report combines and brings up-to-date the information previously given in Trace Elements Investigations Report 325, "Numerical list of U.S. Geological Survey Trace Elements Reports to April 30, 1953," and Trace Elements Investigations Report 301, "Topical index and bibliography of U.S. Geological Survey Trace Elements and related reports." Part I is a numerical list of U.S. Geological Survey Trace Elements Investigations and Memorandum reports. It supersedes TEI-325. This part lists not only reports (followed by a date) that have been transmitted to the U.S. Atomic Energy Commission, but also reports in preparation (followed by an asterisk) for which tentative titles were available on December 31, 1953. Reports that have been published are indicated by the abbreviation of the medium of publication. (See also part II.) Part II is a reference guide to Trace Elements and related reports that are available to the public; this part supersedes Part 2 of the TEI-301 (published as Geological survey Circular 281). These reports are grouped according to the type of publication or release. Abstracts published in Nuclear Science Abstracts are not included in Part II, although certain TEI and TEM reports, the abstracts of which have been published in NSA, are so indicated in Part I. Publications in process on December 31, 1953, are designated by an asterisk. Part III is a finding list of states, areas, and subjects. It is based on information derived mostly from the titles of reports and, where titles are of a general nature, from a cursory review of the reports. This list is not a complete index of the information given in Trace Elements and related reports, but is designed to find subjects of major interest. Because of the numerous entries for Colorado and Utah, information has been listed by counties and, where possible, by subject under these states. Other states have county listings only if a county is included in the title of a report; otherwise, areas may be listed separately under the state. Major subjects are listed separately in the index and also where appropriate under states. Analytical methods and subjects related to analytical research are listed under Analytical Methods and Research, but not separately throughout the index. Most mineralogic studies are included under the heading Mineralogy, but are not necessarily listed according to location. Part IV is a finding list of authors. The words “with” and “and” are used to indicate seniority of authorship. For example, a listing of Jones and Brown indicates that Jones is the senior author. A listing of Jones with Brown indicates that Brown is the senior author. In both parts III and IV all Trace Elements reports are listed, as well as other related reports that have not been issued as Trace Elements reports. The following standard abbreviations have been used: TEI, Trace Elements Investigations report; TEM, Trace Elements Memorandum report; P, Professional Paper; B, Bulletin; C, Circular; J, Journal; OF, open file; TIS, Technical Information Service release; NSA, Nuclear Science Abstracts; QM, Quadrangle Map Series; and OM, Oil and Gas map or Mineral Investigations map or report.

Trace Elements Investigations↗

Geologic mapping of Europa

Galileo data enable the major geological units, structures, and surface features to be identified on Europa. These include five primary units (plains, chaos, band, ridge, and crater materials) and their subunits, along with various tectonic structures such as faults. Plains units are the most widespread. Ridged plains material spans a wide range of geological ages, including the oldest recognizable features on Europa, and appears to represent a style of tectonic resurfacing, rather than cryovolcanism. Smooth plains material typically embays other terrains and units, possibly as a type of fluid emplacement, and is among the youngest material units observed. At global scales, plains are typically mapped as undifferentiated plains material, although in some areas differences can be discerned in the near infrared which might be related to differences in ice grain size. Chaos material is composed of plains and other preexisting materials that have been severely disrupted by inferred internal activity; chaos is characterized by blocks of icy material set in a hummocky matrix. Band material is arrayed in linear, curvilinear, wedge-shaped, or cuspate zones with contrasting albedo and surface textures with respect to the surrounding terrain. Bilateral symmetry observed in some bands and the relationships with the surrounding units suggest that band material forms by the lithosphere fracturing, spreading apart, and infilling with material derived from the subsurface. Ridge material is mapped as a unit on local and some regional maps but shown with symbols at global scales. Ridge material includes single ridges, doublet ridges, and ridge complexes. Ridge materials are considered to represent tectonic processes, possibly accompanied by the extrusion or intrusion of subsurface materials, such as diapirs. The tectonic processes might be related to tidal flexing of the icy lithosphere on diurnal or longer timescales. Crater materials include various interior (smooth central, rough inner, and annular massif) and exterior (continuous ejecta) subunits. Structural features and landforms are shown with conventional symbols. Type localities for the units are identified, along with suggestions for portraying the features on geological maps, including colors and letter abbreviations for material units. Implementing these suggestions by the planetary mapping community would facilitate comparisons of maps for different parts of Europa and contribute to an eventual global synthesis of its complex geology. On the basis of initial mapping results, a stratigraphic sequence is suggested in which ridged plains form the oldest unit on Europa, followed by development of band material and individual ridges. Band materials tend to be somewhat older than ridges, but in many areas the two units formed simultaneously. Similarly, the formation of most chaos follows the development of ridged plains; although chaos is among the youngest materials on Europa, some chaos units might have formed contemporaneously with ridged plains. Smooth plains generally embay all other units and are late-stage in the evolution of the surface. C1 craters are superposed on ridged plains but are crosscut by other materials, including bands and ridges. Most c2 craters postdate all other units, but a few c2 craters are cut by ridge material. C3 craters constitute the youngest recognizable material on Europa. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Geology and origin of Europa's "Mitten" feature (Murias Chaos)

The “Mitten” (provisionally named Murias Chaos by the International Astronomical Union) is a region of elevated chaos-like terrain in the leading hemisphere of Europa. Its origin had been explained under the currently debated theories of melting through a thin lithosphere or convection within a thick one. Galileo observations reveal several characteristics that suggest that the Mitten is distinct from typical chaos terrain and point to a different formational process. Photoclinometric elevation estimates suggest that the Mitten is slightly elevated with respect to the surrounding terrain; geologic relations indicate that it must have raised significantly from the plains in its past, resembling disrupted domes on Europa's trailing hemisphere. Moreover, the Mitten material appears to have extruded onto the plains and flowed for tens of kilometers. The area subsequently subsided as a result of isostatic adjustment, viscous relaxation, and/or plains loading. Using plate flexure models, we estimated the elastic lithosphere in the area to be several kilometers thick. We propose that the Mitten originated by the ascent and extrusion of a large thermal diapir. Thermal-mechanical modeling shows that a Mitten-sized plume would remain sufficiently warm and buoyant to pierce through the crust and flow unconfined on the surface. Such a diapir probably had an initial radius between 5 and 8 km and an initial depth of 20–40 km, consistent with a thick-lithosphere model. In this scenario the Mitten appears to represent the surface expression of the rare ascent of a large diapir, in contrast to lenticulae and chaos terrain, which may form by isolated and clustered small diapirs, respectively.

Journal of Geophysical Research E: Planets↗

Spatial distribution and diet of Lake Michigan juvenile lake trout

Most studies of Lake Michigan lake trout ( Salvelinus namaycush ) have focused on adults, with scant attention to juveniles (<400 mm). We explored the spatial distribution and diet of juvenile lake trout using U.S. Geological Survey September bottom trawl data (2015–2022) and stomach content information opportunistically collected since 2012 by various agencies using multiple gear types. Most juvenile lake trout in the September bottom trawl survey were caught at 37–64 m bottom depths. Length frequency data from the bottom trawl survey identified three size classes likely associated with wild juvenile lake trout age: < 85 mm (∼age-0), 85–170 mm (∼age-1) and > 170 mm (∼age-2+). Largest catches of wild lake trout < 170 mm occurred along a northeastern transect (near Frankfort, Michigan), whereas most > 170 mm were collected along southern transects. Mysis diluviana was the dominant prey for juvenile lake trout < 170 mm, and > 250 mm were primarily piscivorous, while 170–250 mm appeared to be a transitional period of switching from Mysis to fish. Species composition of prey fishes consumed by lake trout varied spatially and we found evidence of seasonal and annual diet variation within Grand Traverse Bay. Diporeia , once an important component of juvenile lake trout diet, appears to no longer be consumed by juvenile lake trout in Lake Michigan to any measurable degree. Continued research on the ecology of juvenile lake trout may provide insight into the effects of a changing ecosystem on juvenile lake trout diet and growth, thereby contributing to the effort to rehabilitate the Lake Michigan lake trout population.

Lake Michigan↗

A Web-Based Decision Support System for Assessing Regional Water-Quality Conditions and Management Actions

The U.S. Geological Survey National Water Quality Assessment Program has completed a number of water-quality prediction models for nitrogen and phosphorus for the conterminous United States as well as for regional areas of the nation. In addition to estimating water-quality conditions at unmonitored streams, the calibrated SPAtially Referenced Regressions On Watershed attributes (SPARROW) models can be used to produce estimates of yield, flow-weighted concentration, or load of constituents in water under various land-use condition, change, or resource management scenarios. A web-based decision support infrastructure has been developed to provide access to SPARROW simulation results on stream water-quality conditions and to offer sophisticated scenario testing capabilities for research and water-quality planning via a graphical user interface with familiar controls. The SPARROW decision support system (DSS) is delivered through a web browser over an Internet connection, making it widely accessible to the public in a format that allows users to easily display water-quality conditions and to describe, test, and share modeled scenarios of future conditions. SPARROW models currently supported by the DSS are based on the modified digital versions of the 1:500,000-scale River Reach File (RF1) and 1:100,000-scale National Hydrography Dataset (medium-resolution, NHDPlus) stream networks. ?? 2011 American Water Resources Association. This article is a U.S. Government work and is in the public domain in the USA.

Journal of the American Water Resources Associatio↗

Joint 3-D tomographic imaging of Vp, Vs and Vp/Vs and hypocenter relocation at Sinabung volcano, Indonesia from November to December 2013

We conducted travel time tomography using P- and S-wave arrival times of volcanic-tectonic (VT) events that occurred between November and December 2013 to determine the three-dimensional (3D) seismic velocity structure (Vp, Vs, and Vp/Vs) beneath Sinabung volcano, Indonesia in order to delineate geological subsurface structure and to enhance our understanding of the volcanism itself. This was a time period when phreatic explosions became phreatomagmatic and then magma migrated to the surface forming a summit lava dome. We used 4846 VT events with 16,138 P- and 16,138 S-wave arrival time phases recorded by 6 stations for the tomographic inversion. The relocated VTs collapse into three clusters at depths from the surface to sea level, from 2 to 4 km below sea level, and from 5 to 8.5 km below sea level. The tomographic inversion results show three prominent regions of high Vp/Vs (~ 1.8) beneath Sinabung volcano at depths consistent with the relocated earthquake clusters. We interpret these anomalies as intrusives associated with previous eruptions and possibly surrounding the magma conduit, which we cannot resolve with this study. One anomalous region might contain partial melt, at sea level and below the eventual eruption site at the summit. Our results are important for the interpretation of a conceptual model of the “plumbing system” of this hazardous volcano.

Sinabung Volcano↗

Active faulting and deformation of the Coalinga anticline as interpreted from three-dimensional velocity structure and seismicity

This work gives a clear picture of the geometry of aftershock seismicity in a large thrust earthquake. Interpretation of hypocenters and fault plane solutions, from the 1983 Coalinga, Coast Range California, earthquake sequence, in combination with the three-dimensional velocity structure shows that the active faulting beneath the fold primarily consists of a set of southwest dipping thrusts uplifting blocks of higher-velocity material. Above the main listric blind thrust there is a conjugate fault, steeply northeast dipping, that provides the western limit of the aftershocks within the Coalinga Anticline and that corresponds in location and spatial extent with the adjacent Pleasant Valley syncline. The character of the seismicity varies with the degree of previous deformation on each section of the anticline. Where the previous uplift was largest, the shallow seismicity shows secondary faulting on either side of the fold with orientations that correspond to the preexisting geologic structure. Diffuse seismicity characterizes the area with the least previous deformation. The mainshock rupture terminated where the fold trend was no longer uniform but had competing north and west trending features. The upward extent of the mainshock rupture ended at the approximate boundary between Franciscan and Great Valley Sequence rocks. Above that depth the main thrust appears to splay into a steeper segment and a near-horizontal segment. Thus the extent of rupture area is limited by the area of uniform structural orientation and by the variation in the type of material. With the three-dimensional velocity model each individual hypocenter moved slightly (0–2 km) in accord with the details of the surrounding velocity structure, so that secondary features in the seismicity pattern are more detailed than with a local one-dimensional model and station corrections. The overall character of the fault plane solutions was not altered by the three-dimensional model, but the more accurate ray paths did result in distinct changes. In particular, the mainshock has a fault plane dipping 30° southwest instead of the 23° obtained with the one-dimensional model.

Journal of Geophysical Research Solid Earth↗

Coastal eolian sand-ramp development related to paleo-sea-level changes during the Latest Pleistocene and Holocene (21–0 ka) in San Miguel Island, California, U.S.A.

Coastal eolian sand ramps (5–130 m elevation) on the northern slope (windward) side of the small San Miguel Island (13 km in W-E length) range in age from late Pleistocene to modern time, though a major hiatus in sand-ramp growth occurred during the early Holocene marine transgression (16–9 ka). The Holocene sand ramps (1–5 m measured thicknesses) currently lack large dune forms, thereby representing deflated erosional remnants, locally covering thicker late Pleistocene sand-ramp deposits. The ramp sand was initially supplied from the adjacent island-shelf platform, extending about 20 km north of the present coastline. The sand-ramp deposits and interbedded loess soils were 14 C dated using 112 samples from 32 archaeological sites and other geologic sections. Latest Pleistocene sand ramps (66–18 ka) were derived from across-shelf eolian sand transport during marine low stands. Shoreward wave transport supplied remobilized late Pleistocene sand from the inner shelf to Holocene beaches, where dominant NW winds supplied sand to the sand ramps. The onset dates of the sand-ramp deposition in San Miguel are 7.2 ± 1.5 ka (sample n = 14). The internal strata dates in the vertically accreting sand ramps are 3.4 ± 1.7 ka ( n = 34). The sand ramps in San Miguel show wide-scale termination of sand supply in the latest Holocene time. The sand-ramp top dates or burial dates are 1.7 ± 0.9 ka ( n = 28). The latest Holocene sand ramps are truncated along most of the island's northern coastline, indicating recent losses of nearshore sand reserves to onshore, alongshore, and, possibly, offshore sand sinks. The truncated sand ramps in San Miguel Island and in other sand-depleted marine coastlines provide warnings about future beach erosion and/or shoreline retreat from accelerated sea-level rise accompanying predicted global warming.

California↗

ShakeAlert® and schools: Incorporating earthquake early warning in school districts in Alaska, California, Oregon, and Washington

The U.S. Geological Survey-managed ShakeAlert® earthquake early warning system is the first public alerting system in the United States to provide rapid mass notification when an earthquake is detected. Although public alert delivery via mobile phones began in California in 2019 followed by Oregon and Washington in 2021, little is known about what might drive widespread implementation in at-risk institutional settings such as schools. For example, there has been limited research on how to best integrate earthquake early warning into existing emergency plans, alert systems, and drills to keep school children and personnel safe in an earthquake. To address this gap, in the spring of 2022, every school district superintendent in Alaska, California, Oregon, and Washington was sent a 15-min online survey. The survey assessed superintendent knowledge of ShakeAlert, preferences for alert messaging, and perceived opportunities and barriers to incorporating the system in schools. The results showed that superintendents had low awareness of ShakeAlert but held positive perceptions of the system's potential to enable life-saving protective actions. A major barrier to adoption included the perceived financial cost of implementing and maintaining the system. There were some statistically significant differences in state responses, and future research could investigate the specific needs of each state based on school district size and composition, hazard exposure, and earthquake experience. Together these findings can help inform targeted strategies to increase ShakeAlert adoption in schools and ultimately improve the safety of school children and staff during earthquakes.

Alaska, California, Oregon, Washington↗

A review of the geologic framework of the Long Island Sound Basin, with some observations relating to postglacial sedimentation

Most of the papers in this thematic section present regional perspectives that build on more than 100 years of geologic investigation in Long Island Sound. When viewed collectively, a common theme emerges in these works. The major geologic components of the Long Island Sound basin (bedrock, buried coastal-plain strata, recessional moraines, glacial-lake deposits, and the remains of a large marine delta) interact with the water body to affect the way the modern sedimentary system functions. Previous work, along with our present knowledge of the geologic framework of the Long Island Sound basin, is comprehensively reviewed with this theme in mind. Aspects of the crystalline bedrock, and the deltaic deposits associated with glacial Lake Connecticut, are examined with respect to their influence on sedimentation along the Connecticut coast and in the northern and western Sound. We also discuss the influence of the glacial drift that mantles the coastal-plain remnant along the north shore of Long Island and in the southern Sound. A total of approximately 22.7 billion m 3 of marine sediment has accumulated in the Long Island Sound basin. A significant portion (44%) of the fine-grained marine section in the central and western basins was redistributed there from the eastern Sound, as tidal scour removed slightly over 5 billion m 3 (5.3 X 10 12 kg) of fine material from glacial lake and early-marine deposits east of the Connecticut River. The remainder of the estimated 1.2 X 10 13 kg of fine-grained marine sediment that now resides in the central and western Sound can be accounted for by riverine input over the past 13.5 ka.

Connecticut, New York↗

Great expectations for earthquake early warnings on the United States West Coast

In October 2019, California became the first state in the United States to fully activate a public earthquake early warning system—ShakeAlert®—managed by the U.S. Geological Survey. The system was subsequently rolled out in March 2021 in Oregon and May 2021 in Washington. Earthquake early warning (EEW) systems can provide seconds of notice to people and technological systems that shaking is imminent, but their effectiveness depends on recipients’ expectations and actions as well as technical performance. To better understand these dependencies, we surveyed representative samples of adults in California (N = 1219), Oregon (N = 1020), and Washington (N = 1037) in February 2021. Most respondents had experienced earthquakes, but few had lived through violent shaking; most had not followed protective action guidance to Drop, Cover, and Hold On (DCHO) in earthquakes; and most reported no personal or social harm from prior earthquakes. Nevertheless, expectations and perceived usefulness of EEW were high, and higher still for those who expected alerts to be accurate and easy to use, expressed tolerance of missed and erroneous warnings, and expected to be affected by a damaging earthquake in their lifetime. Results suggest opportunities to better align public preferences and expectations with ShakeAlert operations. For example, some respondents preferred lower alerting thresholds than those proposed by government and scientists. Moreover, reported tolerance of warning errors was widespread, but respondents wanted explanations quickly, suggesting a need to further develop post-alert messaging. Findings from this study should be informative for future research on the co-evolution of experiences and expectations with EEW systems.

California, Oregon, Washington↗

Novel insights about petroleum systems from source and reservoir rock characterization, Cretaceous Deposits, Babouri-Figuil Basin, Northern Cameroon

The Babouri-Figuil Basin (BFB) is a frontier basin for petroleum in Cameroon. It belongs to the series of Cretaceous rift basins of the West and Central Rift System (WCARS), the origin of which is related to the opening of the South Atlantic. Within the same rift system, commercial hydrocarbon accumulations have been discovered in Chad, Sudan, Niger and, more recently, in Nigeria (Gongola Basin). The study of the geology of the BFB just recently received considerable attention, mainly because of its presumed hydrocarbon potential. In the pursuit of researching possible petroleum systems in the BFB, the current study provides a first look into the characterization of source and reservoir rock and its integration into a 2D lithostratigraphic model. The study was solely based on outcrop samples. Black shale and massive claystone are good to excellent hydrocarbon source rocks [e.g., up to 38 wt% total organic carbon (TOC), up to 943 mg/g hydrogen index, up to 85 m thickness, up to 20–30 km lateral extension], with moderate to high values of extractable organic matter (e.g., >10,000 ppm). Calcareous claystone, on the other hand, are poor source rocks [e.g., <0.20 wt% TOC]. The samples are thermally immature, except for those located close to volcanic intrusion at Golombe that have reached the threshold for oil generation (Tmax >435 °C, production index >0.1). The petrographic analysis of sandstone revealed that they are fine-grained to coarse-grained, poorly to moderately sorted, texturally and compositionally immature to submature, subarkosic to arkosic arenites. The main diagenetic processes that affected sandstones are as follows: moderate to intense compaction characterized by the development of long, concavo-convex, and sutured contacts between grains; cementation through calcite, iron oxide, and quartz cements; alteration of mica and feldspar grains; partial to complete dissolution of feldspar, mica, amphibole grains, and calcite cement; and the replacement of feldspar and mica grains by clay minerals. Alteration and dissolution increase the porosity of sandstone through the creation of secondary pores. However, mechanical compaction through the development of a pseudomatrix and cementation as pore-filling materials have significantly reduced the quality of sandstone beds as conventional petroleum reservoirs. Hence, the best reservoir-quality sandstones in the basin are generally located in the upper portion of the basin in terms of its lithostratigraphic model. They are the cleanest sandstones with the smallest amount of cement and the lowest ductile grain content (pseudomatrix), with a thickness that varies from 3 m to 120 m and a lateral extension of 20 km. The lithostratigraphic model of the basin is characterized by an extensive lacustrine environment that provided a thick sequence of organic-rich formations; sand deposited as extensive reservoirs sandwiched between shale/claystone beds; the development of stratigraphic traps through lateral facies change; and the widespread deposition of lacustrine and floodplain claystone that provide regional seals. The similarities between the Babouri-Figuil Basin and proven petroleum systems in other WCARS rift basins suggest that the basin may host at least one petroleum system where actively generating source rocks are present.

Babouri-Figuil Basin↗

A surrogate regression approach for computing continuous loads for the tributary nutrient and sediment monitoring program on the Great Lakes

Water quality (WQ) in many Great Lake tributaries has been degraded (increased nutrient and sediment concentrations) due to changes in their watersheds, resulting in downstream eutrophication. As part of the Great Lakes Water Quality Agreement, specific goals were established for loading of specific constituents (e.g., phosphorus). In 2010, the Great Lakes Restoration Initiative was launched to identify problem areas, accelerate restoration efforts, and track their progress. In 2011, the U.S. Geological Survey established a monitoring program on 30 tributaries to the lakes, representing ~ 46% of the U.S. draining area and the spectrum of land uses. Discrete measurements of nutrients and suspended sediment, and continuous measurements of flow and WQ surrogates (turbidity, temperature, specific conductance, pH, and dissolved oxygen) are being collected in these tributaries to document their WQ and estimate continuous (5-min) loading. To estimate loadings, two regression models were developed for each constituent for each site: one using continuous flow and a seasonality factor; and one using flow, seasonality, and continuous surrogates. Variables included in the final models for each constituent were chosen from the explanatory variables that worked “best” for all sites. In computing loads, when continuous surrogate data were unavailable for short periods, loads were computed using the flow and seasonality models. Prediction intervals for all loads were calculated using results from both models. These results provide a better understanding of short-term variability and long-term changes in loading affecting the environmental health of the Great Lakes than traditional regression techniques that employ only flow and seasonality parameters.

Great Lake↗

Frequency and initiation of debris flows in Grand Canyon, Arizona

Debris flows from 740 tributaries transport sediment into the Colorado River in Grand Canyon, Arizona, creating rapids that control its longitudinal profile. Debris flows mostly occur when runoff triggers failures in colluvium by a process termed “the fire hose effect.” Debris flows originate from a limited number of geologic strata, almost exclusively shales or other clay‐rich, fine‐grained formations. Observations from 1984 through 2003 provide a 20 year record of all debris flows that reached the Colorado River in Grand Canyon, and repeat photography provides a 100 year record of debris flows from 147 tributaries. Observed frequencies are 5.1 events/year from 1984 to 2003, and historic frequencies are 5.0 events/year from 1890 to 1983. Logistic regression is used to model historic frequencies based on drainage basin parameters observed to control debris flow initiation and transport. From 5 to 7 of the 16 parameters evaluated are statistically significant, including drainage area, basin relief, and the height of and gradient below debris flow source areas, variables which reflect transport distance and potential energy. The aspect of the river channel, which at least partially reflects storm movement within the canyon, is also significant. Model results are used to calculate the probability of debris flow occurrence at the river over a century for all 740 tributaries. Owing to the variability of underlying geomorphic controls, the distribution of this probability is not uniform among tributaries of the Colorado River in Grand Canyon.

Arizona↗

Geophysical constraints on Washington convergent margin structure

Gravity and magnetic maps of western Washington reveal the lateral structure and fabric of the Washington Coast Range, Puget Basin, and southern Washington Cascade Range. The magnetic and gravity maps show large amplitude positive anomalies associated with the shallow but largely buried section of Washington Coast Range mafic rocks which are separated by negative anomalies over deep sedimentary basins. The positive anomalies indicate that the Coast Range mafic basement extends farther east than previously thought, at least as far east as the longitude of Seattle. Linear and steep gravity and magnetic gradients indicate many unmapped, often buried faults in the Washington Coast Range Province. Magnetic highs are also associated with mapped batholiths in the Cascade arc. Several magnetic highs observed east of the Coast Range rocks and west of these batholiths may be associated with buried Tertiary plutons or ophiolites. Two-dimensional gravity and magnetic modeling constrained with geological and other geophysical data indicate that the Coast Range Province rocks are about 1 km thick at the coast, thickening to as much as 30 km near their postulated eastern edge. A maximum boundary on the average density of the upper 15–20 km of the rocks that compose the Coast Range Province of 2920 kg/m 3 was established by the modeling, suggesting a composition largely of basalt and gabbro with little interbedded sediments. Under these rocks may be mantle or a subduction complex composed of dense mafic, ultramafic, and sedimentary rocks like that proposed to underlie Vancouver Island. Previous gravity models of the Washington margin include lower densities for the proposed subduction complex than for Vancouver Island, suggesting a lower component of mafic and ultramafic rocks than the rocks underlying Vancouver Island. However, my Washington model requires that the proposed subduction complex be more dense than the trench sediments and, therefore, that material denser than sediments be incorporated within it. The absence of continental mantle and the modeled wedge shape of the Coast Range Province upper crust suggest that erosion of the bottom of the overriding plate by subduction processes may have occurred.

Washington↗

Modeling tritium transport through a deep unsaturated zone in an arid environment

Understanding transport of tritium ( 3 H) in unsaturated zones is critical to evaluating options for waste isolation. Tritium typically is a large component of low-level radioactive waste (LLRW). Studies at the U.S. Geological Survey's Amargosa Desert Research Site (ADRS) in Nevada investigate 3 H transport from a closed LLRW facility. Two boreholes are 100 and 160 m from the nearest waste trench and extend to the water table at 110 m. Soil-water vapor samples from the deep boreholes show elevated levels of 3 H at all depths. The objectives of this study were to (i) test source thermal and gas-advection mechanisms driving 3 H transport and (ii) evaluate model sensitivity to these mechanisms and to selected physical and hydraulic properties including porosity, tortuosity, and anisotropy. A two-dimensional numerical model incorporated a non-isothermal, heterogeneous domain of the unsaturated zone and instantaneous isotopic equilibrium. The TOUGH2 code was used; however, it required modification to account for temperature dependence of both the Henry's law equilibrium constant and isotopic fractionation with respect to tritiated water. Increases in source temperature, pressure, and porosity enhanced 3 H migration, but failed to match measured 3 H distributions. All anisotropic simulations with a source pressure component resembled, in shape, the upper portion of the 3 H distribution of the nearest borehole. Isotopic equilibrium limited migration of 3 H, while effects of radioactive decay were negligible. A 500 Pa pressure increase above ambient pressure in conjunction with a high degree of anisotropy (1:100) was necessary for simulated 3 H transport to reach the nearest borehole.

Vadose Zone Journal↗

Transport of tritium contamination to the atmosphere in an arid environment

Soil–plant–atmosphere interactions strongly influence water movement in desert unsaturated zones, but little is known about how such interactions affect atmospheric release of subsurface water-borne contaminants. This 2-yr study, performed at the U.S. Geological Survey's Amargosa Desert Research Site in southern Nevada, quantified the magnitude and spatiotemporal variability of tritium (3H) transport from the shallow unsaturated zone to the atmosphere adjacent to a low-level radioactive waste (LLRW) facility. Tritium fluxes were calculated as the product of 3H concentrations in water vapor and respective evaporation and transpiration water-vapor fluxes. Quarterly measured 3H concentrations in soil water vapor and in leaf water of the dominant creosote-bush [ Larrea tridentat a (DC.) Coville] were spatially extrapolated and temporally interpolated to develop daily maps of contamination across the 0.76-km2 study area. Maximum plant and root-zone soil concentrations (4200 and 8700 Bq L−1, respectively) were measured 25 m from the LLRW facility boundary. Continuous evaporation was estimated using a Priestley–Taylor model and transpiration was computed as the difference between measured eddy-covariance evapotranspiration and estimated evaporation. The mean evaporation/transpiration ratio was 3:1. Tritium released from the study area ranged from 0.12 to 12 μg d−1 and totaled 1.5 mg (8.2 × 1010 Bq) over 2 yr. Tritium flux variability was driven spatially by proximity to 3H source areas and temporally by changes in 3H concentrations and in the partitioning between evaporation and transpiration. Evapotranspiration removed and limited penetration of precipitation beneath native vegetation and fostered upward movement and release of 3H from below the root zone.

Vadose Zone Journal↗