Geology ReportsSearch

SEARCH · Geology Reports

Results for “Journal of Interpretation 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 379 records · Page 21Linked to original sources

Stratigraphic architecture and fluvial interpretations of the Upper Cretaceous (Turonian?) Middendorf Formation, Chesterfield County, South Carolina, U.S.A.

The Upper Cretaceous (Turonian?) Middendorf Formation is a sand-rich stratigraphic unit of fluvial origin that forms a large aquifer in the U.S. Atlantic Coastal Plain. In Chesterfield County (South Carolina), which is the site of the type locality, the formation ranges in thickness from 66.5 to > 119.7 meters. The base of the formation is an unconformity above Paleozoic metasiltstone, and the upper contact is an unconformity above which lies sand of the Quaternary Pinehurst Formation. Outcrops display the following five facies assemblages: 1) sandstone to conglomeratic sandstone (fluvial bar and channel deposits), 2) beds of alternating laminae of sandstone and mudstone (fluvial overbank or floodplain deposits), 3) ≥ 1 m-thick beds of clay (swamp deposits, floodplain deposits, and/or sediment that accumulated in abandoned fluvial channels), 4) 0.2–0.5 m-thick planar to slightly undulatory beds of framework-supported sandstone with a mud matrix (traction-dominated current deposits at the top of fluvial bars, upper-flow-regime bedform deposits in subsidiary fluvial channels, or coarse-grained overbank deposits), and 5) sandstone to conglomeratic sandstone cemented by iron (interpreted as fluvial bar and channel deposits, with the iron cement being a diagenetic “groundwater ferricrete” that formed via the circulation of shallow groundwater and the oxidation of iron-bearing minerals). Kaolinite in various forms is pervasive throughout the formation and is interpreted as an early diagenetic phenomenon that formed by prolonged postdepositional weathering and flushing by meteoric water under a warm and humid paleoclimate. The fluvial system that formed the Middendorf Formation prograded into the area from the west or northwest from uplifted margins of Mesozoic rift basins and/or the Appalachian Mountains. This progradation was a response to a base-level fall and the sediment accumulated during base-level lowstand and subsequent early transgression. In Chesterfield County, the Middendorf Formation can be subdivided into three fining-upward sequences. Each sequence consists predominantly of medium to coarse sand with a greater abundance of gravel in the lower part of the sequence and a greater abundance of clay and silt beds in the upper part. Each sequence is interpreted as either a response to autogenic processes or a response to allogenic sea-level changes, specifically a higher-order (higher-frequency) progression from relative lowstand conditions to early transgression whereby coarse sand and gravel (e.g., fluvial bar and channel deposits) were preserved during initial lowstand conditions and a greater proportion of mud and finer-grained sand (floodplain deposits) were preserved during subsequent early transgression. The Middendorf Formation is correlative with several other kaolinite-rich fluvial sandstones in North America including the Raritan Formation in New Jersey, the Tuscaloosa Formation of the eastern Gulf of Mexico (Alabama, Mississippi, Louisiana), the Woodbine Formation of the central Gulf of Mexico (Texas), and the Frontier Formation of Wyoming. The accumulation and preservation of these formations occurred in response to a Turonian eustatic sea-level fall and subsequent transgression, and the early diagenetic kaolinite in these formations is attributed to similar warm and humid paleoclimate conditions.

South Carolina

Applications of remote sensing to structural interpretations in the southern Appalachians

Remote sensing is the technology of studying distant objects by measuring and recording energy from one or more segments of the electromagnetic spectrum. Imaging sensors which operate from medium- and high-altitude aircraft or from spacecraft can provide a synoptic view of large areas and of surface phenomena not evident in the field. Image-acquiring systems and instruments have been designed to partially automate data collection and to reduce the time devoted to analysis, information extraction, and detection of changes of surface phenomena. Among these phenomena are the surface distribution of heat, moisture, snow, water, vegetation, and cultural features. When coupled with ancillary data, including field surveys, sensor data provide useful information for the recognition and mapping of regional structure, jointing patterns, drainage patterns, fault and fracture traces, and rock types. The recognition of several major linear surface features, two of which proved to be traces of previously unrecognized faults (the Canebrake and the Coeburn faults) in the Appalachian Plateaus, demonstrates the pragmatic application of aircraft and spacecraft remote sensing to geological investigations in the Appalachians.

