SEARCH · Geology Reports

Results for “Forest Science”

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.

At least 667 records · Page 37Linked to original sources

USGS2024

Environmental drivers and spatial patterns of antibiotic-resistant, enteric coliforms across a forest–urban riverscape

Antibiotic resistant bacteria are prevalent environmental contaminants in freshwaters, and antibiotic resistance genes circulate throughout the urban water cycle. The increase of antibiotic resistant pathogens threatens public health through direct and indirect exposure, and natural resource managers need information on the spatial patterns of antibiotic resistant bacteria and environmental factors associated with their distribution to improve water quality monitoring and to better assess human, animal, and environmental health risks. We collected water and epilithic biofilm samples and measured physicochemical environmental variables at 29 sites distributed longitudinally in the Green-Duwamish River basin, Washington, USA. We characterized catchment-wide patterns of gram-negative fecal indicator bacteria and hypothesized that the presence of antibiotic resistance would be associated with environmental heterogeneity, bacterial primary ecology, stream compartment, and stream type. Antibiotic resistance was determined by microbial growth on selective media supplemented with 3 different antibiotics (ampicillin, chloramphenicol, or tetracycline). Phenotypic antibiotic resistance was positively associated with disturbance, but resistance to at least 1 antibiotic was also detected in undeveloped river segments, with an 83% overall detection rate (i.e., 24 out of 29 sites, 17 in the mainstem and 7 in tributaries). The most probable number of Escherichia coli was associated with higher levels of antibiotic resistance of non- E. coli coliforms across the basin (ρ = 0.38, p < 0.01) but was not associated with antibiotic resistance of E. coli . Phenotypic resistance was highest among non- E. coli coliforms in the water column of tributaries draining moderately to extensively developed subcatchments. Generalized linear mixed-effects model results showed that 18% of the variance in presence of antibiotic resistance was explained by the fixed effects (summed CV across environmental variables, stream type, primary ecology, and stream compartment), and when a spatial random effect was included, the model explained 27% of the variance. Our study provides new evidence that environmental factors and bacterial primary ecology are important underlying factors associated with spatial patterns of antibiotic resistant enteric coliforms. We used macroecological concepts and a riverscape approach to characterize the distribution of antibiotic resistance with methods applicable to municipalities.

Freshwater Science↗

USGS2023

Dating the penultimate great earthquake in south-central Alaska using tree-ring crossdating and radiocarbon wiggle-matching

A forest bed of tree stumps currently in the intertidal zone at Girdwood, south-central Alaska, records coseismic submergence during the penultimate great earthquake. Tree-ring samples from ten spruce stumps were crossdated to develop a 149-year-long ring-width chronology. Radiocarbon wiggle-matching found that single-ring ages from the chronology were offset 28 ± 7 years older than the IntCal20 calibration curve and that the last ring of the chronology dated as 1169 to 1189 CE (781–761 cal. yr. BP) at the 95% confidence level. Bark was observed on some stumps, six samples had the same year for the last growth ring, and so this wiggle-match date is also the best estimate of the date of the penultimate great earthquake. This date is in good agreement with a date for this event in a seismo-turbidite record from Skilak Lake but not with previous dates from Bayesian models of maximum- and minimum-limiting ages from coastal salt marshes. Reanalysis of the coastal salt marsh ages with the data grouped by area, context and material found that outer wood samples from stumps at coseismic submergence sites and a Bayesian limiting age model based on just herbaceous plant ages from Turnagain Arm and the Copper River area are both consistent with our wiggle-match date. Furthermore, coseismic emergence ages from Cape Suckling and Yakataga are older than the penultimate earthquake and so likely relate to an earlier uplift event in this eastern area. The rupture extent during the penultimate great earthquake appears to have been less than in the 1964 great earthquake and the interseismic interval between these two events was 785 ± 10 years.

Alaska↗

USGS2023

The Factors Affecting Female Black Bear Harvest Rates in Pennsylvania

Pennsylvania’s black bear ( Ursus americanus ) population increased in abundance and distribution during the latter third of the 20th century, leading to an increase in human-bear conflicts. Increases in harvest opportunities from 2002–2018 aimed to stabilize black bear population growth but did not substantially increase harvest, and annual harvest was often below the desired goal of 20% removal. Consequently, additional changes to Pennsylvania’s black bear hunting seasons occurred from 2019–2021, including starting harvest earlier, expanding the length of seasons, and adding additional seasons (i.e., a muzzleloader and special rifle season). Understanding how earlier harvest seasons and new methods of take (i.e., muzzleloader) influence female black bear harvest vulnerability is important to informing harvest management. We trapped and GPS-collared adult female bears in the Sproul State Forest in northcentral Pennsylvania from 2019–2021 to determine home range sizes, patterns of resource selection, and sources of mortality during fall harvest seasons. We assessed annual variability in relative abundance of fall hard mast. We evaluated temporal and spatial variation in hunter activity with road-side surveys and remote camera traps, respectively. We estimated fall and weekly home range size with utilization distributions through an autocorrelated kernel density estimation and evaluated the influence of predictors hypothesized to influence third-order resource selection using generalized linear mixed models. We investigated factors hypothesized to influence female black bear survival during hunting seasons with known-fate models. Mean fall home range size was 248.7 km 2 (range = 6.1–2636.1 km 2 ). Home range sizes varied by year and were generally smaller during archery harvest season than other periods. Patterns of weekly resource selection indicated bears selected steeper slopes and higher elevations outside of harvest seasons but shifted to less-steep areas in the week before harvest and the first week of harvest, and to lower elevations during harvest. Bears selected for areas containing oak ( Quercus spp.) trees throughout the fall. Survival was lower in older age bears, greater relative mast abundance conditions, steeper slopes, and areas of greater hunter space use during the general firearms season. Survival was higher in areas of greater hunter space use during archery season. Harvest rate of adult female bears was 0.345 in 2019, 0.321 in 2020, and 0.150 in 2021, and averaged 0.272 across all three years. The probability of an adult female black bear surviving all harvest seasons each year was 0.611 (SE = 0.086, 95% CI = 0.436, 0.761). The high harvest rate and low predicted survival may lead to population reduction.

