Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Names”

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 559 records · Page 31Linked to original sources
USGS2018

Growth-suppressing and algicidal properties of an extract from Arundo donax, an invasive riparian plant, against Prymnesium parvum, an invasive harmful alga

This study examined the ability of acidic and neutral/alkaline fractions of a methanolic extract from giant reed ( Arundo donax ) and of two of its constituents, gramine and skatole, to inhibit growth of the ichthyotoxic golden alga (Prymnesium parvum ) in batch culture. For this study, growth suppression was defined as inhibition of maximum cell density, algicidal activity as early occurrence of negative growth, and algistatic activity as lack of net growth. The acidic fraction did not affect algal growth. The neutral/alkaline fraction showed growth-suppressing and algicidal activities but no signs of algistatic activity – namely, cells in cultures surviving a partial-algicidal exposure concentration (causing transient negative growth) were later able to initiate positive growth but at higher concentrations, algicidal activity was full and irreversible. Gramine suppressed growth more effectively than skatole and at the highest concentration tested, gramine also showed partial-algicidal and algistatic activity. While the partial-algicidal activities of the neutral/alkaline fraction and of gramine were short-lived (≤6 days) and thus may share similar mechanisms, algistatic activity was unique to gramine and persisted for >3 weeks. Given gramine’s reported concentration in the neutral/alkaline fraction, its corresponding level of algicidal activity is much lower than the fraction’s suggesting the latter contains additional potent algicides. Inhibition of maximum cell density by all test compounds was associated with reductions in exponential growth rate, and in the case of the neutral/alkaline fraction and gramine also reductions in early (pre-exponential) growth. These results indicate that giant reed is a potential source of natural products to control golden alga blooms. Giant reed is an invasive species in North America, thus also providing incentive for research into strategies to couple management efforts for both species.

Harmful Algae↗
USGS2005

The rayed crater Zunil and interpretations of small impact craters on Mars

A 10-km diameter crater named Zunil in the Cerberus Plains of Mars created ∼ 10 7 "> ∼107 secondary craters 10 to 200 m in diameter. Many of these secondary craters are concentrated in radial streaks that extend up to 1600 km from the primary crater, identical to lunar rays. Most of the larger Zunil secondaries are distinctive in both visible and thermal infrared imaging. MOC images of the secondary craters show sharp rims and bright ejecta and rays, but the craters are shallow and often noncircular, as expected for relatively low-velocity impacts. About 80% of the impact craters superimposed over the youngest surfaces in the Cerberus Plains, such as Athabasca Valles, have the distinctive characteristics of Zunil secondaries. We have not identified any other large (⩾10 km diameter) impact crater on Mars with such distinctive rays of young secondary craters, so the age of the crater may be less than a few Ma. Zunil formed in the apparently youngest (least cratered) large-scale lava plains on Mars, and may be an excellent example of how spallation of a competent surface layer can produce high-velocity ejecta (Melosh, 1984, Impact ejection, spallation, and the origin of meteorites, Icarus 59, 234–260). It could be the source crater for some of the basaltic shergottites, consistent with their crystallization and ejection ages, composition, and the fact that Zunil produced abundant high-velocity ejecta fragments. A 3D hydrodynamic simulation of the impact event produced 10 10 rock fragments ⩾10 cm diameter, leading to up to 10 9 secondary craters ⩾10 m diameter. Nearly all of the simulated secondary craters larger than 50 m are within 800 km of the impact site but the more abundant smaller (10–50 m) craters extend out to 3500 km. If Zunil is representative of large impact events on Mars, then secondaries should be more abundant than primaries at diameters a factor of ∼1000 smaller than that of the largest primary crater that contributed secondaries. As a result, most small craters on Mars could be secondaries. Depth/diameter ratios of 1300 small craters (10–500 m diameter) in Isidis Planitia and Gusev crater have a mean value of 0.08; the freshest of these craters give a ratio of 0.11, identical to that of fresh secondary craters on the Moon (Pike and Wilhelms, 1978, Secondary-impact craters on the Moon: topographic form and geologic process, Lunar Planet. Sci. IX, 907–909) and significantly less than the value of ∼0.2 or more expected for fresh primary craters of this size range. Several observations suggest that the production functions of Hartmann and Neukum (2001, Cratering chronology and the evolution of Mars, Space Sci. Rev. 96, 165–194) predict too many primary craters smaller than a few hundred meters in diameter. Fewer small, high-velocity impacts may explain why there appears to be little impact regolith over Amazonian terrains. Martian terrains dated by small craters could be older than reported in recent publications.