Virginia, West Virginia

Fluid flow in the resurgent dome of Long Valley Caldera: Implications from thermal data and deep electrical sounding

Temperatures of 100??C are measured at 3 km depth in a well located on the resurgent dome in the center of Long Valley Caldera, California, despite an assumed >800??C magma chamber at 6-8 km depth. Local downflow of cold meteoric water as a process for cooling the resurgent dome is ruled out by a Pecle??t-number analysis of temperature logs. These analyses reveal zones with fluid circulation at the upper and lower boundaries of the Bishop Tuff, and an upflow zone in the metasedimentary rocks. Vertical Darcy velocities range from 10 to 70 cm a-1. A 21-km-long geoelectrical profile across the caldera provides resistivity values to the order of 100 to >103 ??m down to a depth of 6 km, as well as variations of self-potential. Interpretation of the electrical data with respect to hydrothermal fluid movement confirms that there is no downflow beneath the resurgent dome. To explain the unexpectedly low temperatures in the resurgent dome, we challenge the common view that the caldera as a whole is a regime of high temperatures and the resurgent dome is a local cold anomaly. Instead, we suggest that the caldera was cooled to normal thermal conditions by vigorous hydrothermal activity in the past, and that a present-day hot water flow system is responsible for local hot anomalies, such as Hot Creek and the area of the Casa Diablo geothermal power plant. The source of hot water has been associated with recent shallow intrusions into the West Moat. The focus of planning for future power plants should be to locate this present-day flow system instead of relying on heat from the old magma chamber. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research

Aftershock triggering by complete Coulomb stress changes

We examine the correlation between seismicity rate change following the 1992, M 7.3, Landers, California, earthquake and characteristics of the complete Coulomb failure stress (CFS) changes (ΔCFS( t )) that this earthquake generated. At close distances the time-varying “dynamic” portion of the stress change depends on how the rupture develops temporally and spatially and arises from radiated seismic waves and from permanent coseismic fault displacement. The permanent “static” portion (ΔCFS) depends only on the final coseismic displacement. ΔCFS diminishes much more rapidly with distance than the transient, dynamic stress changes. A common interpretation of the strong correlation between ΔCFS and aftershocks is that load changes can advance or delay failure. Stress changes may also promote failure by physically altering properties of the fault or its environs. Because it is transient, ΔCFS( t ) can alter the failure rate only by the latter means. We calculate both ΔCFS and the maximum positive value of ΔCFS( t ) (peak ΔCFS( t )) using a reflectivity program. Input parameters are constrained by modeling Landers displacement seismograms. We quantify the correlation between maps of seismicity rate changes and maps of modeled ΔCFS and peak ΔCFS( t ) and find agreement for both models. However, rupture directivity, which does not affect ΔCFS, creates larger peak ΔCFS( t ) values northwest of the main shock. This asymmetry is also observed in seismicity rate changes but not in ΔCFS. This result implies that dynamic stress changes are as effective as static stress changes in triggering aftershocks and may trigger earthquakes long after the waves have passed.

Journal of Geophysical Research B: Solid Earth

Light-toned salty soils and co-existing Si-rich species discovered by the Mars Exploration Rover Spirit in Columbia Hills