Pennsylvania↗

USGS2021

Estimating blue carbon sequestration under coastal management scenarios

Restoring and protecting “blue carbon” ecosystems - mangrove forests, tidal marshes, and seagrass meadows - are actions considered for increasing global carbon sequestration. To improve understanding of which management actions produce the greatest gains in sequestration, we used a spatially explicit model to compare carbon sequestration and its economic value over a broad spatial scale (2500 km of coastline in southeastern Australia) for four management scenarios: (1) Managed Retreat, (2) Managed Retreat Plus Levee Removal, (3) Erosion of High Risk Areas, (4) Erosion of Moderate to High Risk Areas. We found that carbon sequestration from avoiding erosion-related emissions (abatement) would far exceed sequestration from coastal restoration. If erosion were limited only to the areas with highest erosion risk, sequestration in the non-eroded area exceeded emissions by 4.2 million Mg CO 2 by 2100. However, losing blue carbon ecosystems in both moderate and high erosion risk areas would result in net emissions of 23.0 million Mg CO 2 by 2100. The removal of levees combined with managed retreat was the strategy that sequestered the most carbon. Across all time points, removal of levees increased sequestration by only an additional 1 to 3% compared to managed retreat alone. Compared to the baseline erosion scenario, the managed retreat scenario increased sequestration by 7.40 million Mg CO 2 by 2030, 8.69 million Mg CO 2 by 2050, and 16.6 million Mg CO 2 by 2100. Associated economic value followed the same patterns, with large potential value loss from erosion greater than potential gains from conserving or restoring ecosystems. This study quantifies the potential benefits of managed retreat and preventing erosion in existing blue carbon ecosystems to help meet climate change mitigation goals by reducing carbon emissions.

Victoria↗

USGS2006

Cenozoic tectonic evolution of the basin and range province in northwestern Nevada

A regional synthesis of new and existing geologic and thermochronologic data document late Cretaceous - early Cenozoic regional erosion, Oligocene - Miocene volcanism, and subsequent late Miocene extension of the Basin and Range Province in northwestern Nevada and northeastern California. Across an ???220-km-wide region between the Santa Rosa and Warner Ranges, conformable sequences of 35 to 15 Ma volcanic rocks are cut by only a single generation of high-angle normal faults that accommodated ???23 km of total east-west extension (???12%). Fission-track, (U-Th)/He, geologic, and structural data from the Pine Forest Range show that faulting there began at 11 to 12 Ma, progressed at a relatively constant rate until at least 3 Ma, and has continued until near the present time. Extension in the Santa Rosa Range to the east took place during the same interval, although the post-6 Ma part of this history is less well constrained. Less complete constraints from adjacent ranges permit a similar timing for faulting, and we infer that extensional faulting in northwestern Nevada began everywhere at 12 Ma and has continued up to the present. Faulting in the Warner Range in northeastern California can only be constrained to have begun between 14 and 3 Ma, but may represent westward migration of Basin and Range extension during the Pliocene. Compared to the many parts of the Basin and Range in central and southern Nevada, extension in northwestern Nevada began more recently, is of lesser total magnitude, and was accommodated entirely by high-angle normal faults. Fission-track data document Late Cretaceous unroofing of Cretaceous (115 - 100 Ma) granitic basement rocks in northwestern Nevada, followed by a long period of relative tectonic quiescence that persisted through Oligocene and Miocene volcanism until the onset of Basin and Range extension at ???12 Ma. The low magnitude of extension (12%) and early Tertiary stability suggest that the modern ???31 km thick crust in northwestern Nevada was only slightly thicker (???35 km) prior to extension at 12 Ma, and was no thicker than ???38 km in the Late Cretaceous. This stands in contrast to other parts of the Basin and Range, where the crust was thickened to at least 45 to 50 km by Cretaceous thrusting and subsequently thinned to ???30 km by large magnitude (>50%) extension.

American Journal of Science↗

USGS2014

Sources and sinks of carbon in boreal ecosystems of interior Alaska: a review

Boreal regions store large quantities of carbon but are increasingly vulnerable to carbon loss due to disturbance and climate warming. The boreal region, underlain by discontinuous permafrost, presents a challenging landscape for itemizing current and potential carbon sources and sinks in the boreal soil and vegetation. The roles of fire, forest succession, and the presence (or absence) of permafrost on carbon cycle, vegetation, and hydrologic processes have been the focus of multidisciplinary research in this area for the past 20 years. However, projections of a warming future climate, an increase in fire severity and extent, and the potential degradation of permafrost could lead to major landscape process changes over the next 20 to 50 years. This provides a major challenge for predicting how the interplay between land management activities and impacts of climate warming will affect carbon sources and sinks in Interior Alaska. To assist land managers in adapting and managing for potential changes in the Interior Alaska carbon cycle we developed this review paper incorporating an overview of the climate, ecosystem processes, vegetation types, and soil regimes in Interior Alaska with a focus on ramifications for the carbon cycle. Our objective is to provide a synthesis of the most current carbon storage estimates and measurements to support policy and land management decisions on how to best manage carbon sources and sinks in Interior Alaska. To support this we have surveyed relevant peer reviewed estimates of carbon stocks in aboveground and belowground biomass for Interior Alaska boreal ecosystems. We have also summarized methane and carbon dioxide fluxes from the same ecosystems. These data have been converted into the same units to facilitate comparison across ecosystem compartments. We identify potential changes in the carbon cycle with climate change and human disturbance including how compounding disturbances can affect the boreal system. Finally, we provide recommendations to address the challenges facing land managers in efforts to manage carbon cycle processes. The results of this study can be used for carbon cycle management in other locations within the boreal biome which encompass a broad distribution from 45° to 83° north.

Alaska↗

USGS2012

Soil-water dynamics and unsaturated storage during snowmelt following wildfire