Icarus↗
USGS2007

Saturn's dynamic D ring

The Cassini spacecraft has provided the first clear images of the D ring since the Voyager missions. These observations show that the structure of the D ring has undergone significant changes over the last 25 years. The brightest of the three ringlets seen in the Voyager images (named D72), has transformed from a narrow, <40-km wide ringlet to a much broader and more diffuse 250-km wide feature. In addition, its center of light has shifted inwards by over 200 km relative to other features in the D ring. Cassini also finds that the locations of other narrow features in the D ring and the structure of the diffuse material in the D ring differ from those measured by Voyager. Furthermore, Cassini has detected additional ringlets and structures in the D ring that were not observed by Voyager. These include a sheet of material just interior to the inner edge of the C ring that is only observable at phase angles below about 60°. New photometric and spectroscopic data from the ISS (Imaging Science Subsystem) and VIMS (Visual and Infrared Mapping Spectrometer) instruments onboard Cassini show the D ring contains a variety of different particle populations with typical particle sizes ranging from 1 to 100 microns. High-resolution images reveal fine-scale structures in the D ring that appear to be variable in time and/or longitude. Particularly interesting is a remarkably regular, periodic structure with a wavelength of ~30 km extending between orbital radii of 73,200 and 74,000 km. A similar structure was previously observed in 1995 during the occultation of the star GSC5249-01240, at which time it had a wavelength of ~60 km. We interpret this structure as a periodic vertical corrugation in the D ring produced by differential nodal regression of an initially inclined ring. We speculate that this structure may have formed in response to an impact with a comet or meteoroid in early 1984.

Icarus↗
USGS2014

Mars Global Digital Dune Database (MGD 3 ): Global dune distribution and wind pattern observations

The Mars Global Digital Dune Database (MGD 3 ) is complete and now extends from 90°N to 90°S latitude. The recently released south pole (SP) portion (MC-30) of MGD 3 adds ∼60,000 km 2 of medium to large-size dark dune fields and ∼15,000 km 2 of sand deposits and smaller dune fields to the previously released equatorial (EQ, ∼70,000 km 2 ), and north pole (NP, ∼845,000 km 2 ) portions of the database, bringing the global total to ∼975,000 km2. Nearly all NP dunes are part of large sand seas, while the majority of EQ and SP dune fields are individual dune fields located in craters. Despite the differences between Mars and Earth, their dune and dune field morphologies are strikingly similar. Bullseye dune fields, named for their concentric ring pattern, are the exception, possibly owing their distinctive appearance to winds that are unique to the crater environment. Ground-based wind directions are derived from slipface (SF) orientation and dune centroid azimuth (DCA), a measure of the relative location of a dune field inside a crater. SF and DCA often preserve evidence of different wind directions, suggesting the importance of local, topographically influenced winds. In general however, ground-based wind directions are broadly consistent with expected global patterns, such as polar easterlies. Intriguingly, between 40°S and 80°S latitude both SF and DCA preserve their strongest, though different, dominant wind direction, with transport toward the west and east for SF-derived winds and toward the north and west for DCA-derived winds.

Icarus↗
USGS2013

Crater-based dating of geological units on Mars: methods and application for the new global geological map

The new, post-Viking generation of Mars orbital imaging and topographical data provide significant higher-resolution details of surface morphologies, which induced a new effort to photo-geologically map the surface of Mars at 1:20,000,000 scale. Although from unit superposition relations a relative stratigraphical framework can be compiled, it was the ambition of this mapping project to provide absolute unit age constraints through crater statistics. In this study, the crater counting method is described in detail, starting with the selection of image data, type locations (both from the mapper’s and crater counter’s perspectives) and the identification of impact craters. We describe the criteria used to validate and analyse measured crater populations, and to derive and interpret crater model ages. We provide examples of how geological information about the unit’s resurfacing history can be retrieved from crater size–frequency distributions. Three cases illustrate short-, intermediate, and long-term resurfacing histories. In addition, we introduce an interpretation-independent visualisation of the crater resurfacing history that uses the reduction of the crater population in a given size range relative to the expected population given the observed crater density at larger sizes. From a set of potential type locations, 48 areas from 22 globally mapped units were deemed suitable for crater counting. Because resurfacing ages were derived from crater statistics, these secondary ages were used to define the unit age rather than the base age. Using the methods described herein, we modelled ages that are consistent with the interpreted stratigraphy. Our derived model ages allow age assignments to be included in unit names. We discuss the limitations of using the crater dating technique for global-scale geological mapping. Finally, we present recommendations for the documentation and presentation of crater statistics in publications.

