Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research: Biogeosciences”

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 19 recordsLinked to original sources
USGS2013

Seasonal changes in peatland surface elevation recorded at GPS stations in the Red Lake Peatlands, northern Minnesota, USA

Northern peatlands appear to hold large volumes of free‐phase gas (e.g., CH 4 and CO 2 ), which has been detected by surface deformations, pore pressure profiles, and electromagnetic surveys. Determining the gas content and its impact in peat is challenging because gas storage depends on both the elastic properties of the peat matrix and the buoyant forces exerted by pore fluids. We therefore used a viscoelastic deformation model to estimate these variables by adjusting model runs to reproduce observed changes in peat surface elevation within a 1300 km 2 peatland. A local GPS network documented significant changes in surface elevations throughout the year with the greatest vertical displacements associated with rapid changes in peat water content and unloadings due to melting of the winter snowpack. These changes were coherent with changes in water table elevation and also abnormal pore pressure changes measured by nests of instrumented piezometers. The deformation model reproduced these changes when the gas content was adjusted to 10% of peat volume, and Young's modulus was varied between 5 and 100 kPa as the peat profile shifted from tension to compression. In contrast, the model predicted little peat deformation when the gas content was 3% or lower. These model simulations are consistent with previous estimates of gas volume in northern peatlands and suggest an upper limit of gas storage controlled by the elastic moduli of the peat fabric.

Minnesota↗
USGS2014

Controls on methane released through ebullition in peatlands affected by permafrost degradation

Permafrost thaw in peat plateaus leads to the flooding of surface soils and the formation of collapse scar bogs, which have the potential to be large emitters of methane (CH 4 ) from surface peat as well as deeper, previously frozen, permafrost carbon (C). We used a network of bubble traps, permanently installed 20 cm and 60 cm beneath the moss surface, to examine controls on ebullition from three collapse bogs in interior Alaska. Overall, ebullition was dominated by episodic events that were associated with changes in atmospheric pressure, and ebullition was mainly a surface process regulated by both seasonal ice dynamics and plant phenology. The majority (>90%) of ebullition occurred in surface peat layers, with little bubble production in deeper peat. During periods of peak plant biomass, bubbles contained acetate-derived CH 4 dominated (>90%) by modern C fixed from the atmosphere following permafrost thaw. Post-senescence, the contribution of CH 4 derived from thawing permafrost C was more variable and accounted for up to 22% (on average 7%), in the most recently thawed site. Thus, the formation of thermokarst features resulting from permafrost thaw in peatlands stimulates ebullition and CH 4 release both by creating flooded surface conditions conducive to CH 4 production and bubbling as well as by exposing thawing permafrost C to mineralization.

Alaska↗
USGS2014

Floodplain biogeochemical processing of floodwaters in the Atchafalaya River Basin during the Mississippi River flood of 2011

The 2011 flood in the Lower Mississippi resulted in the second highest recorded river flow diverted into the Atchafalaya River Basin (ARB). The higher water levels during the flood peak resulted in high hydrologic connectivity between the Atchafalaya River and floodplain, with up to 50% of the Atchafalaya River water moving off channel. Water quality samples were collected throughout the ARB over the course of the flood event. Significant nitrate (NO 3 -) reduction (75%) occurred within the floodplain, resulting in a total NO 3 - reduction of 16.6% over the flood. The floodplain was a small but measurable source of dissolved reactive phosphorus (SRP) and ammonium (NH 4 +). Collectively, these results from this large flood event suggest that enhancing river-floodplain connectivity through freshwater diversions will reduce NO 3 - loads to the Gulf of Mexico during large annual floods.

Louisiana↗
USGS2014

The impact of climate and reservoirs on longitudinal riverine carbon fluxes from two major watersheds in the Central and Intermontane West