Many forested watersheds with a substantial fraction of precipitation delivered as snow have the potential for landscape disturbance by wildfire. Little is known about the immediate effects of wildfire on snowmelt and near-surface hydrologic responses, including soil-water storage. Montane systems at the rain-snow transition have soil-water dynamics that are further complicated during the snowmelt period by strong aspect controls on snowmelt and soil thawing. Here we present data from field measurements of snow hydrology and subsurface hydrologic and temperature responses during the first winter and spring after the September 2010 Fourmile Canyon Fire in Colorado, USA. Our observations of soil-water content and soil temperature show sharp contrasts in hydrologic and thermal conditions between north- and south-facing slopes. South-facing burned soils were ∼1–2 °C warmer on average than north-facing burned soils and ∼1.5 °C warmer than south-facing unburned soils, which affected soil thawing during the snowmelt period. Soil-water dynamics also differed by aspect: in response to soil thawing, soil-water content increased approximately one month earlier on south-facing burned slopes than on north-facing burned slopes. While aspect and wildfire affect soil-water dynamics during snowmelt, soil-water storage at the end of the snowmelt period reached the value at field capacity for each plot, suggesting that post-snowmelt unsaturated storage was not substantially influenced by aspect in wildfire-affected areas. Our data and analysis indicate that the amount of snowmelt-driven groundwater recharge may be larger in wildfire-impacted areas, especially on south-facing slopes, because of earlier soil thaw and longer durations of soil-water contents above field capacity in those areas.

Hydrology and Earth System Sciences↗

USGS2016

Widespread kelp-derived carbon in pelagic and benthic nearshore fishes

Kelp forests provide habitat for diverse and abundant fish assemblages, but the extent to which kelp provides a source of energy to fish and other predators is unclear. To examine the use of kelp-derived energy by fishes we estimated the contribution of kelp- and phytoplankton-derived carbon using carbon (δ 13 C) and nitrogen (δ 15 N) isotopes measured in muscle tissue. Benthic-foraging kelp greenling (Hexagrammos decagrammus) and pelagic-foraging black rockfish (Sebastes melanops) were collected at eight sites spanning ∼35 to 60°N from the California Current (upwelling) to Alaska Coastal Current (downwelling) in the northeast Pacific Ocean. Muscle δ 13 C values were expected to be higher for fish tissue primarily derived from kelp, a benthic macroalgae, and lower for tissue primarily derived from phytoplankton, pelagic microalgae. Muscle δ 13 C values were higher in benthic-feeding kelp greenling than in pelagic-feeding black rockfish at seven of eight sites, indicating more kelp-derived carbon in greenling as expected. Estimates of kelp carbon contributions ranged from 36 to 89% in kelp greenling and 32 to 65% in black rockfish using carbon isotope mixing models. Isotopic evidence suggests that these two nearshore fishes routinely derive energy from kelp and phytoplankton, across coastal upwelling and downwelling systems. Thus, the foraging mode of nearshore predators has a small influence on their ultimate energy source as energy produced by benthic macroalgae and pelagic microalgae were incorporated in fish tissue regardless of feeding mode and suggest strong and widespread benthic-pelagic coupling. Widespread kelp contributions to benthic- and pelagic-feeding fishes suggests that kelp energy provides a benefit to nearshore fishes and highlights the potential for kelp and fish production to be linked.

Estuarine, Coastal and Shelf Science↗

USGS2017

National Park Service Vegetation Mapping Inventory Program: Appalachian National Scenic Trail vegetation mapping project

The National Park Service (NPS) Vegetation Mapping Inventory (VMI) Program classifies, describes, and maps existing vegetation of national park units for the NPS Natural Resource Inventory and Monitoring (I&M) Program. The NPS VMI Program is managed by the NPS I&M Division and provides baseline vegetation information to the NPS Natural Resource I&M Program. The U.S. Geological Survey Upper Midwest Environmental Sciences Center, NatureServe, NPS Northeast Temperate Network, and NPS Appalachian National Scenic Trail (APPA) have completed vegetation classification and mapping of APPA for the NPS VMI Program. Mappers, ecologists, and botanists collaborated to affirm vegetation types within the U.S. National Vegetation Classification (USNVC) of APPA and to determine how best to map the vegetation types by using aerial imagery. Analyses of data from 1,618 vegetation plots were used to describe USNVC associations of APPA. Data from 289 verification sites were collected to test the field key to vegetation associations and the application of vegetation associations to a sample set of map polygons. Data from 269 validation sites were collected to assess vegetation mapping prior to submitting the vegetation map for accuracy assessment (AA). Data from 3,265 AA sites were collected, of which 3,204 were used to test accuracy of the vegetation map layer. The collective of these datasets affirmed 280 USNVC associations for the APPA vegetation mapping project. To map the vegetation and land cover of APPA, 169 map classes were developed. The 169 map classes consist of 150 that represent natural (including ruderal) vegetation types in the USNVC, 11 that represent cultural (agricultural and developed) vegetation types in the USNVC, 5 that represent natural landscapes with catastrophic disturbance or some other modification to natural vegetation preventing accurate classification in the USNVC, and 3 that represent nonvegetated water (non-USNVC). Features were interpreted from viewing 4-band digital aerial imagery using digital onscreen three-dimensional stereoscopic workflow systems in geographic information systems (GIS). (Digital aerial imagery was collected each fall during 2009–11 to capture leaf-phenology change of hardwood trees across the latitudinal range of APPA.) The interpreted data were digitally and spatially referenced, thus making the spatial-database layers usable in GIS. Polygon units were mapped to either a 0.5-hectare (ha) or 0.25-ha minimum mapping unit, depending on vegetation type or scenario; however, polygon units were mapped to 0.1 ha for alpine vegetation. A geodatabase containing various feature-class layers and tables provide locations and support data to USNVC vegetation types (vegetation map layer), vegetation plots, verification sites, validation sites, AA sites, project boundary extent and zones, and aerial image centers and flight lines. The feature-class layer and related tables of the vegetation map layer provide 30,395 polygons of detailed attribute data covering 110,919.7 ha, with an average polygon size of 3.6 ha; the vegetation map coincides closely with the administrative boundary for APPA. Summary reports generated from the vegetation map layer of the map classes representing USNVC natural (including ruderal) vegetation types apply to 28,242 polygons (92.9% of polygons) and cover 106,413.0 ha (95.9%) of the map extent for APPA. The map layer indicates APPA to be 92.4% forest and woodland (102,480.8 ha), 1.7% shrubland (1866.3 ha), and 1.8% herbaceous cover (2,065.9 ha). Map classes representing park-special vegetation (undefined in the USNVC) apply to 58 polygons (0.2% of polygons) and cover 404.3 ha (0.4%) of the map extent. Map classes representing USNVC cultural types apply to 1,777 polygons (5.8% of polygons) and cover 2,516.3 ha (2.3%) of the map extent. Map classes representing nonvegetated water (non-USNVC) apply to 332 polygons (1.1% of polygons) and cover 1,586.2 ha (1.4%) of the map extent.

