USGS Science⌕ Search

SEARCH · USGS Science

Results for “Fast Times”

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 55 records · Page 3Linked to original sources

Higher nest predation favors rapid fledging at the cost of plumage quality in nestling birds

Life-history theory predicts that rapid growth comes at a cost to offspring quality and adult longevity. However, trade-offs have been examined primarily based on proximate variation within species rather than evolved differences across species. Evolved differences are important to examine because species may co-evolve mechanisms to reduce long-term costs of rapid growth. For example, trade-offs associated with fast growth and short development time might be expressed in ephemeral traits because costs associated with these traits are paid for only a short duration. We explored this idea using nestling body feathers, an ephemeral trait, and examined whether shorter development time yielded lower quality feathers. We found a strong trade-off whereby nestlings that spend more time in the nest produced higher quality plumage across 123 temperate and tropical species. Using a subset of these species (n=67), we found that plumage quality was lower in species with higher rates of nest predation, which favors rapid nestling growth and short development times. Our results suggest that evolution of more rapid growth may co-evolve with other mechanisms whereby long-term costs are reduced by shifting costs of allocational trade-offs to ephemeral traits that are short-lived.

The American Naturalist↗

Stress- and structure-controlled anisotropy in a region of complex faulting—Yuha Desert, California

We examine shear velocity anisotropy in the Yuha Desert, California using aftershocks of the 2010 M7.2 El Mayor-Cucapah earthquake. The Yuha Desert is underlain by a complex network of right- and left-lateral conjugate faults, some of which experienced triggered slip during the El Mayor-Cucapah earthquake. An automated method that implements multiple measurement windows and a range of bandpass filters is used to estimate the fast direction ( ϕ ) and delay time ( δt ) of the split shear waves. We find an average ϕ oriented approximately north–south suggesting it is primarily controlled by the regional maximum compressive stress direction. However, the spatial variability in ϕ reveals that the fault structures that underlie the Yuha Desert also influence the measured splitting parameters. We infer that the northeast- and northwest-oriented ϕ reflect shear fabric subparallel to the conjugate fault structures. We do not observe a simple correlation between δt and hypocentral distance. Instead, the observed spatial variation in δt suggests that near-source variation in anisotropic strength may be equal to or more important than effects local to the station. No temporal variation in splitting parameters is observed during the 70-day period following the main shock. In this region of complex faulting, we observe a spatially variable pattern of anisotropy that is both stress- and structure-controlled. This study suggests that shear fabric can form even along short, discontinuous fault strands with minimal offset.

California↗

Species variation in osmotic, cryoprotectant, and cooling rate tolerance in poultry, eagle, and Peregrine Falcon spermatozoa

Potential factors influencing spermatozoa survival to cryopreservation and thawing were analyzed across a range of the following avian species: domestic chicken ( Gallus domesticus ), domestic turkey ( Meleagris gallopavo ), golden eagle ( Aquila chrysaetos ), Bonelli’s eagle ( Hieraaetus fasciatus ), imperial eagle ( Aquila adalberti ), and peregrine falcon ( Falco peregrinus ). Studies focused on spermatozoa tolerance to the following: 1) osmotic stress, 2) different extracellular concentrations of the cryoprotectant dimethylacetamide (DMA), 3) equilibration times of 1 versus 4 h, 4) equilibration temperature of 4 versus 21°C, and 5) rapid versus slow cooling before cryopreservation and standard thawing. Sperm viability was assessed with the live/dead stain (SYBR-14/propidium iodine). Sperm viability at osmolalities ≥800 mOsm was higher ( P < 0.05) in raptor than poultry semen. Return to isotonicity after exposure to hypertonicity (3000 mOsm) decreased ( P < 0.05) number of viable spermatozoa in chicken, turkey, and golden and Bonelli’s eagle spermatozoa but not in imperial eagle or peregrine falcon spermatozoa. Differences were found in spermatozoa resistance to hypotonic conditions, with eagle species demonstrating the most tolerance. Semen, equilibrated for 1 h (4°C) in diluent containing DMA (≥2.06 M), experienced decreased ( P < 0.05) spermatozoa survival in all species, except the golden eagle and peregrine falcon. Number of surviving spermatozoa diminished progressively with increasing DMA concentrations in all species. Increased equilibration temperature (from 4 to 21°C) markedly reduced ( P < 0.05) spermatozoa survival in all species except the Bonelli’s eagle and turkey. Rapid cooling was detrimental ( P < 0.05) to spermatozoa from all species except the imperial eagle and the chicken. These results demonstrate that avian spermatozoa differ remarkably in response to osmotic changes, DMA concentrations, equilibration time, temperature, and survival after fast or slow freezing. These differences emphasize the need for species-specific studies in the development and enhancement of assisted breeding for poultry and endangered species.