Light-toned soils were exposed, through serendipitous excavations by Spirit Rover wheels, at eight locations in the Columbia Hills. Their occurrences were grouped into four types on the basis of geomorphic settings. At three major exposures, the light-toned soils are hydrous and sulfate-rich. The spatial distributions of distinct types of salty soils vary substantially: with centimeter-scaled heterogeneities at Paso Robles, Dead Sea, Shredded, and Champagne-Penny, a well-mixed nature for light-toned soils occurring near and at the summit of Husband Hill, and relatively homogeneous distributions in the two layers at the Tyrone site. Aeolian, fumarolic, and hydrothermal fluid processes are suggested to be responsible for the deposition, transportation, and accumulation of these light-toned soils. In addition, a change in Pancam spectra of Tyrone yellowish soils was observed after being exposed to current Martian surface conditions for 175 sols. This change is interpreted to be caused by the dehydration of ferric sulfates on the basis of laboratory simulations and suggests a relative humidity gradient beneath the surface. Si-rich nodules and soils were observed near the major exposures of S-rich soils. They possess a characteristic feature in Pancam visible near-infrared (Vis-NIR) spectra that may be diagnostic of hydrated species, and this spectral feature can be used to search for additional Si-rich species. The exposures of hydrated salty soils within various geomorphic settings imply the potential existence of hydrous minerals in similar settings over a much wider area. Hydrous sulfates represent one of the candidates that may contribute the high level of water equivalent hydrogen in equatorial regions detected by the Neutron Spectrometer on Mars Odyssey.

Journal of Geophysical Research

Gravity survey of the Mount Toondina impact structure, South Australia

Gravity and seismic reflection data, together with geologic mapping, indicate that the Mount Toondina feature in South Australia is best interpreted as an eroded 4-km-diameter impact structure consisting of a ring structural depression surrounding a pronounced central uplift. Beds at the center of the structure within the central uplift have been raised as much as 200 m from depth and deformed by convergent flow. Seismic reflection data indicate that deformation extends to depths of only ≈800 m; at greater depths the reflectors are nearly flat lying, indicating little or no deformation. Gravity data show residual anomalies of +1.0 mGal coincident with the central uplift and a −0.5 Mgal low associated with the ring structural depression. Modeling of the gravity data indicates that relatively high-density material occurs within the central uplift, whereas the ring depression is filled with low-density material. The deformation at Mount Toondina is typical of a complex impact crater; the 4-km diameter is consistent with the expected threshold size for complex craters formed in weak to moderate strength sedimentary rocks.

South Australia

The use of MAGSAT data to determine secular variation

A combined spatial and secular variation model of the geomagnetic field, labeled M061581, is derived from a selection of Magsat data. Secular variation (SV) data computed from linear fits to midnight hourly values from 19 magnetic observatories were also included in the analysis but were seen to have little effect on the model. The SV patterns from this new model are compared with those from the 1980 IGRF and with those for 1970 computed by the DGRF and with the 1960 patterns computed using the GSFC(12/66) model. Most of the features of the M061581 are identical in location and level with those of the 1980 IGRF. Together they confirm that the reversals in sign of field change seen over Asia and North America between 1965 and 1975 are reverting to the pre-1965 states. The M061581 model gives −32 nT/yr for the dipole decay rate, larger than the 70% increase already reported since 1965. This abnormally high value is interpreted as being a defect of the model because it appears to result from a much larger (−100 nT/yr) drop in field over the polar regions not indicated by the 1980 IGRF. This north polar decrease is shown to be of external origin as the result of a combination of the seasonal effect of the north polar ionospheric (counterclockwise) afternoon Sq p cell increasing in intensity from the beginning (November 1979) to the end (June 1980) of the Magsat data collection period, coupled with an enhancement of its effect as the orbit lowered from the 350– 550 km initial altitudes to near 200 km just prior to burnup. This experiment indicates that secular variation can be obtained from satellite data for intervals of less than a full year if corrections can be made for seasonal effects and that ‘annual snapshots’ of the field by a satellite would allow easy and accurate models of secular change without the use of any surface data.

Journal of Geophysical Research Solid Earth

Seismicity and detection/location threshold in the southern Great Basin seismic network