Appalachian National Scenic Trail↗

USGS2018

Predicting biological conditions for small headwater streams in the Chesapeake Bay watershed

A primary goal for Chesapeake Bay watershed restoration is to improve stream health and function in 10% of stream miles by 2025. Predictive spatial modeling of stream conditions, when accurate, is one method to fill gaps in monitoring coverage and estimate baseline conditions for restoration goals. Predictive modeling can also monitor progress as additional data become available. We developed a random forests model to predict biological condition of small streams (<200 km 2 in drainage) in the Chesapeake Bay watershed. Biological condition was measured with the Chesapeake Bay Basin-wide Index of Biotic Integrity (Chessie BIBI), a stream macroinvertebrate index. Our goal was to predict biological condition in all unsurveyed small streams present in a 1:24,000 scale catchment layer as a 2004–2008 baseline. We reclassified the 5-category Chessie BIBI ratings into two categories, poor and fair/good, to align with management goals of the Chesapeake Bay Program. The model included 12 geospatial predictor variables including measures on spatial location, bioregion, land cover, soil, precipitation, and number of dams in local catchments. We trained the model with a random 75% subset of Chessie BIBI data ( n = 1449), and used the remaining 25% of Chessie BIBI data ( n = 484) as test data. The model performed well, correctly predicting 72% of samples in training data and 73% of samples in test data, but model accuracy varied among bioregions. We performed uncertainty analyses by adding bands of either ±0.05 or ±0.10 BIBI units to the cutoff between poor and fair/good. These uncertainty analyses resulted in 14.5% (±0.05 band) and 24.8% (±0.10 band) of samples in test data being classified as in uncertain condition. For 95,877 small stream reaches in the Chesapeake Bay watershed, the model predicted 64% in fair/good condition, the ±0.05 uncertainty analyses predicted 57% in fair/good condition, and the ±0.10 uncertainty analysis predicted 50% in fair/good condition. These reported values have different implications for the number of improved stream miles required to meet the goal of improving 10%. Incorporating uncertainty provides an assessment of model strength as well as confidence in predictions. We, therefore, suggest increased reporting of uncertainty in studies that spatially predict stream conditions.

Chesapeake Bay watershed↗

USGS2024

Predicting redox conditions in groundwater at a national scale using random forest classification

Redox conditions in groundwater may markedly affect the fate and transport of nutrients, volatile organic compounds, and trace metals, with significant implications for human health. While many local assessments of redox conditions have been made, the spatial variability of redox reaction rates makes the determination of redox conditions at regional or national scales problematic. In this study, redox conditions in groundwater were predicted for the contiguous United States using random forest classification by relating measured water quality data from over 30,000 wells to natural and anthropogenic factors. The model correctly predicted the oxic/suboxic classification for 78 and 79% of the samples in the out-of-bag and hold-out data sets, respectively. Variables describing geology, hydrology, soil properties, and hydrologic position were among the most important factors affecting the likelihood of oxic conditions in groundwater. Important model variables tended to relate to aquifer recharge, groundwater travel time, or prevalence of electron donors, which are key drivers of redox conditions in groundwater. Partial dependence plots suggested that the likelihood of oxic conditions in groundwater decreased sharply as streams were approached and gradually as the depth below the water table increased. The probability of oxic groundwater increased as base flow index values increased, likely due to the prevalence of well-drained soils and geologic materials in high base flow index areas. The likelihood of oxic conditions increased as topographic wetness index (TWI) values decreased. High topographic wetness index values occur in areas with a propensity for standing water and overland flow, conditions that limit the delivery of dissolved oxygen to groundwater by recharge; higher TWI values also tend to occur in discharge areas, which may contain groundwater with long travel times. A second model was developed to predict the probability of elevated manganese (Mn) concentrations in groundwater (i.e., ≥50 μg/L). The Mn model relied on many of the same variables as the oxic/suboxic model and may be used to identify areas where Mn-reducing conditions occur and where there is an increased risk to domestic water supplies due to high Mn concentrations. Model predictions of redox conditions in groundwater produced in this study may help identify regions of the country with elevated groundwater vulnerability and stream vulnerability to groundwater-derived contaminants.

contiguous United States↗

USGS2015

An integrated approach to modeling changes in land use, land cover, and disturbance and their impact on ecosystem carbon dynamics: a case study in the Sierra Nevada Mountains of California

Increased land-use intensity (e.g. clearing of forests for cultivation, urbanization), often results in the loss of ecosystem carbon storage, while changes in productivity resulting from climate change may either help offset or exacerbate losses. However, there are large uncertainties in how land and climate systems will evolve and interact to shape future ecosystem carbon dynamics. To address this we developed the Land Use and Carbon Scenario Simulator (LUCAS) to track changes in land use, land cover, land management, and disturbance, and their impact on ecosystem carbon storage and flux within a scenario-based framework. We have combined a state-and-transition simulation model (STSM) of land change with a stock and flow model of carbon dynamics. Land-change projections downscaled from the Intergovernmental Panel on Climate Change’s (IPCC) Special Report on Emission Scenarios (SRES) were used to drive changes within the STSM, while the Integrated Biosphere Simulator (IBIS) ecosystem model was used to derive input parameters for the carbon stock and flow model. The model was applied to the Sierra Nevada Mountains ecoregion in California, USA, a region prone to large wildfires and a forestry sector projected to intensify over the next century. Three scenario simulations were conducted, including a calibration scenario, a climate-change scenario, and an integrated climate- and land-change scenario. Based on results from the calibration scenario, the LUCAS age-structured carbon accounting model was able to accurately reproduce results obtained from the process-based biogeochemical model. Under the climate-only scenario, the ecoregion was projected to be a reliable net sink of carbon, however, when land use and disturbance were introduced, the ecoregion switched to become a net source. This research demonstrates how an integrated approach to carbon accounting can be used to evaluate various drivers of ecosystem carbon change in a robust, yet transparent modeling environment.