Icarus↗
USGS2018

A new generation of the United States National Land Cover Database: Requirements, research priorities, design, and implementation strategies

The U.S. Geological Survey (USGS), in partnership with several federal agencies, has developed and released four National Land Cover Database (NLCD) products over the past two decades: NLCD 1992, 2001, 2006, and 2011. These products provide spatially explicit and reliable information on the Nation’s land cover and land cover change. To continue the legacy of NLCD and further establish a long-term monitoring capability for the Nation’s land resources, the USGS has designed a new generation of NLCD products named NLCD 2016. The NLCD 2016 design aims to provide innovative, consistent, and robust methodologies for production of a multi-temporal land cover and land cover change database from 2001 to 2016 at 2–3-year intervals. Comprehensive research was conducted and resulted in developed strategies for NLCD 2016: a streamlined process for assembling and preprocessing Landsat imagery and geospatial ancillary datasets; a multi-source integrated training data development and decision-tree based land cover classifications; a temporally, spectrally, and spatially integrated land cover change analysis strategy; a hierarchical theme-based post-classification and integration protocol for generating land cover and change products; a continuous fields biophysical parameters modeling method; and an automated scripted operational system for the NLCD 2016 production. The performance of the developed strategies and methods were tested in twenty World Reference System-2 path/row throughout the conterminous U.S. An overall agreement ranging from 71% to 97% between land cover classification and reference data was achieved for all tested area and all years. Results from this study confirm the robustness of this comprehensive and highly automated procedure for NLCD 2016 operational mapping.

Conterminous United States↗
USGS2022

A novel regression method for harmonic analysis of time series

Harmonic analysis of time series is an important technique in remote sensing to reveal seasonal land surface dynamics. However, frequency selection in the harmonic analysis is often difficult because high-frequency components are useful for delineating seasonal dynamics but sensitive to noise and gaps in time series. On the other hand, it is challenging to obtain temporally continuous satellite data with high quality because of atmospheric contamination. We developed a novel regression method named Harmonic Adaptive Penalty Operator (HAPO) for harmonic analysis of unevenly distributed time series. We introduced a new penalty function to minimize unexpected fluctuations in the model, which can substantially reduce the overfitting issue of regression in time series with temporal gaps. Specifically, the new penalty function minimizes the length of the model curve and the value range difference between the model and the time series observations. We compared HAPO with three widely used regression methods (OLS: Ordinary Least Squares; LASSO: Least Absolute Shrinkage and Selection Operator; and Ridge) in different scenarios using Landsat time series data across the United States. First, we evaluated methods using the Landsat surface reflectance time series within a single year. HAPO showed low and consistent monthly Root Mean Square Deviation (RMSD) values, in which most of the time RMSD of predicted reflectance were less than 0.04. More importantly, HAPO showed consistent and less bias given varying density and irregularity of time series. Second, we evaluated methods using multi-year time series. HAPO was a better predictor of relatively short time series (< 4 years) with steady low RMSD values. When a longer time series ( 4 years) was used, all four methods showed similar RMSD values, but HAPO outperformed the other methods if there were temporal gaps. Therefore, for places with large seasonal observation gaps or for time series that are relatively short (less than 4 years), HAPO can provide more consistent and accurate results in harmonic analysis of time series.

ISPRS Journal of Photogrammetry and Remote Sensing↗
USGS2013

Bactericidal efficacy of elevated pH on fish pathogenic and environmental bacteria

