Search USGS⌕ Search

SEARCH · Search USGS

Results for “Research Square”

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 307 records · Page 17Linked to original sources
USGS2005

An efficient algorithm for double-difference tomography and location in heterogeneous media, with an application to the Kilauea volcano

Improving our understanding of crustal processes requires a better knowledge of the geometry and the position of geological bodies. In this study we have designed a method based upon double-difference relocation and tomography to image, as accurately as possible, a heterogeneous medium containing seismogenic objects. Our approach consisted not only of incorporating double difference in tomography but also partly in revisiting tomographic schemes for choosing accurate and stable numerical strategies, adapted to the use of cross-spectral time delays. We used a finite difference solution to the eikonal equation for travel time computation and a Tarantola-Valette approach for both the classical and double-difference three-dimensional tomographic inversion to find accurate earthquake locations and seismic velocity estimates. We estimated efficiently the square root of the inverse model's covariance matrix in the case of a Gaussian correlation function. It allows the use of correlation length and a priori model variance criteria to determine the optimal solution. Double-difference relocation of similar earthquakes is performed in the optimal velocity model, making absolute and relative locations less biased by the velocity model. Double-difference tomography is achieved by using high-accuracy time delay measurements. These algorithms have been applied to earthquake data recorded in the vicinity of Kilauea and Mauna Loa volcanoes for imaging the volcanic structures. Stable and detailed velocity models are obtained: the regional tomography unambiguously highlights the structure of the island of Hawaii and the double-difference tomography shows a detailed image of the southern Kilauea caldera-upper east rift zone magmatic complex. Copyright 2005 by the American Geophysical Union.

Hawaii↗
USGS2009

Hydroecological factors governing surface water flow on a low-gradient floodplain

"Interrelationships between hydrology and aquatic ecosystems are better understood in streams and rivers compared to their surrounding floodplains. Our goal was to characterize the hydrology of the Everglades ridge and slough floodplain ecosystem, which is valued for the comparatively high biodiversity and connectivity of its parallel-drainage features but which has been degraded over the past century in response to flow reductions associated with flood control. We measured flow velocity, water depth, and wind velocity continuously for 3 years in an area of the Everglades with well-preserved parallel-drainage features (i.e., 200-m wide sloughs interspersed with slightly higher elevation and more densely vegetated ridges). Mean daily flow velocity averaged 0.32 cm s1 and ranged between 0.02 and 0.79 cm s1. Highest sustained velocities were associated with flow pulses caused by water releases from upstream hydraulic control structures that increased flow velocity by a factor of 2–3 on the floodplain for weeks at a time. The highest instantaneous measurements of flow velocity were associated with the passage of Hurricane Wilma in 2005 when the inverse barometric pressure effect increased flow velocity up to 5 cm s1 for several hours. Time-averaged flow velocities were 29% greater in sloughs compared to ridges because of marginally higher vegetative drag in ridges compared to sloughs, which contributed modestly (relative to greater water depth and flow duration in sloughs compared to ridges) to the predominant fraction (86%) of total discharge through the landscape occurring in sloughs. Univariate scaling relationships developed from theory of flow through vegetation, and our field data indicated that flow velocity increases with the square of water surface slope and the fourth power of stem diameter, decreases in direct proportion with increasing frontal area of vegetation, and is unrelated to water depth except for the influence that water depth has in controlling the submergence height of vegetation that varies vertically in its architectural characteristics. In the Everglades the result of interactions among controlling variables was that flow velocity was dominantly controlled by water surface slope variations responding to flow pulses more than spatial variation in vegetation characteristics or fluctuating water depth. Our findings indicate that floodplain managers could, in addition to managing water depth, manipulate the frequency and duration of inflow pulses to manage water surface slope, which would add further control over flow velocities, water residence times, sediment settling, biogeochemical transformations, and other processes that are important to floodplain function."

Water Resources Research↗
USGS2009

Global distribution of bedrock exposures on Mars using THEMIS high-resolution thermal inertia