California↗

USGS2012

Estimated probability of postwildfire debris flows in the 2012 Whitewater-Baldy Fire burn area, southwestern New Mexico

In May and June 2012, the Whitewater-Baldy Fire burned approximately 1,200 square kilometers (300,000 acres) of the Gila National Forest, in southwestern New Mexico. The burned landscape is now at risk of damage from postwildfire erosion, such as that caused by debris flows and flash floods. This report presents a preliminary hazard assessment of the debris-flow potential from 128 basins burned by the Whitewater-Baldy Fire. A pair of empirical hazard-assessment models developed by using data from recently burned basins throughout the intermountain Western United States was used to estimate the probability of debris-flow occurrence and volume of debris flows along the burned area drainage network and for selected drainage basins within the burned area. The models incorporate measures of areal burned extent and severity, topography, soils, and storm rainfall intensity to estimate the probability and volume of debris flows following the fire. In response to the 2-year-recurrence, 30-minute-duration rainfall, modeling indicated that four basins have high probabilities of debris-flow occurrence (greater than or equal to 80 percent). For the 10-year-recurrence, 30-minute-duration rainfall, an additional 14 basins are included, and for the 25-year-recurrence, 30-minute-duration rainfall, an additional eight basins, 20 percent of the total, have high probabilities of debris-flow occurrence. In addition, probability analysis along the stream segments can identify specific reaches of greatest concern for debris flows within a basin. Basins with a high probability of debris-flow occurrence were concentrated in the west and central parts of the burned area, including tributaries to Whitewater Creek, Mineral Creek, and Willow Creek. Estimated debris-flow volumes ranged from about 3,000-4,000 cubic meters (m 3 ) to greater than 500,000 m 3 for all design storms modeled. Drainage basins with estimated volumes greater than 500,000 m 3 included tributaries to Whitewater Creek, Willow Creek, Iron Creek, and West Fork Mogollon Creek. Drainage basins with estimated debris-flow volumes greater than 100,000 m 3 for the 25-year-recurrence event, 24 percent of the basins modeled, also include tributaries to Deep Creek, Mineral Creek, Gilita Creek, West Fork Gila River, Mogollon Creek, and Turkey Creek, among others. Basins with the highest combined probability and volume relative hazard rankings for the 25-year-recurrence rainfall include tributaries to Whitewater Creek, Mineral Creek, Willow Creek, West Fork Gila River, West Fork Mogollon Creek, and Turkey Creek. Debris flows from Whitewater, Mineral, and Willow Creeks could affect the southwestern New Mexico communities of Glenwood, Alma, and Willow Creek. The maps presented herein may be used to prioritize areas where emergency erosion mitigation or other protective measures may be necessary within a 2- to 3-year period of vulnerability following the Whitewater-Baldy Fire. This work is preliminary and is subject to revision. It is being provided because of the need for timely "best science" information. The assessment herein is provided on the condition that neither the U.S. Geological Survey nor the U.S. Government may be held liable for any damages resulting from the authorized or unauthorized use of the assessment.

New Mexico↗

USGS2016

National Park Service Vegetation Mapping Inventory Program: Natchez Trace Parkway vegetation mapping project report

The National Park Service (NPS) Vegetation Mapping Inventory (VMI) Program is an effort to classify, describe, and map existing vegetation of national park units for the NPS Natural Resource Inventory and Monitoring (I&M) Program. The NPS VMI Program is managed by the NPS I&M Division and provides baseline vegetation information to the NPS Natural Resource I&M Program. The U.S. Geological Survey Upper Midwest Environmental Sciences Center, NatureServe, NPS Gulf Coast Network, and NPS Natchez Trace Parkway (NATR; also referred to as Parkway) have completed vegetation classification and mapping of NATR for the NPS VMI Program. Mappers, ecologists, and botanists collaborated to affirm vegetation types within the U.S. National Vegetation Classification (USNVC) of NATR and to determine how best to map them by using aerial imagery. Analyses of data from 589 vegetation plots had been used to describe an initial 99 USNVC associations in the Parkway; this classification work was completed prior to beginning this NATR vegetation mapping project. Data were collected during this project from another eight quick plots to support new vegetation types not previously identified at the Parkway. Data from 120 verification sites were collected to test the field key to vegetation associations and the application of vegetation associations to a sample set of map polygons. Furthermore, data from 900 accuracy assessment (AA) sites were collected (of which 894 were used to test accuracy of the vegetation map layer). The collective of all these datasets resulted in affirming 122 USNVC associations at NATR. To map the vegetation and open water of NATR, 63 map classes were developed. including the following: 54 map classes represent natural (including ruderal) vegetation types in the USNVC, 5 map classes represent cultural (agricultural and developed) vegetation types in the USNVC, 3 map classes represent nonvegetation open-water bodies (non-USNVC), and 1 map class represents landscapes that had received tornado damage a few months prior to the time of aerial imagery collection. Features were interpreted from viewing 4-band digital aerial imagery by means of digital onscreen three-dimensional stereoscopic workflow systems in geographic information systems. (The aerial imagery was collected during mid-October 2011 for the northern reach of the Parkway and mid-November 2011 for the southern reach of the Parkway to capture peak leaf-phenology of trees.) The interpreted data were digitally and spatially referenced, thus making the spatial-database layers usable in geographic information systems. Polygon units were mapped to either a 0.5 hectare (ha) or 0.25 ha minimum mapping unit, depending on vegetation type or scenario. A geodatabase containing various feature-class layers and tables present the locations of USNVC vegetation types (vegetation map), vegetation plot samples, verification sites, AA sites, project boundary extent, and aerial image centers. The feature-class layer and related tables for the vegetation map provide 13,529 polygons of detailed attribute data covering 21,655.5 ha, with an average polygon size of 1.6 ha; the vegetation map coincides closely with the administrative boundary for NATR. Summary reports generated from the vegetation map layer of the map classes representing USNVC natural (including ruderal) vegetation types apply to 12,648 polygons (93.5% of polygons) and cover 18,542.7 ha (85.6%) of the map extent for NATR. The map layer indicates the Parkway to be 70.5% forest and woodland (15,258.7 ha), 0.3% shrubland (63.0 ha), and 14.9% herbaceous cover (3,221.0 ha). Map classes representing USNVC cultural types apply to 678 polygons (5.0% of polygons) and cover 2,413.9 ha (11.1%) of the map extent.

