Geology Reports⌕ Search

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.

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

At least 649 records · Page 36Linked to original sources

Predicting the spatiotemporal exposure of aquatic species to intrusions of fire retardant in streams with limited data

Because fire retardant can enter streams and harm aquatic species including endangered fish, agencies such as the U.S. Forest Service (USFS) must estimate the downstream extent of toxic effects every time fire retardant enters streams (denoted as an “intrusion”). A challenge in estimating the length of stream affected by the intrusion and the exposure time of species in the affected reach is the lack of data typically available on the stream's geometry and flow characteristics. Previously, the USFS estimated the affected reach length assuming instantaneous mixing of the retardant over the reach; however, this approach neglects key river mixing processes. An approach is described that accounts for advection and dispersion of the retardant as well as the downstream growth of the stream. Applied to 13 intrusions documented by the USFS, the new approach shows affected reach lengths range between 8.0 and 362 km; all 13 cases exceeded previous estimates from an instantaneous mixing model. The time that a stationary individual in the affected reach is exposed to concentrations above a pre-defined toxicity threshold (10% of 96-hour LC 50 , for example) ranges from 0.17 to 2.73 h, with all but one case having a maximum exposure time less than 1.5 h. Results from 1152 hypothetical intrusions provided by the USFS confirm that exposure times rarely exceed 5 h. This result suggests that 96-hour tests to determine toxicity (LC 50 ) to various species should be reconsidered. Although the approach described can be improved in several ways, it provides a first estimate of the effects of fire retardant intrusions.

Science of the Total Environment↗

Sources and temporal dynamics of arsenic in a New Jersey watershed, USA

We examined potential sources and the temporal dynamics of arsenic (As) in the slightly alkaline waters of the Wallkill River, northwestern New Jersey, where violations of water-quality standards have occurred. The study design included synoptic sampling of stream water and bed sediments in tributaries and the mainstem, hyporheic-zone/ground water on the mainstem, and seasonal and diurnal sampling of water at selected mainstem sites. The river valley is bordered by gneiss and granite highlands and shale lowlands and underlain by glacial deposits over faulted dolomites and the Franklin Marble. Ore bodies in the Marble, which have been mined for rare Zn ore minerals, also contain As minerals. Tributaries, which drain predominantly forested and agricultural land, contributed relatively little As to the river. The highest concentrations of As (up to 34 μg/L) emanated from the outlet of man-made Lake Mohawk at the river's headwaters; these inputs varied substantially with season—high during warm months, low during cold months, apparently because of biological activity in the lake. Dissolved As concentrations were lower (3.3 μg/L) in river water than those in ground water discharging into the riverbed (22 μg/L) near the now-closed Franklin Mine. High total As concentrations (100–190 mg/kg) on the < 0.63 μm fraction of bed sediments near the mine apparently result from sorption of the As in the ground-water discharge as well as from the As minerals in the streambed. As concentrations in river water were diluted during high stream flow in fall, winter and spring, and concentrated during low flow in summer. In unfiltered samples from a wetlands site, diurnal cycles in trace-element concentrations occurred; As concentrations appeared to peak during late afternoon as pH increased, but Fe, Mn, and Zn concentrations peaked shortly after midnight. The temporal variability of As and its presence at elevated concentrations in ground water and sediments as well as streamwater demonstrate the importance of (1) sampling a variety of media and (2) determining the time scales of As variability to fully characterize its passage through a river system.

New Jersey↗

Late-Glacial paleoecology of the Middle Susitna Valley, Alaska: Environmental context for human dispersal

We present here the results of multi-proxy analyses (sediment geochemistry, diatoms, and pollen) from sediment cores collected at four lakes in the middle Susitna Valley, Alaska. These lakes form a transect from the tundra to the boreal forest. The retrieved cores span from ~12,000 cal yr BP to the present, with age control provided by radiometric dates and tephra deposits, some of which are newly identified. Results indicate that deglaciation occurred before 12,000 cal yr BP and that by that time, the lakes were deep, productive, and surrounded by shrub tundra. The lake with the highest sampling resolution indicates a brief climatic reversal ~11,500 cal yr BP with decreased diatom-inferred lake level and lowered lake productivity, and reduced shrub presence. During the early to middle Holocene, all of the sedimentary records provide evidence of climatic amelioration with tree expansion and productive lakes. A middle to late Holocene climatic deterioration with reduced trees and a shallower, less productive lake is also indicated. In addition, the prominent Watana tephra at ~4,000 cal yr BP likely reduced lake productivity and affected the vegetation. Even though the region was relatively productive soon after deglaciation, people did not occupy the region until ~11,000 cal yr BP, about 1000 years later, and then only sparsely. By the middle and late Holocene, the region was more densely populated and this shift in human occupancy presumably reflects changes in resource abundance, especially caribou. Whether the Watana ashfall influenced caribou abundance and thus people, is still under investigation, but given the tephra’s effect on vegetation and lake productivity, it seems likely.