Biology of Reproduction↗

Solution of the one-dimensional consolidation theory equation with a pseudospectral method

The one-dimensional consolidation theory equation is solved for an aquifer system using a pseudospectral method. The spatial derivatives are computed using Fast Fourier Transforms and the time derivative is solved using a fourth-order Runge-Kutta scheme. The computer model calculates compaction based on the void ratio changes accumulated during the simulated periods of time. Compactions and expansions resulting from groundwater withdrawals and recharges are simulated for two observation wells in Santa Clara Valley and two in San Joaquin Valley, California. Field data previously published are used to obtain mean values for the soil grain density and the compression index and to generate depth-dependent profiles for hydraulic conductivity and initial void ratio. The water-level plots for the wells studied were digitized and used to obtain the time dependent profiles of effective stress.

Conference Paper↗

Hydrographic characterization of two tidal creeks with implications for watershed land use, flushing times, and benthic production

Many coastal ecosystems are undergoing anthropogenic stress from large increases in population and urbanization. In many regions changes in freshwater and material inputs to the coastal zone are altering the biogeochemical and biological capacities of ecosystems. Despite increased watershed inputs, large tidal volumes and flushing indicative of macrotidal estuaries can modulate the fate of introduced materials masking some of the symptoms of eutrophication. The Land Use Coastal Ecosystem Study (LU-CES) examined linkages between land use and environmental properties of Malind and Okatee Creeks in South Carolina from 2001 to 2004. The objectives of this particular study were to assess the hydrography of the two macrotidal creek ecosystems, explore differences in dissolved oxygen (DO), and develop a better understanding of the variations in primary and benthic secondary production in southeastern creek ecosystems. Depth, pH, salinity, and DO were reduced and more variable in Malind Creek than in Okatee Creek, although both creeks had strong semidiurnal frequencies in salinity time signatures. While time series analyses of DO saturation in Malind Creek revealed a dominant semidiurnal pattern, Okatee Creek had a distinctly diel DO pattern. The strongly semidiurnal fluctuations in DO and reduced flushing time indicated that biological processes were not fast enough to influence DO in Malind Creek. The Okatee Creek system had a much greater storage volume, a wider marsh, and a dominant 25-h DO frequency. These attributes contributed to an estimated 8-10 times more phytoplankton-based carbon in Okatee Creek and twice the annual benthic production. As expected from their proximity to the upland, low surface area, and high organic content, both ecosystems were net heterotrophic. This fundamental understanding of tidal creek hydrography is being used to help define linkages among differential watershed land uses, flushing characteristics, and levels of biological production in coastal ecosystems of the southeastern United States.

South Carolina↗

Characterization of subsurface groundwater pathways from treated wastewater effluent infiltration along the South Rim of Grand Canyon National Park, Arizona, USA

Recharge of treated wastewater to the environment is increasingly becoming a component of water resources management, especially in water-limited areas such as the Grand Canyon. Treated wastewater from the Grand Canyon South Rim Wastewater Treatment Plant has been discharged to a dry wash located along the surface expression of a fault for decades, but the fate of the wastewater is unclear. Previous research suggested possible connections between the treated wastewater discharge sites and springs along the South Rim of Grand Canyon, particularly Monument spring. Additional investigations presented in this report add further evidence to the connection between the wastewater and Monument spring, utilizing per- and polyfluoroalkyl substances, artificial sweeteners, and stable isotopes of water. A dye trace established the first direct connection between wastewater and Monument spring where the dye added to treated wastewater took about 4–6 months to move through the subsurface and was detectable for 5 months. The persistence of wastewater tracers over time and the timing of dye movement suggests a fast subsurface flow pathway. On the basis of geologic observations of faults and jointing, groundwater flow follows a possible combination of karst and other structural pathways. The dye and wastewater tracers were monitored at other springs along the South Rim of Grand Canyon and the combination of repeated dye detections and multiple wastewater tracers were only found at Monument spring. Multiple lines of evidence provide a robust tracer suite to support evidence of structural/karst controls on preferential flowpaths; faults alone are unlikely to explain the flow pathway.

