SEARCH · USGS Science

Results for “Forests”

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 1,441 records · Page 80Linked to original sources

USGS2017

Mapping tree density in forests of the southwestern USA using Landsat 8 data

The increase of tree density in forests of the American Southwest promotes extreme fire events, understory biodiversity losses, and degraded habitat conditions for many wildlife species. To ameliorate these changes, managers and scientists have begun planning treatments aimed at reducing fuels and increasing understory biodiversity. However, spatial variability in tree density across the landscape is not well-characterized, and if better known, could greatly influence planning efforts. We used reflectance values from individual Landsat 8 bands (bands 2, 3, 4, 5, 6, and 7) and calculated vegetation indices (difference vegetation index, simple ratios, and normalized vegetation indices) to estimate tree density in an area planned for treatment in the Jemez Mountains, New Mexico, characterized by multiple vegetation types and a complex topography. Because different vegetation types have different spectral signatures, we derived models with multiple predictor variables for each vegetation type, rather than using a single model for the entire project area, and compared the model-derived values to values collected from on-the-ground transects. Among conifer-dominated areas (73% of the project area), the best models (as determined by corrected Akaike Information Criteria (AICc)) included Landsat bands 2, 3, 4, and 7 along with simple ratios, normalized vegetation indices, and the difference vegetation index ( R 2 values for ponderosa: 0.47, piñon-juniper: 0.52, and spruce-fir: 0.66). On the other hand, in aspen-dominated areas (9% of the project area), the best model included individual bands 4 and 2, simple ratio, and normalized vegetation index ( R 2 value: 0.97). Most areas dominated by ponderosa, pinyon-juniper, or spruce-fir had more than 100 trees per hectare. About 54% of the study area has medium to high density of trees (100–1000 trees/hectare), and a small fraction (4.5%) of the area has very high density (>1000 trees/hectare). Our results provide a better understanding of tree density for identifying areas in need of treatment and planning for more effective treatment. Our analysis also provides an integrated method of estimating tree density across complex landscapes that could be useful for further restoration planning.

New Mexico↗

USGS2017

Differences in human versus lightning fires between urban and rural areas of the boreal forest in interior Alaska

In western North America, the carbon-rich boreal forest is experiencing warmer temperatures, drier conditions and larger and more frequent wildfires. However, the fire regime is also affected by direct human activities through suppression, ignition, and land use changes. Models are important predictive tools for understanding future conditions but they are based on regional generalizations of wildfire behavior and weather that do not adequately account for the complexity of human–fire interactions. To achieve a better understanding of the intensity of human influence on fires in this sparsely populated area and to quantify differences between human and lightning fires, we analyzed fires by both ignition types in regard to human proximity in urban (the Fairbanks subregion) and rural areas of interior Alaska using spatial (Geographic Information Systems) and quantitative analysis methods. We found substantial differences in drivers of wildfire: while increases in fire ignitions and area burned were caused by lightning in rural interior Alaska, in the Fairbanks subregion these increases were due to human fires, especially in the wildland urban interface. Lightning fires are starting earlier and fires are burning longer, which is much more pronounced in the Fairbanks subregion than in rural areas. Human fires differed from lightning fires in several ways: they started closer to settlements and highways, burned for a shorter duration, were concentrated in the Fairbanks subregion, and often occurred outside the brief seasonal window for lightning fires. This study provides important insights that improve our understanding of the direct human influence on recently observed changes in wildfire regime with implications for both fire modeling and fire management.

Alaska↗

USGS1998

Mixed conifer forest mortality and establishment before and after prescribed fire in Sequoia National Park, California

Pre-and post-burn tree mortality rates, size structure, basal area, and ingrowth were determined for four 1.0 ha mixed conifer forest stands in the Log Creek and Tharp's Creek watersheds of Sequoia National Park. Mean annual mortality between 1986 and 1990 was 0.8% for both watersheds. In the fall of 1990, the Tharp's Creek watershed was treated with a prescribed burn. Between 1991 and 1995, mean annual mortality was 1.4% in the unburned Log Creek watershed and 17.2% in the burned Tharp's Creek watershed. A drought from 1987 to 1992 likely contributed to the mortality increase in the Log Creek watershed. The high mortality in the Tharp's Creek watershed was primarily related to crown scorch from the 1990 fire and was modeled with logistic regression for white fir (Abies concolor [Gord. and Glend.]) and sugar pine (Pinus lambertiana [Dougl.]). From 1989 to 1994, basal area declined an average of 5% per year in the burned Tharp's Creek watershed, compared to average annual increases of less than 1% per year in the unburned Log Creek watershed and in the Tharp's watershed prior to burning. Post-burn size structure was dramatically changed in the Tharp's Creek stands: 75% of trees ???50 cm and 25% of trees >50 cm were killed by the fire.