Natural Resource Report↗

USGS2023

Dissolved organic carbon dynamics and fluxes in Mississippi-Atchafalaya deltaic system impacted by an extreme flood event and hurricanes: A multi-satellite approach using Sentinel-2/3 and Landsat-8/9 data

Transport of riverine and wetland-derived dissolved organic carbon (DOC) spanning tidal wetlands, estuaries, and continental shelf waters functionally connects terrestrial and aquatic carbon reservoirs, yet the magnitude and ecological significance of this variable and its spatiotemporal linkage remains uncertain for coastal deltaic regions, such as Mississippi River Delta Plain, which includes Mississippi (MR) and Atchafalaya (AR) rivers and estuaries with vast expanses of wetlands and coastal forests. We examined DOC dynamics and fluxes in this large river-dominated wetland-estuarine system for the period between 2019 and 2021 that included an extreme river flood event in 2019, two major hurricanes (Barry in 2019 and Ida in 2021), and cold front passage using an improved adaptive quasi-analytical algorithm (QAA-AD) applied to multi-satellite sensors (Sentinel 3A/B OLCI, Landsat-8/OLI and Sentinel-2A/B MSI) with varying spectral and spatial (10/30/300 m) resolutions. The DOC estimates from multi-satellite sensors in combination with water fluxes were used to assess DOC fluxes from two large rivers (MR and AR) and small channels across the delta plain. Overall, this system delivered a total of 6.7 Tg C yr -1 (1 Tg = 10 12 g) into the estuarine zone and the northern Gulf of Mexico (nGoM) during 2019. High DOC fluxes from the AR (1.3 Tg C yr -1 ) and MR (4.5 Tg C yr -1 ) were associated with the extreme flood event in 2019. Hurricanes that occurred in the study period also contributed to the wetland and estuarine DOC fluxes into continental shelf waters; for example, the passage of Hurricane Barry in July 2019, delivered over a 3-day period ~1.33 ×10 9 g DOC from Barataria Basin into the nGoM. Sentinel 2-MSI land and water classification revealed that Hurricane Ida eroded a total of 1.34×10 8 m 2 of marshes in middle Barataria Basin, converting those habitats into open water with 3.0 m inundation depth and high DOC concentrations (16.4 mg L -1 ), a potentially large DOC source to the coastal waters. Overall, storms and flood events are major sources of DOC flux that facilitate transport of upstream carbon as well as transformation of carbon in the wetlands, through the conversion of vegetated wetland to open water.

Louisiana↗

USGS2021

National Park Service Vegetation Mapping Inventory Program: Great Smoky Mountains National Park vegetation mapping project

The National Park Service (NPS) Vegetation Mapping Inventory (VMI) Program is an effort to classify, describe, and map existing vegetation communities in national park units throughout the United States. The NPS VMI Program is managed by the NPS Natural Resource Stewardship and Science Inventory and Monitoring Program and provides baseline vegetation information to natural resource managers, researchers, and ecologists. The U.S. Geological Survey Upper Midwest Environmental Sciences Center, NatureServe, and NPS Great Smoky Mountains National Park (GRSM, also referred to as the “Park”) have completed vegetation classification and mapping of GRSM, including the Foothills Parkway, for the NPS VMI Program. Mappers, ecologists, and botanists collaborated to affirm vegetation types of GRSM and to determine how best to map the vegetation types by using aerial imagery. A vegetation classification developed in 2003 by NatureServe and the NPS served as a foundation to further classify and map the vegetation types of the Park. Data from an additional 10 vegetation plots supported vegetation types either rare or not documented in the 2003 classification. Data from 203 verification sites were collected to test the field key to vegetation types and the application of vegetation types to a sample set of map polygons. Furthermore, data from 972 accuracy assessment (AA) sites were collected (of which 966 were used to test accuracy of the vegetation map layer). This GRSM vegetation mapping project identified 112 vegetation types consisting of 105 association types in the U.S. National Vegetation Classification (USNVC), 2 “park-special” types, 1 “map-special” type, and 4 cultural types in the USNVC. To map the vegetation and land cover of GRSM, 52 map classes were developed. Of these 52 map classes, 46 represent natural (including ruderal) vegetation types, most of which types are recognized in the USNVC. For the remaining 6 of the 52 map classes, 4 represent USNVC cultural types for agricultural and developed areas, and 2 represent non-USNVC types for nonvegetated open water and nonvegetated rock. Features were interpreted from viewing four-band digital aerial imagery using digital onscreen three-dimensional stereoscopic workflow systems in geographic information systems; digital aerial imagery was collected during September 23–October 30, 2015. The interpreted data were digitally and spatially referenced, thus making the spatial-database layers usable in a geographic information system. Polygon units were mapped to either a 0.5- or 0.25- hectare (ha) minimum mapping unit, depending on vegetation type. A geodatabase containing several feature-class layers and tables provides the locations and data of USNVC vegetation types (vegetation map layer), vegetation plots, verification sites, AA sites, project boundary extent, and aerial image centers and flight lines. Covering 210,875 ha, the feature-class layer and related tables for the vegetation map layer provide 34,084 polygons of detailed attribute data when special modifiers are not considered (average polygon size of 6.2 ha) and 36,589 polygons of detailed attribute data when special modifiers are considered (average polygon size of 5.8 ha). Each map polygon is assigned a map-class code and name and, when applicable, are linked to USNVC classification tables within the geodatabase. The vegetation map extent includes the administrative boundary for GRSM and the Foothills Parkway. A summary report, generated from the vegetation map layer, concludes that the 46 map classes representing natural (including ruderal) vegetation types apply to 99.2% of polygons (33,797 polygons; average size of 6.2 ha) and cover 98.6% of the Park (207,971.4 ha). Further broken down, map classes representing natural vegetation types indicate that the Park is 97.7% forest and woodland (205,882.5 ha), 0.6% shrubland (1,174.6 ha), and 0.4% herbaceous (914.3 ha). Map classes representing cultural vegetation types apply to 0.8% of polygons (259 polygons; average size of 4.9 ha) and cover 0.6% of the Park (1,277.4 ha). Map classes representing nonvegetation open and flowing water and unvegetated rock apply to 0.08% of polygons (28 polygons; average size of 58.1 ha) and cover 0.8% of the Park (1,625.9 ha). A thematic AA study was completed of map classes representing the natural (including ruderal) vegetation types of the Park. Initial AA results were discussed with NPS staff from the Park. Following input from NPS staff on how to handle map classes that fell below accuracy standards, adjustments were made to the vegetation map layer. Final results indicate an overall accuracy of 80.64% (kappa index of 79.96% for chance agreements) based on data from 966 of the 972 AA sites. Most individual map-class themes exceed the NPS VMI Program standard of 80% with a 90% confidence interval. The GRSM vegetation mapping project delivers many geospatial and vegetation data products, including an in-depth project report discussing methods and results, which includes map classification and map-class descriptions. This suite of products also includes descriptions and a field key to vegetation types; a database of vegetation plots, verification sites, and AA sites; digital images of field sites; field data sheets; digital aerial imagery; hardcopy and digital maps; a geodatabase of vegetation and land cover (map layer), field sites (vegetation plots, verification sites, and AA sites), aerial imagery index, project boundary, and metadata; and a contingency table listing AA results. Geospatial products are projected in the Universal Transverse Mercator, Zone 17 North, by using the North American Datum of 1983. Information on the NPS VMI Program and completed mapping projects are on the internet at https://www.nps.gov/im/vegetation-inventory.htm.

