USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computer contribution”

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 253 records · Page 14Linked to original sources
USGS2015

Methods for estimating the magnitude and frequency of peak streamflows at ungaged sites in and near the Oklahoma Panhandle

This report presents the results of a cooperative study by the U.S. Geological Survey and the Oklahoma Department of Transportation to estimate the magnitude and frequency of peak streamflows from regional regression equations for ungaged stream sites in and near the Oklahoma Panhandle. These methods relate basin characteristics (physiographic and climatic attributes) to selected peak streamflow frequency statistics with the 50-, 20-, 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probabilities. These relations were developed based on data from 32 selected streamflow-gaging stations in the Oklahoma Panhandle and in neighboring parts of Colorado, Kansas, New Mexico, and Texas. The basin characteristics for the selected streamflow-gaging stations were determined by using a geographic information system and the Oklahoma StreamStats application. Peak-streamflow frequency statistics were computed from annual peak-streamflow records from the irrigated period of record from water year 1978 through water year 2014. Generalized-least-squares multiple-linear regression analysis was used to formulate regression relations between peak-streamflow frequency statistics and basin characteristics. Contributing drainage area was the only basin characteristic determined to be statistically significant for all percentage of annual exceedance probabilities and was the only basin characteristic used in regional regression equations for estimating peak-streamflow frequency statistics on unregulated streams in and near the Oklahoma Panhandle. The regression model pseudo-coefficient of determination, converted to percent, for the Oklahoma Panhandle regional regression equations ranged from about 38 to 63 percent. The standard errors of prediction and the standard model errors for the Oklahoma Panhandle regional regression equations ranged from about 84 to 148 percent and from about 76 to 138 percent, respectively. These errors were comparable to those reported for regional peak-streamflow frequency regression equations for the High Plains areas of Texas and Colorado. The root mean square errors for the Oklahoma Panhandle regional regression equations (ranging from 3,170 to 92,000 cubic feet per second) were less than the root mean square errors for the Oklahoma statewide regression equations (ranging from 18,900 to 412,000 cubic feet per second); therefore, the Oklahoma Panhandle regional regression equations produce more accurate peak-streamflow statistic estimates for the irrigated period of record in the Oklahoma Panhandle than do the Oklahoma statewide regression equations. The regression equations developed in this report are applicable to streams that are not substantially affected by regulation, impoundment, or surface-water withdrawals. These regression equations are intended for use for stream sites with contributing drainage areas less than or equal to about 2,060 square miles, the maximum value for the independent variable used in the regression analysis.

Colorado, Kansas, New Mexico, Oklahoma, Texas↗
USGS2000

Factors affecting nutrient trends in major rivers of the Chesapeake Bay Watershed

Trends in nutrient loads and flow-adjusted concentrations in the major rivers entering Chesapeake Bay were computed on the basis of water-quality data collected between 1985 and 1998 at 29 monitoring stations in the Susquehanna, Potomac, James, Rappahannock, York, Patuxent, and Choptank River Basins. Two computer models-the Chesapeake Bay Watershed Model (WSM) and the U.S. Geological Survey's 'Spatially Referenced Regressions on Watershed attributes' (SPARROW) Model-were used to help explain the major factors affecting the trends. Results from WSM simulations provided information on temporal changes in contributions from major nutrient sources, and results from SPARROW model simulations provided spatial detail on the distribution of nutrient yields in these basins. Additional data on nutrient sources, basin characteristics, implementation of management practices, and ground-water inputs to surface water were analyzed to help explain the trends. The major factors affecting the trends were changes in nutrient sources and natural variations in streamflow. The dominant source of nitrogen and phosphorus from 1985 to 1998 in six of the seven tributary basins to Chesapeake Bay was determined to be agriculture. Because of the predominance of agricultural inputs, changes in agricultural nutrient sources such as manure and fertilizer, combined with decreases in agricultural acreage and implementation of best management practices (BMPs), had the greatest impact on the trends in flow-adjusted nutrient concentrations. Urban acreage and population, however, were noted to be increasing throughout the Chesapeake Bay Watershed, and as a result, delivered loads of nutrients from urban areas increased during the study period. Overall, agricultural nutrient management, in combination with load decreases from point sources due to facility upgrades and the phosphate detergent ban, led to downward trends in flow-adjusted nutrient concentrations atmany of the monitoring stations in the watershed. The loads of nutrients, however, were not reduced significantly at most of the monitoring stations. This is due primarily to higher streamflow in the latter years of the monitoring period, which led to higher loading in those years. Results of this study indicate a need for more detailed information on BMP effectiveness under a full range of hydrologic conditions and in different areas of the watershed; an internally consistent fertilizer data set; greater consideration of the effects of watershed processes on nutrient transport; a refinement of current modeling efforts; and an expansion of the non-tidal monitoring network in the Chesapeake Bay Watershed.