Ship ballast water is a recognized medium for transfer and introductions of nonindigenous species. There is a need for new ballast water treatment methods that effectively and safely eliminate or greatly minimize movements of these species. The present study employed laboratory methods to evaluate the bactericidal efficacy of increased pH (pH 10.0–12.0) for exposure durations of up to 72 h to kill a variety of Gram-negative and Gram-positive bacteria including fish pathogens (Aeromonas spp., Yersinia ruckeri, Edwardsiella ictaluri, Serratia liquefaciens, Carnobacterium sp.), other common aquatic-inhabitant bacteria (Serratia marcescens, Pseudomonas fluorescens, Staphylococcus sp., Bacillus sp.) and indicators listed in International Maritime Organization D2 Standards; namely, Vibrio cholera (an environmental isolate from fish), Escherichia coli and Enterococcus faecalis. Volumes of 5 N NaOH were added to tryptic soy broth to obtain desired pH adjustments. Viable cells were determined after 0, 4, 12, 24, 48, and 72 h. Initial (0 h) cell numbers ranged from 3.40 × 10 4 cfu/mL for Bacillus sp. to 2.44 × 10 7 cfu/mL for E. faecalis. The effective endpoints of pH and treatment duration necessary to realize 100% bactericidal effect varied; however, all bacteria tested were killed within 72 h at pH 12.0 or lower. The lowest parameters examined, 4 h at pH 10.0, were bactericidal to V. cholera, E. ictaluri, three of four isolates of E. coli, and (three of four) Aeromonas salmonicida subsp. salmonicida. Bactericidal effect was attained at pH 10.0 within 12 h for the other A. salmonicida subsp. salmonicida, and within 24 h for P. fluorescens, and the remaining E. coli.

Journal of Advanced Research↗
USGS2008

Biomedical and veterinary science can increase our understanding of coral disease

A balanced approach to coral disease investigation is critical for understanding the global decline of corals. Such an approach should involve the proper use of biomedical concepts, tools, and terminology to address confusion and promote clarity in the coral disease literature. Investigating disease in corals should follow a logical series of steps including identification of disease, systematic morphologic descriptions of lesions at the gross and cellular levels, measurement of health indices, and experiments to understand disease pathogenesis and the complex interactions between host, pathogen, and the environment. This model for disease investigation is widely accepted in the medical, veterinary and invertebrate pathology disciplines. We present standard biomedical rationale behind the detection, description, and naming of diseases and offer examples of the application of Koch's postulates to elucidate the etiology of some infectious diseases. Basic epidemiologic concepts are introduced to help investigators think systematically about the cause(s) of complex diseases. A major goal of disease investigation in corals and other organisms is to gather data that will enable the establishment of standardized case definitions to distinguish among diseases. Concepts and facts amassed from empirical studies over the centuries by medical and veterinary pathologists have standardized disease investigation and are invaluable to coral researchers because of the robust comparisons they enable; examples of these are given throughout this paper. Arguments over whether coral diseases are caused by primary versus opportunistic pathogens reflect the lack of data available to prove or refute such hypotheses and emphasize the need for coral disease investigations that focus on: characterizing the normal microbiota and physiology of the healthy host; defining ecological interactions within the microbial community associated with the host; and investigating host immunity, host-agent interactions, pathology, pathogenesis, and factors that promote the pathogenicity of the causative agent(s) of disease.

Journal of Experimental Marine Biology and Ecology↗
USGS2016

Eurytemora carolleeae in the Laurentian Great Lakes revealed by phylogenetic and morphological analysis

In the Laurentian Great Lakes, specimens of Eurytemora have been reported as Eurytemora affinis since its invasion in the late 1950s. During an intensive collection of aquatic invertebrates for morphological and molecular identification in Western Lake Erie in 2012-2013, several specimens of Eurytemora were collected. Analysis of these specimens identified them as the recently described species Eurytemora carolleeae Alekseev and Souissi 2011. This result led us to assess E. carolleeae’s identifying features, geographic distribution and historical presence in the Laurentian Great Lakes in view of its recent description in 2011. Cytochrome oxidase I (COI) DNA sequences of Eurytemora specimens were identified as closer (2 - 4% different) to recently described E. carolleeae than to most E. affinis sequences (14% different). Eurytemora from other areas of the Great Lakes and from North American rivers as far west as South Dakota (Missouri River) and east to Delaware (Christina River) also keyed to E. carolleeae . Morphological analysis of archival specimens from 1962 and from all the Great Lakes was identified as E. carolleeae . Additionally, Eurytemora drawings in previous publications were reassessed to determine if the species was E. carolleeae and are reported here. Additional morphological characters that may distinguish North American E. carolleeae from other taxa are also described. We conclude that E. carolleeae is the correct name for the species of Eurytemora that has inhabited the Great Lakes since its invasion, as established by both morphological and COI sequence comparisons to reference keys and sequence databases in present and archival specimens.

Great Lakes↗
USGS2025

Visioning and conceptual framework for coordinating Great Lakes connecting waters research and monitoring