Arizona↗

Movement patterns, habitat use, and survival of Lahontan cutthroat trout in the Truckee River

Habitat fragmentation, hybridization, and competition with nonnative salmonids are viewed as major threats to Lahontan cutthroat trout Oncorhynchus clarkii henshawi . Understanding Lahontan cutthroat trout behavior and survival is a necessary step in the reintroduction and establishment of naturally reproducing populations of Lahontan cutthroat trout. We used weekly radiotelemetry monitoring to examine movement patterns, habitat use, and apparent survival of 42 hatchery-reared Lahontan cutthroat trout in a 16.5-km stretch of the Truckee River, Nevada, across three reaches separated by barriers to upstream movement. We found differences in total movement distances and home range sizes of fish in different reaches within our study area. Fish used pool habitats more than fast water habitats in all reaches. Time of year, stream temperature, and fish standard length covariates had the strongest relationship with apparent survival. Monthly apparent survival was lowest in January, which coincided with the lowest flows and temperatures during the study period. Our results verify the mobility of Lahontan cutthroat trout and indicate that conditions during winter may limit the survival and reintroduction success in the portions of the Truckee River evaluated in this study.

California, Nevada↗

Supergene iron ores of minas Gerais, Brazil

The iron ores of Minas Gerais , Brazil , fall into two categories: (1) hypogene hematite ore averaging 66 percent or more Fe, and (2) lower-grade supergene ores . Most ore now extracted is high-grade hypogene ore ; lower-grade supergene ores will be of much future value. All supergene ores formed by weathering of itabirite, a metamorphosed oxide-facies iron formation averaging about 38 percent Fe and 44 percent Si02. The Caue Itabirite crops out for about 540 linear kilometers in central Minas Gerais . Supergene ores fall into three intergradational categories: (1) enriched itabirite, averaging 49 percent Fe, easily con-centratable, with reserves about 25,000 million tons; (2) intermediate grade ores , averaging perhaps 63 percent Fe, with indicated and inferred reserves more than 600 million tons; and (3) canga, averaging between 57 and 62 percent Fe, with reserves in the hundreds of millions of tons. Disaggregation of hard and brittle itabirite by solution primarily of quartz and secondarily of other soluble constituents causes residual enrichment in iron with minor hydration of hematite. As weathering continues, most of the quartz is removed and more of the hematite is hydrated, producing intermediate grade ore . Secondary enrichment by limonite is important. The final weathering product is canga. In canga, almost all the iron is hydrated, and the rock is still further impoverished in Si02 and residually enriched in A1203 and P. Canga also forms by cementation of iron -rich detritus by limonite. Four factors control the supergene ore -forming process: (1) physiography, for ores occur on ridges and plateaus; (2) climate, for seasonal rainfall is apparently needed for the formation of extensive canga blankets; (3) grain size of the original itabirite, for quartz with a grain size much greater than 0.1 millimeter is not readily soluble; and (4) composition of the iron formation. The chemically inert and physically resistant canga blanket is essential to supergene ore formation because soft weathering products would otherwise be removed as fast as formed. Thus, canga gives time for the formation of other ore types. Iron fixed as cementing limonite in canga and as enriching limonite in intermediate grade ore was derived by leaching and hydration of hematite from itabirite. It probably moved in the ferrous state and was precipitated as insoluble hydrous ferric oxide primarily by oxygenation of the solutions, to a lesser extent by their evaporation at or near the surface, and to a still smaller degree by pH changes. The data show that the high-grade hematite ore , 66 percent Fe or higher, cannot have formed directly by supergene action. Geochemical processes resulting in supergene concentration of iron also concentrated alumina and phosphorus. The high-grade ore contains about the same low percentages of these materials as unweathered itabirite.

Minas Gerais↗

Population dynamics of Yellowstone Cutthroat Trout in Henrys Lake, Idaho