A spatially varying model of the detection/location capabilities of the Southern Great Basin seismic network (SGBSN) has been derived that is based on simple empirical relations and statistics. This permits use of almost all the catalog data gathered; instead of ignoring data that are below the threshold of completeness, a spatially varying threshold model is developed so that subregions having lower completeness levels than the network as a whole can be outlined and the completeness level of each sub-region determined. Such a model is required to unambiguously identify regions that are aseismic due to natural processes rather than to limited detection and/or location capabilities. Accounting for spatial variations in detection/location threshold is also important for studies in which magnitude-frequency distributions are interpreted in terms of source scaling properties. The characteristics of the spatial distribution of earthquakes, where earthquake clusters and aseismic regions locate, appear to be stable at all magnitude levels so that inferences about where strain is being accommodated will be the same whether numbers of earthquakes or strain estimated from seismic moments are examined. For the southern Great Basin region these principal characteristics include clusters at the northern end of the Furnace Creek fault and in the Pahranagat Shear Zone, and a relatively large number of earthquakes in the northern and southeastern portions of the Nevada Test Site. These clusters cover regions much larger than the surface projections of any of the mapped faults. The extent to which seismicity is induced by nuclear testing is unclear. The predominantly aseismic regions include the area west of the Death Valley/Furnace Creek fault system and an almost complete absence of events at Yucca Mountain. Finally, a considerable number of isolated events in the SGBSN catalog cannot be correlated with mapped faults.

Journal of Geophysical Research

Three-dimensional seismic velocity structure of Mauna Loa and Kilauea volcanoes in Hawaii from local seismic tomography

We present a new three-dimensional seismic velocity model of the crustal and upper mantle structure for Mauna Loa and Kilauea volcanoes in Hawaii. Our model is derived from the first-arrival times of the compressional and shear waves from about 53,000 events on and near the Island of Hawaii between 1992 and 2009 recorded by the Hawaiian Volcano Observatory stations. The V p model generally agrees with previous studies, showing high-velocity anomalies near the calderas and rift zones and low-velocity anomalies in the fault systems. The most significant difference from previous models is in V p / V s structure. The high- V p and high- V p / V s anomalies below Mauna Loa caldera are interpreted as mafic magmatic cumulates. The observed low- V p and high- V p / V s bodies in the Kaoiki seismic zone between 5 and 15 km depth are attributed to the underlying volcaniclastic sediments. The high- V p and moderate- to low- V p / V s anomalies beneath Kilauea caldera can be explained by a combination of different mafic compositions, likely to be olivine-rich gabbro and dunite. The systematically low- V p and low- V p / V s bodies in the southeast flank of Kilauea may be caused by the presence of volatiles. Another difference between this study and previous ones is the improved V p model resolution in deeper layers, owing to the inclusion of events with large epicentral distances. The new velocity model is used to relocate the seismicity of Mauna Loa and Kilauea for improved absolute locations and ultimately to develop a high-precision earthquake catalog using waveform cross-correlation data.

Hawaii

Tracking magma volume recovery at Okmok Volcano using GPS and an unscented kalman filter

Changes beneath a volcano can be observed through position changes in a GPS network, but distinguishing the source of site motion is not always straightforward. The records of continuous GPS sites provide a favorable data set for tracking magma migration. Dense campaign observations usually provide a better spatial picture of the overall deformation field, at the expense of an episodic temporal record. Combining these observations provides the best of both worlds. A Kalman filter provides a means for integrating discrete and continuous measurements and for interpreting subtle signals. The unscented Kalman filter (UKF) is a nonlinear method for time‐dependent observations. We demonstrate the application of this technique to deformation data by applying it to GPS data collected at Okmok volcano. Seven years of GPS observations at Okmok are analyzed using a Mogi source model and the UKF. The deformation source at Okmok is relatively stable at 2.5 km depth below sea level, located beneath the center of the caldera, which means the surface deformation is caused by changes in the strength of the source. During the 7 years of GPS observations more than 0.5 m of uplift has occurred, a majority of that during the time period January 2003 to July 2004. The total volume recovery at Okmok since the last eruption in 1997 is ∼60–80%. The UKF allows us to solve simultaneously for the time‐dependence of the source strength and for the location without a priori information about the source. 2009 by the American Geophysical Union.

Alaska

Viscoelastic coupling model of the San Andreas fault along the Big Bend, southern California