A nested sampling network on the Colorado (CR) and Missouri Rivers (MR) provided data to assess impacts of large-scale reservoir systems and climate on carbon export. The Load Estimator (LOADEST) model was used to estimate both dissolved inorganic and organic carbon (DIC and DOC) fluxes for a total of 22 sites along the main stems of the CR and MR. Both the upper CR and MR DIC and DOC fluxes increased longitudinally, but the lower CR fluxes decreased while the lower MRs continued to increase. We examined multiple factors through space and time that help explain these flux patterns. Seasonal variability in precipitation and temperature, along with site-level concentration versus discharge relationships proved to be significant factors explaining much of the difference among sites located below reservoirs as compared to sites located in more free-flowing segments of the river. The characterization of variability in carbon exports over space and time provides a basis for understanding carbon cycling and transport within river basins affected by large reservoir systems, particular in arid-to semi-arid ecosystems.

Journal of Geophysical Research: Biogeosciences↗
USGS2014

Spatial variability and landscape controls of near-surface permafrost within the Alaskan Yukon River Basin

The distribution of permafrost is important to understand because of permafrost's influence on high-latitude ecosystem structure and functions. Moreover, near-surface (defined here as within 1 m of the Earth's surface) permafrost is particularly susceptible to a warming climate and is generally poorly mapped at regional scales. Subsequently, our objectives were to (1) develop the first-known binary and probabilistic maps of near-surface permafrost distributions at a 30 m resolution in the Alaskan Yukon River Basin by employing decision tree models, field measurements, and remotely sensed and mapped biophysical data; (2) evaluate the relative contribution of 39 biophysical variables used in the models; and (3) assess the landscape-scale factors controlling spatial variations in permafrost extent. Areas estimated to be present and absent of near-surface permafrost occupy approximately 46% and 45% of the Alaskan Yukon River Basin, respectively; masked areas (e.g., water and developed) account for the remaining 9% of the landscape. Strong predictors of near-surface permafrost include climatic indices, land cover, topography, and Landsat 7 Enhanced Thematic Mapper Plus spectral information. Our quantitative modeling approach enabled us to generate regional near-surface permafrost maps and provide essential information for resource managers and modelers to better understand near-surface permafrost distribution and how it relates to environmental factors and conditions.

Alaska↗
USGS2014

Differential response of carbon fluxes to climate in three peatland ecosystems that vary in the presence and stability of permafrost

Changes in vegetation and soil properties following permafrost degradation and thermokarst development in peatlands may cause changes in net carbon storage. To better understand these dynamics, we established three sites in Alaska that vary in permafrost regime, including a black spruce peat plateau forest with stable permafrost, an internal collapse scar bog formed as a result of thermokarst, and a rich fen without permafrost. Measurements include year-round eddy covariance estimates of carbon dioxide (CO 2 ), water, and energy fluxes, associated environmental variables, and methane (CH 4 ) fluxes at the collapse scar bog. The ecosystems all acted as net sinks of CO 2 in 2011 and 2012, when air temperature and precipitation remained near long-term means. In 2013, under a late snowmelt and late leaf out followed by a hot, dry summer, the permafrost forest and collapse scar bog were sources of CO 2 . In this same year, CO 2 uptake in the fen increased, largely because summer inundation from groundwater inputs suppressed ecosystem respiration. CO 2 exchange in the permafrost forest and collapse scar bog was sensitive to warm air temperatures, with 0.5 g C m −2 lost each day when maximum air temperature was very warm (≥29°C). The bog lost 4981 ± 300 mg CH 4 m −2 between April and September 2013, indicating that this ecosystem acted as a significant source of both CO 2 and CH 4 to the atmosphere in 2013. These results suggest that boreal peatland responses to warming and drying, both of which are expected to occur in a changing climate, will depend on permafrost regime.

Alaska↗
USGS2014

Radiocarbon age-offsets in an arctic lake reveal the long-term response of permafrost carbon to climate change