The Laurentian Great Lakes are connected via naturally occurring straits and rivers: St. Marys River, Straits of Mackinac, St. Clair-Detroit River System, Niagara River, and the St. Lawrence River. Despite the historical ecological and economic importance of these waters, international agreements (e.g., Great Lakes Water Quality Agreement) only recently explicitly named the Great Lakes Connecting Waters (GLCWs), requiring governments to address the challenges of adequate restoration and protection from historical use and degradation. Standardized research and monitoring activities are needed; however, there is no established mechanism for coordination across the GLCWs. A three-day summit in 2023 convened experts to form the initial framework for a GLCWs Collaborative to increase standardizations and knowledge transfer. Participants drafted a governance structure and priorities following the principles of collective impact but allowed for place-based specificity for local connecting water organizations. Priorities and suggestions for success included: 1) co-development of the collaborative with all rights holders, stakeholders, and surrounding communities; 2) investment in research and technology specific to GLCWs; 3) investment in information transfer and training; 4) increased communication; and 5) better integration into existing Great Lakes research, monitoring, and funding programs. Expanding participation in all local GLCWs organizations with principles of inclusivity was identified as a larger collaborative goal. Next steps in the development of a GLCWs Collaborative include increased communication and formation of working groups and obtaining funding for a dedicated organization to begin supporting activities (communication, facilitation, logistics). By using a deliberate process for establishment, the potentially slower time frame for establishment may result in increased participation and success.

Laurentian Great Lakes↗
USGS2026

Predatory impact of freshwater drum on dreissenid populations in western Lake Erie

Impacts of dreissenid mussels ( Dreissena spp.) on Great Lakes ecosystems are well documented, but the mechanisms driving variation in their abundance remain poorly understood. Dreissenid mussels have been incorporated into fish diets throughout the Great Lakes; however, studies quantifying the amount of dreissenid mussels consumed by fish predators are limited. To date, attention has mainly focused on invasive round goby ( Neogobius melanostomus ) predation of dreissenid mussels. Biomass of native molluscivores, namely the freshwater drum ( Aplodinotus grunniens ), may exceed round goby biomass by an order of magnitude in some areas. Thus, the role of predation on dreissenid mussel population dynamics may be greater than currently assumed. Here, we combine estimates of diet composition and fish biomass to estimate kg/ha of dreissenid mussels consumed by freshwater drum in the West Basin of Lake Erie. Annual consumption estimates of dreissenids by freshwater drum were large (averaging 23.79 kg/ha shell-free mass), and generally exceeded existing dreissenid consumption estimates for round goby. Our results support evaluation of ecological mechanisms, such as predation, to improve our knowledge of factors that may influence dreissenid mussel abundance.

western Lake Erie↗
USGS2026

Evolution of arsenic speciation during thermal treatment simulating wildfire heating in arsenic-rich sediments

Understanding arsenic transformations during wildfires can help better constrain the environmental impacts of increased wildfire intensity, frequency, and burned area. Previous studies have monitored the evolution of metal(loid) speciation at specific time/temperature endpoints, hampering the comprehension of the dynamic transformations of metal(loid)s during wildfires. Here, a novel approach, namely in situ time-resolved X-ray absorption near edge structure (TR-XANES) spectroscopy, was used for the first time to monitor the heat-induced redox transformations of arsenic (As) in two As-rich sediments with different characteristics. No arsenic transformations were observed in either sediment at room temperature or when the sediments were heated to a maximum temperature of 400 ºC. When heated to a maximum temperature of 600 or 700 ºC, As underwent complex, dynamic, and partially reversible redox transformations, the extent of which varied with sediment properties, initial As speciation, and heating temperature and duration. In the case of sediments initially containing As(V) and a low sulfide/sulfate ratio, the As(V) was reduced to As(III) over a short period of time, followed by immediate reoxidation of As(III) to As(V). The extent of the transient As reduction increased with increasing maximum temperature. The final As speciation at the end of the heating experiment was not substantially different from that at the beginning of the experiment. In contrast, in sediments that were best described by a more complex combination of reference compounds (As(V) coprecipitated with or adsorbed on iron oxide, arseniosiderite, arsenopyrite) with a high sulfide/sulfate ratio, As(V) phases were reduced to As(III) and As(III) was further reduced to As(-I) (at 700 ºC, only), followed by delayed reoxidation of As(-I) to As(III) (at 700 ºC, only) and reoxidation of As(III) to As(V). Approximately 9% and 26% of all As remained as As(III) ( i.e., was not reoxidized to As(V)) at the end of the experiment when heated to a maximum temperature of 600 and 700 ºC, respectively, which is higher than the relative abundance of As(III) in the native sediment. These fire-induced As transformations can potentially enhance As mobilization to surface water and groundwater, posing risks to environmental and human health.