Delaware, Maryland, New York, Pennsylvania, Virgin↗
USGS2002

Streamflow and nutrient data for the Yazoo River below Steele Bayou near Long Lake, Mississippi, 1996-2000

Increased nutrient loading to the Gulf of Mexico from off-continent flux has been identified as contributing to the increase in the areal extent of the low dissolved-oxygen zone that develops annually off the Louisiana and Texas coast. The proximity of the Yazoo River Basin in northwestern Mississippi to the Gulf of Mexico, and the intensive agricultural activities in the basin have led to speculation that the Yazoo River Basin contributes a disproportionate amount of nitrogen and phosphorus to the Mississippi River and ultimately to the Gulf of Mexico. An empirical measurement of the flux of nitrogen and phosphorus from the Yazoo Basin has not been possible due to the hydrology of the lower Yazoo River Basin. Streamflow for the Yazoo River below Steele Bayou is affected by backwater from the Mississippi River. Flow at the gage is non-uniform and varying, with bi-directional and reverse flows possible. Streamflow was computed by using remote sensing and acoustic and conventional discharge and velocity measurement techniques. Streamflow from the Yazoo River for the 1996-2000 period accounted for 2.8 percent of the flow of the Mississippi River for the same period. Water samples from the Yazoo River were collected from February 1996 through December 2000 and were analyzed for total nitrogen, nitrate, total phosphorus, and orthophosphorus as part of the U.S. Geological Survey National Water-Quality Assessment Program. These data were used to compute annual loads of nitrogen and phosphorus discharged from the Yazoo River for the period 1996-2000. Annual loads of nitrogen and phosphorus were calculated by two methods. The first method used multivariate regression and the second method multiplied the mean annual concentration by the total annual flow. Load estimates based on the product of the mean annual concentration and the total annual flow were within the 95 percent confidence interval for the load calculated by multivariate regression in 10 of 20 cases. The Yazoo River loads, compared to average annual loads in the Mississippi River, indicated that the Yazoo River was contributing 1.4 percent of the total nitrogen load, 0.7 percent of the nitrate load, 3.4 percent of the total phosphorus load, and 1.6 percent of the orthophosphorus load during 1996 - 2000. The total nitrogen, nitrate, and orthophosphorus loads in the Yazoo River Basin were less than expected, whereas the total phosphorus load was slightly higher than expected based on discharge.

Water-Resources Investigations Report↗
USGS2024

Skill assessment of a total water level and coastal change forecast during the landfall of a hurricane