The big bend segment of the San Andreas fault is the 300-km-long segment in southern California that strikes about N65°W, roughly 25° counterclockwise from the local tangent to the small circle about the Pacific-North America pole of rotation. The broad distribution of deformation of trilateration networks along this segment implies a locking depth of at least 25 km as interpreted by the conventional model of strain accumulation (continuous slip on the fault below the locking depth at the rate of relative plate motion), whereas the observed seismicity and laboratory data on fault strength suggest that the locking depth should be no greater than 10 to 15 km. The discrepancy is explained by the viscoelastic coupling model which accounts for the viscoelastic response of the lower crust. Thus the broad distribution of deformation observed across the big bend segment can be largely associated with the San Andreas fault itself, not subsidiary faults distributed throughout the region. The Working Group on California Earthquake Probabilities [1995] in using geodetic data to estimate the seismic risk in southern California has assumed that strain accumulated off the San Andreas fault is released by earthquakes located off the San Andreas fault. Thus they count the San Andreas contribution to total seismic moment accumulation more than once, leading to an overestimate of the seismicity for magnitude 6 and greater earthquakes in their Type C zones.

Journal of Geophysical Research B: Solid Earth

Chemistry of rock-forming minerals of the Cretaceous-Paleocene batholith in southwestern Japan and implications for magma genesis

Petrographic descriptions and electron microprobe analyses of minerals are presented for 35 specimens from seven suites chosen to examine the transition from magnetite series to ilmenite series granitoids along two transects across the Cretaceous-Paleocene Inner Zone batholith of southwestern Japan. Regularities in chemical compositions of amphiboles, biotites, and feldspars suggest that fundamentally similar processes produced the magmas that formed the two series. Constant or decreasing Fe/(Fe+Mg) for biotites and amphiboles with increasing host rock silica content, coupled with the absence of early formed magnetite and sphene, suggest that magnetite series rocks may have become oxidized during crystallization near the level of intrusion, through the processes of second boiling and differential loss of hydrogen. For the Daito-Yokota, magnetite series suite, Fe/(Fe+Mg) for biotites decreased from 0.48 to 0.37 as SiO 2 content of the host rock increased from 55.3 to 75.5 wt %; for an ilmenite series suite from the Takanawa Peninsula, Fe/(Fe+Mg) for biotites increased from 0.51 to 0.77 with an increase in host rock SiO 2 from 53.4 to 75.5. Detailed consideration of amphibole chemistry shows predominance of edenitic and tschermakitic substitution schemes as well as coupling between substitutions of Ti in octahedral sites and Al IV . Interrelations between amphibole and biotite chemistry show that Fe/(Fe+Mg) and Mn contents can be interpreted in terms of equilibration, whereas Ti content cannot. The chemistry of chlorites correlates well with that of biotites; primary and secondary muscovites are distinct in composition. Plagioclase in all studied suites shows igneous zoning appropriate to host rock composition; perthitic alkali feldspars in all samples have lost albite component, and temperatures based on the two-feldspar geothermometer are low. The biotite-apatite geothermometer is also inoperative for this group of samples because fully fluorinated apatites typically occur in biotites of modest F content. Whereas magnetites have reequilibrated, analyses of ilmenites for the representative Daito-Yokota and Takanawa suites corroborate biotite compositional data and suggest that f O2 probably differed by 2 to 3 orders of magnitude during crystallization of silica-rich magnetite and ilmenite series granites. Whole-rock chemistry supports mineral chemistry in suggesting that the studied granitoids have crystallized from magmas generated in a lower crustal environment in which mantle-derived magmas partially melted source materials with igneous characteristics.

Journal of Geophysical Research Solid Earth

Combined interpretation of radar, hydraulic, and tracer data from a fractured-rock aquifer near Mirror Lake, New Hampshire, USA

An integrated interpretation of field experimental cross-hole radar, tracer, and hydraulic data demonstrates the value of combining time-lapse geophysical monitoring with conventional hydrologic measurements for improved characterization of a fractured-rock aquifer. Time-lapse difference-attenuation radar tomography was conducted during saline tracer experiments at the US Geological Survey Fractured Rock Hydrology Research Site near Mirror Lake, Grafton County, New Hampshire, USA. The presence of electrically conductive saline tracer effectively illuminates permeable fractures or pathways for geophysical imaging. The geophysical results guide the construction of three-dimensional numerical models of ground-water flow and solute transport. In an effort to explore alternative explanations for the tracer and tomographic data, a suite of conceptual models involving heterogeneous hydraulic conductivity fields and rate-limited mass transfer are considered. Calibration data include tracer concentrations, the arrival time of peak concentration at the outlet, and steady-state hydraulic head. Results from the coupled inversion procedure suggest that much of the tracer mass migrated outside the three tomographic image planes, and that solute is likely transported by two pathways through the system. This work provides basic and site-specific insights into the control of permeability heterogeneity on ground-water flow and solute transport in fractured rock.