California↗
USGS2005

Repeated surveys by acoustic Doppler current profiler for flow and sediment dynamics in a tidal river

A strategy of repeated surveys by acoustic Doppler current profiler (ADCP) was applied in a tidal river to map velocity vectors and suspended-sediment indicators. The Sacramento River at the junction with the Delta Cross Channel at Walnut Grove, California, was surveyed over several tidal cycles in the Fall of 2000 and 2001 with a vessel-mounted ADCP. Velocity profiles were recorded along flow-defining survey paths, with surveys repeated every 27 min through a diurnal tidal cycle. Velocity vectors along each survey path were interpolated to a three-dimensional Cartesian grid that conformed to local bathymetry. A separate array of vectors was interpolated onto a grid from each survey. By displaying interpolated vector grids sequentially with computer animation, flow dynamics of the reach could be studied in three-dimensions as flow responded to the tidal cycle. Velocity streamtraces in the grid showed the upwelling of flow from the bottom of the Sacramento River channel into the Delta Cross Channel. The sequential display of vector grids showed that water in the canal briefly returned into the Sacramento River after peak flood tides, which had not been known previously. In addition to velocity vectors, ADCP data were processed to derive channel bathymetry and a spatial indicator for suspended-sediment concentration. Individual beam distances to bed, recorded by the ADCP, were transformed to yield bathymetry accurate enough to resolve small bedforms within the study reach. While recording velocity, ADCPs also record the intensity of acoustic backscatter from particles suspended in the flow. Sequential surveys of backscatter intensity were interpolated to grids and animated to indicate the spatial movement of suspended sediment through the study reach. Calculation of backscatter flux through cross-sectional grids provided a first step for computation of suspended-sediment discharge, the second step being a calibrated relation between backscatter intensity and sediment concentration. Spatial analyses of ADCP data showed that a strategy of repeated surveys and flow-field interpolation has the potential to simplify computation of flow and sediment discharge through complex waterways. The use of trade, product, industry, or firm names in this report is for descriptive purposes only and does not constitute endorsement of products by the US Government. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology↗
USGS2010

Comment on “Two statistics for evaluating parameter identifiability and error reduction” by John Doherty and Randall J. Hunt

Doherty and Hunt (2009) present important ideas for first-order-second moment sensitivity analysis, but five issues are discussed in this comment. First, considering the composite-scaled sensitivity (CSS) jointly with parameter correlation coefficients (PCC) in a CSS/PCC analysis addresses the difficulties with CSS mentioned in the introduction. Second, their new parameter identifiability statistic actually is likely to do a poor job of parameter identifiability in common situations. The statistic instead performs the very useful role of showing how model parameters are included in the estimated singular value decomposition (SVD) parameters. Its close relation to CSS is shown. Third, the idea from p. 125 that a suitable truncation point for SVD parameters can be identified using the prediction variance is challenged using results from Moore and Doherty (2005) . Fourth, the relative error reduction statistic of Doherty and Hunt is shown to belong to an emerging set of statistics here named perturbed calculated variance statistics. Finally, the perturbed calculated variance statistics OPR and PPR mentioned on p. 121 are shown to explicitly include the parameter null-space component of uncertainty. Indeed, OPR and PPR results that account for null-space uncertainty have appeared in the literature since 2000.

Journal of Hydrology↗
USGS2013

Modeling unsaturated zone flow and runoff processes by integrating MODFLOW-LGR and VSF, and creating the new CFL package

In this paper two modifications to the MODFLOW code are presented. One concerns an extension of Local Grid Refinement (LGR) to Variable Saturated Flow process (VSF) capability. This modification allows the user to solve the 3D Richards’ equation only in selected parts of the model domain. The second modification introduces a new package, named CFL (Cascading Flow), which improves the computation of overland flow when ground surface saturation is simulated using either VSF or the Unsaturated Zone Flow (UZF) package. The modeling concepts are presented and demonstrated. Programmer documentation is included in appendices.

Journal of Hydrology↗
USGS2013

Appraising options to reduce shallow groundwater tables and enhance flow conditions over regional scales in an irrigated alluvial aquifer system