The Total Water Level and Coastal Change Forecast (TWL&CC Forecast) provides coastal communities with 6-day notice of potential elevated water levels and coastal change (i.e., dune erosion, overwash, or inundation) on sandy beaches that threatens safety, infrastructure, or resources. This continuously operating model provides hourly information for select regions along U.S. Gulf of Mexico and Atlantic Ocean coastlines. The objective of this work is to assess the skill of forecasts during a period of elevated water levels along the coasts of North Carolina (NC) and South Carolina, USA caused by Hurricane Isaias in August 2020, using a combination of observations and model hindcasts. Water levels and waves were observed throughout the storm at three locations near Wrightsville Beach, NC, which provided information to assess forecast skill; a wave buoy offshore, a tide gage at a local pier, and a pressure sensor deployed at the pier. In addition to observations, the non-hydrostatic phase-resolving model SWASH (Simulating WAves till SHore) was forced with hourly wave energy spectra derived from a coupled Delft3D-SWAN simulation during the peak of Isaias, to complement observations by computing nearshore wave height and wave-induced setup and runup at the shoreline. During the storm peak, SWASH-simulated water levels at the sensor position were comparable to those at the maximum landward extent (bias = −0.05 m; gain = 0.26; r 2 = 0.99), suggesting that observations at the USGS sensor location were a useful proxy for total water level (TWL; sum of tide, surge and wave runup) at the shoreline that are predicted by the TWL&CC Forecast. The TWL forecast at Wrightsville Beach was consistent with observations from the USGS sensor (bias = −0.38 m and −0.74 m, scatter index = 0.22 and 0.28 for the two forecast model grids considered, respectively; weighted regression considering model uncertainty explained 95 percent of variability in observed TWL). Observed TWL was within the confidence interval of the TWL&CC Forecast for the 5 h at the storm peak. Forecast mean water levels (MWL; sum of tide, surge and wave setup) and tide gage observations were also consistent (bias = 0.07 m and 0.02 m for the forecast model grids; scatter index = 0.46; r 2 = 0.80). Forecast MWL at the storm peak was within 0.06 m of the observed MWL from the tide gage for both sites. In the region where Isaias made landfall, eight additional pressure sensors were compared to the peak TWL forecast (bias = 0.14 m; scatter index = 0.18). Forecast TWL explained 90 percent of observed variability in TWL when considering uncertainty of the forecast with a weighted regression. The results demonstrate that wave-driven water levels contributed a significant portion of the forecast TWL during Isaias (52 percent during the three peak hours of the storm), and that TWL were represented using the forecast model. Mean absolute error of the coastal change forecast and observed overwash is 0.4 and 0.14 for the two forecast model grids considered. The skill demonstrated by this computationally efficient method indicates that the forecasting system can provide fast and reliable predictions of TWL across hundreds of km of coastline at sub-km resolution, days to hours in advance of when storms threaten coastal regions.

Coastal Engineering↗
USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗
USGS2005

Effects of tag loss on direct estimates of population growth rate

The temporal symmetry approach of R. Pradel can be used with capture– recapture data to produce retrospective estimates of a population's growth rate, λ i , and the relative contributions to λ i from different components of the population. Direct estimation of λ i provides an alternative to using population projection matrices to estimate asymptotic λ and is seeing increased use. However, the robustness of direct estimates of λ i to violations of several key assumptions has not yet been investigated. Here, we consider tag loss as a possible source of bias for scenarios in which the rate of tag loss is (1) the same for all marked animals in the population and (2) a function of tag age. We computed analytic approximations of the expected values for each of the parameter estimators involved in direct estimation and used those values to calculate bias and precision for each parameter estimator. Estimates of λ i were robust to homogeneous rates of tag loss. When tag loss rates varied by tag age, bias occurred for some of the sampling situations evaluated, especially those with low capture probability, a high rate of tag loss, or both. For situations with low rates of tag loss and high capture probability, bias was low and often negligible. Estimates of contributions of demographic components to λ i were not robust to tag loss. Tag loss reduced the precision of all estimates because tag loss results in fewer marked animals remaining available for estimation. Clearly tag loss should be prevented if possible, and should be considered in analyses of λ i , but tag loss does not necessarily preclude unbiased estimation of λ i .

Ecology↗
USGS1990

Sediment movement along the U.S. east coast continental shelf-I. Estimates of bottom stress using the Grant-Madsen model and near-bottom wave and current measurements