We investigate high thermal inertia surfaces using the Mars Odyssey Thermal Emission Imaging System (THEMIS) nighttime temperature images (100 m/pixel spatial sampling). For this study, we interpret any pixel in a THEMIS image with a thermal inertia over 1200 J m -2 K -1 s -1/2 as "bedrock" which represents either in situ rock exposures or rock-dominated surfaces. Three distinct morphologies, ranked from most to least common, are associated with these high thermal inertia surfaces: (1) valley and crater walls associated with mass wasting and high surface slope angles; (2) floors of craters with diameters >25 km and containing melt or volcanics associated with larger, high-energy impacts; and (3) intercrater surfaces with compositions significantly more mafic than the surrounding regolith. In general, bedrock instances on Mars occur as small exposures (less than several square kilometers) situated in lower-albedo (<0.18), moderate to high thermal inertia (>350 J m -2 K -1 s -1/2 ), and relatively dust-free (dust cover index <0.95) regions; however, there are instances that do not follow these generalizations. Most instances are concentrated in the southern highlands, with very few located at high latitudes (poleward of 45 o N and 58 o S), suggesting enhanced mechanical breakdown probably associated with permafrost. Overall, Mars has very little exposed bedrock with only 960 instances identified from 75 o S to 75 o N with likely <3500 km 2 exposed, representing???1% of the total surface area. These data indicate that Mars has likely undergone large-scale surface processing and reworking, both chemically and mechanically, either destroying or masking a majority of the bedrock exposures on the planet. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2017

Geologic field trip guide to Mount Mazama and Crater Lake Caldera, Oregon

Crater Lake partly fills one of the most spectacular calderas of the world—an 8 by 10 kilometer (km) basin more than 1 km deep formed by collapse of the Mount Mazama volcano during a rapid series of explosive eruptions ~7,700 years ago. Having a maximum depth of 594 meters (m), Crater Lake is the deepest lake in the United States. Crater Lake National Park, dedicated in 1902, encompasses 645 square kilometers (km 2 ) of pristine forested and alpine terrain, including the lake itself, and virtually all of Mount Mazama. The geology of the area was first described in detail by Diller and Patton (1902) and later by Williams (1942), whose vivid account led to international recognition of Crater Lake as the classic collapse caldera. Because of excellent preservation and access, Mount Mazama, Crater Lake caldera, and the deposits formed by the climactic eruption constitute a natural laboratory for study of volcanic and magmatic processes. For example, the climactic ejecta are renowned among volcanologists as evidence for systematic compositional zonation within a subterranean magma chamber. Mount Mazama’s climactic eruption also is important as the source of the widespread Mazama ash, a useful Holocene stratigraphic marker throughout the Pacific Northwest United States, adjacent Canada, and offshore. A detailed bathymetric survey of the floor of Crater Lake in 2000 (Bacon and others, 2002) provides a unique record of postcaldera eruptions, the interplay between volcanism and filling of the lake, and sediment transport within this closed basin. Knowledge of the geology and eruptive history of the Mount Mazama edifice, enhanced by the caldera wall exposures, gives exceptional insight into how large volcanoes of magmatic arcs grow and evolve. In addition, many smaller volcanoes of the High Cascades beyond the limits of Mount Mazama provide information on the flux of mantle-derived magma through the region. General principles of magmatic and eruptive processes revealed by geologic research at Crater Lake have been incorporated not only in scientific investigations elsewhere, but also in the practical evaluation of local hazards (Bacon and others, 1997b) and geothermal resources (Bacon and Nathenson, 1996). The 1:24,000-scale geologic map of Mount Mazama and Crater Lake caldera (Bacon, 2008) is unusual because it portrays bedrock (outcrop), surficial, and lake floor geology. Caldera wall geology is depicted in detail on the accompanying geologic panoramas, and bedrock geology is shown in a 1:50,000-scale geologic map. This field guide supersedes earlier geology guides of Crater Lake (Bacon, 1987, 1989).

Oregon↗
USGS1978

Heavy-mineral variability in the Baltimore Canyon trough area