Alaska↗

National Park Service Vegetation Inventory Program, Cuyahoga Valley National Park, Ohio

The National Park Service (NPS) Vegetation Inventory Program (VIP) 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 VIP is managed by the NPS Biological Resources Management Division and provides baseline vegetation information to the NPS Natural Resource I&M Program. The U.S. Geological Survey (USGS) Vegetation Characterization Program lends a cooperative role in the NPS VIP. The USGS Upper Midwest Environmental Sciences Center, NatureServe, and NPS Cuyahoga Valley National Park (CUVA) have completed vegetation classification and mapping of CUVA. Mappers, ecologists, and botanists collaborated to identify and describe vegetation types within the National Vegetation Classification Standard (NVCS) and to determine how best to map them by using aerial imagery. The team collected data from 221 vegetation plots within CUVA to develop detailed descriptions of vegetation types. Data from 50 verification sites were also collected to test both the key to vegetation types and the application of vegetation types to a sample set of map polygons. Furthermore, data from 647 accuracy assessment (AA) sites were collected (of which 643 were used to test accuracy of the vegetation map layer). These data sets led to the identification of 45 vegetation types at the association level in the NVCS at CUVA. A total of 44 map classes were developed to map the vegetation and general land cover of CUVA, including the following: 29 map classes represent natural/semi-natural vegetation types in the NVCS, 12 map classes represent cultural vegetation (agricultural and developed) in the NVCS, and 3 map classes represent non-vegetation features (open-water bodies). Features were interpreted from viewing color-infrared digital aerial imagery dated October 2010 (during peak leaf-phenology change of trees) via digital onscreen three-dimensional stereoscopic workflow systems in geographic information systems (GIS). 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 ha or 0.25 ha minimum mapping unit, depending on vegetation type. A geodatabase containing various feature-class layers and tables shows the locations of vegetation types and general land cover (vegetation map), vegetation plot samples, verification sites, AA sites, project boundary extent, and aerial photographic centers. The feature-class layer and relate tables for the CUVA vegetation map provides 4,640 polygons of detailed attribute data covering 13,288.4 ha, with an average polygon size of 2.9 ha. Summary reports generated from the vegetation map layer show map classes representing natural/semi-natural types in the NVCS apply to 4,151 polygons (89.4% of polygons) and cover 11,225.0 ha (84.5%) of the map extent. Of these polygons, the map layer shows CUVA to be 74.4% forest (9,888.8 ha), 2.5% shrubland (329.7 ha), and 7.6% herbaceous vegetation cover (1,006.5 ha). Map classes representing cultural types in the NVCS apply to 435 polygons (9.4% of polygons) and cover 1,825.7 ha (13.7%) of the map extent. Map classes representing non-NVCS units (open water) apply to 54 polygons (1.2% of polygons) and cover 237.7 ha (1.8%) of the map extent. A thematic AA study was conducted of map classes representing natural/semi-natural types in the NVCS. Results present an overall accuracy of 80.7% (kappa index of 79.5%) based on data from 643 of the 647 AA sites. Most individual map-class themes exceed the NPS VIP standard of 80% with a 90% confidence interval. The CUVA vegetation mapping project delivers many geospatial and vegetation data products in hardcopy and/or digital formats. These products consist of an in-depth project report discussing methods and results, which include descriptions and a dichotomous key to vegetation types, map classification and map-class descriptions, and a contingency table showing AA results. The suite of products also includes a database of vegetation plots, verification sites, and AA sites; digital pictures of field sites; field data sheets; aerial photographic imagery; hardcopy and digital maps; and a geodatabase of vegetation types and land cover (map layer), fieldwork locations (vegetation plots, verification sites, and AA sites), aerial photographic index, project boundary, and metadata. All geospatial products are projected in Universal Transverse Mercator, Zone 17, by using the North American Datum of 1983. Information on the NPS VIP and completed park mapping projects are located on the Internet at <http://science.nature.nps.gov/im/inventory/veg> and <http://www.usgs.gov/core_science_systems/csas/vip/index.html>.

Natural Resource Technical Report↗