Continued warming of the Arctic may cause permafrost to thaw and speed the decomposition of large stores of soil organic carbon (OC), thereby accentuating global warming. However, it is unclear if recent warming has raised the current rates of permafrost OC release to anomalous levels or to what extent soil carbon release is sensitive to climate forcing. Here we use a time series of radiocarbon age-offsets ( 14 C) between the bulk lake sediment and plant macrofossils deposited in an arctic lake as an archive for soil and permafrost OC release over the last 14,500 years. The lake traps and archives OC imported from the watershed and allows us to test whether prior warming events stimulated old carbon release and heightened age-offsets. Today, the age-offset (2 ka; thousand of calibrated years before A.D. 1950) and the depositional rate of ancient OC from the watershed into the lake are relatively low and similar to those during the Younger Dryas cold interval (occurring 12.9–11.7 ka). In contrast, age-offsets were higher (3.0–5.0 ka) when summer air temperatures were warmer than present during the Holocene Thermal Maximum (11.7–9.0 ka) and Bølling-Allerød periods (14.5–12.9 ka). During these warm times, permafrost thaw contributed to ancient OC depositional rates that were ~10 times greater than today. Although permafrost OC was vulnerable to climate warming in the past, we suggest surface soil organic horizons and peat are presently limiting summer thaw and carbon release. As a result, the temperature threshold to trigger widespread permafrost OC release is higher than during previous warming events.

Alaska↗
USGS2014

Using dissolved organic matter age and composition to detect permafrost thaw in boreal watersheds of interior Alaska

Recent warming at high latitudes has accelerated permafrost thaw, which can modify soil carbon dynamics and watershed hydrology. The flux and composition of dissolved organic matter (DOM) from soils to rivers are sensitive to permafrost configuration and its impact on subsurface hydrology and groundwater discharge. Here, we evaluate the utility of DOM composition and age as a tool for detecting permafrost thaw in three rivers (Beaver, Birch, and Hess Creeks) within the discontinuous permafrost zone of interior Alaska. We observed strong temporal controls on Δ 14 C content of hydrophobic acid isolates (Δ 14 C-HPOA) across all rivers, with the most enriched values occurring during spring snowmelt (75 ± 8‰) and most depleted during winter flow (−21 ± 8‰). Radiocarbon ages of winter flow samples ranged from 35 to 445 yr BP, closely tracking estimated median base flow travel times for this region (335 years). During spring snowmelt, young DOM was composed of highly aromatic, high molecular-weight compounds, whereas older DOM of winter flow had lower aromaticity and molecular weight. We observed a significant correlation between Δ 14 C-HPOA and UV absorbance coefficient at 254 nm ( α 254 ) across all study rivers. Using α 254 as an optical indicator for Δ 14 C-HPOA, we also observed a long-term decline in α 254 during maximum annual thaw depth over the last decade at the Hess Creek study site. These findings suggest a shift in watershed hydrology associated with increasing active layer thickness. Further development of DOM optical indicators may serve as a novel and inexpensive tool for detecting permafrost degradation in northern watersheds.

Alaska↗
USGS2014

Nitrogen cycling processes and microbial community composition in bed sediments in the Yukon River at Pilot Station

Information on the contribution of nitrogen (N)-cycling processes in bed sediments to river nutrient fluxes in large northern latitude river systems is limited. This study examined the relationship between N-cycling processes in bed sediments and N speciation and loading in the Yukon River near its mouth at the Bering Sea. We conducted laboratory bioassays to measure N-cycling processes in sediment samples collected over distinct water cycle seasons. In conjunction, the microbial community composition in the bed sediments using genes involved in N-cycling ( narG , napA , nosZ , and amoA ) and 16S rRNA gene pyrosequences was examined. Temporal variation was observed in net N mineralization, nitrate uptake, and denitrification rate potentials and correlated strongly with sediment carbon (C) and extractable N content and microbial community composition rather than with river water nutrient concentrations. The C content of the bed sediment was notably impacted by the spring flood, ranging from 1.1% in the midst of an ice-jam to 0.1% immediately after ice-out, suggesting a buildup of organic material (OM) prior to scouring of the bed sediments during ice break up. The dominant members of the microbial community that explained differences in N-processing rates belonged to the genera Crenothrix , Flavobacterium , and the family of Comamonadaceae. Our results suggest that biogeochemical processing rates in the bed sediments appear to be more coupled to hydrology, nutrient availability in the sediments, and microbial community composition rather than river nutrient concentrations at Pilot Station.

Alaska↗
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↗
USGS2015

Complex terrain alters temperature and moisture limitations of forest soil respiration across a semiarid to subalpine gradient