New Hampshire

Surface processes recorded by rocks and soils on Meridiani Planum, Mars: Microscopic Imager observations during Opportunity's first three extended missions

The Microscopic Imager (MI) on the Mars Exploration Rover Opportunity has returned images of Mars with higher resolution than any previous camera system, allowing detailed petrographic and sedimentological studies of the rocks and soils at the Meridiani Planum landing site. Designed to simulate a geologist's hand lens, the MI is mounted on Opportunity's instrument arm and can resolve objects 0.1 mm across or larger. This paper provides an overview of MI operations, data calibration, and analysis of MI data returned during the first 900 sols (Mars days) of the Opportunity landed mission. Analyses of Opportunity MI data have helped to resolve major questions about the origin of observed textures and features. These studies support eolian sediment transport, rather than impact surge processes, as the dominant depositional mechanism for Burns formation strata. MI stereo observations of a rock outcrop near the rim of Erebus Crater support the previous interpretation of similar sedimentary structures in Eagle Crater as being formed by surficial flow of liquid water. Well‐sorted spherules dominate ripple surfaces on the Meridiani plains, and the size of spherules between ripples decreases by about 1 mm from north to south along Opportunity's traverse between Endurance and Erebus craters.

Journal of Geophysical Research E: Planets

Three-dimensional P velocity structure of the crust below Newberry Volcano, Oregon

Utilizing teleseismic P residuals, we have detected a column of high P velocity material extending from within 10 km of the surface below the summit of Newberry Volcano, Oregon, to midcrustal depths near 25 km. We interpret this column to be the expression of a swarm of predominantly subsolidus gabbroic sills and dikes which were intruded as the volcano was built. The high P velocities observed below the volcano severely limit the size of magma chambers which could presently exist in the crust below Newberry Volcano. Those possible include a few percent of partial melt distributed through large volumes of a mafic intrusion zone in the midcrust; a few, smaller, higher melt fraction zones in the midcrust with dimensions less than 6 km and whose aggregate volume is only a few percent of enclosing volumes of 200 km 3 ; small magma bodies with dimensions of a few kilometers located within the upper 10 km of the crust; or a mafic, crystal-rich magma of arbitrary dimensions located in the upper few km. The low P velocities detected in the upper 4 km beneath the center of the summit caldera may be partially caused by a magma chamber in the second of these catagories.

Journal of Geophysical Research Solid Earth

Algal blooms and "Marine snow": Mechanisms that enhance preservation of organic carbon in ancient fine-grained sediments

Combined petrographic and geochemical methods are used to investigate the microfabrics present in thin sections prepared from representative organic carbon-rich mudstones collected from three successions (the Kimmeridge Clay Formation, the Jet Rock Member of the Whitby Mudstone Formation, and the pebble shale and Hue Shale). This study was initiated to determine how organic carbon-rich materials were being delivered to the sediment–water interface, and what happened to them after deposition, prior to deep burial. Analyses of the fabrics present shows that they exhibit many common attributes. In particular they are all: (1) highly heterogeneous on the scale of a thin section, (2) organized into thin beds (< 10 mm thick) composed mainly of mineral mixtures of fine-grained siliciclastic detritus and carbonate materials, and (3) contain significant concentrations of organic carbon, much of which is organized into laminasets that contain abundant organomineralic aggregates and pellets. In addition, framboidal pyrite (range of sizes from < 20 μm to < 1 μm) and abundant agglutinated foraminifers are present in some units. The individual beds are commonly sharp-based and overlain by thin, silt lags. The tops of many of the beds have been homogenized and some regions of the pelleted laminasets contain small horizontal burrows. The organomineralic aggregates present in these mudstones are interpreted to be ancient examples of marine snow. This marine snow likely formed in the water column, particularly during phytoplankton blooms, and was then transported rapidly to the seafloor. The existence of the thin beds with homogenized tops and an in-situ infauna indicates that between blooms there was sufficient oxygen and time for a mixed layer to develop as a result of sediment colonization by diminutive organisms using either aerobic or dysaerobic metabolic pathways. These textures suggest that the constituents of these mudstones were delivered neither as a continuous rain of sediment nor were the bottom waters persistently anoxic. In addition, the presence of thin lags and sharp-based beds suggests that the seafloor was being episodically reworked during deposition. These fabrics indicate that conditions in the water columns and at the seafloors while these rocks were being deposited were very dynamic, and episodic fluxes of high concentrations of organic carbon to the seafloor, during phytoplankton blooms, likely enhanced preservation of organic carbon.