Bottom stress is calculated for several long-term time-series observations, made on the U.S. east coast continental shelf during winter, using the wave-current interaction and moveable bed models of Grant and Madsen (1979, Journal of Geophysical Research, 84, 1797-1808; 1982, Journal of Geophysical Research, 87, 469-482). The wave and current measurements were obtained by means of a bottom tripod system which measured current using a Savonius rotor and vane and waves by means of a pressure sensor. The variables were burst sampled about 10% of the time. Wave energy was reasonably resolved, although aliased by wave groupiness, and wave period was accurate to 1-2 s during large storms. Errors in current speed and direction depend on the speed of the mean current relative to the wave current. In general, errors in bottom stress caused by uncertainties in measured current speed and wave characteristics were 10-20%. During storms, the bottom stress calculated using the Grant-Madsen models exceeded stress computed from conventional drag laws by a factor of about 1.5 on average and 3 or more during storm peaks. Thus, even in water as deep as 80 m, oscillatory near-bottom currents associated with surface gravity waves of period 12 s or longer will contribute substantially to bottom stress. Given that the Grant-Madsen model is correct, parameterizations of bottom stress that do not incorporate wave effects will substantially underestimate stress and sediment transport in this region of the continental shelf.

Georges Bank, Mid-Atlantic Bight↗
USGS2001

Nutrient and suspended-sediment concentrations and loads and benthic-invertebrate data for tributaries to the St. Croix River, Wisconsin and Minnesota, 1997-99

Nutrient and suspended-sediment data were collected on major tributaries to the St. Croix River during 1997-99 as part of three studies. The first study, done in 1997 as part of the U.S. Geological Survey?s National Water-Quality Assessment Program Upper Mississippi Study Unit, was a widespread synoptic survey of nutrient and suspended-sediment concentrations, loads, and yields during snowmelt. Runoff from snowmelt in agricultural areas and other areas with low permeability soils had significantly greater nutrient concentrations than forested areas, whereas differences in suspended-sediment loading were not detected. In 1998, synoptic samplings of 11 tributaries were done during snowmelt, base-flow, and storm-runoff periods. These studies showed that the Apple, Willow, and Kinnickinnic Rivers were major contributors of suspended sediments and nutrients to the St. Croix River during base flow and storm-runoff events. Nitrate concentrations were high during base flow in the agricultural tributaries?specifically, the Kinnickinnic (4.83 mg/L), Willow (1.53 mg/L), and Apple (0.79 mg/L) Rivers?possibly from ground-water recharge or point-source contributions. Extensive water-quality sampling was done monthly and during high-flow events in water year 1999 (October 1, 1998 to September 30, 1999) in coordination with continuous streamflow monitoring at 12 sites in the St. Croix River Basin. These data were used to compute annual nutrient and suspended-sediment loads and yields at the monitored sites for water year 1999. Relations among environmental characteristics and calculated annual nutrient and suspended-sediment yields were used to estimate loading from unmonitored parts of the basin. The environmental characteristics found to best estimate annual yields were soil characteristics (clay, permeability of soil, and erodibility), basin slope and area, and the percentages of wetland and urban areas in the basins. Variability in 1999 rainfall intensity resulted in annual yields from several northern, forested basins being higher than those from the southern, agricultural basins. The Sunrise River had the highest annual suspended-sediment and nutrient yields in the basin in 1999. Concentrations and instantaneous loading rates varied as much among various flow conditions at individual sites as among sites during the three years of study. Benthic invertebrates were sampled and indices of water quality were calculated at 16 tributaries in fall 1999. Benthic invertebrate indices indicated excellent to good water quality at all tributaries except Valley Creek, Willow River, and Kettle River. No relations were found between benthic invertebrate indices and the calculated and estimated 1999 annual tributary loads and yields.

Minnesota, Wisconsin↗
USGS1988

Simulation of three lahars in the Mount St Helens area, Washington using a one-dimensional, unsteady-state streamflow model