Some of the world’s key agricultural production systems face big challenges to both water quantity and quality due to shallow groundwater that results from long-term intensive irrigation, namely waterlogging and salinity, water losses, and environmental problems. This paper focuses on water quantity issues, presenting finite-difference groundwater models developed to describe shallow water table levels, non-beneficial groundwater consumptive use, and return flows to streams across two regions within an irrigated alluvial river valley in southeastern Colorado, USA. The models are calibrated and applied to simulate current baseline conditions in the alluvial aquifer system and to examine actions for potentially improving these conditions. The models provide a detailed description of regional-scale subsurface unsaturated and saturated flow processes, thereby enabling detailed spatiotemporal description of groundwater levels, recharge to infiltration ratios, partitioning of ET originating from the unsaturated and saturated zones, and groundwater flows, among other variables. Hybrid automated and manual calibration of the models is achieved using extensive observations of groundwater hydraulic head, groundwater return flow to streams, aquifer stratigraphy, canal seepage, total evapotranspiration, the portion of evapotranspiration supplied by upflux from the shallow water table, and irrigation flows. Baseline results from the two regional-scale models are compared to model predictions under variations of four alternative management schemes: (1) reduced seepage from earthen canals, (2) reduced irrigation applications, (3) rotational lease fallowing (irrigation water leased to municipalities, resulting in temporary dry-up of fields), and (4) combinations of these. The potential for increasing the average water table depth by up to 1.1 and 0.7 m in the two respective modeled regions, thereby reducing the threat of waterlogging and lowering non-beneficial consumptive use from adjacent fallow and naturally-vegetated lands, is demonstrated for the alternative management intervention scenarios considered. Net annual average savings of up to about 9.9 million m 3 (8000 ac ft) and 2.3 million m 3 (1900 ac ft) of non-beneficial groundwater consumptive use is demonstrated for the study periods in each of the two respective study regions. Alternative water management interventions achieve varying degrees of benefits in each of the two regions, suggesting a need to adopt region-specific interventions and avoid a ‘one-size-fits-all’ approach. Impacts of the considered interventions on return flows to the river were predicted to be significant, highlighting the need for flow augmentation to comply with an interstate river compact and portending beneficial impacts on solute loading.

Colorado↗
USGS2018

Satellite remote sensing estimation of river discharge: Application to the Yukon River Alaska

A methodology based on general hydraulic relations for rivers has been developed to estimate the discharge (flow rate) of rivers using satellite remote sensing observations. The estimates of discharge, flow depth, and flow velocity are derived from remotely observed water surface area, water surface slope, and water surface height, and demonstrated for two reaches of the Yukon River in Alaska, at Eagle (reach length 34.7 km) and near Stevens Village (reach length 38.3 km). The method is based on fundamental equations of hydraulic flow resistance in rivers, including the Manning equation and the Prandtl-von Karman universal velocity distribution equation. The method employs some new hydraulic relations to help define flow resistance and height of the zero flow boundary in the channel. Estimates are made both with and without calibration. The water surface area of the river reach is measured by using a provisional version of the U.S. Geological Survey (USGS) Landsat based product named Dynamic Surface Water Extent (DSWE). The water surface height and slope measurements require a self-consistent datum, and are derived from observations from the Jason-2 satellite altimeter mission. At both reach locations, the Jason-2 radar altimeter non-winter heights consistently tracked the stage recorded at USGS streamgages with a standard deviation of differences (error) during the non-winter periods of less than 7%. Part of the error may be due to differences in the gage and altimeter crossing locations with respect to the range of stage change and the response to changes in discharge at the upstream and downstream locations. For the non-winter periods, the radar derived slope estimates (mean = 0.0003) were constant over the mission lifetime, and in agreement with previously measured USGS water surface slopes and slopes determined from USGS topographic maps. The accuracy of the mean of the uncalibrated daily estimates of discharge varied between reaches, ranging from 13% near Stevens Village (N = 90) to −21% at Eagle (N = 246) based on the absolute error, and 5% to −6% based on the error of the log of the estimates. Calibrating to the mean of USGS daily discharge estimates from the streamflow rating for the same period of record at each streamgage resulted in mean absolute errors ranging from 1% to 2%, and log errors ranging from 1% or less. The error pattern of the estimates shows that without calibration, even though the mean is well simulated, the high and low end values over the range of estimates may have significant bias.

Alaska↗