Forest soil respiration is a major carbon (C) flux that is characterized by significant variability in space and time. We quantified growing season soil respiration during both a drought year and a nondrought year across a complex landscape to identify how landscape and climate interact to control soil respiration. We asked the following questions: (1) How does soil respiration vary across the catchments due to terrain-induced variability in moisture availability and temperature? (2) Does the relative importance of moisture versus temperature limitation of respiration vary across space and time? And (3) what terrain elements are important for dictating the pattern of soil respiration and its controls? Moisture superseded temperature in explaining watershed respiration patterns, with wetter yet cooler areas higher up and on north facing slopes yielding greater soil respiration than lower and south facing areas. Wetter subalpine forests had reduced moisture limitation in favor of greater seasonal temperature limitation, and the reverse was true for low-elevation semiarid forests. Coincident climate poorly predicted soil respiration in the montane transition zone; however, antecedent precipitation from the prior 10 days provided additional explanatory power. A seasonal trend in respiration remained after accounting for microclimate effects, suggesting that local climate alone may not adequately predict seasonal variability in soil respiration in montane forests. Soil respiration climate controls were more strongly related to topography during the drought year highlighting the importance of landscape complexity in ecosystem response to drought.

Journal of Geophysical Research: Biogeosciences↗
USGS2015

Soil greenhouse gas emissions and carbon budgeting in a short-hydroperiod floodplain wetland

Understanding the controls on floodplain carbon (C) cycling is important for assessing greenhouse gas emissions and the potential for C sequestration in river-floodplain ecosystems. We hypothesized that greater hydrologic connectivity would increase C inputs to floodplains that would not only stimulate soil C gas emissions but also sequester more C in soils. In an urban Piedmont river (151 km 2 watershed) with a floodplain that is dry most of the year, we quantified soil CO 2 , CH 4 , and N 2 O net emissions along gradients of floodplain hydrologic connectivity, identified controls on soil aerobic and anaerobic respiration, and developed a floodplain soil C budget. Sites were chosen along a longitudinal river gradient and across lateral floodplain geomorphic units (levee, backswamp, and toe slope). CO 2 emissions decreased downstream in backswamps and toe slopes and were high on the levees. CH 4 and N 2 O fluxes were near zero; however, CH 4 emissions were highest in the backswamp. Annual CO 2 emissions correlated negatively with soil water-filled pore space and positively with variables related to drier, coarser soil. Conversely, annual CH 4 emissions had the opposite pattern of CO 2 . Spatial variation in aerobic and anaerobic respiration was thus controlled by oxygen availability but was not related to C inputs from sedimentation or vegetation. The annual mean soil CO 2 emission rate was 1091 g C m −2  yr −1 , the net sedimentation rate was 111 g C m −2  yr −1 , and the vegetation production rate was 240 g C m −2  yr −1 , with a soil C balance (loss) of −338 g C m −2  yr −1 . This floodplain is losing C likely due to long-term drying from watershed urbanization.

Virginia↗
USGS2015

Soil surface organic layers in Arctic Alaska: spatial distribution, rates of formation, and microclimatic effects

Organic layers of living and dead vegetation cover the ground surface in many permafrost landscapes and play important roles in ecosystem processes. These soil surface organic layers (SSOLs) store large amounts of carbon and buffer the underlying permafrost and its contained carbon from changes in aboveground climate. Understanding the dynamics of SSOLs is a prerequisite for predicting how permafrost and carbon stocks will respond to warming climate. Here we ask three questions about SSOLs in a representative area of the Arctic Foothills region of northern Alaska: (1) What environmental factors control the thickness of SSOLs and the carbon they store? (2) How long do SSOLs take to develop on newly stabilized point bars? (3) How do SSOLs affect temperature in the underlying ground? Results show that SSOL thickness and distribution correlate with elevation, drainage area, vegetation productivity, and incoming solar radiation. A multiple regression model based on these correlations can simulate spatial distribution of SSOLs and estimate the organic carbon stored there. SSOLs develop within a few decades after a new, sandy, geomorphic surface stabilizes but require 500–700 years to reach steady state thickness. Mature SSOLs lower the growing season temperature and mean annual temperature of the underlying mineral soil by 8 and 3°C, respectively. We suggest that the proximate effects of warming climate on permafrost landscapes now covered by SSOLs will occur indirectly via climate's effects on the frequency, extent, and severity of disturbances like fires and landslides that disrupt the SSOLs and interfere with their protection of the underlying permafrost.