Short- and long-term responses of riparian cottonwoods (Populus spp.) to flow diversion: Analysis of tree-ring radial growth and stable carbon isotopes

Long duration tree-ring records with annual precision allow for the reconstruction of past growing conditions. Investigations limited to the most common tree-ring proxy of ring width can be difficult to interpret, however, because radial growth is affected by multiple environmental processes. Furthermore, studies of living trees may miss important effects of drought on tree survival and forest changes. Stable carbon isotopes can help distinguish drought from other environmental factors that influence tree-ring width and forest stand condition. We quantified tree-ring radial expansion and stable carbon isotope ratios (δ 13 C) in riparian cottonwoods ( Populus angustifolia and P. angustifolia x P. trichocarpa) along Snake Creek in Nevada, USA. We investigated how hydrological drought affected tree growth and death at annual to half-century scales in a partially dewatered reach (DW) compared to reference reaches immediately upstream and downstream. A gradual decline in tree-ring basal area increment (BAI) began at DW concurrent to streamflow diversion in 1961. BAI at DW diverged from one reference reach immediately but not from the other until nearly 50 years later. In contrast, tree-ring δ 13 C had a rapid and sustained increase following diversion at DW only, providing the stronger and clearer drought signal. BAI and δ 13 C were not significantly correlated prior to diversion; after diversion they both reflected drought and were correlated for DW trees only. Cluster analyses distinguished all trees in DW from those in reference reaches based on δ 13 C, but BAI patterns left trees intermixed across reaches. Branch and tree mortality were also highest and canopy vigor was lowest in DW. Results indicate that water scarcity strongly limited cottonwood photosynthesis following flow diversion, thus reducing carbon assimilation, basal growth and survival. The dieback was not sudden, but occurred over decades as carbon deficits mounted and depleted streamflow left trees increasingly vulnerable to local meteorological drought.

Nevada↗

Mapping global coral vulnerability to stony coral tissue loss disease: Implications for biosecurity and conservation

Stony Coral Tissue Loss Disease (SCTLD) has devastated Caribbean coral reefs since 2014, but its potential for global impact remains uncertain. We developed predictive models to assess the worldwide vulnerability of coral reefs to SCTLD under different origin and spread hypotheses. Using random forest regression models incorporating coral taxonomy and zooxanthellae clade associations from 52 taxa, we projected SCTLD susceptibility and mortality patterns globally using six indices: Mean susceptibility per genus per location, Summed susceptibilities across genera per location, Summed susceptibilities across genera per realm, Mean mortality per genus per location, Summed mortalities across genera per location, and Summed mortalities across genera per realm. Models demonstrated strong predictive performance (R² = 0.57 for susceptibility; R² = 0.73 for mortality) and revealed that about 7% of coral genera per location are potentially susceptible to SCTLD. While mean susceptibility and mortality per genus were highest in the Tropical Atlantic, the summed susceptibility and mortality across genera were much higher in the biodiverse Central Indo-Pacific. Natural barriers could limit SCTLD’s spread, including the mid-Atlantic gap and the low diversity of the Tropical Eastern Pacific, supporting the contained disease hypothesis. However, the widespread distribution of susceptible genera across coral reef realms indicates significant vulnerability should SCTLD circumvent these barriers through human-mediated transport, particularly via ballast water or the aquarium trade. If SCTLD is an invasive pathogen originating in the Pacific, as shipping patterns for the aquarium trade suggest, mortality in its native range would likely be lower than our projections. These findings point to targeted intervention strategies, including enhanced monitoring at key locations, assessment of biosecurity needs in high-risk areas, and prioritized conservation efforts in vulnerable high-diversity regions to prevent SCTLD from spreading globally.

Frontiers in Marine Science↗

Patterns and predictability in the intra-annual organic carbon variability across the boreal and hemiboreal landscape