Forest Science↗

USGS1993

Epidemic pox and malaria in native forest birds

Studies by Warner in the 1950’s and van Riper in the 1970’s identified disease as a potential limiting factor in the distribution and abundance of Hawaii’s native forest birds. Mosquito-transmitted protozoan and viral infections caused by malarial parasites and pox virus were especially significant. Both organisms were introduced to the islands after the arrival of Europeans and are thought to have affected avian communities the same way that measles devastated native Hawaiian peoples.

Hawai'i↗

USGS1997

How wide is a road? The association of roads and mass-wasting in a forested montane environment

A spatial data base of 1609 landslides was analysed using a geographic information system to determine landslide frequency in relation to highways. A 126 km long transportation network in a 201km2 area of humid-tropical, mountainous, forested terrain in Puerto Rico was used in conjunction with a series of 20 buffer (disturbance) zones varying from 5 to 400m in length, measured perpendicular to the highways. Average landslide frequency in the study area at distances greater than 85m from roads was about six landslides per square kilometre. At distances of 85m or less on either side of a highway, landslide frequency was about 30 landslides per square kilometre. On average, this elevated disturbance rate affected 330m2km-2a-1 within the 170m swath. The mass-wasting rate outside of the disturbance zone affected 40m2km-2 a-1. These results indicate that the rate of mass-wasting disturbance is increased from five to eight times in a 170m wide swath along road corridors. The lateral extent of the environmental impact of roads in the study area is greater than is commonly perceived. The approach described herein demonstrates a simple method to assess the spatial association of mass-wasting with highways. ?? 1997 by John Wiley & Sons, Ltd.

Earth Surface Processes and Landforms↗

USGS2014

Parameter optimization, sensitivity, and uncertainty analysis of an ecosystem model at a forest flux tower site in the United States

Ecosystem models are useful tools for understanding ecological processes and for sustainable management of resources. In biogeochemical field, numerical models have been widely used for investigating carbon dynamics under global changes from site to regional and global scales. However, it is still challenging to optimize parameters and estimate parameterization uncertainty for complex process-based models such as the Erosion Deposition Carbon Model (EDCM), a modified version of CENTURY, that consider carbon, water, and nutrient cycles of ecosystems. This study was designed to conduct the parameter identifiability, optimization, sensitivity, and uncertainty analysis of EDCM using our developed EDCM-Auto, which incorporated a comprehensive R package—Flexible Modeling Framework (FME) and the Shuffled Complex Evolution (SCE) algorithm. Using a forest flux tower site as a case study, we implemented a comprehensive modeling analysis involving nine parameters and four target variables (carbon and water fluxes) with their corresponding measurements based on the eddy covariance technique. The local sensitivity analysis shows that the plant production-related parameters (e.g., PPDF1 and PRDX) are most sensitive to the model cost function. Both SCE and FME are comparable and performed well in deriving the optimal parameter set with satisfactory simulations of target variables. Global sensitivity and uncertainty analysis indicate that the parameter uncertainty and the resulting output uncertainty can be quantified, and that the magnitude of parameter-uncertainty effects depends on variables and seasons. This study also demonstrates that using the cutting-edge R functions such as FME can be feasible and attractive for conducting comprehensive parameter analysis for ecosystem modeling.

Massachusetts↗

USGS2014

Temporal variability of carbon and nutrient burial, sediment accretion, and mass accumulation over the past century in a carbonate platform mangrove forest of the Florida Everglades.

The objective of this research was to measure temporal variability in accretion and mass sedimentation rates (including organic carbon (OC), total nitrogen (TN), and total phosphorous (TP)) from the past century in a mangrove forest on the Shark River in Everglades National Park, USA. The 210 Pb Constant Rate of Supply model was applied to six soil cores to calculate annual rates over the most recent 10, 50, and 100 year time spans. Our results show that rates integrated over longer timeframes are lower than those for shorter, recent periods of observation. Additionally, the substantial spatial variability between cores over the 10 year period is diminished over the 100 year record, raising two important implications. First, a multiple-decade assessment of soil accretion and OC burial provides a more conservative estimate and is likely to be most relevant for forecasting these rates relative to long-term processes of sea level rise and climate change mitigation. Second, a small number of sampling locations are better able to account for spatial variability over the longer periods than for the shorter periods. The site average 100 year OC burial rate, 123 ± 19 (standard deviation) g m -2 yr -1 , is low compared with global mangrove values. High TN and TP burial rates in recent decades may lead to increased soil carbon remineralization, contributing to the low carbon burial rates. Finally, the strong correlation between OC burial and accretion across this site signals the substantial contribution of OC to soil building in addition to the ecosystem service of CO 2 sequestration.