Alaska↗
USGS2016

Effects of permafrost aggradation on peat properties as determined from a pan-Arctic synthesis of plant macrofossils

Permafrost dynamics play an important role in high-latitude peatland carbon balance and are key to understanding the future response of soil carbon stocks. Permafrost aggradation can control the magnitude of the carbon feedback in peatlands through effects on peat properties. We compiled peatland plant macrofossil records for the northern permafrost zone (515 cores from 280 sites) and classified samples by vegetation type and environmental class (fen, bog, tundra and boreal permafrost, and thawed permafrost). We examined differences in peat properties (bulk density, carbon (C), nitrogen (N) and organic matter content, and C/N ratio) and C accumulation rates among vegetation types and environmental classes. Consequences of permafrost aggradation differed between boreal and tundra biomes, including differences in vegetation composition, C/N ratios, and N content. The vegetation composition of tundra permafrost peatlands was similar to permafrost-free fens, while boreal permafrost peatlands more closely resembled permafrost-free bogs. Nitrogen content in boreal permafrost and thawed permafrost peatlands was significantly lower than in permafrost-free bogs despite similar vegetation types (0.9% versus 1.5% N). Median long-term C accumulation rates were higher in fens (23 g C m −2   yr −1 ) than in permafrost-free bogs (18 g C m −2  yr −1 ) and were lowest in boreal permafrost peatlands (14 g C m −2  yr −1 ). The plant macrofossil record demonstrated transitions from fens to bogs to permafrost peatlands, bogs to fens, permafrost aggradation within fens, and permafrost thaw and reaggradation. Using data synthesis, we have identified predominant peatland successional pathways, changes in vegetation type, peat properties, and C accumulation rates associated with permafrost aggradation.

Journal of Geophysical Research: Biogeosciences↗
USGS2016

Variation of energy and carbon fluxes from a restored temperate freshwater wetland and implications for carbon market verification protocols

Temperate freshwater wetlands are among the most productive terrestrial ecosystems, stimulating interest in using restored wetlands as biological carbon sequestration projects for greenhouse gas reduction programs. In this study, we used the eddy covariance technique to measure surface energy carbon fluxes from a constructed, impounded freshwater wetland during two annual periods that were 8 years apart: 2002–2003 and 2010–2011. During 2010–2011, we measured methane (CH 4 ) fluxes to quantify the annual atmospheric carbon mass balance and its concomitant influence on global warming potential (GWP). Peak growing season fluxes of latent heat and carbon dioxide (CO 2 ) were greater in 2002–2003 compared to 2010–2011. In 2002, the daily net ecosystem exchange reached as low as −10.6 g C m −2  d −1 , which was greater than 3 times the magnitude observed in 2010 (−2.9 g C m −2  d −1 ). CH 4 fluxes during 2010–2011 were positive throughout the year and followed a strong seasonal pattern, ranging from 38.1 mg C m −2  d −1 in the winter to 375.9 mg C m −2  d −1 during the summer. The results of this study suggest that the wetland had reduced gross ecosystem productivity in 2010–2011, likely due to the increase in dead plant biomass (standing litter) that inhibited the generation of new vegetation growth. In 2010–2011, there was a net positive GWP (675.3 g C m −2  yr −1 ), and when these values are evaluated as a sustained flux, the wetland will not reach radiative balance even after 500 years.

California↗
USGS2015

Abundance, stable isotopic composition, and export fluxes of DOC, POC, and DIC from the Lower Mississippi River during 2006–2008