Factors affecting total organic carbon (TOC) concentrations in 215 watercourses across Sweden were investigated using parameter parsimonious regression approaches to explain spatial and temporal variabilities of the TOC water quality responses. We systematically quantified the effects of discharge, seasonality, and long-term trend as factors controlling intra-annual (among year) and inter-annual (within year) variabilities of TOC by evaluating the spatial variability in model coefficients and catchment characteristics (e.g. land cover, retention time, soil type). Catchment area (0.18–47,000 km 2 ) and land cover types (forests, agriculture and alpine terrain) are typical for the boreal and hemiboreal zones across Fennoscandia. Watercourses had at least 6 years of monthly water quality observations between 1990 and 2010. Statistically significant models (p < 0.05) describing variation of TOC in streamflow were identified in 209 of 215 watercourses with a mean Nash-Sutcliffe efficiency index of 0.44. Increasing long-term trends were observed in 149 (70%) of the watercourses, and intra-annual variation in TOC far exceeded inter-annual variation. The average influences of the discharge and seasonality terms on intra-annual variations in daily TOC concentration were 1.4 and 1.3 mg l − 1 (13 and 12% of the mean annual TOC), respectively. The average increase in TOC was 0.17 mg l − 1 year − 1 (1.6% year − 1 ). Multivariate regression with over 90 different catchment characteristics explained 21% of the spatial variation in the linear trend coefficient, less than 20% of the variation in the discharge coefficient and 73% of the spatial variation in mean TOC. Specific discharge, water residence time, the variance of daily precipitation, and lake area, explained 45% of the spatial variation in the amplitude of the TOC seasonality. Because the main drivers of temporal variability in TOC are seasonality and discharge, first-order estimates of the influences of climatic variability and change on TOC concentration should be predictable if the studied catchments continue to respond similarly.

Science of the Total Environment↗

Temperature response of soil respiration largely unaltered with experimental warming

The respiratory release of carbon dioxide (CO 2 ) from soil is a major yet poorly understood flux in the global carbon cycle. Climatic warming is hypothesized to increase rates of soil respiration, potentially fueling further increases in global temperatures. However, despite considerable scientific attention in recent decades, the overall response of soil respiration to anticipated climatic warming remains unclear. We synthesize the largest global dataset to date of soil respiration, moisture, and temperature measurements, totaling >3,800 observations representing 27 temperature manipulation studies, spanning nine biomes and over 2 decades of warming. Our analysis reveals no significant differences in the temperature sensitivity of soil respiration between control and warmed plots in all biomes, with the exception of deserts and boreal forests. Thus, our data provide limited evidence of acclimation of soil respiration to experimental warming in several major biome types, contrary to the results from multiple single-site studies. Moreover, across all nondesert biomes, respiration rates with and without experimental warming follow a Gaussian response, increasing with soil temperature up to a threshold of ∼25 °C, above which respiration rates decrease with further increases in temperature. This consistent decrease in temperature sensitivity at higher temperatures demonstrates that rising global temperatures may result in regionally variable responses in soil respiration, with colder climates being considerably more responsive to increased ambient temperatures compared with warmer regions. Our analysis adds a unique cross-biome perspective on the temperature response of soil respiration, information critical to improving our mechanistic understanding of how soil carbon dynamics change with climatic warming.

Proceedings of the National Academy of Sciences of↗

Forest soil chemistry and terrain attributes in a Catskills watershed

Knowledge of soil chemistry is useful in assessing the sensitivity of forested areas to natural and anthropogenic disturbances, but characterizing large areas is expensive because of the large sample numbers required and the cost of soil chemical analyses. We collected and chemically analyzed soil samples from 72 sites within a 214-ha watershed in the Catskill Mountains of New York to evaluate factors that influence soil chemistry and whether terrain features could be used to predict soil chemical properties. Using geographic information system (GIS) techniques, we determined five terrain attributes at each sampling location: (i) slope, (ii) aspect, (iii) elevation, (iv) topographic index, and (v) flow accumulation. These attributes were ineffective in predicting the chemical properties of organic and mineral soil samples; together they explained only 4 to 25% of the variance in pH(w), effective cation-exchange capacity (CEC(e)), exchangeable bases, exchangeable acidity, total C, total N, and C/N ratio. Regressions among soil properties were much better; total C and pH(w) together explained 33 to 66% of the variation in exchangeable bases and CEC(e). Total C was positively correlated with N (r = 0.91 and 0.96 in Oa horizons and mineral soil, respectively), exchangeable bases (r = 0.65, 0.76), and CEC(e) (r = 0.54, 0.44), indicating the importance of organic matter to the chemistry of these acidic soils. The fraction of CEC(e) occupied by H explained 44% of the variation in pH(w). Soil chemical properties at this site vary on spatial scales finer than typical GIS analyses, resulting in relationships with poor predictive power. Thus, interrelationships among soil properties are more reliable for prediction.Knowledge of soil chemistry is useful in assessing the sensitivity of forested areas to natural and anthropogenic disturbances, but characterizing large areas is expensive because of the large sample numbers required and the cost of soil chemical analyses. We collected and chemically analyzed soil samples from 72 sites within a 214-ha watershed in the Catskill Mountains of New York to evaluate factors that influence soil chemistry and whether terrain features could be used to predict soil chemical properties. Using geographic information system (GIS) techniques, we determined five terrain attributes at each sampling location: (i) slope, (ii) aspect, (iii) elevation, (iv) topographic index, and (v) flow accumulation. These attributes were ineffective in predicting the chemical properties of organic and mineral soil samples; together they explained only 4 to 25% of the variance in pHw, effective cation-exchange capacity (CECe), exchangeable bases, exchangeable acidity, total C, total N, and C/N ratio. Regressions among soil properties were much better; total C and pHw together explained 33 to 66% of the variation in exchangeable bases and CECe. Total C was positively correlated with N (r = 0.91 and 0.96 in Oa horizons and mineral soil, respectively), exchangeable bases (r = 0.65, 0.76), and CECe (r = 0.54, 0.44), indicating the importance of organic matter to the chemistry of these acidic soils. The fraction of CECe occupied by H explained 44% of the variation in pHw. Soil chemical properties at this site vary on spatial scales finer than typical GIS analyses, resulting in relationships with poor predictive power. Thus, interrelationships among soil properties are more reliable for prediction.