Florida↗

USGS2017

Deposition of mercury in forests across a montane elevation gradient: Elevational and seasonal patterns in methylmercury inputs and production

Global mercury contamination largely results from direct primary atmospheric and secondary legacy emissions, which can be deposited to ecosystems, converted to methylmercury, and bioaccumulated along food chains. We examined organic horizon soil samples collected across an elevational gradient on Whiteface Mountain in the Adirondack region of New York State, USA to determine spatial patterns in methylmercury concentrations across a forested montane landscape. We found that soil methylmercury concentrations were highest in the midelevation coniferous zone (0.39 ± 0.07 ng/g) compared to the higher elevation alpine zone (0.28 ± 0.04 ng/g) and particularly the lower elevation deciduous zone (0.17 ± 0.02 ng/g), while the percent of total mercury as methylmercury in soils decreased with elevation. We also found a seasonal pattern in soil methylmercury concentrations, with peak methylmercury values occurring in July. Given elevational patterns in temperature and bioavailable total mercury (derived from mineralization of soil organic matter), soil methylmercury concentrations appear to be driven by soil processing of ionic Hg, as opposed to atmospheric deposition of methylmercury. These methylmercury results are consistent with spatial patterns of mercury concentrations in songbird species observed from other studies, suggesting that future declines in mercury emissions could be important for reducing exposure of mercury to montane avian species.

New York↗

USGS2017

Estimation and uncertainty of recent carbon accumulation and vertical accretion in drained and undrained forested peatlands of the southeastern USA

The purpose of this study was to determine how drainage impacts carbon densities and recent rates (past 50 years) of vertical accretion and carbon accumulation in southeastern forested peatlands. We compared these parameters in drained maple-gum (MAPL), Atlantic white cedar (CDR), and pocosin (POC) communities in the Great Dismal Swamp National Wildlife Refuge (GDS) of Virginia/North Carolina and in an intact (undrained) CDR swamp in the Alligator River National Wildlife Refuge (AR) of North Carolina. Peat cores were analyzed for bulk density, percent organic carbon, and 137 Cs and 210 Pb. An uncertainty analysis of both 137 Cs and 210 Pb approaches was used to constrain error at least partially related to mobility of both radioisotopes. GDS peats had lower porosities (89.6% (SD = 1.71) versus 95.3% (0.18)) and higher carbon densities (0.082 (0.021) versus 0.037 (0.009) g C cm −3 ) than AR. Vertical accretion rates (0.10–0.56 cm yr −1 ) were used to estimate a time period of ~84–362 years for reestablishment of peat lost during the 2011 Lateral West fire at the GDS. Carbon accumulation rates ranged from 51 to 389 g C m −2 yr −1 for all sites. In the drained (GDS) versus intact (AR) CDR sites, carbon accumulation rates were similar with 137 Cs (87 GDS versus 92 AR g C m −2 yr −1 ) and somewhat less at the GDS than AR as determined with 210 Pb (111 GDS versus 159 AR g C m −2 yr −1 ). Heightened productivity and high polyphenol content of peat may be responsible for similar rates of carbon accumulation in both drained and intact CDR peatlands.

Alligator Rivers National Wildlife Refuge, Great D↗

USGS2018

Quantifying climate-related interactions in shallow and deep storage and evapotranspiration in a forested, seasonally water-limited watershed in the Southeastern United States