Petrographic analyses of bottom sediments from 87 stations within a relatively large subarea (1700 square kilometers) define the local variability and the distributional processes of heavy minerals in the Baltimore Canyon Trough area (13500km 2 ). Of the 29 mineral groups that were identified, those most diagnostic of differences between stations were opaque minerals, garnet, hornblende, orthopyroxene, other amphiboles, and staurolite. Some of the common components cannot be used to define areal trends. The association of minerals by their specific gravities reflects the modern reworking and sorting of sediments in this area. The poor correlation between mineral abundance and the ridge-and trough topography suggests either that our analyses and samples did not resolve bathymetrically induced changes or that smaller bed forms (such as ripples) largely control the heavy-mineral distribution. Data from this study provide a baseline for evaluating manmade disturbances of the bottom sediments and show that relatively large changes in the heavy mineral composition can take place over small distances on this part of the Continental Shelf.

Baltimore Canyon Trough↗
USGS2012

Dynamic rupture modeling with laboratory-derived constitutive relations

A laboratory-derived state variable friction constitutive relation is used in the numerical simulation of the dynamic growth of an in-plane or mode II shear crack. According to this formulation, originally presented by J. H. Dieterich, frictional resistance varies with the logarithm of the slip rate and with the logarithm of the frictional state variable as identified by A. L. Ruina. Under conditions of steady sliding, the state variable is proportional to (slip rate) −1 . Following suddenly introduced increases in slip rate, the rate and state dependencies combine to produce behavior which resembles slip weakening. When rupture nucleation is artificially forced at fixed rupture velocity, rupture models calculated with the state variable friction in a uniformly distributed initial stress field closely resemble earlier rupture models calculated with a slip weakening fault constitutive relation. Additional rupture models are calculated in which rupture nucleation is achieved naturally, with numerical simulations of the quasi-static response of the fault leading to the onset of unstable, dynamic rupture. When rupture nucleation with the state variable friction law takes place naturally, a large fraction of the fault accelerates before accelerating slip is concentrated in what ultimately becomes the rupture nucleation patch. The state evolution accompanying this accelerating slip leads to higher average rupture speeds or a more rapid rupture acceleration to near P wave rupture speeds. Rupture models are also calculated for the seismological asperity problem, that is, the failure of a highly stressed fault patch surrounded by a region of zero stress drop. Dynamic overshoot of slip into the region of zero stress drop roughly agrees with a simple energy balance analysis; the final size of the rupture is proportional to the square of the size of the high stress patch. Earlier frictional stability analyses have led to the definition of a critical fault patch size for rupture nucleation. This critical patch size is generally different from critical crack lengths determined from crack tip energy balance considerations applied to a simpler slip weakening law. In the model calculations, dynamic rupture does not nucleate if the starting patch size is less than the critical patch size. This is consistent with the frictional stability analyses. Thus these model calculations suggest that dynamic rupture following a state variable friction relation is similar to that following a simpler fault slip weakening law. However, when modeling the full cycle of fault motions, rate-dependent frictional responses included in the state variable formulation are important at low slip rates associated with rupture nucleation. The critical rupture nucleation dimension appropriate for a slip weakening fault does not predict the critical nucleation dimension for a state variable fault.

Journal of Geophysical Research Solid Earth↗
USGS2002

Flow distribution in selected branches of St. Clair and Detroit rivers

St. Clair and Detroit rivers, which are connecting channels between Lake Huron and Lake Erie in the Great Lakes basin, form part of the boundary between the state of Michigan and the province of Ontario. In 13 reaches, this flow divides locally around islands and dikes to form 31 branches. This study develops a set of simple linear regression equations for computing expected flow proportions in branches, generally as a function of the total flow within the reach. The equations are based on 533 acoustic Doppler current profiler measurements of flow obtained between 1996 and 2000. Root-mean-square errors of these regressions range from 0.00323 to 0.0895. In seven upstream reaches where flow is known because of flow specifications at the boundaries of the waterway and continuity constraints, the uncertainties of the flow proportions can be used to directly infer the uncertainties of the corresponding flows. In six downstream reaches, the uncertainties of flows are determined by both the uncertainties of the flow proportions and the uncertainties of the total flow in the reach. For these reaches, Monte Carlo simulations quantify the ratios of total uncertainty to flow proportion uncertainty, which range from 1.0026 to 13.984. To facilitate routine calculation, polynomial regression equations are developed to approximate these ratios as a function of flow. Results provide a mechanism for computing the magnitudes and uncertainties of steady-state flows within selected branches of the connecting channels by specifying inflows at the headwaters of St. Clair River, seven intervening tributaries, and Lake St. Clair.