New York↗

A comparison of litter production in young and old baldcypress (Taxodium distichum L.) stands at Caddo Lake, Texas

Aboveground primary productivity for cypress forests was assessed from measurements of litter production in two age groups and in two hydrological regimes (standing water and free-flowing). Caddo Lake, located in northeast Texas on the Texas-Louisiana border, offered a unique study site since it is dominated by extensive stands composed entirely of Taxodium distichum (L.) Rich, (baldcypress) in different age groups. Young stands (approximately 100 years old) are found along the shoreline and on shallow flooded islands. Old stands (-150 to 300 years old) are found in deeper water where they were continuously flooded. Litter production over three years from October 1998 to September 2001 was measured. Litter consisting of leaves, twigs, bark, reproductive parts, and Tillandsia usneoides (L.) L. (Spanish moss) was collected monthly using 0.5 m 2 floating traps. Tree diameters were measured within 200 m 2 circular plots in each stand. The young stands supported densities greater than 2,000 stems/ha and a mean stand basal area of 72.3 m 2 /ha, whereas old stands supported lower densities of about 500 stems/ha but with a similar mean stand basal area of 73.3 m 2 /ha. There was a significant difference between old and young stands for overall yearly litter production, averaging about 670 g/m 2 /yr in the young stands and 460 g/m 2 /yr in the old stands. Leaves and twigs were significantly greater in the young stands, while reproductive parts were higher in old stands. Litter collections between years or hydrological regimes were not significantly different.

Texas↗

Predicting regional fluoride concentrations at public and domestic supply depths in basin-fill aquifers of the western United States using a random forest model

A random forest regression (RFR) model was applied to over 12,000 wells with measured fluoride (F) concentrations in untreated groundwater to predict F concentrations at depths used for domestic and public supply in basin-fill aquifers of the western United States. The model relied on twenty-two regional-scale environmental and surficial predictor variables selected to represent factors known to control F concentrations in groundwater. The testing model fit R 2 and RMSE were 0.52 and 0.78 mg/L. Comparisons of measured to predicted proportions of four F-concentrations categories (<0.7 mg/L, 0.7–2 mg/L, >2 mg/L – 4 mg/L, and > 4 mg/L) indicate that the model performed well at making regional-scale predictions. Differences between measured and predicted proportions indicate underprediction of measured F at values by between 4 and 20 mg/L, representing less than 1% of the regional scale predicted values. These residuals most often map to geographic regions where local-scale processes including evaporative discharge in closed basins or intermittent streams concentrate fluoride in shallow groundwater. Despite this, the RFR model provides spatially continuous F predictions across the basin-fill aquifers where discrete samples are missing. Further, the predictions capture documented areas that exceed the F maximum contaminant level for drinking water of 4 mg/L and areas that are below the oral-health benchmark of 0.7 mg/L. These predictions can be used to estimate fluoride concentrations in unmonitored areas and to aid in identifying geographic areas that may require further investigation at localized scales.

Arizona, California, Colorado, New Mexico, Nevada,↗

Wildfire-driven changes in hydrology mobilize arsenic and metals from legacy mine waste