Yellowstone Cutthroat Trout (YCT) Oncorhynchus clarkii bouvieri is a species with significant ecological and recreational value. In many YCT fisheries, managers are tasked with balancing angler expectations and fish conservation. Henrys Lake supports a popular trophy trout fishery, but the increase of nonnative Utah Chub Gila atraria has caused concern for YCT. We summarized long-term trends in abundance, length structure, body condition, and growth of YCT to evaluate the effect of Utah Chub. Additionally, we investigated abiotic and biotic factors influencing YCT. We examined archived hard structures to provide a comprehensive evaluation of changes in age and growth of YCT in the system. We used stocking records and catch rates of Utah Chub and trout in Henrys Lake as covariates to explain changes in YCT catch rates and growth. Catch rates varied from 1.5 to 15.4 YCT per net night during the 2002 to 2020 sampling period, but we did not identify consistent patterns. Length structure was consistently dominated by stock- to quality-length fish, and we captured few fish >600 mm in total length. Relative weight of YCT was decreased from a mean ± standard deviation (SD) of 115.9 ± 16.5 in 2004 to 93.2 ± 8.2 in 2020. The age of YCT varied between 1 and 11 years; fish that we captured during 2010 to 2020 were the oldest. The majority of fish that we sampled were age 4 and younger. Total annual mortality of age-2 and older YCT was higher than other Cutthroat Trout populations (i.e., 0.70 during 2002 to 2010 and 0.60 during 2011 to 2020). Based on regression models, we identified positive relationships between catch rates of YCT, Brook Trout Salvelinus fontinalis, and Rainbow Trout Oncorhynchus mykiss × YCT hybrid trout. We observed negative relationships between growth of YCT and abundance of Utah Chub and Brook Trout. Although we identified negative relationships, YCT growth in recent decades is as fast as or faster than earlier time periods. Results from this research suggest that major changes in YCT population dynamics are not evident over the last 20 years. This study provides insight into the factors influencing an adfluvial trout population. In particular, results from this research may be useful for managers of systems where Utah Chub have been introduced.

Idaho↗

Distributed Processing of Projections of Large Datasets: A Preliminary Study

Modern information needs have resulted in very large amounts of data being used in geographic information systems. Problems arise when trying to project these data in a reasonable amount of time and accuracy, however. Current single-threaded methods can suffer from two problems: fast projection with poor accuracy, or accurate projection with long processing time. A possible solution may be to combine accurate interpolation methods and distributed processing algorithms to quickly and accurately convert digital geospatial data between coordinate systems. Modern technology has made it possible to construct systems, such as Beowulf clusters, for a low cost and provide access to supercomputer-class technology. Combining these techniques may result in the ability to use large amounts of geographic data in time-critical situations.

Open-File Report↗

Grassification and fast-evolving fire connectivity and risk in the Sonoran Desert, United States

In the southwestern United States, non-native grass invasions have increased wildfire occurrence in deserts and the likelihood of fire spread to and from other biomes with disparate fire regimes. The elevational transition between desertscrub and montane grasslands, woodlands, and forests generally occurs at ∼1,200 masl and has experienced fast suburbanization and an expanding wildland-urban interface (WUI). In summer 2020, the Bighorn Fire in the Santa Catalina Mountains burned 486 km 2 and prompted alerts and evacuations along a 40-km stretch of WUI below 1,200 masl on the outskirts of Tucson, Arizona, a metropolitan area of >1M people. To better understand the changing nature of the WUI here and elsewhere in the region, we took a multidimensional and timely approach to assess fire dynamics along the Desertscrub-Semi-desert Grassland ecotone in the Catalina foothills, which is in various stages of non-native grass invasion. The Bighorn Fire was principally a forest fire driven by a long-history of fire suppression, accumulation of fine fuels following a wet winter and spring, and two decades of hotter droughts, culminating in the hottest and second driest summer in the 125-yr Tucson weather record. Saguaro ( Carnegia gigantea ), a giant columnar cactus, experienced high mortality. Resprouting by several desert shrub species may confer some post-fire resiliency in desertscrub. Buffelgrass and other non-native species played a minor role in carrying the fire due to the patchiness of infestation at the upper edge of the Desertscrub biome. Coupled state-and-transition fire-spread simulation models suggest a marked increase in both burned area and fire frequency if buffelgrass patches continue to expand and coalesce at the Desertscrub/Semi-desert Grassland interface. A survey of area residents six months after the fire showed awareness of buffelgrass was significantly higher among residents that were evacuated or lost recreation access, with higher awareness of fire risk, saguaro loss and declining property values, in that order. Sustained and timely efforts to document and assess fast-evolving fire connectivity due to grass invasions, and social awareness and perceptions, are needed to understand and motivate mitigation of an increasingly fire-prone future in the region.