Michigan, Ontario↗
USGS2020

Exploring the exceptional performance of a deep learning stream temperature model and the value of streamflow data

Stream water temperature ( T s ) is a variable of critical importance for aquatic ecosystem health. T s is strongly affected by groundwater-surface water interactions which can be learned from streamflow records, but previously such information was challenging to effectively absorb with process-based models due to parameter equifinality. Based on the long short-term memory (LSTM) deep learning architecture, we developed a basin-centric lumped daily mean T s model, which was trained over 118 data-rich basins with no major dams in the conterminous United States, and showed strong results. At a national scale, we obtained a median root-mean-square error of 0.69°C, Nash–Sutcliffe model efficiency coefficient of 0.985, and correlation of 0.994, which are marked improvements over previous values reported in literature. The addition of streamflow observations as a model input strongly elevated the performance of this model. In the absence of measured streamflow, we showed that a two-stage model could be used, where simulated streamflow from a pre-trained LSTM model ( Q sim ) still benefited the T s model even though no new information was brought directly into the inputs of the T s model. The model indirectly used information learned from streamflow observations provided during the training of Q sim , potentially to improve internal representation of physically meaningful variables. Our results indicate that strong relationships exist between basin-averaged forcing variables, catchment attributes, and T s that can be simulated by a single model trained by data on the continental scale.

Environmental Research Letters↗
USGS2009

Methodology for Estimation of Flood Magnitude and Frequency for New Jersey Streams

Methodologies were developed for estimating flood magnitudes at the 2-, 5-, 10-, 25-, 50-, 100-, and 500-year recurrence intervals for unregulated or slightly regulated streams in New Jersey. Regression equations that incorporate basin characteristics were developed to estimate flood magnitude and frequency for streams throughout the State by use of a generalized least squares regression analysis. Relations between flood-frequency estimates based on streamflow-gaging-station discharge and basin characteristics were determined by multiple regression analysis, and weighted by effective years of record. The State was divided into five hydrologically similar regions to refine the regression equations. The regression analysis indicated that flood discharge, as determined by the streamflow-gaging-station annual peak flows, is related to the drainage area, main channel slope, percentage of lake and wetland areas in the basin, population density, and the flood-frequency region, at the 95-percent confidence level. The standard errors of estimate for the various recurrence-interval floods ranged from 48.1 to 62.7 percent. Annual-maximum peak flows observed at streamflow-gaging stations through water year 2007 and basin characteristics determined using geographic information system techniques for 254 streamflow-gaging stations were used for the regression analysis. Drainage areas of the streamflow-gaging stations range from 0.18 to 779 mi2. Peak-flow data and basin characteristics for 191 streamflow-gaging stations located in New Jersey were used, along with peak-flow data for stations located in adjoining States, including 25 stations in Pennsylvania, 17 stations in New York, 16 stations in Delaware, and 5 stations in Maryland. Streamflow records for selected stations outside of New Jersey were included in the present study because hydrologic, physiographic, and geologic boundaries commonly extend beyond political boundaries. The StreamStats web application was developed cooperatively by the U.S. Geological Survey and the Environmental Systems Research Institute, Inc., and was designed for national implementation. This web application has been recently implemented for use in New Jersey. This program used in conjunction with a geographic information system provides the computation of values for selected basin characteristics, estimates of flood magnitudes and frequencies, and statistics for stream locations in New Jersey chosen by the user, whether the site is gaged or ungaged.

Scientific Investigations Report↗
USGS2023

Assessing global elevation models for mapping the low elevation coastal zone