North Carolina, Tennessee↗

USGS1996

Chemical evaluation of soil-solution in acid forest soils

Soil-solution chemistry is commonly studied in forests through the use of soil lysimeters.This approach is impractical for regional survey studies, however, because lysimeter installation and operation is expensive and time consuming. To address these problems, a new technique was developed to compare soil-solution chemistry among red spruce stands in New York, Vermont, New Hampshire, Maine. Soil solutions were expelled by positive air pressure from soil that had been placed in a sealed cylinder. Before the air pressure was applied, a solution chemically similar to throughfall was added to the soil to bring it to approximate field capacity. After the solution sample was expelled, the soil was removed from the cylinder and chemically analyzed. The method was tested with homogenized Oa and Bs horizon soils collected from a red spruce stand in the Adirondack Mountains of New York, a red spruce stand in east-central Vermont, and a mixed hardwood stand in the Catskill Mountains of New York. Reproducibility, effects of varying the reaction time between adding throughfall and expelling soil solution (5-65 minutes) and effects of varying the chemical composition of added throughfall, were evaluated. In general, results showed that (i) the method was reproducible (coefficients of variation were generally < 15%), (ii) variations in the length of reaction-time did not affect expelled solution concentrations, and (iii) adding and expelling solution did not cause detectable changes in soil exchange chemistry. Concentrations of expelled solutions varied with the concentrations of added throughfall; the lower the CEC, the more sensitive expelled solution concentrations were to the chemical concentrations of added throughfall. Addition of a tracer (NaBr) showed that the expelled solution was a mixture of added solution and solution that preexisted in the soil. Comparisons of expelled solution concentrations with concentrations of soil solutions collected by zero-tension and tension lysimetry indicated that expelled solution concentrations were higher than those obtained with either type of lysimeter, although there was less difference with tension lysimeters than zero-tension lysimeters. The method used for collection of soil solution should be taken into consideration whenever soil solution data are being interpreted.

New York, Vermont↗

USGS2013

Sources and characteristics of organic matter in the Clackamas River, Oregon, related to the formation of disinfection by-products in treated drinking water