Arizona↗

Rapid Source Characterization of the 2011 Mw 9.0 off the Pacific coast of Tohoku Earthquake

On March 11th, 2011, a moment magnitude 9.0 earthquake struck off the coast of northeast Honshu, Japan, generating what may well turn out to be the most costly natural disaster ever. In the hours following the event, the U.S. Geological Survey National Earthquake Information Center led a rapid response to characterize the earthquake in terms of its location, size, faulting source, shaking and slip distributions, and population exposure, in order to place the disaster in a framework necessary for timely humanitarian response. As part of this effort, fast finite-fault inversions using globally distributed body- and surface-wave data were used to estimate the slip distribution of the earthquake rupture. Models generated within 7 hours of the earthquake origin time indicated that the event ruptured a fault up to 300 km long, roughly centered on the earthquake hypocenter, and involved peak slips of 20 m or more. Updates since this preliminary solution improve the details of this inversion solution and thus our understanding of the rupture process. However, significant observations such as the up-dip nature of rupture propagation and the along-strike length of faulting did not significantly change, demonstrating the usefulness of rapid source characterization for understanding the first order characteristics of major earthquakes.

Tohoku, Honshu↗

Near real-time monitoring of volcanic surface deformation from GPS measurements at Long Valley Caldera, California

Long Valley Caldera in eastern California is an active volcanic area and has shown continued unrest in the last three decades. We have monitored surface deformation from Global Positioning System (GPS) data by using a projection method that we call Targeted Projection Operator (TPO). TPO projects residual time series with secular rates and periodic terms removed onto a predefined spatial pattern. We used the 2009–2010 slow deflation as a target spatial pattern. The resulting TPO time series shows a detailed deformation history including the 2007–2009 inflation, the 2009–2010 deflation, and a recent inflation that started in late-2011 and is continuing at the present time (November 2012). The recent inflation event is about four times faster than the previous 2007–2009 event. A Mogi source of the recent event is located beneath the resurgent dome at about 6.6 km depth at a rate of 0.009 km3/yr volume change. TPO is simple and fast and can provide a near real-time continuous monitoring tool without directly looking at all the data from many GPS sites in this potentially eruptive volcanic system.

California↗

The perception of volcanic risk in Kona communities from Mauna Loa and Hualālai volcanoes, Hawai'i