Elevation data are critical for assessments of coastal hazards, including sea-level rise (SLR), flooding, storm surge, tsunami impacts, and wave run-up. Previous research has demonstrated that the quality of data used in elevation-based hazard assessments must be well documented and applied properly to assess potential impacts. Global digital elevation models (DEMs), at 30- to 90-meter resolution, have been used extensively to map and characterize coastal environments and the at-risk resources (population and built structures) contained therein. The inherent absolute vertical accuracy of global DEMs precludes their usefulness for assessing exposure to fine increments (< 1 meter) of coastal inundation at high confidence levels. However, global DEMs are highly suitable for delineation of the global low elevation coastal zone (LECZ) (elevation < 10 meters). An accuracy evaluation of global DEMs over the United States has been conducted to quantify their performance in correctly mapping the LECZ, namely in terms of vertical uncertainty and corresponding confidence levels for several representations of the coastal zone. The evaluation approach includes comparison of the DEMs with an extensive set of high-accuracy geodetic control points as the independent reference data covering a variety of coastal relief settings. The 1-arc-second (30-meter) global DEMs evaluated include ALOS World 3D, ASTER GDEM, Copernicus, FABDEM, and NASADEM, and the 3-arc-second (90-meter) global DEMs include CoastalDEM, Copernicus, MERIT, and TanDEM-X. Additionally, lower resolution (1-kilometer) global DEMs were also assessed, namely the Global Lidar Lowland DTM (derived from ICESat-2) and the GEDI 1-km DEM. The results of the accuracy characterization show that FABDEM performs the best (minimal vertical bias and lowest vertical root mean square error) for high-confidence mapping of the LECZ. Among 90-m DEMs, CoastalDEM performs best, although the differences across datasets are minimal. The results also demonstrate the importance of rigorously accounting for elevation uncertainty when applying global DEMs for coastal mapping applications.

Conference Paper↗
USGS1970

Annotated bibliography on artificial recharge of ground water, 1955-67

Artificial ground-water recharge has become more important as water use by agriculture, industry, and municipalities increases. Water management agencies are increasingly interested in potential use of recharge for pollution abatement, waste-water disposal, and re-use and reclamation of locally available supplies. Research projects and theoretical analyses of operational recharge systems show increased scientific emphasis on the practice. Overall ground-water basin management systems generally now contain considerations of artificial recharge, whether by direct or indirect methods. Artificial ground-water recharge is a means of conserving surface runoff for future use in places where it would otherwise be lost, of protecting ground-water basins from salt-water encroachment along coastal areas, and of storing and distributing imported water. The biblio-graphy emphasizes technology; however, annotations of articles on waste-water reclamation, ground-water management and ground-water basin management are included. Subjects closely related to artificial recharge, including colloidal flow through porous media, field or laboratory instrumentation, and waste disposal by deep well injection are included where they specifically relate to potential recharge problems. Where almost the same material has been published in several journals, all references are included on the assumption that some publications may be more readily available to interested persons than others. Other publications, especially those of foreign literature, provided abstracts that were used freely as time limitations precluded obtaining and annotating all materials. Abstracts taken from published sources are noted. These are: "Abstracts of North American Geology," U.S. Department of the Interior, Geological Survey; "Abstracts of Recent Published Material on Foil and Water Conservation," ARS-41 series, Agricultural F.esearch Service, U.S. Department of Agriculture; "Water and1 Water Engineering," published by Fuel and Metallurgical Journals, Ltd., London, England; "Journal of Geophysical Research," American Geophysical Union, Washington, D.C.; "American Society of Civil Engineers Transactions," New York; "Selected Bibliography of Hydrology, United Kingdom, for the Years 1955-59," International Association of Scientific Hydrology; "Water Wells, an Annotated Bibliography," California University Water Resources Center Archives Report 13; "Re-use of Effluent in the Future With an Annotated Bibliography," by G. A. Whetstone, Texas Water Development Board Report 8, Austin, Tex.; "Journal of Water Pollution Control Federation," Washington, D.C.; and "A List of Selected Technical References on Artificial Recharge of Ground-Water Reservoirs," compiled by Roy W. Graves, Tulsa University, Information Services Department, Tulsa, Okla. Other notations are self-explanatory, and initials are those of the authors (DCS, DJG, WK). An unpublished compilation of recharge references by Arnon Arad sponsored by the United Nations Educational, Scientific, and Cultural Organization during a training period with the U.S. Geological Survey was also used. The bibliography is arranged alphabetically by author. Where an author has more than one publication, the arrangement is chronological; where an author has more than one publication in a given year, a, b, c, . . . are added. The indexing is by subject and geographic location. Each article was assigned the key words or phrases to best characterize its contents. Units of measure are as they were in the original article; abbreviations retained are generally those in common use such as mg/1 (milligrams per liter), ppm (parts per million), gpm (gallons per minute), km (kilometers), m (meters), cu m per hr (cubic meters p^r hour), cfs (cubic feet per second), me/1 (milliequivalents per liter), psi (pounds per square inch), BOD (biochemical oxygen demand), sq m (square meters), gpd (gallons per day), and mgd (million gallons per day). The bibliography was prepared because of the worldwide interest in the field of artificial recharge and the need for a single source of references to the literature published since 1954. The work is a sequel to the "Annotated Bibliography on Artificial Recharge of Ground Water Through 1954," by D. K. Todd, U.S. Geological Survey Water-Supply Paper 1477, published in 1959.

