SEARCH · Geology Reports

Results for “Applied Sciences”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,068 records · Page 60Linked to original sources

USGS2020

An open source database for the synthesis of soil radiocarbon data: ISRaD version 1.0

Radiocarbon is a critical constraint on our estimates of the timescales of soil carbon cycling that can aid in identifying mechanisms of carbon stabilization and destabilization and improve the forecast of soil carbon response to management or environmental change. Despite the wealth of soil radiocarbon data that have been reported over the past 75 years, the ability to apply these data to global-scale questions is limited by our capacity to synthesize and compare measurements generated using a variety of methods. Here, we present the International Soil Radiocarbon Database (ISRaD; http://soilradiocarbon.org , last access: 16 December 2019), an open-source archive of soil data that include reported measurements from bulk soils, distinct soil carbon pools isolated in the laboratory by a variety of soil fractionation methods, samples of soil gas or water collected interstitially from within an intact soil profile, CO 2 gas isolated from laboratory soil incubations, and fluxes collected in situ from a soil profile. The core of ISRaD is a relational database structured around individual datasets (entries) and organized hierarchically to report soil radiocarbon data, measured at different physical and temporal scales as well as other soil or environmental properties that may also be measured and may assist with interpretation and context. Anyone may contribute their own data to the database by entering it into the ISRaD template and subjecting it to quality assurance protocols. ISRaD can be accessed through (1) a web-based interface, (2) an R package (ISRaD), or (3) direct access to code and data through the GitHub repository, which hosts both code and data. The design of ISRaD allows for participants to become directly involved in the management, design, and application of ISRaD data. The synthesized dataset is available in two forms: the original data as reported by the authors of the datasets and an enhanced dataset that includes ancillary geospatial data calculated within the ISRaD framework. ISRaD also provides data management tools in the ISRaD-R package that provide a starting point for data analysis; as an open-source project, the broader soil community is invited and encouraged to add data, tools, and ideas for improvement. As a whole, ISRaD provides resources to aid our evaluation of soil dynamics across a range of spatial and temporal scales. The ISRaD v1.0 dataset is archived and freely available at https://doi.org/10.5281/zenodo.2613911 (Lawrence et al., 2019).

Earth System Science Data↗

USGS2015

SHRIMP U–Pb and REE data pertaining to the origins of xenotime in Belt Supergroup rocks: evidence for ages of deposition, hydrothermal alteration, and metamorphism

The Belt–Purcell Supergroup, northern Idaho, western Montana, and southern British Columbia, is a thick succession of Mesoproterozoic sedimentary rocks with an age range of about 1470–1400 Ma. Stratigraphic layers within several sedimentary units were sampled to apply the new technique of U–Pb dating of xenotime that sometimes forms as rims on detrital zircon during burial diagenesis; xenotime also can form epitaxial overgrowths on zircon during hydrothermal and metamorphic events. Belt Supergroup units sampled are the Prichard and Revett Formations in the lower Belt, and the McNamara and Garnet Range Formations and Pilcher Quartzite in the upper Belt. Additionally, all samples that yielded xenotime were also processed for detrital zircon to provide maximum age constraints for the time of deposition and information about provenances; the sample of Prichard Formation yielded monazite that was also analyzed. Ten xenotime overgrowths from the Prichard Formation yielded a U–Pb age of 1458 ± 4 Ma. However, because scanning electron microscope – backscattered electrons (SEM–BSE) imagery suggests complications due to possible analysis of multiple age zones, we prefer a slightly older age of 1462 ± 6 Ma derived from the three oldest samples, within error of a previous U–Pb zircon age on the syn-sedimentary Plains sill. We interpret the Prichard xenotime as diagenetic in origin. Monazite from the Prichard Formation, originally thought to be detrital, yielded Cretaceous metamorphic ages. Xenotime from the McNamara and Garnet Range Formations and Pilcher Quartzite formed at about 1160– 1050 Ma, several hundred million years after deposition, and probably also experienced Early Cretaceous growth. These xenotime overgrowths are interpreted as metamorphic–diagenetic in origin (i.e., derived during greenschist facies metamorphism elsewhere in the basin, but deposited in sub-greenschist facies rocks). Several xenotime grains are older detrital grains of igneous derivation. A previous study on the Revett Formation at the Spar Lake Ag–Cu deposit provides data for xenotime overgrowths in several ore zones formed by hydrothermal processes; herein, those results are compared with data from newly analyzed diagenetic, metamorphic, and magmatic xenotime overgrowths. The origin of a xenotime overgrowth is reflected in its rareearth element (REE) pattern. Detrital (i.e., igneous) xenotime has a large negative Eu anomaly and is heavy rare-earth element (HREE)-enriched (similar to REE in igneous zircon). Diagenetic xenotime has a small negative Eu anomaly and flat HREE (Tb to Lu). Hydrothermal xenotime is depleted in light rare-earth element (LREE), has a small negative Eu anomaly, and decreasing HREE. Metamorphic xenotime is very LREE-depleted, has a very small negative Eu anomaly, and is strongly depleted in HREE (from Gd to Lu). Because these characteristics seem to be process related, they may be useful for interpretation of xenotime of unknown origin. The occurrence of 1.16–1.05 Ga metamorphic xenotime, in the apparent absence of pervasive deformation structures, suggests that the heating may be related to poorly understood regional heating due to broad regional underplating of mafic magma. These results may be additional evidence (together with published ages from metamorphic titanite, zircon, monazite, and garnet) for an enigmatic, Grenville-age metamorphic event that is more widely recognized in the southwestern and eastern United States

Alberta, British Columbia, Idaho, Montana, Oregon,↗

USGS2016

A review of single-sample-based models and other approaches for radiocarbon dating of dissolved inorganic carbon in groundwater

Numerous methods have been proposed to estimate the pre-nuclear-detonation 14 C content of dissolved inorganic carbon (DIC) recharged to groundwater that has been corrected/adjusted for geochemical processes in the absence of radioactive decay ( 14 C 0 ) - a quantity that is essential for estimation of radiocarbon age of DIC in groundwater. The models/approaches most commonly used are grouped as follows: (1) single-sample-based models, (2) a statistical approach based on the observed (curved) relationship between 14 C and δ 13 C data for the aquifer, and (3) the geochemical mass-balance approach that constructs adjustment models accounting for all the geochemical reactions known to occur along a groundwater flow path. This review discusses first the geochemical processes behind each of the single-sample-based models, followed by discussions of the statistical approach and the geochemical mass-balance approach. Finally, the applications, advantages and limitations of the three groups of models/approaches are discussed. The single-sample-based models constitute the prevailing use of 14 C data in hydrogeology and hydrological studies. This is in part because the models are applied to an individual water sample to estimate the 14 C age, therefore the measurement data are easily available. These models have been shown to provide realistic radiocarbon ages in many studies. However, they usually are limited to simple carbonate aquifers and selection of model may have significant effects on 14 C 0 often resulting in a wide range of estimates of 14 C ages. Of the single-sample-based models, four are recommended for the estimation of 14 C 0 of DIC in groundwater: Pearson's model, (Ingerson and Pearson, 1964; Pearson and White, 1967), Han & Plummer's model (Han and Plummer, 2013), the IAEA model (Gonfiantini, 1972; Salem et al., 1980), and Oeschger's model (Geyh, 2000). These four models include all processes considered in single-sample-based models, and can be used in different ranges of 13 C values. In contrast to the single-sample-based models, the extended Gonfiantini & Zuppi model (Gonfiantini and Zuppi, 2003; Han et al., 2014) is a statistical approach. This approach can be used to estimate 14 C ages when a curved relationship between the 14 C and 13 C values of the DIC data is observed. In addition to estimation of groundwater ages, the relationship between 14 C and δ 13 C data can be used to interpret hydrogeological characteristics of the aquifer, e.g. estimating apparent rates of geochemical reactions and revealing the complexity of the geochemical environment, and identify samples that are not affected by the same set of reactions/processes as the rest of the dataset. The investigated water samples may have a wide range of ages, and for waters with very low values of 14 C, the model based on statistics may give more reliable age estimates than those obtained from single-sample-based models. In the extended Gonfiantini & Zuppi model, a representative system-wide value of the initial 14 C content is derived from the 14 C and δ 13 C data of DIC and can differ from that used in single-sample-based models. Therefore, the extended Gonfiantini & Zuppi model usually avoids the effect of modern water components which might retain ‘bomb’ pulse signatures. The geochemical mass-balance approach constructs an adjustment model that accounts for all the geochemical reactions known to occur along an aquifer flow path (Plummer et al., 1983; Wigley et al., 1978; Plummer et al., 1994; Plummer and Glynn, 2013), and includes, in addition to DIC, dissolved organic carbon (DOC) and methane (CH 4 ). If sufficient chemical, mineralogical and isotopic data are available, the geochemical mass-balance method can yield the most accurate estimates of the adjusted radiocarbon age. The main limitation of this approach is that complete information is necessary on chemical, mineralogical and isotopic data and these data are often limited. Failure to recognize the limitations and underlying assumptions on which the various models and approaches are based can result in a wide range of estimates of 14 C 0 and limit the usefulness of radiocarbon as a dating tool for groundwater. In each of the three generalized approaches (single-sample-based models, statistical approach, and geochemical mass-balance approach), successful application depends on scrutiny of the isotopic ( 14 C and 13 C) and chemical data to conceptualize the reactions and processes that affect the 14 C content of DIC in aquifers. The recently developed graphical analysis method is shown to aid in determining which approach is most appropriate for the isotopic and chemical data from a groundwater system.

Earth-Science Reviews↗

USGS2016

Quantifying resilience

The biosphere is under unprecedented pressure, reflected in rapid changes in our global ecological, social, technological and economic systems. In many cases, ecological and social systems can adapt to these changes over time, but when a critical threshold is surpassed, a system under stress can undergo catastrophic change and reorganize into a different state. The concept of resilience, introduced more than 40 years ago in the ecological sciences, captures the behaviour of systems that can occur in alternative states. The original definition of resilience forwarded by Holling ( 1973 ) is still the most useful. It defines resilience as the amount of disturbance that a system can withstand before it shifts into an alternative stable state. The idea of alternative stable states has clear and profound implications for ecological management. Coral reefs, for example, are high-diversity systems that provide key ecosystem services such as fisheries and coastal protection. Human impacts are causing significant, ongoing reef degradation, and many reefs have shifted from coral- to algal-dominated states in response to anthropogenic pressures such as elevated water temperatures and overfishing. Understanding and differentiating between the factors that help maintain reefs in coral-dominated states vs. those that facilitate a shift to an undesired algal-dominated state is a critical step towards sound management and conservation of these, and other, important social–ecological systems. Resilience has gained popularity among both academicians and laypeople, as a term meant to describe a systems’ ability to withstand disturbance. Resilience has become a buzzword in the last decade, as shown by its increasing appearance in calls for research proposals and scientific citation data bases. The term resilience has in many cases lost the clarity of the original definition and in fact is frequently used in a manner in direct opposition to the original definition. Many current uses of the concept are loose and incorrect. The term is becoming increasingly used in a normative sense (Brand & Jax 2007 ), as if resilience were a desirable quality of systems. However, even systems in highly undesirable states, such as macro-algae dominated reefs, or city cores in poverty traps, may be highly resilient, which is to say they withstand attempts to transform them into different (desirable) states. Operationalizing the concept of resilience for application and management has been difficult. Misuse of the term can have significant negative impacts, because resilience is being used to help guide responses to natural disasters and to assess the sustainability of ecosystems and urban systems and has been driving international research priorities. Resilience has been argued to be a basic emergent property of systems, a process or a rate. We focus on the original concept as described by Holling, which is that of an emergent system property; when a system is in a desirable state and managers wish to enhance resilience, or when the system is in an undesirable state and managers wish to erode resilience and foster a transformation to an alternative state. Fostering or eroding resilience is a process. When a system is perturbed but resilience is not exceeded, then the recovery can be measured as a rate. Several frameworks to operationalize resilience have been proposed. A decade ago, a special feature focused on quantifying resilience was published in the journal Ecosystems (Carpenter, Westley & Turner 2005 ). The approach there was towards identifying surrogates of resilience, but few of the papers proposed quantifiable metrics. Consequently, many ecological resilience frameworks remain vague and difficult to quantify, a problem that this special feature aims to address. However, considerable progress has been made during the last decade (e.g. Pope, Allen & Angeler 2014 ). Although some argue that resilience is best kept as an unquantifiable, vague concept (Quinlan et al . 2016 ), to be useful for managers, there must be concrete guidance regarding how and what to manage and how to measure success (Garmestani, Allen & Benson 2013 ; Spears et al . 2015 ). Ideas such as ‘resilience thinking’ have utility in helping stakeholders conceptualize their systems, but provide little guidance on how to make resilience useful for ecosystem management, other than suggesting an ambiguous, Goldilocks approach of being just right (e.g. diverse, but not too diverse; connected, but not too connected). Here, we clarify some prominent resilience terms and concepts, introduce and synthesize the papers in this special feature on quantifying resilience and identify core unanswered questions related to resilience.

Journal of Applied Ecology↗

USGS2020

The Missoula and Bonneville floods—A review of ice-age megafloods in the Columbia River basin

The Channeled Scabland of eastern Washington State, USA, brought megafloods to the scientific forefront. A 30,000-km2 landscape of coulees and cataracts carved into the region’s loess-covered basalt attests to overwhelming volumes of energetic water. The scarred landscape, garnished by huge boulder bars and far-travelled ice-rafted erratics, spurred J Harlen Bretz’s vigorously disputed flood hypothesis in the 1920s. First known as the Spokane flood, it was rebranded the Missoula flood once understood that the water came from glacial Lake Missoula, formed when the Purcell Trench lobe of the last-glacial Cordilleran ice sheet dammed the Clark Fork valley in northwestern Idaho with ice a kilometer thick. Bretz’s flood evidence in the once-remote Channeled Scabland, widely seen and elaborated by the 1950s, eventually swayed consensus for cataclysmic flooding. Missoula flood questions then turned to some that continue today: how many? when? how big? what routes? what processes? The Missoula floods passed through eastern Washington by a multitude of valleys, coulees and scabland tracts, some contemporaneously, some sequentially. Which routings and their timing depended on the positions of various lobes of the multi-pronged Cordilleran ice sheet and the erosional development of the channels themselves. The first floods mostly followed the big bend of Columbia valley looping through north-central Washington. But the south-advancing Okanogan ice lobe soon blocked that path, forming long-lasting glacial Lake Columbia in the impounded Columbia valley. Missoula floods into this lake were diverted south out of the Columbia valley and into eastern Washington coulees and scabland tracts. At least four floods entered Moses Coulee, but then as the Okanogan lobe advanced over and blocked the head of that coulee, more eastern paths took the water, including Grand Coulee and the Telford-Crab-Creek and Cheney-Palouse scabland tracts. Flood routing also depended on the erosion of the coulees. At some point, headward erosion of upper Grand Coulee lowered the divide saddle between the west-running Columbia valley and the deep and wide Grand Coulee heading southwest. Still uncertain is when this happened and the consequences with respect to the stage and extent of glacial Lake Columbia and to flood access to the other, higher, flood routes. Downstream, all flood routes converged onto Pasco Basin, flowed through Wallula Gap and the Columbia River Gorge into the Pacific Ocean, following submarine canyons and depositing sediment layers on abyssal plains. Stratigraphic studies indicate dozens—likely more than a hundred—separate Missoula floods during the last glacial period. Over the length of the flood route, backwater areas and depositional basins preserve multiple flood beds, many of which are separated by signs of time, including volcanic ash layers and soil development in subaerial environments; and varve-like beds and pelagic mud layers in lacustrine and marine settings. Evidence also comes from the glacial Lake Missoula basin, where stratigraphy indicates dozens of filling and emptying cycles. Varve counts in conjunction of radiocarbon dating and paleomagnetic secular variation show the repeated filling-and-release cycles of glacial Lake Missoula had intervals possibly as long as 100 years early in the lake’s history but diminished to just one or two years for the last few floods. This behavior accords with jökulhlaup-style floods released by subglacial drainage from a self-dumping ice-dammed lake. But not yet clear is whether such a mechanism applies to all the floods or if some emptied more cataclysmically as hypothesized by some. Radiocarbon dating of sparse organic materials remains key to defining flood chronology but has been lately bolstered by analyses of terrestrial cosmogenic nuclides and optically stimulated luminescence. Varve counts and paleomagnetic secular variation studies help to define durations and intervals represented by sequences of flood beds. The ~16 ka Mount St. Helens Set S tephra is commonly interbedded within flood deposits, enabling correlation of deposits among sites. Tephra from the 13.7–13.4 ka eruption of Glacier Peak overlies all glacial Lake Missoula and Missoula flood deposits, defining an end time. Overall conclusions are that glacial Lake Missoula was extant and producing floods for at least 3–4 ky during 20–14 ka. At least ~75 floods preceded Mount St Helens Set S, followed by 30 or more after the tephra fall. Most floods entered glacial Lake Columbia, impounded by the Okanogan lobe, for 2–5 ky between about 18.5 and 15 ka. Glacial Lake Columbia outlived Lake Missoula by >200–400 yr but may have been born later since at least one flood came down the Columbia valley before the Okanogan ice lobe blocked the Columbia valley at 18.5–18 ka. The maximum extent of the Okanogan and Purcell Trench lobes, many Missoula floods, substantial erosion of upper Grand Coulee, and the widespread tephra falls from Mount St. Helens eruptions all happened about 17–15 ka. People, in the area since 16.6–15.3 ka, almost certainly witnessed the last of the Missoula floods and later large floods from other ice-dammed lakes in the Columbia River basin. Quantitative flow analyses give peak discharge estimates and support understanding of erosional and depositional processes. The first flow assessments were simple cross-section calculations but recent assessments employ two-dimensional hydrodynamic models. The general finding is that emplacement of the maximum stage evidence requires about 20 million m3/s near the Lake Missoula outlet and about 5–15 million m3/s through Wallula Gap and downstream in the Columbia River Gorge. These hydraulic analyses raise still-unresolved questions regarding canyon erosion and possible additional water sources. The large Pleistocene Bonneville flood entered the Columbia River system from the southeast from pluvial Lake Bonneville, the Pleistocene predecessor to Great Salt Lake in the eastern Great Basin. During the last glacial, the lake basin filled, covering >50,000 km2 with 10,400 km3 of water before reaching its maximum possible stage governed by Red Rock Pass, the lowest divide separating the basin from the Snake River basin to the north. The overtopping lake rapidly incised 108–125 m into the Red Rock Pass outlet, spilling half of its total lake volume. G.K. Gilbert described the essential sequence in the 1870s, but the flood was mostly forgotten until the late 1950s when Harold Malde linked the spectacular scabland topography and bouldery “melon gravel” on the Snake River Plain to the Lake Bonneville overflow. The Bonneville flood appears to have been a singular event at about 18 ka. No evidence of multiple or pre-last-glacial spillovers has yet been found. Its total volume was about twice that of a maximum Lake Missoula flood yet its peak discharge was ~1 million m3/s, less than a tenth of the largest Missoula floods. Its comparatively simple flow path and much steadier flow make the Bonneville flood ideal for new studies of erosional and depositional processes. At least two floods seem to have passed down the Columbia valley after the last of the Missoula floods, including a large flood about ~14 ka likely from cataclysmic demise of the thinning Okanogan ice lobe dam impounding glacial Lake Columbia. Floods from earlier glacial ages left scant yet clear evidence in the Channeled Scabland and Columbia valley. But their source, timing, and magnitudes are little understood. Some deposits are paleomagnetically reversed, thus older than ~800 ka. Last-glacial floods and perhaps older ones affected the Snake River Plain, some likely sourced in lakes dammed by alpine glaciers in central Idaho.

Idaho, Oregon, Washington↗

USGS2025

Calibration of the Stream Salmonid Simulator (S3) model to estimate annual survival, movement, and food consumption by juvenile Chinook salmon (Oncorhynchus tshawytscha) in the restoration reach of the Trinity River, California, 2006–18

Executive Summary The Trinity River is managed in two sections: (1) from the upper 64-kilometer “restoration reach” downstream from Lewiston Dam to the confluence with the North Fork Trinity River, and (2) the 120-kilometer lower Trinity River downstream from the restoration reach. The Stream Salmonid Simulator (S3) has been previously applied to these reaches and the Klamath River. To estimate fish growth, past S3 calibration efforts in the Trinity and Klamath Rivers used maximum likelihood methods that considered only the abundance of juvenile Chinook salmon ( Oncorhynchus tshawytscha ) passing a fish trap to estimate survival and movement parameters, but not fish consumption. Previous calibrations did not estimate the average proportion of maximum consumption ( C y ) when estimating survival ( S y ) and movement ( M 0 y ) parameters across years ( y ) of data, but because no other information was available in the literature a fixed value of C y =0.66 was assumed. Therefore, the goal of this report is to present an alternative approach that calibrates the S3 model to multivariate data (that is, abundance and size), enabling the estimation of the average proportion of maximum consumption, in conjunction with survival and movement parameters for a particular migration year. We fit the S3 model to individual years of weekly trap abundance estimates and mean fish sizes (fork length) at the Pear Tree Gulch (hereafter referred to as Pear Tree) fish trap representing the restoration reach. We used the Earth Mover’s Distance (EMD) as the objective value to be minimized in parameter optimization. This approach estimated survival, movement, and consumption parameters for each migration year. Because we had information on the abundance of natural and hatchery produced juvenile salmon at the fish traps, we estimated survival and movement for natural and hatchery fish. S3 is a deterministic life-stage-structured population model that tracks daily growth, movement, and survival of juvenile Chinook Salmon. A key theme of the model is that river discharge affects habitat availability and capacity, which in turn drives density-dependent population dynamics. To explicitly link population dynamics to habitat quality and quantity, the river environment is constructed as a one-dimensional series of linked habitat units, each of which has an associated daily timeseries of discharge, water temperature, and useable habitat area or carrying capacity. In turn, the physical characteristics of each habitat unit and the number of fish occupying each unit drive survival and growth within each habitat unit and movement of fish among habitat units. The physical template of the restoration reach of the Trinity River was classified into 356 meso-habitat units comprised of runs, riffles, and pools. For each habitat unit, we developed a timeseries of daily discharge, water temperature, amount of available spawning habitat, and fry and parr carrying capacity. Capacity time series were constructed using state-of-the-art models of spatially explicit hydrodynamics and quantitative fish habitat relationships developed for the Trinity River. These variables were then used to drive population dynamics such as egg maturation and survival, and in turn, juvenile movement, growth, and survival. We estimated movement, survival, and consumption parameters by calibrating the model to 12 years of weekly juvenile abundance estimates and fish sizes at the Pear Tree fish trap near the downstream end of the restoration reach. We estimated parameters for 12 years that included a wide range of female spawner abundances (1,414–11,494) and water year types (critically dry–extremely wet). We contrast the estimated parameters to the corresponding number of female spawners and the total annual volume of water discharged for the Trinity River (Trinity River Restoration Program; https://www.trrp.net/restoration/flows/summary/ ). The calibration consisted of replicating historical conditions as closely as possible (for example, discharge; temperature; spawner abundance, spawning location and timing, and hatchery releases), and then running the model to predict weekly abundance passing the trap location from each brood year of adults and subsequent migration year of their juvenile progeny. Because density-dependent movement was favored in past evaluations, we estimated S3 parameters based on density-independent survival and density-dependent movement. Likewise, each year’s estimated survival parameter for natural ( S N y ) and hatchery ( S H y ) fish may be interpreted as the mean daily survival probability from emergence or hatchery release to the Pear Tree fish trap. Under density dependence, the estimated movement parameter for natural ( M 0N y ) and hatchery ( M 0H y ) fish represents the intercept of the Beverton-Holt model; the probability of remaining in a habitat at near-zero abundance. We estimated C y by using EMD and incorporating abundance and fish size into model calibration. Average daily proportions of maximum consumption, , across the years were generally high (=0.640; standard deviation (SD) SD=0.176), suggesting that fish were feeding at about two-thirds of expected maximum consumption rates. This average proportion of maximum consumption,is very similar to what has been assumed (=0.66) in previous Trinity and Klamath River S3 calibration and simulation efforts. In 2017, we estimated the lowest C y , suggesting lower average consumption for juvenile salmon in high-discharge water years. When this high discharge year was excluded, there was no apparent trend in C y with annual water volume. Estimates of survival showed little trend over the range in spawner abundances, but a trend towards higher natural and hatchery fish survival with higher annual volumes of water was apparent. Over the 12 years, the average survival of hatchery fish was =0.888 (SD=0.079) and the average survival natural fish was=0.969 (SD=0.01). With respect to fish movement, we estimated higher M 0N y and M 0H y with higher annual volumes of water in the Trinity River. Higher M N0 y or M H0 y suggest greater probability of remaining in a habitat at low fish densities, with potential for density-dependent processes in movement to occur. The highest M 0N y = 0.676 was estimated during brood year 2012, and the overall average for natural fish was =0.276 (SD=0.188) and for hatchery fish was=0.467 (SD=0.235). Under the Beverton-Holt model, as M 0N y or M 0H y approach zero, there is less capacity for change in fish movement as fish density increases. The S3 model was initialized with only the spatiotemporal distribution of spawners, so it performed well at capturing the essential outmigration features that are ultimately governed by rates of growth, movement, and mortality. We used a new optimization method that could accommodate multivariate data on abundance and fish size collected at the Pear Tree fish trap, enabling the calibration of S3 to estimate five parameters for 12 separate years of data. Incorporating weekly fish size data for each year in our parameter optimization process made the estimation of C y possible and represents a step forward in the fitting of the S3 model to fish trap data for the purposes of parameter calibration and the estimation of growth parameters with respect to annual conditions. We identified lack of fit and adding important effects into the S3 model may improve the S3 estimation and simulation of water scenarios. The Trinity River Restoration Program (TRRP) Science Advisory Board recommended that the TRRP focus on developing core elements of a decision support system (DSS; Buffington and others, 2014). Toward that end, the habitat and S3 models described in this report are both core elements of the DSS. The structure of S3 makes it a particularly useful fish production model for the DSS because population dynamics are sensitive to (1) water temperature, (2) daily discharge management, and (3) habitat quality and quantity. Each of these variables are key management parameters under consideration in the TRRP. As such, the S3 model may provide valuable insights into the potentially variable effects of different management decisions on the Trinity River.

California↗