The Southeastern United States experiences recurring hydrological droughts, which can reduce water availability and can result in water-limiting conditions. Long-term monitoring at Panola Mountain Research Watershed, a small, forested, seasonally water-limited watershed near Atlanta, Georgia, was used to quantify the interactions of climatic variability with shallow and deep storage and evapotranspiration. Watershed storage (WS) and actual evapotranspiration (AET) were estimated monthly from 1985 through 2015 using a water-budget approach combined with a WS-baseflow relationship. Shallow storage (SS) was assessed from a soil moisture profile. Soil moisture transitioned from recharge to surplus as SS increased from its field capacity to a nearly saturated state during the dormant season, and transitioned from utilization to climatic water deficits as SS declined from its field capacity to its wilting point during the growing season. Deeper storage was unavailable to AET during dry conditions. The majority of deeper storage recharge occurred during the dormant season and required SS to be wet. WS was an effective drought indicator. Growing season droughts typically occurred when WS was below normal at the end of the dormant season and growing season precipitation (P) was below or near normal. A hydrologic persistence analysis found that monthly-standardized WS was significantly correlated (p-value <0.05) with past monthly-standardized WS for the previous 19 months and with past monthly P for the previous 11 months, indicating the importance of past hydrologic conditions on WS. Expected climatic changes affected recharge during the dormant season and deficits during the growing season.

Water Resources Research↗

USGS2017

Similarity of plant functional traits and aggregation pattern in a subtropical forest

The distribution of species and communities in relation to environmental heterogeneity is a central focus in ecology. Co-occurrence of species with similar functional traits is an indication that communities are determined in part by environmental filters. However, few studies have been designed to test how functional traits are selectively filtered by environmental conditions at local scales. Exploring the relationship between soil characteristics and plant traits is a step toward understanding the filtering hypothesis in determining plant distribution at local scale. Toward this end, we mapped all individual trees (diameter >1 cm) in a one-ha subtropical forest of China in 2007 and 2015. We measured topographic and detailed soil properties within the field site, as well as plant leaf functional traits and demographic rates of the seven most common tree species. A second one-ha study plot was established in 2015, to test and validate the general patterns that were drawn from first plot. We found that variation in species distribution at local scale can be explained by soil heterogeneity and plant functional traits. (From first plot). (1) Species dominant in habitats with high soil ammonium nitrogen and total phosphorus tended to have high specific leaf area (SLA) and relative growth rate (RGR). (2) Species dominant in low-fertility habitats tended to have high leaf dry matter content (LDMC), ratio of chlorophyll a and b (ratioab), and leaf thickness (LT). The hypothesis that functional traits are selected in part by environmental filters and determine plant distribution at local scale was confirmed by the data of the first plot and a second regional site showed similar species distribution patterns.

Zhejiang Province↗

USGS2019

Impact of prey occupancy and other ecological and anthropogenic factors on Tiger distribution in Thailand’s Western Forest Complex

Despite conservation efforts, large mammals such as tigers (Panthera tigris) and their main prey, gaur (Bos gaurus), banteng (Bos javanicus), and sambar (Rusa unicolor), are highly threatened and declining across their entire range. The only large viable source population of tigers in mainland Southeast Asia occurs in Thailand's Western Forest Complex (WEFCOM), an approximately 19,000 km 2 landscape of 17 contiguous protected areas. We used an occupancy modeling framework, which accounts for imperfect detection, to identify the factors that affect tiger distribution at the approximate scale of a female tiger's home range, 64 km 2 , and site use at a scale of 1-km 2 . At the larger scale, we estimated the proportion of sites at WEFCOM that were occupied by tigers; at the finer scale, we identified the key variables that influence site-use and developed a predictive distribution map. At both scales, we examined key anthropogenic and ecological factors that help explain tiger distribution and habitat use, including probabilities of gaur, banteng, and sambar occurrence from a companion study. Occupancy estimated at the 64-km 2 scale was primarily influenced by the combined presence of all three large prey species, and 37% or 5,858 km 2 of the landscape was predicted to be occupied by tigers. In contrast, site use estimated at the scale of 1 km 2 was most strongly influenced by the presence of sambar. By modeling occupancy while accounting for imperfect probability of detection, we established reliable benchmark data on the distribution of tigers in WEFCOM. This study also identified factors that limit tiger distributions; which managers can then target to expand tiger distribution and guide recovery elsewhere in Southeast Asia.

Ecology and Evolution↗

USGS2015

Proximity to encroaching coconut palm limits native forest water use and persistence on a Pacific atoll