This study characterized the amount and quality of organic matter in the Clackamas River, Oregon, to gain an understanding of sources that contribute to the formation of chlorinated and brominated disinfection by-products (DBPs), focusing on regulated DBPs in treated drinking water from two direct-filtration treatment plants that together serve approximately 100,000 customers. The central hypothesis guiding this study was that natural organic matter leaching out of the forested watershed, in-stream growth of benthic algae, and phytoplankton blooms in the reservoirs contribute different and varying proportions of organic carbon to the river. Differences in the amount and composition of carbon derived from each source affects the types and concentrations of DBP precursors entering the treatment plants and, as a result, yield varying DBP concentrations and species in finished water. The two classes of DBPs analyzed in this study-trihalomethanes (THMs) and haloacetic acids (HAAs)-form from precursors within the dissolved and particulate pools of organic matter present in source water. The five principal objectives of the study were to (1) describe the seasonal quantity and character of organic matter in the Clackamas River; (2) relate the amount and composition of organic matter to the formation of DBPs; (3) evaluate sources of DBP precursors in the watershed; (4) assess the use of optical measurements, including in-situ fluorescence, for estimating dissolved organic carbon (DOC) concentrations and DBP formation; and (5) assess the removal of DBP precursors during treatment by conducting treatability "jar-test" experiments at one of the treatment plants. Data collection consisted of (1) monthly sampling of source and finished water at two drinking-water treatment plants; (2) event-based sampling in the mainstem, tributaries, and North Fork Reservoir; and (3) in-situ continuous monitoring of fluorescent dissolved organic matter (FDOM), turbidity, chlorophyll- a , and other constituents to continuously track source-water conditions in near real-time. Treatability tests were conducted during the four event-based surveys to determine the effectiveness of coagulant and powdered activated carbon (PAC) on the removal of DBP precursors. Sample analyses included DOC, total particulate carbon (TPC), total and dissolved nutrients, absorbance and fluorescence spectroscopy, and, for regulated DBPs, concentrations of THMs and HAAs in finished water and laboratory-based THM and HAA formation potentials (THMFP and HAAFP, respectively) for source water and selected locations throughout the watershed. The results of this study may not be typical given the record and near record amounts of precipitation that occurred during spring that produced streamflow much higher than average in 2010-11. Although there were algal blooms, lower concentrations of chlorophyll- a were observed in the water column during the study period compared to historical data. Concentrations of DBPs in finished (treated) water averaged 0.024 milligrams per liter (mg/L) for THMs and 0.022 mg/L for HAAs; maximum values were about 0.040 mg/L for both classes of DBPs. Although DBP concentrations were somewhat higher within the distribution system, none of the samples collected for this study or for the quarterly compliance monitoring by the water utilities exceeded levels permissible under existing U.S. Environmental Protection Agency (USEPA) regulations: 0.080 mg/L for THMs and 0.060 mg/L for HAAs. DOC concentrations were generally low in the Clackamas River, typically about 1.0-1.5 mg/L. Concentrations in the mainstem occasionally increased to nearly 2.5 mg/L during storms; DOC concentrations in tributaries were sometimes much higher (up to 7.8 mg/L). The continuous in-situ FDOM measurements indicated sharp rises in DOC concentrations in the mainstem following rainfall events; concentrations were relatively stable during summer base flow. Even though the first autumn storm mobilized appreciable quantities of carbon, higher concentrations of DBPs in finished water were observed 3-weeks later, after the ground was saturated from additional rainfall. The majority of the DOC in the lower Clackamas River appears to originate from the upper basin, suggesting terrestrial carbon was commonly the dominant source. Lower-basin tributaries typically contained the highest concentrations of DOC and DBP precursors and contributed substantially to the overall loads in the mainstem during storms. During low-flow periods, tributaries were not major sources of DOC or DBP precursors to the Clackamas River. Although the dissolved fraction of organic carbon contributed the majority of DBP precursors, at times the particulate fraction (inorganic sediment and organic particles including detritus and algal material) contributed a substantial fraction of DBP precursors. Considering just the main-stem sites, on average, 10 percent of THMFP and 32 percent of HAAFP were attributed to particulate carbon. This finding suggests water-treatment methods that remove particles prior to chlorination would reduce finished-water DBP concentrations to some degree. Overall, concentrations of THM and HAA precursors were closely linked to DOC concentrations; laboratory DBP formation potentials (DBPFPs) clearly showed that THMFP and HAAFP were greatest in the downstream tributaries that contained elevated carbon concentrations. However, carbon-normalized "specific" formation potentials for THMs and HAAs (STHMFP and SHAAFP, respectively) revealed changes in carbon character over time that affected the two types of DBP classes differently. HAA precursors were elevated in waters containing aromatic-rich soil-derived material arising from forested areas. In contrast, THM precursors were associated with carbon having a lower aromatic content; highest STHMFP occurred in autumn 2011 in the mainstem from North Fork Reservoir downstream to LO DWTP. This pattern suggests the potential for a link between THM precursors and algal-derived carbon. The highest STHMFP value was measured within North Fork Reservoir, indicating reservoir derived carbon may be important for this class of DBPs. Weak correlations between STHMFP and SHAAFP emphasize that precursor sources for these types of DBPs may be different. This highlights not only that different locations within the watershed produce carbon with different reactivity (specific DBPFP), but also that different management approaches for each class of DBP precursors could be required for control. Treatability tests conducted on source water during four basin-wide surveys demonstrated that an average of about 40 percent of DOC can be removed by coagulation. While the decrease in THMFP following coagulation was similar to DOC, the decrease in HAAFP was much greater (approximately 70 percent), indicating coagulation is particularly effective at removing HAA precursors'likely because of the aromatic nature of the carbon associated with HAA precursors. Several findings from this study have direct implications for managing drinking-water resources and for providing useful information that may help improve treatment-plant operations. For example, the use of in-situ fluorometers that measure FDOM provided an excellent proxy for DOC concentration in this system and revealed short-term, rapid changes in DOC concentration during storm events. In addition, the strong correlation between FDOM values measured in-situ and HAA5 concentrations in finished water may permit estimation of continuous HAA concentrations, as was done here. As part of this study, multiple in-situ FDOM sensors were deployed continuously and in real-time to characterize the composition of dissolved organic matter. Although the initial results were promising, additional research and engineering developments will be needed to demonstrate the full utility of these sensors for this purpose. In conclusion, although DBPFPs were strongly correlated to DOC concentration, some DBPs formed from particulate carbon, including terrestrial leaf material and algal material such as planktonic species of blue-green algae and sloughed filaments, stalks, and cells of benthic algae. Different precursor sources in the watershed were evident from the data, suggesting specific actions may be available to address some of these sources. In-situ measurements of FDOM proved to be an excellent proxy for DOC concentration as well as HAA formation during treatment, which suggests further development and refinement of these sensors have the potential to provide real-time information about complex watershed processes to operators at the drinking-water treatment plants. Follow-up studies could examine the relative roles that terrestrial and algal sources have on the DBP precursor pool to better understand how watershed-management activities may be affecting the transport of these compounds to Clackamas River drinking-water intakes. Given the low concentrations of algae in the water column during this study, additional surveys during more typical river conditions could provide a more complete understanding of how algae contribute DBP precursors. Further development of FDOM-sensor technology can improve our understanding of carbon dynamics in the river and how concentrations may be trending over time. This study was conducted in collaboration with Clackamas River Water and the City of Lake Oswego water utilities. Other research partners included Oregon Health and Science University in Hillsboro, Oregon, Alexin Laboratory in Tigard, Oregon, U.S. Geological Survey National Research Program Laboratory in Denver, Colorado, and the U.S. Geological Survey Water Science Centers in Portland, Oregon, and Sacramento, California. This project was supported with funding from Clackamas River Water, City of Lake Oswego, the U.S. Geological Survey, and the Water Research Foundation.

Oregon↗