Sources, abundance, isotopic compositions, and export fluxes of dissolved inorganic carbon (DIC), dissolved and colloidal organic carbon (DOC and COC), and particulate organic carbon (POC), and their response to hydrologic regimes were examined through monthly sampling from the Lower Mississippi River during 2006–2008. DIC was the most abundant carbon species, followed by POC and DOC. Concentration and δ 13 C of DIC decreased with increasing river discharge, while those of DOC remained fairly stable. COC comprised 61 ± 3% of the bulk DOC with similar δ 13 C abundances but higher percentages of hydrophobic organic acids than DOC, suggesting its aromatic and diagenetically younger status. POC showed peak concentrations during medium flooding events and at the rising limb of large flooding events. While δ 13 C-POC increased, δ 15 N of particulate nitrogen decreased with increasing discharge. Overall, the differences in δ 13 C between DOC or DIC and POC show an inverse correlation with river discharge. The higher input of soil organic matter and respired CO 2 during wet seasons was likely the main driver for the convergence of δ 13 C between DIC and DOC or POC, whereas enhanced in situ primary production and respiration during dry seasons might be responsible for their isotopic divergence. Carbon export fluxes from the Mississippi River were estimated to be 13.6 Tg C yr −1 for DIC, 1.88 Tg C yr −1 for DOC, and 2.30 Tg C yr −1 for POC during 2006–2008. The discharge-normalized DIC yield decreased during wet seasons, while those of POC and DOC increased and remained constant, respectively, implying variable responses in carbon export to the increasing discharge.

Louisiana↗
USGS2016

Component greenhouse gas fluxes and radiative balance from two deltaic marshes in Louisiana: Pairing chamber techniques and eddy covariance

Coastal marshes take up atmospheric CO 2 while emitting CO 2 , CH 4 , and N 2 O. This ability to sequester carbon (C) is much greater for wetlands on a per-area basis than from most ecosystems, facilitating scientific, political, and economic interest in their value as greenhouse gas sinks. However, the greenhouse gas balance of Gulf of Mexico wetlands is particularly understudied. We describe the net ecosystem exchange (NEE c ) of CO 2 and CH 4 using eddy covariance (EC) in comparison with fluxes of CO 2 , CH 4 , and N 2 O using chambers from brackish and freshwater marshes in Louisiana, USA. From EC, we found that 182 g C m -2 y -1 was lost through NEE c from the brackish marsh. Of this, 11 g C m -2 y -1 resulted from net CH 4 emissions and the remaining 171 g C m -2 y -1 resulted from net CO 2 emissions. In contrast, -290 g C m 2 y -1 was taken up through NEE c by the freshwater marsh, with 47 g C m -2 y -1 emitted as CH 4 and -337 g C m -2 y -1 taken up as CO 2 . From chambers, we discovered that neither site had large fluxes of N 2 O. Sustained-flux greenhouse gas accounting metrics indicated that both marshes had a positive (warming) radiative balance, with the brackish marsh having a substantially greater warming effect than the freshwater marsh. That net respiratory emissions of CO 2 and CH 4 as estimated through chamber techniques were 2-4 times different from emissions estimated through EC requires additional understanding of the artifacts created by different spatial and temporal sampling footprints between techniques.

Journal of Geophysical Research: Biogeosciences↗
USGS2016

Climate change and dissolved organic carbon export to the Gulf of Maine

Ongoing climate change is affecting the concentration, export (flux), and timing of dissolved organic carbon (DOC) exported to the Gulf of Maine (GoM) through changes in hydrologic regime. DOC export was calculated for water years 1950 through 2013 for 20 rivers and for water years 1930 through 2013 for 14 rivers draining to the GoM. DOC export was also estimated for the 21st century based on climate and hydrologic modeling in a previously published study. DOC export was calculated by using the regression model LOADEST to fit seasonally adjusted concentration discharge (C-Q) relations. Our results are an analysis of the sensitivity of DOC export to changes in hydrologic conditions over time since land cover and vegetation were held constant over time. Despite large interannual variability, all rivers had increasing DOC export during winter and these trends were significant ( p < 0.05) in 10 out of 20 rivers for 1950 to 2013 and in 13 out of 14 rivers for 1930 to 2013. All rivers also had increasing annual export of DOC although fewer trends were statistically significant than for winter export. Projections for DOC export during the 21st century were variable depending on the climate model and greenhouse gas emission scenario that affected future river discharge through effects on precipitation and evapotranspiration. The most consistent result was a significant increase in DOC export in winter in all model-by-emission scenarios. DOC export was projected to decrease during the summer in all model-by-emission scenarios, with statistically significant decreases in half of the scenarios.

Maine↗