Wildfires burning in watersheds that have been mined and since revegetated pose unique risks to downstream water supplies. A wildfire near Boulder, Colorado that burned a forested watershed recovering from mining disturbance that occurred 80-160 years ago allowed us to 1) assess arsenic and metal contamination in streams draining the burned area for a five-year period after the wildfire and 2) determine the fire-affected hydrologic drivers that convey arsenic and metals to surface water. Most metal concentrations were low in the circumneutral waters draining the burned area. Water and sediment collected from streams downstream of the burned area had elevated arsenic concentrations during and after post-fire storms. Mining-related deposits were the main source of arsenic to streams. An increased proportion of overland flow relative to infiltration after the fire mobilized arsenic- and metal- rich surface deposits and wildfire ash into streams within and downstream of the burned area. The deposition of this sediment into stream channels resulted in the remobilization of arsenic for the five-year post-fire study period. It is also possible that enhanced subsurface flow after the fire increased contact of water with arsenic-bearing minerals exposed in underground mine workings. Other studies have reported that wildfire ash can be an important source of arsenic and metals to surface waters, but wildfire ash was not an important source of arsenic in this study. Predicted increases in frequency, size, and intensity of wildfires in the western U.S., a region with widely dispersed historical mines, suggest that the intersection of legacy mining and post-wildfire hydrologic response poses an increasing risk for water supplies.

Colorado↗

Whitebark pine in the national parks of the Pacific states: An assessment of population vulnerability

Whitebark pine ( Pinus albicaulis ) is a long-lived tree found in high-elevation forests of western North America that is declining due to the non-native white pine blister rust ( Cronartium ribicola ) and climate-driven outbreaks of mountain pine beetle ( Dendroctonus ponderosae ; MPB). The National Park Service established a monitoring program for whitebark pine in seven parks, including Sequoia & Kings Canyon, Yosemite, Lassen Volcanic, Crater Lake, Mount Rainier, Olympic, and North Cascades National Parks. Using these data, we summarized stand structure, presence of blister rust, and MPB prevalence to provide a baseline for future monitoring. Next, we used a stochastic, size-structured population model to speculate on future trends in the seven national park populations under conditions of increased MPB activity and ongoing blister rust infection observed in Crater Lake. We found that blister rust infected 29 to 54% of whitebark pine in all the parks except the two southernmost, Sequoia & Kings Canyon and Yosemite, where infections rates were 0.3% and 0.2%, respectively. The proportion of dead trees in Sequoia & Kings Canyon and Yosemite was low (0 to 1%), while they ranged from 10 to 43% in the other parks. Model projections suggested an average population decline of 25% in the parks over the next century using Crater Lake conditions, declines which are possible if blister rust continues to spread and climate change results in a significant increase in the frequency or severity of MPB outbreaks. Overall, our study describes conditions at seven western parks and illustrates potential rates of whitebark pine decline if pest outbreaks and/or blister rust infections worsen.

California, Oregon, Washington↗

Relations among geology, physiography, land use, and stream habitat conditions in the Buffalo and Current River systems, Missouri and Arkansas

This study investigated links between drainage-basin characteristics and stream habitat conditions in the Buffalo National River, Arkansas and the Ozark National Scenic Riverways, Missouri. It was designed as an associative study - the two parks were divided into their principle tributary drainage basins and then basin-scale and stream-habitat data sets were gathered and compared between them. Analyses explored the relative influence of different drainage-basin characteristics on stream habitat conditions. They also investigated whether a relation between land use and stream characteristics could be detected after accounting for geologic and physiographic differences among drainage basins. Data were collected for three spatial scales: tributary drainage basins, tributary stream reaches, and main-stem river segments of the Current and Buffalo Rivers. Tributary drainage-basin characteristics were inventoried using a Geographic Information System (GIS) and included aspects of drainage-basin physiography, geology, and land use. Reach-scale habitat surveys measured channel longitudinal and cross-sectional geometry, substrate particle size and embeddedness, and indicators of channel stability. Segment-scale aerial-photo based inventories measured gravel-bar area, an indicator of coarse sediment load, along main-stem rivers. Relations within and among data sets from each spatial scale were investigated using correlation analysis and multiple linear regression. Study basins encompassed physiographically distinct regions of the Ozarks. The Buffalo River system drains parts of the sandstone-dominated Boston Mountains and of the carbonate-dominated Springfield and Salem Plateaus. The Current River system is within the Salem Plateau. Analyses of drainage-basin variables highlighted the importance of these physiographic differences and demonstrated links among geology, physiography, and land-use patterns. Buffalo River tributaries have greater relief, steeper slopes, and more streamside bluffs than the Current River tributaries. Land use patterns in both river systems correlate with physiography - cleared land area is negatively associated with drainage-basin average slope. Both river systems are dominantly forested (0-35 per-cent cleared land), however, the potential for landscape disturbance may be greater in the Buffalo River system where a larger proportion of cleared land occurs on steep slopes (>15 degrees). When all drainage basins are grouped together, reach-scale channel characteristics show the strongest relations with drainage-basin physiography. Bankfull channel geometry and residual pool dimensions are positively correlated with drainage area and topographic relief variables. After accounting for differences in drainage area, channel dimensions in Buffalo River tributaries tend to be larger than in Current River tributaries. This trend is consistent with the flashy runoff and large storm flows that can be generated in rugged, sandstone-dominate terrain. Substrate particle size is also most strongly associated with physiography; particle size is positively correlated with topographic relief variables. When tributaries are subset by river system, relations with geology and land use variables become apparent. Buffalo River tributaries with larger proportions of carbonate bedrock and cleared land area have shallower channels, better-sorted, gravel-rich substrate, and more eroding banks than those with little cleared land and abundant sandstone bedrock. Gravel-bar area on the Buffalo River main stem was also larger within 1-km of carbonate-rich tributary junctions. Because geology and cleared land are themselves correlated, relations with anthropogenic and natural factors could often not be separated. Channel characteristics in the Current River system show stronger associations with physiography than with land use. Channels are shallower and have finer substrates in the less rugged, karst-rich, western basins than in the