Journal of Sedimentary Research

Since “Groundwater and surface water–A single resource”: some U.S. Geological Survey advances in modeling groundwater/surface-water interactions

Field and interpretive studies conducted by T.C. Winter and U.S. Geological Survey colleagues, and summarized in the 1998 publication “Groundwater and Surface Water – A Single Resource”, inspired a new generation of research centered on extensions of the groundwater-flow code MODFLOW to more sophisticated simulation of coupled groundwater and surface-water systems. Guided by emerging concerns with water availability, safe yields from wells, health of aquatic habitat, and climate change, the changes to MODFLOW involve: 1) the ability to more precisely and accurately represent the interface between surface and subsurface flows and 2) the consideration of a variety of mechanisms that influence their interaction. A review of the most important changes to the code is supplemented in this article by selected case studies in an effort to show the scope of the advances. The updates discussed include the Streamflow Routing (SFR), Lake (LAK), and Unsaturated-Zone Flow (UZF) Packages in MODFLOW-2005 and the Groundwater Management (GWM), Local Grid Refinement (LGR), and Newton (NWT) formulation versions of MODFLOW-2005. New developments feature the integration of rainfall-runoff modeling with MODFLOW in GSFLOW, coupling of GFLOW and MODFLOW in a hybrid code, and the forthcoming unstructured grid version of MODFLOW. They promise continued advances in the ability to use science to protect the single water resource.

Acque Sotterranee: Italian Journal of Groundwater

Field measurement of basal forces generated by erosive debris flows

It has been proposed that debris flows cut bedrock valleys in steeplands worldwide, but field measurements needed to constrain mechanistic models of this process remain sparse due to the difficulty of instrumenting natural flows. Here we present and analyze measurements made using an automated sensor network, erosion bolts, and a 15.24 cm by 15.24 cm force plate installed in the bedrock channel floor of a steep catchment. These measurements allow us to quantify the distribution of basal forces from natural debris‒flow events that incised bedrock. Over the 4 year monitoring period, 11 debris‒flow events scoured the bedrock channel floor. No clear water flows were observed. Measurements of erosion bolts at the beginning and end of the study indicated that the bedrock channel floor was lowered by 36 to 64 mm. The basal force during these erosive debris‒flow events had a large‒magnitude (up to 21 kN, which was approximately 50 times larger than the concurrent time‒averaged mean force), high‒frequency (greater than 1 Hz) fluctuating component. We interpret these fluctuations as flow particles impacting the bed. The resulting variability in force magnitude increased linearly with the time‒averaged mean basal force. Probability density functions of basal normal forces were consistent with a generalized Pareto distribution, rather than the exponential distribution that is commonly found in experimental and simulated monodispersed granular flows and which has a lower probability of large forces. When the bed sediment thickness covering the force plate was greater than ~ 20 times the median bed sediment grain size, no significant fluctuations about the time‒averaged mean force were measured, indicating that a thin layer of sediment (~ 5 cm in the monitored cases) can effectively shield the subjacent bed from erosive impacts. Coarse‒grained granular surges and water‒rich, intersurge flow had very similar basal force distributions despite differences in appearance and bulk‒flow density. These results demonstrate that debris flows can have strong control on rates of steepland evolution and contribute to a foundation needed for modeling debris‒flow incision stochastically.

Journal of Geophysical Research F: Earth Surface