Water Supply Paper↗
USGS2012

Fitting degradation of shoreline scarps by a nonlinear diffusion model

The diffusion model of degradation of topographic features is a promising means by which vertical offsets on Holocene faults might be dated. In order to calibrate the method, we have examined present-day profiles of wave-cut shoreline scarps of late Pleistocene lakes Bonneville and Lahontan. It may be assumed that these scarps were initially at least as steep as the angle of repose. Offsets range from 1 to 12 m, and present slope angles range from 9° to 29°. A parameter called apparent diffusion age, defined as half the mean square horizontal extent of the slope function of each profile, is plotted as a function of scarp offset. The points show a clear trend of apparent age increasing nearly linearly with offset. If linear diffusion held and scarps were initially vertical, apparent diffusion age would be the same for all the shoreline profiles. The increasing trend can only partly be explained by nonvertical initial scarp slope, and therefore the rate of transport of material downslope must increase significantly faster than a linear law in the range of slopes spanned by the data. The transport law must become linear at small slope to reduce scatter between profiles with varying ambient slopes. The transport law adopted for the purpose of dating is k 0 (1 + 5 s 2 ), where s is local slope. The transport coefficient k 0 is correlated inversely with fan slope, suggesting that there is a dependence on the particle size distribution. A table is included that allows easy application of the model to scarps with simple initial shape.

Journal of Geophysical Research Solid Earth↗
USGS2022

Subsampling large-scale digital elevation models to expedite geospatial analyses in coastal regions

Large-area, high-resolution digital elevation models (DEMs) created from light detection and ranging (LIDAR) and/or multibeam echosounder data sets are commonly used in many scientific disciplines. These DEMs can span thousands of square kilometers, typically with a spatial resolution of 1 m or finer, and can be difficult to process and analyze without specialized computers and software. Such DEMs often can be subsampled to expedite analysis with negligible impact on results for large-scale geospatial analyses. Subsampling can be achieved by creating a grid of points that specify the locations from which to extract elevation values from the DEM. This paper presents a method that can be used to accurately perform subsampling of large-scale, high-resolution DEMs using GIS software. This subsampling method was applied to two LIDAR-derived DEMs encompassing 242 km 2 of the northern Florida Reef Tract as an example application and to test subsampling accuracy. Results indicate that subsampling 1-m-resolution DEMs using a 2-m-spaced grid results in no significant difference in mean elevation or other basic statistics for analyses performed over multiple spatial scales ranging from 1 km 2 to 242 km 2 .

Journal of Coastal Research↗
USGS2018

Estimating soil respiration in a subalpine landscape using point, terrain, climate and greenness data