A one-dimensional, unsteady-state, open-channel model was used to analytically reproduce three lahar events. Factors contributing to the success of the modeling were: (1) the lahars were confined to a channel, (2) channel roughness was defined by field information, and (3) the volume of the flow remained relatively unchanged for the duration of the peak. Manning 's 'n ' values used in computing conveyance in the model were subject to the changing rheology of the debris flow and were calculated from field cross-section information (velocities used in these calculations were derived from super-elevation or run-up formulas). For the events modeled in this exercise, Manning 's 'n ' calculations ranged from 0.020 to 0.099. In all lahar simulations, the rheology of the flow changed in a downstream direction during the course of the event. Chen's 'U ', the mudflow consistency index, changed approximately an order of magnitude for each event. The ' u ' values ranged from 5-2,260 kg/m for three events modeled. The empirical approach adopted in this paper is useful as a tool to help predict debris-flow behavior, but does not lead to understanding the physical processes of debris flows.

Washington↗
USGS2011

Changes in low-flow frequency from 1976-2006 at selected streamgages in New York, excluding Long Island

Many Federal, State, and local agencies use low-flow data to establish water-use policy and help determine the total maximum daily loads and effluent limits of point and nonpoint sources of contamination of surface water during periods of decreased streamflow. Low-flow magnitude and frequency are used often by water-supply planners, reservoir managers, and hydroelectric facilities to manage water availability for supply and power generation. Low-flow statistics for eight selected U.S. Geological Survey streamgages in New York State were calculated for the period from 1976 through 2006 and for the entire period of continuous streamflow record. The 7-day, 2-year and 10-year low flows were computed and compared with those low flows published in the1979 U.S. Geological Survey report, Low-flow frequency analysis of streams in New York, Bulletin 74. Observed changes in low-flow frequency at each gage were then examined and compared to changes in precipitation and land use to determine whether a relation between similar patterns could be identified. A statewide U.S. Geological Survey study has not been done to develop equations for estimating low flows on rural unregulated streams in New York. Currently (2010) only one regional study developed for parts of the lower Hudson River Basin in 1986 is available to assist in estimating low flows on rural streams with unregulated streamflow in New York. Low-flow statistics published in the 1979 report need to be updated by using additional data collected since 1976 to determine current low-flow conditions across New York State. At-site low-flow statistics were updated for eight streamgages in New York by using continuous daily streamflow data through 2006 for the future development of a statewide research study. Selection of the eight streamgages used in this study identified a major deficiency in the number of available unregulated long-term U.S. Geological Survey streamgages needed for the development of regional low-flow equations in New York. A limited analysis of the changes in land use for the contributing drainage areas for each streamgage, changes in precipitation, and trends in the annual 7-day minimum flow also are presented. The 7-day, 2-year low flow showed increases of 14 to 35 percent and the 7-day 10-year low flow showed zero to 19 percent increases at rural streamgages with unregulated streamflows when statistics were computed by using data from 1976 through 2006 and compared with published data in Bulletin 74. When the entire period of record was used to compute low flow frequencies, the 7-day, 2-year low flows increased from about 6 to 15 percent whereas the 7-day 10-year low flows showed zero to 5 percent increases. Streamgages affected by urbanization and regulation for water supply showed the most significant changes in the 7-day, 2-year and 10-year low-flow frequencies. These streamgages are included to help identify the effects of urbanization and regulation on streamflow at these locations. The 7-day 10-year low flow increased by 65 percent at the U.S. Geological Survey streamgage Hackensack River at West Nyack, N.Y., and increased 120 percent at the U.S. Geological Survey streamgage Neversink River at Godeffroy, N.Y., when statistics were computed by using data from 1976 through 2006 and compared with the statistics for the regulated period computed in Bulletin 74.

New York↗
USGS2000

Creating a standardized watersheds database for the lower Rio Grande/Rio Bravo, Texas