Competition for fresh water between native and introduced plants is one important challenge facing native forests as rainfall variability increases. Competition can be especially acute for vegetation on Pacific atolls, which depend upon consistent rainfall to replenish shallow groundwater stores. Patterns of sap flow, water use, and diameter growth of Pisonia grandis trees were investigated on Sand Islet, Palmyra Atoll, Line Islands, during a period of low rainfall. Sap flow in the outer sapwood was reduced by 53% for P . grandis trees growing within coconut palm ( Cocos nucifera ) stands ( n  = 9) versus away from coconut palm ( n  = 9). This suggested that water uptake was being limited by coconut palm. Radial patterns of sap flow into the sapwood of P . grandis also differed between stands with and without coconut palm, such that individual tree water use for P . grandis ranged from 14 to 67 L day −1 , averaging 47·8 L day −1 without coconut palm and 23·6 L day −1 with coconut palm. Diameter growth of P . grandis was measured from nine islets. In contrast to sap flow, competition with coconut palm increased diameter growth by 89%, equating to an individual tree basal area increment of 5·4 versus 10·3 mm 2  day −1 . Greater diameter growth countered by lower rates of water use by P . grandis trees growing in competition with coconut palm suggests that stem swell may be associated with water storage when positioned in the understory of coconut palm, and may facilitate survival when water becomes limiting until too much shading overwhelms P . grandis .

Palmyra Atoll↗

USGS2017

Effects of flood inundation and invasion by Phalaris arundinacea on nitrogen cycling in an Upper Mississippi River floodplain forest

Although floodplains are thought to serve as important buffers against nitrogen (N) transport to aquatic systems, frequent flooding and high levels of nutrient availability also make these systems prone to invasion by exotic plant species. Invasive plants could modify the cycling and availability of nutrients within floodplains, with effects that could feedback to promote the persistence of the invasive species and impact N export to riverine and coastal areas. We examined the effect of flooding on soil properties and N cycling at a floodplain site in Pool 8 of the Upper Mississippi River with 2 plant communities: mature native forest ( Acer saccharinum ) and patches of an invasive grass ( Phalaris arundinacea ). Plots were established within each vegetation type along an elevation gradient and sampled throughout the summers of 2013 and 2014. Spatial trends in flooding resulted in higher soil organic matter, porosity, and total nitrogen and carbon in low elevations. Nutrient processes and NH 4 + and NO 3 − availability, however, were best explained by vegetation type and time after flooding. Phalaris plots maintained higher rates of nitrification and higher concentrations of available NH 4 + and NO 3 − . These results suggest that invasion by Phalaris may make nitrogen more readily available and could help to reinforce this species' persistence in floodplain wetlands. They also raise the possibility that Phalaris may decrease floodplain N storage capacity and influence downstream transport of N to coastal zones.

Wisconsin↗

USGS2016

Static and dynamic controls on fire activity at moderate spatial and temporal scales in the Alaskan boreal forest

Wildfire, a dominant disturbance in boreal forests, is highly variable in occurrence and behavior at multiple spatiotemporal scales. New data sets provide more detailed spatial and temporal observations of active fires and the post-burn environment in Alaska. In this study, we employ some of these new data to analyze variations in fire activity by developing three explanatory models to examine the occurrence of (1) seasonal periods of elevated fire activity using the number of MODIS active fire detections data set (MCD14DL) within an 11-day moving window, (2) unburned patches within a burned area using the Monitoring Trends in Burn Severity fire severity product, and (3) short-to-moderate interval (<60 yr) fires using areas of burned area overlap in the Alaska Large Fire Database. Explanatory variables for these three models included dynamic variables that can change over the course of the fire season, such as weather and burn date, as well as static variables that remain constant over a fire season, such as topography, drainage, vegetation cover, and fire history. We found that seasonal periods of high fire activity are associated with both seasonal timing and aggregated weather conditions, as well as the landscape composition of areas that are burning. Important static inputs to the model of seasonal fire activity indicate that when fire weather conditions are suitable, areas that typically resist fire (e.g., deciduous stands) may become more vulnerable to burning and therefore less effective as fire breaks. The occurrence of short-to-moderate interval fires appears to be primarily driven by weather conditions, as these were the only relevant explanatory variables in the model. The unique importance of weather in explaining short-to-moderate interval fires implies that fire return intervals (FRIs) will be sensitive to projected climate changes in the region. Unburned patches occur most often in younger stands, which may be related to a greater deciduous fraction of vegetation as well as lower fuel loads compared with mature stands. The fraction of unburned patches may therefore increase in response to decreasing FRIs and increased deciduousness in the region, or these may decrease if fire weather conditions become more severe.

Alaska↗

USGS2016

Topographic and fire weather controls of fire refugia in forested ecosystems of northwestern North America