Arkansas, Missouri↗

Lessons learned from 20 y of monitoring suburban development with distributed stormwater management in Clarksburg, Maryland, USA

Urban development is a well-known stressor for stream ecosystems, presenting a challenge to managers tasked with mitigating its effects. For the past 20 y, streamflow, water quality, geomorphology, and benthic communities were monitored in 5 watersheds in Montgomery County, Maryland, USA. This study presents a synthesis of multiple studies of monitoring efforts in the study area and new analysis of more recent monitoring data to document the primary lessons learned from monitoring. The monitored watersheds include a forested control, an urban control with centralized stormwater management, and 3 suburban treatment watersheds featuring low-impact development and a high density of infiltration-focused stormwater facilities distributed across the watershed. Treatment watersheds were monitored before development, during construction, and after development. Monitoring was initiated to inform adaptive management of stormwater and impervious cover limits within the study area, with a focus on the impacts of distributed stormwater management. Results from our synthesis indicate that distributed stormwater management is advantageous compared with centralized stormwater management in numerous ways. Hydrologic benefits were greater with distributed stormwater infrastructure, demonstrating the ability to mitigate runoff volumes and peak flows and, for small storms, replicate predevelopment conditions. Baseflow temporarily increased during the construction phase in the treatment watersheds. Water-quality benefits were mixed, with declines in baseflow nitrate concentrations but limited changes to nitrate export and increases in specific conductance after development. Substantial topographic changes occurred during construction in the treatment watersheds, including changes within the riparian zone, despite riparian buffer protections. Ecological monitoring indicated that even though index of biotic integrity scores rebounded in some cases, sensitive benthic macroinvertebrate families did not fully recover in the treatment watersheds. Lessons learned from this synthesis highlight the importance of tracking multiple indicators of stream health and considering past land use and that more stormwater facilities distributed across the watershed is beneficial but cannot mitigate the effects of all urban stressors on aquatic ecosystems.

Maryland↗

The response of stream ecosystems in the Adirondack region of New York to historical and future changes in atmospheric deposition of sulfur and nitrogen

The present-day acid-base chemistry of surface waters can be directly linked to contemporary observations of acid deposition; however, pre-industrial conditions are key to predicting the potential future recovery of stream ecosystems under decreasing loads of atmospheric sulfur (S) and nitrogen (N) deposition. The integrated biogeochemical model PnET-BGC was applied to 25 forest watersheds that represent a range of acid sensitivity in the Adirondack region of New York, USA to simulate the response of streams to past and future changes in atmospheric S and N deposition, and calculate the target loads of acidity for protecting and restoring stream water quality and ecosystem health. Using measured data, the model was calibrated and applied to simulate soil and stream chemistry at all study sites. Model hindcasts indicate that historically stream water chemistry in the Adirondacks was variable, but inherently sensitive to acid deposition. The median model-simulated acid neutralizing capacity (ANC) of the streams was projected to be 55 μeq L −1 before the advent of anthropogenic acid deposition (~1850), decreasing to minimum values of 10 μeq L −1 around the year 2000. The median simulated ANC increased to 13 μeq L −1 by 2015 in response to decreases in acid deposition that have occurred over recent decades. Model projections suggest that simultaneous decreases in sulfate, nitrate and ammonium deposition are more effective in restoring stream ANC than individual decreases in sulfur or nitrogen deposition. However, the increases in stream ANC per unit equivalent decrease in S deposition is greater compared to decreases in N deposition. Using empirical algorithms, fish community density and biomass are projected to increase under several deposition-control scenarios that coincide with increases in stream ANC. Model projections suggest that even under the most aggressive deposition-reduction scenarios, stream chemistry and fisheries will not fully recover from historical acidification by 2200.