This report describes the creation of a large-scale watershed database for the lower Rio Grande/Rio Bravo Basin in Texas. The watershed database includes watersheds delineated to all 1:24,000-scale mapped stream confluences and other hydrologically significant points, selected watershed characteristics, and hydrologic derivative datasets. Computer technology allows generation of preliminary watershed boundaries in a fraction of the time needed for manual methods. This automated process reduces development time and results in quality improvements in watershed boundaries and characteristics. These data can then be compiled in a permanent database, eliminating the time-consuming step of data creation at the beginning of a project and providing a stable base dataset that can give users greater confidence when further subdividing watersheds. A standardized dataset of watershed characteristics is a valuable contribution to the understanding and management of natural resources. Vertical integration of the input datasets used to automatically generate watershed boundaries is crucial to the success of such an effort. The optimum situation would be to use the digital orthophoto quadrangles as the source of all the input datasets. While the hydrographic data from the digital line graphs can be revised to match the digital orthophoto quadrangles, hypsography data cannot be revised to match the digital orthophoto quadrangles. Revised hydrography from the digital orthophoto quadrangle should be used to create an updated digital elevation model that incorporates the stream channels as revised from the digital orthophoto quadrangle. Computer-generated, standardized watersheds that are vertically integrated with existing digital line graph hydrographic data will continue to be difficult to create until revisions can be made to existing source datasets. Until such time, manual editing will be necessary to make adjustments for man-made features and changes in the natural landscape that are not reflected in the digital elevation model data.

Texas↗
USGS2004

Contribution of natural history collection data to biodiversity assessment in national parks

There has been mounting interest in the use of museum and herbaria collections to assess biodiversity; information is often difficult to locate and access, however, and few recommendations are available for effectively using natural history collections. As part of an effort to inventory vertebrates and vascular plants in U.S. national parks, we searched manually and by computer for specimens originating within or adjacent to 14 parks throughout the northeastern United States. We compared the number of specimens located to collection size to determine whether there was any effect on detection rate of specimens. We evaluated the importance of park characteristics (e.g., age since establishment, size, theme [natural vs. cultural]) for influencing the number of specimens found in a collection. We located >31,000 specimens and compiled associated records (hereafter referred to as specimens) from 78 collections; >9000 specimens were park-significant, originating either within park boundaries or in the local township where the park was located. We found >2000 specimens by means of manual searches, which cost $0.001–0.15 per specimen searched and $0.81–151.95 per specimen found. Collection effort appeared relatively uniform between 1890 and 1980, with low periods corresponding to significant sociopolitical events. Detection rates for specimens were inversely related to collection size. Although specimens were most often located in collections within the region of interest, specimens can be found anywhere, particularly in large collections international in scope, suggesting that global searches will be necessary to evaluate historical biodiversity. Park characteristics indicated that more collecting effort occurred within or adjacent to larger parks established for natural resources than in smaller historical sites. Because many institutions have not yet established electronic databases for collections, manual searches can be useful for retrieving specimens. Our results show that thorough, systematic searching of natural history collections for park-significant specimens can provide a historical perspective on biodiversity for park managers.

Conservation Biology↗
USGS2009

Identifying sources of dissolved organic carbon in agriculturally dominated rivers using radiocarbon age dating: Sacramento-San Joaquin River Basin, California

We used radiocarbon measurements of dissolved organic carbon (DOC) to resolve sources of riverine carbon within agriculturally dominated landscapes in California. During 2003 and 2004, average Δ 14 C for DOC was −254‰ in agricultural drains in the Sacramento–San Joaquin Delta, −218‰ in the San Joaquin River, −175‰ in the California State Water Project and −152‰ in the Sacramento River. The age of bulk DOC transiting the rivers of California’s Central Valley is the oldest reported for large rivers and suggests wide-spread loss of soil organic matter caused by agriculture and urbanization. Using DAX 8 adsorbent, we isolated and measured 14 C concentrations in hydrophobic acid fractions (HPOA); river samples showed evidence of bomb-pulse carbon with average Δ 14 C of 91 and 76‰ for the San Joaquin and Sacramento Rivers, respectively, with older HPOA, −204‰, observed in agricultural drains. An operationally defined non-HPOA fraction of DOC was observed in the San Joaquin River with seasonally computed Δ 14 C values of between −275 and −687‰; the source of this aged material was hypothesized to be physically protected organic-matter in high clay-content soils and agrochemicals (i.e., radiocarbon-dead material) applied to farmlands. Mixing models suggest that the Sacramento River contributes about 50% of the DOC load in the California State Water Project, and agricultural drains contribute approximately one-third of the load. In contrast to studies showing stabilization of soil carbon pools within one or two decades following land conversion, sustained loss of soil organic matter, occurring many decades after the initial agricultural-land conversion, was observed in California’s Central Valley.