Landscape carbon (C) flux estimates are necessary for assessing the ability of terrestrial ecosystems to buffer further increases in anthropogenic carbon dioxide (CO2) emissions. Advances in remote sensing have allowed for coarse-scale estimates of gross primary productivity (GPP) (e.g., MODIS 17), yet efforts to assess spatial patterns in respiration lag behind those of GPP. Here, we demonstrate a method to predict growing season soil respiration at a regional scale in a forested ecosystem. We related field measurements (n=144) of growing season soil respiration across subalpine forests in the Southern Rocky Mountains ecoregion to a suite of biophysical predictors with a Random Forest model (30 m pixel size). We found that Landsat Enhanced Vegetation Index (EVI), growing season AI, temperature, precipitation, elevation, and slope aspect explained spatiotemporal variability in soil respiration. Our model had a psuedo-r2 of 0.45 and root mean squared error (RMSE) of roughly one-quarter of the mean value of respiration. Predicted growing season soil respiration across the region was remarkably consistent across 2004, 2005 and 2006 (150-d averages of 542.8, 544.3, and 536.5 g C m-2, respectively). Yet, we observed substantial variability in spatial patterns of soil respiration predictions that varied between years, suggesting that our method is sensitive to changes in respiration drivers. We compared our estimates to MODIS GPP and nocturnal net ecosystem exchange (NEE) derived from eddy covariance towers as a proxy for ecosystem respiration. Averaged across the predictive region, mean predicted growing season soil respiration was 73% of MODIS GPP, while predicted soil respiration was generally within 20% of nocturnal NEE from eddy covariance towers. This study demonstrated that geospatial and remotely-sensed datasets can be used in a statistical modeling framework to estimate soil respiration at landscape scales.

Colorado, Wyoming↗
USGS2026

Towards affordable wetland evapotranspiration monitoring using the Variance-Bowen Ratio method: Insights from three contrasting wetlands

Accurate measurement of evapotranspiration (ET) is essential for sustainable water management. Standard methods such as eddy covariance (EC) are costly, while alternatives such as surface renewal are cheaper but require calibration and complex data processing. This study evaluates the utility of the Variance-Bowen Ratio (VBR) method for estimating ET across three California’s wetlands. Using data from 2023, half-hourly latent heat flux (λE) and daily/monthly ET from VBR were compared with EC at one non-tidal (site A) and two tidal (sites B and C) wetlands. λE VBR consistently underestimated λE EC , with root mean squared errors (RMSE) of 61.2 W m -2 at sites A, 106.1 W m -2 at site B, and 137.2 W m -2 at site C, largely due to storage fluxes across sites. Temporal integration improved VBR’s performance at tidal sites, where compensating water heat storage errors yielded low daily and monthly biases (site B: RMSE = 0.78 mm/d and 12 mm/month; r = 0.93; site C: RMSE = 0.90 mm/d and 13 mm/month; r = 0.93), with reduced major axis (RMA) regression slopes of 0.98 and ~0.91. In contrast, biomass heat storage at site A caused persistent biases (RMSEs = 0.97 mm/d and 23 mm/month; daily and monthly RMA slopes ~0.75; r = 0.85). These results highlight VBR’s limitations in environments with substantial storage fluxes. Despite this, VBR is cost-effective for estimating daily and monthly ET, with sensor costs at least tenfold lower than EC and simpler setup, making it suitable for ET monitoring in resource-limited and hard-to-access regions.

California↗
USGS2017

Rock-weathering rates as functions of time