New York↗

Estimating the value of mangrove leaf litter in sesarmid crab diets: The importance of fractionation factors

Sesarmid crabs play an important role in organic matter and carbon cycling of mangrove forests. Visual observations and gut content studies have verified that sesarmid crabs are feeding on mangrove leaves, yet stable isotopes of carbon and nitrogen ( 13 C and 15 N) have indicated that leaf litter is not assimilated as a food source. Sesarmid crabs tend to be much more enriched in 13 C than leaf litter (0.9‰ – 11.6‰) and have C values that are often more like microphytobenthos (MPB). General 13 C trophic enrichment factors (TEF; 0.1‰ – 0.5‰) suggest crabs feed more heavily on MPB. Field and laboratory-based evidence reveal that general 13 C TEF for crabs feeding on mangrove leaves may be incorrect and much greater than 0.1‰ – 0.5‰. A food web study conducted annually over 2 yrs revealed a shift in the δ 13 C and δ 15 N of Parasesarma sp. crabs similar to mangrove leaves also sampled. This suggested Parasesarma sp. may be feeding more heavily on mangrove leaves than previously reported despite crabs being 4.4‰ – 11.6‰ more enriched in 13 C than mangrove leaves. A laboratory feeding study confirmed that average 13 C TEF between Parasesarma sp. and decayed Rhizophora sp. leaves was 3.3‰ (SE 0.5). The Stable Isotope Analysis in R package (SIAR) used with our TEF and the general 0.5‰ 13 C TEF revealed that published TEFs may underestimate mangrove leaf contributions to sesarmid crab diets on average by 33.3% (SE 0.1) and overestimate MPB and epiphytic algal contributions by 31.3% (SE 0.1). Food web studies in mangroves and other ecosystems will continue to inaccurately identify important food resources or food web structures unless more accurate and species-specific isotope fractionation values are determined.

Bulletin of Marine Science↗

Hysteretic response of solutes and turbidity at the event scale across forested tropical montane watersheds

Concentration-discharge relationships are a key tool for understanding the sourcing and transport of material from watersheds to fluvial networks. Storm events in particular provide insight into variability in the sources of solutes and sediment within watersheds, and the hydrologic pathways that connect hillslope to stream channel. Here we examine high-frequency sensor-based specific conductance and turbidity data from multiple storm events across two watersheds (Quebrada Sonadora and Rio Icacos) with different lithology in the Luquillo Mountains of Puerto Rico, a forested tropical ecosystem. Our analyses include Hurricane Maria, a category 5 hurricane. To analyze hysteresis, we used a recently developed set of metrics to describe and quantify storm events including the hysteresis index (HI), which describes the directionality of hysteresis loops, and the flushing index (FI), which can be used to infer whether the mobilization of material is source or transport limited. We also examine the role of antecedent discharge to predict hysteretic behavior during storms. Overall, specific conductance and turbidity showed contrasting responses to storms. The hysteretic behavior of specific conductance was very similar across sites, displaying clockwise hysteresis and a negative flushing index indicating proximal sources of solutes and consistent source limitation. In contrast, the directionality of turbidity hysteresis was significantly different between watersheds, although both had strong flushing behavior indicative of transport limitation. Overall, models that included antecedent discharge did not perform any better than models with peak discharge alone, suggesting that the magnitude and trajectory of an individual event was the strongest driver of material flux and hysteretic behavior. Hurricane Maria produced unique hysteresis metrics within both watersheds, indicating a distinctive response to this major hydrological event. The similarity in response of specific conductance to storms suggests that solute sources and pathways are similar in the two watersheds. The divergence in behavior for turbidity suggests that sources and pathways of particulate matter vary between the two watersheds. The use of high-frequency sensor data allows the quantification of storm events while index-based metrics of hysteresis allow for the direct comparison of complex storm events across a heterogeneous landscape and variable flow conditions.

Luquillo Mountains↗