California↗
USGS2019

A review of machine learning applications to coastal sediment transport and morphodynamics

A range of computer science methods under the heading of machine learning (ML) enables the extraction of insight and quantitative relationships from multidimensional datasets. Here, we review some common ML methods and their application to studies of coastal morphodynamics and sediment transport. We examine aspects of ‘what’ and ‘why’ ML methods contribute, such as ‘what’ science problems ML tools have been used to address, ‘what’ was learned when using ML, and ‘why’ authors used ML methods. We find a variety of research questions have been addressed, ranging from small-scale predictions of sediment transport to larger-scale sand bar morphodynamics and coastal overwash on a developed island. We find various reasons justify the use of ML, including maximize predictability, emulation of model components, smooth and continuous nonlinear regression through data, and explicit inclusion of uncertainty. Overall the expanding use of ML has allowed for an expanding set of questions to be addressed. After reviewing the studies we outline a set of ‘best practices’ for coastal researchers using machine learning methods. Finally we suggest possible areas for future research, including the use of novel machine learning techniques and exploring ‘open data’ that is becoming increasingly available.

Earth-Science Reviews↗
USGS2000

MODFLOW-2000, the U.S. Geological Survey modular ground-water model; user guide to the observation, sensitivity, and parameter-estimation processes and three post-processing programs

This report documents the Observation, Sensitivity, and Parameter-Estimation Processes of the ground-water modeling computer program MODFLOW-2000. The Observation Process generates model-calculated values for comparison with measured, or observed, quantities. A variety of statistics is calculated to quantify this comparison, including a weighted least-squares objective function. In addition, a number of files are produced that can be used to compare the values graphically. The Sensitivity Process calculates the sensitivity of hydraulic heads throughout the model with respect to specified parameters using the accurate sensitivity-equation method. These are called grid sensitivities. If the Observation Process is active, it uses the grid sensitivities to calculate sensitivities for the simulated values associated with the observations. These are called observation sensitivities. Observation sensitivities are used to calculate a number of statistics that can be used (1) to diagnose inadequate data, (2) to identify parameters that probably cannot be estimated by regression using the available observations, and (3) to evaluate the utility of proposed new data. The Parameter-Estimation Process uses a modified Gauss-Newton method to adjust values of user-selected input parameters in an iterative procedure to minimize the value of the weighted least-squares objective function. Statistics produced by the Parameter-Estimation Process can be used to evaluate estimated parameter values; statistics produced by the Observation Process and post-processing program RESAN-2000 can be used to evaluate how accurately the model represents the actual processes; statistics produced by post-processing program YCINT-2000 can be used to quantify the uncertainty of model simulated values. Parameters are defined in the Ground-Water Flow Process input files and can be used to calculate most model inputs, such as: for explicitly defined model layers, horizontal hydraulic conductivity, horizontal anisotropy, vertical hydraulic conductivity or vertical anisotropy, specific storage, and specific yield; and, for implicitly represented layers, vertical hydraulic conductivity. In addition, parameters can be defined to calculate the hydraulic conductance of the River, General-Head Boundary, and Drain Packages; areal recharge rates of the Recharge Package; maximum evapotranspiration of the Evapotranspiration Package; pumpage or the rate of flow at defined-flux boundaries of the Well Package; and the hydraulic head at constant-head boundaries. The spatial variation of model inputs produced using defined parameters is very flexible, including interpolated distributions that require the summation of contributions from different parameters. Observations can include measured hydraulic heads or temporal changes in hydraulic heads, measured gains and losses along head-dependent boundaries (such as streams), flows through constant-head boundaries, and advective transport through the system, which generally would be inferred from measured concentrations. MODFLOW-2000 is intended for use on any computer operating system. The program consists of algorithms programmed in Fortran 90, which efficiently performs numerical calculations and is fully compatible with the newer Fortran 95. The code is easily modified to be compatible with FORTRAN 77. Coordination for multiple processors is accommodated using Message Passing Interface (MPI) commands. The program is designed in a modular fashion that is intended to support inclusion of new capabilities.