Fire refugia, sometimes referred to as fire islands, shadows, skips, residuals, or fire remnants, are an important element of the burn mosaic, but we lack a quantitative framework that links observations of fire refugia from different environmental contexts. Here, we develop and test a conceptual model for how predictability of fire refugia varies according to topographic complexity and fire weather conditions. Refugia were quantified as areas unburned or burned at comparatively low severity based on remotely sensed burn severity data. We assessed the relationship between refugia and a suite of terrain-related explanatory metrics by fitting a collection of boosted regression tree models. The models were developed for seven study fires that burned in conifer-dominated forested landscapes of the Western Cordillera of Canada between 2001 and 2014. We fit nine models, each for distinct levels of fire weather and terrain ruggedness. Our framework revealed that the predictability and abundance of fire refugia varied among these environmental settings. We observed highest predictability under moderate fire weather conditions and moderate terrain ruggedness (ROC-AUC = 0.77), and lowest predictability in flatter landscapes and under high fire weather conditions (ROC-AUC = 0.63–0.68). Catchment slope, local aspect, relative position, topographic wetness, topographic convergence, and local slope all contributed to discriminating where refugia occur but the relative importance of these topographic controls differed among environments. Our framework allows us to characterize the predictability of contemporary fire refugia across multiple environmental settings and provides important insights for ecosystem resilience, wildfire management, conservation planning, and climate change adaptation.

Ecosphere↗

USGS2023

Patterns and controls of foliar nutrient stoichiometry and flexibility across United States forests

Plant element stoichiometry and stoichiometric flexibility strongly regulate ecosystem responses to global change. Here, we tested three potential mechanistic drivers (climate, soil nutrients, and plant taxonomy) of both using paired foliar and soil nutrient data from terrestrial forested National Ecological Observatory Network sites across the USA. We found that broad patterns of foliar nitrogen (N) and foliar phosphorus (P) are explained by different mechanisms. Plant taxonomy was an important control over all foliar nutrient stoichiometries and concentrations, especially foliar N, which was dominantly related to taxonomy and did not vary across climate or soil gradients. Despite a lack of site-level correlations between N and environment variables, foliar N exhibited intraspecific flexibility, with numerous species-specific correlations between foliar N and various environmental factors, demonstrating the variable spatial and temporal scales on which foliar chemistry and stoichiometric flexibility can manifest. In addition to plant taxonomy, foliar P and N:P ratios were also linked to soil nutrient status (extractable P) and climate, especially actual evapotranspiration rates. Our findings highlight the myriad factors that influence foliar chemistry and show that broad patterns cannot be explained by a single consistent mechanism. Furthermore, differing controls over foliar N versus P suggests that each may be sensitive to global change drivers on distinct spatial and temporal scales, potentially resulting in altered ecosystem N:P ratios that have implications for processes ranging from productivity to carbon sequestration.

Ecology↗

USGS2019

Effects of climate-related variability in storage on streamwater solute concentrations and fluxes in a small forested watershed in the Southeastern United States

Streamwater quality can be affected by climate-related variability in hydrologic state, which controls flow paths and affects biogeochemical processes. Thirty-one years of input/output solute fluxes at Panola Mountain Research Watershed, a small, forested, seasonally water-limited watershed near Atlanta, Georgia, were used to quantify the effects of climatic-related variability in storage on streamwater solute concentrations and fluxes. Streamwater fluxes were estimated for ten solutes from weekly and event sample concentrations using regression-based methods. The most pertinent storage attribute (current or antecedent watershed, shallow, and deep storage) for each solute was determined by fitting separate concentration relationships. The concentration-discharge relationships varied more for reactive solutes such as potassium, sulfate, and DOC and less for weathering products (base cations and dissolved silica) and conservative chloride. Many solutes exhibited higher concentrations when storage levels were lower or wetting up, which was likely the result of the concentrating effects of evapotranspiration and/or the buildup and flushing of weathering products associated with longer residence times. The impacts of storage modeling on annual fluxes varied by solute, ranging from about 5% (magnesium) to 52% (nitrate) as relative standard deviations, and sufficiently removed climate-related patterns observed in streamwater concentrations. Sulfate was particularly mobilized following growing season droughts but only if deep storage was sufficiently recharged, possibly indicating that sulfides in the deep storage pool were oxidized to sulfate during droughts and mobilized when re-wetted. The lack of streamwater sulfate response to 61% declines in atmospheric deposition indicates the importance of watershed biogeochemical processes on controls of streamwater export of sulfate. The approach of explicitly incorporating storage in the streamwater concentration modeling elucidated the effects of climate on streamwater water-quality and may provide insight into the effects of climatic change on future fluxes.

Georgia↗