The scarcity of documented numerical relations between rock weathering and time has led to a common assumption that rates of weathering are linear. This assumption has been strengthened by studies that have calculated long-term average rates. However, little theoretical or empirical evidence exists to support linear rates for most chemical-weathering processes, with the exception of congruent dissolution processes. The few previous studies of rock-weathering rates that contain quantitative documentation of the relation between chemical weathering and time suggest that the rates of most weathering processes decrease with time. Recent studies of weathering rinds on basaltic and andesitic stones in glacial deposits in the western United States also clearly demonstrate that rock-weathering processes slow with time. Some weathering processes appear to conform to exponential functions of time, such as the square-root time function for hydration of volcanic glass, which conforms to the theoretical predictions of diffusion kinetics. However, weathering of mineralogically heterogeneous rocks involves complex physical and chemical processes that generally can be expressed only empirically, commonly by way of logarithmic time functions. Incongruent dissolution and other weathering processes produce residues, which are commonly used as measures of weathering. These residues appear to slow movement of water to unaltered material and impede chemical transport away from it. If weathering residues impede weathering processes then rates of weathering and rates of residue production are inversely proportional to some function of the residue thickness. This results in simple mathematical analogs for weathering that imply nonlinear time functions. The rate of weathering becomes constant only when an equilibrium thickness of the residue is reached. Because weathering residues are relatively stable chemically, and because physical removal of residues below the ground surface is slight, many weathering features require considerable time to reach constant rates of change. For weathering rinds on volcanic stones in the western United States, this time is at least 0.5 my.

Quaternary Research↗
USGS2017

Revised tephra volumes for Cascade Range volcanoes

Isopach maps from tephra eruptions from Mount St. Helens were reported in Carey et al. (1995) and for tephra eruptions from Glacier Peak in Gardner et al. (1998). For exponential thinning, the isopach data only define a single slope on a log thickness versus square root of area plot. Carey et al. (1995) proposed a model that was used to estimate a second slope, and volumes were presented in both studies using this model. A study by Sulpizio (2005) for estimating the second slope and square root of area where the lines intersect involves a systematic analysis of many eruptions to provide correlation equations. The purpose of this paper is to recalculate the volumes of Cascades eruptions and compare results from the two methods. In order to gain some perspective on the methods for estimating the second slope, we use data for thickness versus distance beyond the last isopach that are available for some of the larger eruptions in the Cascades. The thickness versus square root of area method is extended to thickness versus distance by developing an approximate relation between the two assuming elliptical isopachs with the source at one of the foci. Based on the comparisons made between the Carey et al. (1995) and Sulpizio (2005) methods, it is felt that the later method provides a better estimate of the second slope. For Mount St. Helens, the estimates of total volume using the Sulpizio (2005) method are generally smaller than those using the Carey et al. (1995) method. For the volume estimates of Carey et al. (1995), the volume of the May 18, 1980, eruption of Mount St. Helens is smaller than six of the eight previous eruptions. With the new volumes using the Sulpizio (2005) method, the 1980 eruption is smaller in volume than the upper end of the range for only three of the layers (Wn, Ye, and Yn) and is the same size as layer We. Thus the 1980 eruption becomes representative of the mid-range of volumes rather than being in the lower range.

Cascade Range volcanoes↗
USGS2022

Evidence of alternative trophic pathways for fish consumers in a large river system in the face of invasion

Large rivers are susceptible to anthropogenic alteration, which can result in drastic changes to their functional ecology. We evaluated spatial–temporal changes in the functional fish communities of the Upper Mississippi River System (UMRS) using data from six study reaches. Species were classified into one of 14 feeding guilds and mass per unit effort (MPUE) was then calculated for each feeding guild annually per gear type. MPUE was standardized using the multigear mean standardization method (MGMS) and log-transformed. Both ANOSIM and Chi-square tests were used to determine differences in MPUE among reaches. We then estimated functional diversity by calculating the number of functional groups ( N ), Margalef's d , Pielou's J′, Shannon's Diversity, and Simpson's Diversity Index. An AR(1) time series model was used to investigate proportional changes in each guild over 25 years. To evaluate the effect of invasive Carp species in invaded reaches, a Chow test was applied to observations between 2000 and 2005. Analyses revealed differences in the functional fish community among reaches. We found differences in functional diversity metrics among study reaches, but there was little evidence that this differed between invaded and non-invaded reaches. Results determined that invertivore/detritivores have been consistently declining system-wide, with few groups showing a net change. There was also little evidence that invasion altered the proportion of any functional guild. Evaluating the spatial–temporal patterns of functional communities is beneficial to understanding the resilience of a system and can provide further insight into its trophic needs when considering future restoration initiatives.

River Research and Applications↗