Open-File Report↗
USGS2008

Investigation and hazard assessment of the 2003 and 2007 Staircase Falls rock falls, Yosemite National Park, California, USA

Since 1857 more than 600 rock falls, rock slides, debris slides, and debris flows have been documented in Yosemite National Park, with rock falls in Yosemite Valley representing the majority of the events. On 26 December 2003, a rock fall originating from west of Glacier Point sent approximately 200 m 3 of rock debris down a series of joint-controlled ledges to the floor of Yosemite Valley. The debris impacted talus near the base of Staircase Falls, producing fragments of flying rock that struck occupied cabins in Curry Village. Several years later on 9 June 2007, and again on 26 July 2007, smaller rock falls originated from the same source area. The 26 December 2003 event coincided with a severe winter storm and was likely triggered by precipitation and/or frost wedging, but the 9 June and 26 July 2007 events lack recognizable triggering mechanisms. We investigated the geologic and hydrologic factors contributing to the Staircase Falls rock falls, including bedrock lithology, weathering, joint spacing and orientations, and hydrologic processes affecting slope stability. We improved upon previous geomorphic assessment of rock-fall hazards, based on a shadow angle approach, by using STONE, a three-dimensional rock-fall simulation computer program. STONE produced simulated rock-fall runout patterns similar to the mapped extent of the 2003 and 2007 events, allowing us to simulate potential future rock falls from the Staircase Falls detachment area. Observations of recent rock falls, mapping of rock debris, and simulations of rock fall runouts beneath the Staircase Falls detachment area suggest that rock-fall hazard zones extend farther downslope than the extent previously defined by mapped surface talus deposits.

Natural Hazards and Earth System Sciences↗
USGS2006

Numerical modeling of magnetic moments for UXO applications

The surface magnetic anomaly observed in UXO clearance is mainly dipolar and, consequently, the dipole is the only magnetic moment regularly recovered in UXO applications. The dipole moment contains information about intensity of magnetization but lacks information about shape. In contrast, higher-order moments, such as quadrupole and octupole, encode asymmetry properties of the magnetization distribution within the buried targets. In order to improve our understanding of magnetization distribution within UXO and non-UXO objects and its potential utility in UXO clearance, we present a 3D numerical modeling study for highly susceptible metallic objects. The basis for the modeling is the solution of a nonlinear integral equation describing magnetization within isolated objects. A solution for magnetization distribution then allows us to compute magnetic moments of the object, analyze their relationships, and provide a depiction of the surface anomaly produced by different moments within the object. Our modeling results show significant high-order moments for more asymmetric objects situated at depths typical of UXO burial, and suggest that the increased relative contribution to magnetic gradient data from these higher-order moments may provide a practical tool for improved UXO discrimination.

Conference Paper↗
USGS2006

Relative importance of magnetic moments in UXO clearance applications

Surface magnetic anomaly observed in UXO clearance is mainly dipolar and, as a result, the dipole is the only moment used regularly in UXO applications. The dipole moment contains intensity of magnetization information but lacks shape information. Unlike dipole, higher-order moments, such as quadrupole and octupole, encode asymmetry properties of magnetization distribution within buried targets. In order to improve our understanding of magnetization distribution within UXO and non-UXO objects and its potential utility in UXO clearance, we present results of a 3D numerical modeling study for highly susceptible metallic objects. The basis for modeling is the solution of a nonlinear integral equation, describing magnetization within isolated objects, allowing us to compute magnetic moments of the object, analyze their relationships, and provide a depiction of the surface anomaly produced by the different moments within the object. Our modeling results show significant high-order moments for more asymmetric objects situated at typical UXO burial depths, and suggest that the increased relative contribution to magnetic gradient data from these higher-order moments may provide a practical tool for improved UXO discrimination. ?? 2005 Society of Exploration Geophysicists.

SEG Technical Program Expanded Abstracts↗