Volcanic hazards in Kona (i.e. the western side of the island of Hawai'i) stem primarily from Mauna Loa and Hualālai volcanoes. The former has erupted 39 times since 1832. Lava flows were emplaced in Kona during seven of these eruptions and last impacted Kona in 1950. Hualālai last erupted in ca. 1800. Society's proximity to potential eruptive sources and the potential for relatively fast-moving lava flows, coupled with relatively long time intervals since the last eruptions in Kona, are the underlying stimuli for this study of risk perception. Target populations were high-school students and adults ( n =462). Using these data, we discuss threat knowledge as an influence on risk perception, and perception as a driving mechanism for preparedness. Threat knowledge and perception of risk were found to be low to moderate. On average, fewer than two-thirds of the residents were aware of the most recent eruptions that impacted Kona, and a minority felt that Mauna Loa and Hualālai could ever erupt again. Furthermore, only about one-third were aware that lava flows could reach the coast in Kona in less than 3 h. Lava flows and ash fall were perceived to be among the least likely hazards to affect the respondent's community within the next 10 years, whereas vog (volcanic smog) was ranked the most likely. Less than 18% identified volcanic hazards as amongst the most likely hazards to affect them at home, school, or work. Not surprisingly, individual preparedness measures were found on average to be limited to simple tasks of value in frequently occurring domestic emergencies, whereas measures specific to infrequent hazard events such as volcanic eruptions were seldom adopted. Furthermore, our data show that respondents exhibit an 'unrealistic optimism bias' and infer that responsibility for community preparedness for future eruptions primarily rests with officials. We infer that these respondents may be less likely to attend to hazard information, react to warnings as directed, and undertake preparedness measures than other populations who perceive responsibility to lie with themselves. There are significant differences in hazard awareness and risk perception between students and adults, between subpopulations representing local areas, and between varying ethnicities. We conclude that long time intervals since damaging lava flows have occurred in Kona have contributed to lower levels of awareness and risk perceptions of the threat from lava flows, and that the on-going eruption at Kilauea has facilitated greater awareness and perception of risk of vog but not of other volcanic hazards. Low levels of preparedness may be explained by low perceptions of threat and risk and perhaps by the lack of a clear motivation or incentive to seek new modes of adjustment. 2003 Published by Elsevier B.V.

Hawaii↗

Growth rate variation in Brown Treesnakes (Boiga irregularis): An invasive species of conservation concern

Somatic growth rate is a fundamental trait that influences metabolism, lifespan and reproductive maturity and is critical for understanding population dynamics and informing management actions. Brown Treesnakes ( Boiga irregularis ) introduced to Guam are highly invasive and can reproduce year-round without discrete cohorts. We compared snake size trajectories described by the conventionally used von Bertalanffy growth function versus the Gompertz model. Using quantile regression with a regularized effect for individual snakes we modeled growth rates of 270 marked, wild snakes as a function of size. The Gompertz model explained more of the variation in growth and rendered more realistic predictions of asymptotic sizes than did the von Bertalanffy model. With the Gompertz model, growth rates were 1.05–1.16× faster in males than in females. Females reached asymptotic sizes at shorter snout-vent lengths than males. Growth rate was positively correlated with amount of precipitation, and modeling wet-dry seasonality on Guam as a sinusoidal function identified a growth peak in September—October. Effects of seasonality and precipitation, however, were minor compared to individual and sex related differences in size-adjusted growth rates. We estimated that the 50th (and 5th, 95th) growth-rate percentile males in our study population become sexually mature at an age of 33 (∞, 15) months, while females mature at 41 (∞, 18) months, where ∞ indicates that the slowest growing snakes never reach maturity. However, 50% of the snakes mature at a size below the median, and age at maturity may be as low as 10.4 (males) and 13.7 (females) months for average-sized hatchlings that grow fast. Our results have implications for the timing of management options for this species and our approach can be broadly applied to animals where repeated growth data are obtained and age is unknown.

Ecology and Evolution↗

Determination of nonylphenol isomers in landfill leachate and municipal wastewater using steam distillation extraction coupled with comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry

4-Nonylphenols (4-NPs) are known endocrine disruptors and by-products of the microbial degradation of nonylphenol polyethoxylate surfactants. One of the challenges to understanding the toxic effects of nonylphenols is the large number of isomers that may exist in environmental samples. In order to attribute toxic effects to specific compounds, a method is needed for the separation and quantitation of individual nonylphenol isomers. The pre-concentration methods of solvent sublimation, solid-phase extraction or liquid–liquid extraction prior to chromatographic analysis can be problematic because of co-extraction of thousands of compounds typically found in complex matrices such as municipal wastewater or landfill leachate. In the present study, steam distillation extraction (SDE) was found to be an effective pre-concentration method for extraction of 4-NPs from leachate and wastewater, and comprehensive two-dimensional gas chromatography (GC × GC) coupled with fast mass spectral data acquisition by time-of-flight mass spectrometry (ToFMS) enhanced the resolution and identification of 4-NP isomers. Concentrations of eight 4-NP isomers were determined in leachate from landfill cells of different age and wastewater influent and effluent samples. 4-NP isomers were about 3 times more abundant in leachate from the younger cell than the older one, whereas concentrations in wastewater effluent were either below detection limits or <1% of influent concentrations. 4-NP isomer distribution patterns were found to have been altered following release to the environment. This is believed to reflect isomer-specific degradation and accumulation of 4-NPs in the aquatic environment.

Journal of Chromatography A↗

Growth and smolting in lower-mode Atlantic Salmon stocked into the Penobscot River, Maine

Restoration of Atlantic Salmon Salmo salar in Maine has relied on hatchery-produced fry and smolts for critical stocking strategies. Stocking fry minimizes domestication selection, but these fish have poor survival. Conversely, stocked smolts have little freshwater experience but provide higher adult returns. Lower-mode (LM) fish, those not growing fast enough to ensure smolting by the time of stocking, are a by-product of the smolt program and are an intermediate hatchery product. From 2002 to 2009, between 70,000 and 170,000 marked LM Atlantic Salmon were stocked into the Pleasant River (a tributary in the Penobscot River drainage, Maine) in late September to early October. These fish were recaptured as actively migrating smolts (screw trapping), as nonmigrants (electrofishing), and as returning adults to the Penobscot River (Veazie Dam trap). Fork length (FL) was measured and a scale sample was taken to retrospectively estimate FL at winter annulus one (FW1) using the intercept-corrected direct proportion model. The LM fish were observed to migrate as age-1, age-2, and infrequently as age-3 smolts. Those migrating as age-1 smolts had a distinctly larger estimated FL at FW1 (>112 mm) than those that remained in the river for at least one additional year. At the time of migration, age-2 and age-3 smolts were substantially larger than age-1 smolts. Returning adult Atlantic Salmon of LM origin had estimated FLs at FW1 that corresponded to smolt age (greater FL for age 1 than age 2). The LM product produces both age-1 and age-2 smolts that have greater freshwater experience than hatchery smolts and may have growth and fitness advantages. The data from this study will allow managers to better assess the probability of smolting age and manipulate hatchery growth rates to produce a targeted-size LM product.

Maine↗

Uncertainty quantification for environmental models

Environmental models are used to evaluate the fate of fertilizers in agricultural settings (including soil denitrification), the degradation of hydrocarbons at spill sites, and water supply for people and ecosystems in small to large basins and cities—to mention but a few applications of these models. They also play a role in understanding and diagnosing potential environmental impacts of global climate change. The models are typically mildly to extremely nonlinear. The persistent demand for enhanced dynamics and resolution to improve model realism [17] means that lengthy individual model execution times will remain common, notwithstanding continued enhancements in computer power. In addition, high-dimensional parameter spaces are often defined, which increases the number of model runs required to quantify uncertainty [2]. Some environmental modeling projects have access to extensive funding and computational resources; many do not. The many recent studies of uncertainty quantification in environmental model predictions have focused on uncertainties related to data error and sparsity of data, expert judgment expressed mathematically through prior information, poorly known parameter values, and model structure (see, for example, [1,7,9,10,13,18]). Approaches for quantifying uncertainty include frequentist (potentially with prior information [7,9]), Bayesian [13,18,19], and likelihood-based. A few of the numerous methods, including some sensitivity and inverse methods with consequences for understanding and quantifying uncertainty, are as follows: Bayesian hierarchical modeling and Bayesian model averaging; single-objective optimization with error-based weighting [7] and multi-objective optimization [3]; methods based on local derivatives [2,7,10]; screening methods like OAT (one at a time) and the method of Morris [14]; FAST (Fourier amplitude sensitivity testing) [14]; the Sobol' method [14]; randomized maximum likelihood [10]; Markov chain Monte Carlo (MCMC) [10]. There are also bootstrapping and cross-validation approaches.Sometimes analyses are conducted using surrogate models [12]. The availability of so many options can be confusing. Categorizing methods based on fundamental questions assists in communicating the essential results of uncertainty analyses to stakeholders. Such questions can focus on model adequacy (e.g., How well does the model reproduce observed system characteristics and dynamics?) and sensitivity analysis (e.g., What parameters can be estimated with available data? What observations are important to parameters and predictions? What parameters are important to predictions?), as well as on the uncertainty quantification (e.g., How accurate and precise are the predictions?). The methods can also be classified by the number of model runs required: few (10s to 1000s) or many (10,000s to 1,000,000s). Of the methods listed above, the most computationally frugal are generally those based on local derivatives; MCMC methods tend to be among the most computationally demanding. Surrogate models (emulators)do not necessarily produce computational frugality because many runs of the full model are generally needed to create a meaningful surrogate model. With this categorization, we can, in general, address all the fundamental questions mentioned above using either computationally frugal or demanding methods. Model development and analysis can thus be conducted consistently using either computation-ally frugal or demanding methods; alternatively, different fundamental questions can be addressed using methods that require different levels of effort. Based on this perspective, we pose the question: Can computationally frugal methods be useful companions to computationally demanding meth-ods? The reliability of computationally frugal methods generally depends on the model being reasonably linear, which usually means smooth nonlin-earities and the assumption of Gaussian errors; both tend to be more valid with more linear

SIAM News↗