Effects of harvest, fire, and pest/pathogen disturbances on the West Cascades ecoregion carbon balance

Abstract

Background

Disturbance is a key influence on forest carbon dynamics, but the complexity of spatial and temporal patterns in forest disturbance makes it difficult to quantify their impacts on carbon flux over broad spatial domains. Here we used a time series of Landsat remote sensing images and a climate-driven carbon cycle process model to evaluate carbon fluxes at the ecoregion scale in western Oregon.

Results

Thirteen percent of total forest area in the West Cascades ecoregion was disturbed during the reference interval (1991-2010). The disturbance regime was dominated by harvesting (59 % of all area disturbed), with lower levels of fire (23 %), and pest/pathogen mortality (18 %). Ecoregion total Net Ecosystem Production was positive (a carbon sink) in all years, with greater carbon uptake in relatively cool years. Localized carbon source areas were associated with recent harvests and fire. Net Ecosystem Exchange (including direct fire emissions) showed greater interannual variation and became negative (a source) in the highest fire years. Net Ecosystem Carbon Balance (i.e. change in carbon stocks) was more positive on public that private forestland, because of a lower disturbance rate, and more positive in the decade of the 1990s than in the warmer and drier 2000s because of lower net ecosystem production and higher direct fire emissions in the 2000s.

Conclusion

Despite recurrent disturbances, the West Cascades ecoregion has maintained a positive carbon balance in recent decades. The high degree of spatial and temporal resolution in these simulations permits improved attribution of regional carbon sources and sinks.

Turner et al. Carbon Balance and Management (2015) 10:12 DOI 10.1186/s13021-015-0022-9

RESEARCH Open Access

Effects of harvest, fire, and pest/pathogen disturbances on the West Cascades ecoregion carbon balance

David P Turner1*, William D Ritts1, Robert E Kennedy2, Andrew N Gray3 and Zhiqiang Yang1

Abstract

Background: Disturbance is a key influence on forest carbon dynamics, but the complexity of spatial and temporal patterns in forest disturbance makes it difficult to quantify their impacts on carbon flux over broad spatial domains. Here we used a time series of Landsat remote sensing images and a climate-driven carbon cycle process model to evaluate carbon fluxes at the ecoregion scale in western Oregon.

Results: Thirteen percent of total forest area in the West Cascades ecoregion was disturbed during the reference interval (1991-2010). The disturbance regime was dominated by harvesting (59 % of all area disturbed), with lower levels of fire (23 %), and pest/pathogen mortality (18 %). Ecoregion total Net Ecosystem Production was positive (a carbon sink) in all years, with greater carbon uptake in relatively cool years. Localized carbon source areas were associated with recent harvests and fire. Net Ecosystem Exchange (including direct fire emissions) showed greater interannual variation and became negative (a source) in the highest fire years. Net Ecosystem Carbon Balance(i.e. change in carbon stocks) was more positive on public that private forestland, because of a lower disturbance rate, and more positive in the decade of the 1990s than in the warmer and drier 2000s because of lower net ecosystem production and higher direct fire emissions in the 2000s.

Conclusion: Despite recurrent disturbances, the West Cascades ecoregion has maintained a positive carbon balance in recent decades. The high degree of spatial and temporal resolution in these simulations permits improved attribution of regional carbon sources and sinks.

Keywords: Forests; Carbon; Net ecosystem production; Net ecosystem exchange; Net ecosystem carbon balance; Disturbance; West Cascades ecoregion

Background

Net uptake of carbon by forests provides a significant offset to anthropogenic carbon emissions at the global[1], national [2, 3], regional [4], and landscape [5] scales.

However, forest carbon sinks are vulnerable to disturbances in the form of harvesting, fire, and pest/pathogen outbreaks. At the regional scale, we have a poor understanding of the relative contribution of these disturbances to overall carbon budgets [6], but such knowledge is important in understanding how the carbon cycle is responding to on-going management and climate change

[7, 8]. It is also critical for developing policies for greenhouse gas mitigation through altered land use [9].

Forest disturbances are a strong determinant of carbon stocks and fluxes on both managed and unmanaged landscapes [5, 10]. Clear-cut harvesting, as is commonly practiced in coniferous forests of western Oregon, shuts down the photosynthetic carbon sink and increases the carbon source from heterotrophic respiration of harvest residues. Partial harvests for thinning likewise induce a near term reduction in carbon sequestration [11]. Wildfire is similar to harvesting in reducing carbon uptake, but has a longer term impact on heterotrophic respiration because of the slow conversion of snags to more readily decomposed woody debris on the ground [12, 13]. Pest/pathogen outbreaks reduce leaf area and leave slow decomposing snags, thus altering ecosystem carbon flux for decades [14,

* Correspondence: mailto:[email protected]

Web End [email protected]

1Department of Forest Ecosystems and Society, Oregon State University, 97331 Corvallis, OR, USAFull list of author information is available at the end of the article

2015 Turner et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0

Web End =http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Turner et al. Carbon Balance and Management (2015) 10:12 Page 2 of 12

15]. Spatially-explicit carbon cycle assessments in the western U.S. have generally emphasized effects of harvests and fire [1618] and not explicitly captured the impact of slow pest/pathogen disturbances.

Satellite remote sensing, particularly from the Landsat series of sensors (~30 m resolution), offers the opportunity to monitor forest disturbances [19, 20]. In combination with spatially-distributed ecosystem process models that simulate carbon cycle responses to specific disturbances, remote sensing data can be used to map and monitor forest carbon stocks and flux [21, 22]. Here, we take advantage of a new Landsat-based time series analysis of forest disturbance (LandTrendr) [23, 24] and a well-established modeling infrastructure for simulating regional carbon flux based on the Biome-BGC carbon cycle process model [25] to quantify carbon cycle impacts of harvesting, fire, and pest/pathogen outbreaks on forests of the West Cascades (WC) ecoregion in the Northwestern U.S. We make extensive use of plot scale and aggregated U.S.D.A. Forest Service Forest Inventory and Analysis (FIA) data [26, 27] for calibration and validation of both LandTrendr and Biome-BGC.

We quantify several distinct carbon fluxes [28] at the ecoregion scale. Net ecosystem production (NEP) is the balance of net primary production (NPP) and heterotrophic respiration (Rh). It reflects ecosystem metabolism as it responds to variation in weather, and to disturbance events. Net ecosystem exchange (NEE) is the absolute vertical flux of CO2 over a given geographical domain. This is the flux as seen by a continental to global scale inversion, e.g. [29]. In addition to NEP and direct fire emissions, it includes river/stream evasion as well as emissions associated with harvested products [30]. Net ecosystem carbon balance (NECB) refers to the absolute change in carbon stocks, and is affected by NEP as well as removals in the form of harvested products, lateral transfers of dissolved organics, and by direct fire emissions. NECB is the equivalent of carbon sequestration as would be relevant to offsetting fossil fuel emissions. For the purposes of comparisons here, we report NEE using the same convention as with NEP and NECB, i.e. a positive value is a carbon sink. The capacity to isolate these fluxes is required to fully understand the role of forests and forest management in regulating the atmospheric CO2 concentration.

Results

Domain characterization

Our study domain was the West Cascades ecoregion in western Oregon, U.S.A. It is characterized by a strong elevation gradient from west to east, with corresponding gradients in temperature and precipitation (Fig. 1). Land cover is predominantly conifer forest (Fig. 2a). Forest stand age tends to be <60 on the low elevation private lands that are managed for wood production (Figs. 2b,

and 3). On public lands at higher elevations, there is a

broader range of stand ages.

Disturbance patterns

The total proportion of the study area that was disturbed during the 1991-2010 interval was 13 %. The proportion of total disturbed area attributed to harvest was 59 %. The corresponding proportion for fire was 23 %, and for pest/pathogen outbreak was 18 %. The location of the harvests have been predominantly on lower elevation private forestland (Fig. 4a). The time series of annual area harvested shows a decrease on public lands in the late 1980s and an increase on private lands in the most recent decade (Fig. 5). On public lands, there was a shift from stand replacing harvests to partial harvests whereas on private forestland stand replacing harvests were most common over the whole time series (Fig. 6).

The incidence of fire was low in the 1980s and 1990s but increased appreciably in the 2000s (Fig. 7), with fires located at both high and low elevations (Fig. 4b). The overall proportion of fires at high, medium, and low intensity was 28 %, 33 %, and 40 % respectively. These proportions did not change much from year to year.

Areas of pest/pathogen disturbance occurred primarily at high elevations (Fig. 6c). Intensity tended to be higher in the years with relatively large areas disturbed.

Carbon flux

In our stand-level simulations of pest/pathogen disturbance, NPP falls in parallel with the drop in stem and foliar biomass and begins to recover after the year of maximum intensity (e.g. Fig. 8). NEP correspondingly decreases, falling below zero in the case of a short, high intensity slow disturbance.

In the ecoregion-wide simulations, the majority of the surface area of the WC ecoregion had a positive NEP during the study period (e.g. Fig. 9). NEP tends to decrease as elevation increases (Fig. 10) because of decreasing rates of wood productivity associated with a shorter growing season, and a shift towards older, slower growing, age classes on public lands (Fig. 2b). Carbon sources (negative NEP) are indicated in areas recently burned or harvested.

The ecoregion total NEP was also generally positive (Fig. 11). Total direct fire emissions increased in the 2000s relative to the 1990s, but only exceeded NEP in 2003. Over the 1991-2010 interval, harvest removals offset 28 % of NEP and direct fire emissions offset 7 % of NEP.

The interannual variation in NEE is a function of both NEP and direct fire emissions. In 2003, which was a high fire year (Fig. 7), high temperatures and soil drought reduced NPP more than Rh, and hence reduced NEP.

Direct fire emissions were also relatively high, causing a dip in NEE (Fig. 11).

Turner et al. Carbon Balance and Management (2015) 10:12 Page 3 of 12

Fig. 1 Climate in the study area: a annual precipitation, b annual average temperature

The ecoregion total NECB averaged 1.9 (SD = 1.3) TgC yr-1 over the 1991-2010 interval. Only in 2003 did it fall below zero. Thus, the ecoregion has been a sustained sink for atmospheric CO2 in recent decades. The area weighted mean NECB (1991-2010) for public lands was 65 gC m-2 year-1 compared to 8 gC m-2 year-1 on private forestland, reflecting a proportionally lower harvest rate on public lands in recent years.

Discussion

Harvest

Since much of the WC ecoregion is public forestland, the reduction in harvests associated with implementation of the Northwest Forest Plan (NWFP) in the early 1990s [31] had a notable impact on the overall rate of harvesting [24, 32]. The effect has been stabilization of a long-term decline in the proportion of public lands that is in the old-growth condition [33].

Fig. 2 The West Cascades domain: a land cover, b stand age

Turner et al. Carbon Balance and Management (2015) 10:12 Page 4 of 12

Fig. 3 The distribution of public and private forestland

The rate of harvest on private forestland also shifted over the course of the study period, but in this case it was an increase that reflected a period of economic growth and high demand for wood products. The ob-served harvest rate of 1.5 % per year on private forest-land over the 2000-2010 period is consistent with the 4560 year rotation that is common in the area.

The notable shift on public lands from high intensity (clear-cut) harvesting to low and moderate intensities (thinning) reflects policy changes. Beginning in the 1990s, litigation largely prevented further harvesting of old growth stands. However, quotas established by the NWFP did allow for a low level of harvesting. Thinning of young to mature stands thus became the standard practice [34, 35], as reflected in the LandTrendr observations. On private lands that are managed primarily for wood production, the practice of clear-cut harvesting has remained the standard.

Fire

The fire regime in the WC ecoregion has varied over the last several hundred years. Fire scar data and tree age

class distributions indicate a period of high fire about 1500, possibly associated with a relatively warm climate [36, 37]. During the settlement era in the late 19th century, the incidence of fire increased because of anthropogenic factors. That period was followed by a large decrease in the incidence of fire in the 20th century, associated with successful fire suppression. Most recently, an increase in the incidence of forest wildfire has been noted over much of the western U.S. and attributed in part to climate warming [38]. The post-2000 increase in area burned in this study is consistent with this broad pattern.

The sharp increase in the area burned in the WC eco-region in the 2000s (3.6 % per year) compared to the 1990s (0.4 % per year) is associated with a 20 % decrease in mean May through September precipitation (Fig. 12). There was also a minor increase (0.2 C) in mean May through September temperature. This natural experiment in interdecadal climate variation mimics to some degree the summer precipitation trends expected in the region for the 21st Century [39]. The observations here lend support to projections of an increased incidence of fire over the course of the 21st Century in the Cascade Mountains [40, 41]. Note that detecting relationships between temporal trends in burned area and climate is potentially confounded by changes in policy with respect to managing fire [42].

Pests/pathogens

Our trajectory-based LandTrendr algorithm is complimentary to the U.S. Forest Service airborne surveys for pests/pathogens in that it better resolves spatial heterogeneity and is sensitive to severity [43]. The duration of the slow disturbances is also informative with respect to identifying the relevant organisms: western spruce budworm show a more consistent long-term decline in the reference spectral vegetation index compared to mountain pine beetle [44].

The elevated incidence of pest/pathogen disturbances around 1990 is associated with a multiyear outbreak of western spruce budworm at relatively high elevations. This outbreak was widespread over the western U.S. [45,44]. After 1995, the annual area disturbed became stable, presumably at an endemic population level characteristic of the native pests/pathogens. Meigs et al. [44] examined the incidence of area in the Cascade Mountains subject to pest/pathogen impacts and later burned. While that sequence was not uncommon, there did not appear to be a strong relationship. Indeed, mountain pine beetle outbreaks in lodgepole pine forests, may reduce the probability of active crown fire [46].

Carbon flux

NEP. In forest ecosystems, NEP is generally negative after a stand replacing harvest (because NPP is low and

Turner et al. Carbon Balance and Management (2015) 10:12 Page 5 of 12

Fig. 4 Location of disturbances: a harvest, b fire, c pest/pathogen

Rh from decay of residues is high), becomes moderately positive for multiple decades while wood mass accumulates, and then falls to near zero in late succession [47,48]. Stand-level simulations with Biome-BGC are in agreement with that pattern [17, 49].

The stand age for crossover from carbon source to carbon sink varies widely in the WC ecoregion because of differences in the amount of residues present at the time of stand initiation [50]. The post disturbance pulse of Rh and negative NEP is of greater magnitude in the

case of harvest but of greater duration in the case of fire because the proportion of the biomass actually burned in a wildfire is often small [51] and standing dead trees (snags) decay relatively slowly [13].

Moderate to low intensity abrupt disturbances (e.g. thinning) may also introduce periods of negative NEP [52]. Lost leaf area means a smaller fraction of photosynthetically active radiation is captured by the canopy, and heterotrophic respiration is boosted by the metabolism of the pests/ pathogens themselves [53] or by heterotrophic respiration

Fig. 5 Area harvested by year for public and private forestland (1986-2010)

Turner et al. Carbon Balance and Management (2015) 10:12 Page 6 of 12

Fig. 6 Trends in harvest magnitude (1986-2010): a public forestland, b private forestland

of the increased leaf and fine root litter [54]. In our Biome-BGC simulations, the time to recovery of positive NEP is closely tied to the magnitude of the disturbance.

Here we also treat slow disturbances presumed to be caused by pests/pathogens. These disturbances are similar to fires in leaving large amounts of dead wood bio-mass behind that are slowly respired by heterotrophs and hence reduce NEP [14].

The most controversial aspect of the generalized trajectory for NEP over the course of succession in the WC ecoregion is the degree to which relatively old stands remain a carbon sink [55, 56]. Given that maritime Douglas-fir trees can survive for over 1000 years[57], and possibly benefit from increasing temperature and CO2 concentration, a long term carbon sink might be expected. Mean simulated NEP (2001-2009) for old growth stands in the study area was 70 gC m-2 year-1, which compares to and eddy covariance based mean of 49 gC m-2 year-1 (1998-2008) at the Wind River old growth Douglas-fir site [58].

Interannual variation in climate variables drives the large interannual variation in ecoregion-wide NEP in

our simulations. Observations at eddy covariance flux towers in the Cascade Mountains indicate that both gross primary production and ecosystem respiration contribute to interannual variation in NEP [58, 59], with smaller carbon sinks associated with warmer, drier years. Tree ring observations find that relatively warm years increase tree growth at high elevations [60, 61] and reduce it at mid to low elevations [62] in the Cascade Mountains. Thus, spatially distributed simulations are needed to evaluate if projected climate change driven warming and possible summer drying in the Pacific Northwest[39] would act to increase or decrease ecoregion mean NEP [7, 63].

NEE. The declining trend of NEE between the 1990s and 2000s was driven by both a decrease in NEP (dominant factor) and an increase in direct fire emissions. Schwalm et al. [7] also report a drop in regional NEP associated with turn of the century drought in western North America. Previous Biome-BGC simulations in Oregon suggest that both NPP and Rh have declined in recent decades, with a greater decline in NPP [64]. However, there remains considerable uncertainty about the

Fig. 7 Area burned by year for total forestland (1986-2010)

Turner et al. Carbon Balance and Management (2015) 10:12 Page 7 of 12

Fig. 10 Relationship of mean Net Ecosystem Production (2006-2010) to elevation band

Fig. 8 Representative simulation results for a disturbed stand. Conditions are high intensity and short duration pest/pathogen disturbance beginning in 1990 at a mid-elevation site. a Stem mass, b Net Primary Production (NPP) and Net Ecosystem Production (NEP)

impact of changing climate and CO2 concentration on regional tree growth. Long term observations at 21 temperate and boreal zone eddy covariance flux towers support a broadly realized CO2 benefit on forest water use efficiency [65], which may offset the impacts of lower summer precipitation.

By estimating NEE, we are in a position to compare our ecoregion fluxes with independent scaling approaches. Our mean NEE for 2004 was 148 gC m-2 years-1 (a carbon sink), much larger that an estimate from the Carbon-Tracker inversion [29, 66] of 18 gC m-2 years-1. Differences of this magnitude have been observed in other NEE comparisons [67, 68] and point to the need to reconcile results from alternative scaling approaches [69]. Land-Trendr could potentially deliver detailed information on the disturbance regime over much larger domains than the present study, thus opening the possibility for improved regional carbon budgets.

NECB. The positive NECB for the ecoregion using our flux scaling approach is consistent with carbon sinks in western Oregon reported using other modeling approaches [70] as well as inventory based approaches [35, 71, 72]. An additional forest sector sink is associated with harvested wood removals that accumulate in long-lived wood products [73] but is not estimated here. The high carbon accumulation rate on public lands, as driven by a policy-based lower harvest rate, supports the inclusion of carbon sequestration in the suite of ecosystem services to be considered in public forestland management.

Forest carbon stocks and flux are now widely recognized as highly relevant to mitigating the on-going rise in atmospheric CO2 [74] and the combined remote sensing/

modeling approach described here offers the opportunity

Fig. 9 Net ecosystem production over the West Cascades ecoregion (mean for 2006-2010)

Turner et al. Carbon Balance and Management (2015) 10:12 Page 8 of 12

to monitor forest carbon budgets at the fine spatial and temporal resolution that may be needed for quantifying carbon sinks. Ongoing efforts to characterize uncertainty associated with the remote sensing of disturbance and with stand-age-specific carbon fluxes (e.g. [75]) will increase the policy relevance of these carbon flux maps. The end points of our analysis are also the prerequisites for a realistic landscape simulation into the future that accounts for climate change and land use [76].

Conclusions

Harvests are the dominant form of disturbance in the West Cascades ecoregion, followed by fire, and pests/ pathogens. The majority of the WC landscape has a positive NEP most of the time, and annual total NEP has been positive in recent decades. In high fire years, eco-region total NEE can fall below zero (become a source) because of low NEP and high direct fire emissions. Harvest removals offset 4 times more NEP than do direct fire emissions. The sustained carbon accumulation on the West Cascades land base contributes to the regional forest carbon sink.

Methods

Overview

Our NEP/NECB scaling methods and validation results for applications in Oregon, Washington, and California have been reported in Turner et al. [18, 21, 22, 64] and Law et al. [77, 78]. These previous studies considered only stand-replacing disturbances and employed simple Landsat-based change detection algorithms to specify the approximate year of disturbance. Here we added treatment of pest/pathogen outbreaks (slow disturbance), allowed for multiple disturbance intensity classes, and accounted for sequential disturbance events based on an annual time series of Landsat data. To quantify NEE, we added a previously derived estimate of emissions from harvested products (wood and crops) to NEP and fire emissions.

NEP modeling

Our primary tool for scaling carbon flux across the eco-region is the Biome-BGC carbon cycle process model[25]. The model has a daily time step and is run over multiple years to simulate succession. Simulated carbon cycle processes include photosynthesis, autotrophic respiration, heterotrophic respiration, plant C allocation, and mortality. Simulated C pools include stemwood, coarse roots, fine roots, foliage, litter, coarse woody debris, snags, and soil organic matter. A water balance is calculated based on the Penman-Montieth formulation of evapotranspiration. For each model run, there is a spin-up using a 25-year repeating loop of climate data to bring the soil carbon pool into near equilibrium with the

Fig. 11 Annual ecoregion totals for (a) net ecosystem production (NEP), (b) fire emissions, (c) harvest removals, (d) net ecosystem exchange (NEE), (e) net ecosystem carbon balance (NECB)

Fig. 12 Anomaly in May 1 through September 30 precipitation

Turner et al. Carbon Balance and Management (2015) 10:12 Page 9 of 12

climate. At its end, one or two disturbance events are prescribed by year, type, duration, and intensity to bring the simulation up to the current condition.

A look-up table determines the partitioning of the extant carbon stocks existing at the time of disturbance into removals, direct fire emissions, and transfers of necromass from one ecosystem component to another, thus maintaining ecosystem mass balance (Table 1). The disturbance history of a stand is prescribed in terms of one or two disturbances based on the record of Landsat imagery.

The version of Biome-BGC used here was adapted from version 4.1.2 [25] to simulate stand-replacing disturbances [77], dynamic allocation over the course of succession [78], and mixed severity fire [17]. To accommodate partial disturbances, the proportional disturbance intensity from satellite data (see below) is applied directly to the live carbon pools (Table 1).

The initiation year, duration, and maximum magnitude are prescribed by LandTrendr in the case of pest/pathogen disturbances. A linear ramp from initiation year to the maximum magnitude year is used to prescribe an annual mortality amount such that the prescribed maximum mortality is achieved in the maximum mortality year (Fig. 8). The stand recovers from disturbance prognostically. Transfers of dead foliar and fine root C are made to the litter pools and transfers of tree C to a standing dead (snag) pool. There is subsequent transfer to coarse woody debris on the ground [7981]. Partitioning factors at the time of disturbance are based on observations within the region [51, 8284].

Biome-BGC has 20 cover-type specific ecophysiological parameters that must be specified. For the most part, we adopted the recommendations in White et al. [85] and a representative set of parameters is given in Turner et al.[64]. In the case of the evergreen needle leaf cover type, we did an ecoregion specific adjustment on two of the parameters (the fraction of leaf nitrogen as rubisco, and the mortality fraction) that had been identified in sensitivity analyses as strongly impacting wood mass and NEP. Reference data for these parameter selections were site- and age-specific estimates of net stem growth and wood mass from USDA Forest Service FIA plots [35]. We ran the model at all FIA plot locations (approximated within 500 m to preclude potential disclosure of confidential locations) within the ecoregion to the stand age in the plot data, and compared predicted and observed net stem production using a range of possible parameter values. The optimal value for the mortality parameter was 0.0125 (proportion of biomass per year) and for the fraction of leaf nitrogen as rubisco parameter was 0.035 (unitless).

Model inputs

Our base land cover dataset (Fig. 2a) is the 2011 National Land Cover Database (NLCD) [86] which is derived from Landsat data. Areas that had been harvested and were classified by NLCD as shrubland were reclassified as forest. The distribution of public and private land (Fig. 3) was from U.S. Geological Survey [87].

The attributes of the disturbance regime were from the LandTrendr analysis of Landsat Thematic Mapper time series data [23, 24]. It captures both abrupt events, such as fire and harvest, and slow (multiyear)

Table 1 Partitioning of biomass pools at the time of disturbance

Biomass componentLeaf/fine root Stem Coarse root Litter CWD Snag Disturbance DisturbanceType MagnitudeHarvest Low 0.25a 0.25b 0.25c

Medium 0.55a 0.55b 0.55cHigh 1.00a 1.00b 1.00cFired Low 0.125 0.02e 0.02f 0.66 0.17 0.11

Medium 0.50 0.03e 0.03f 0.66 0.22 0.14 High 1.00 0.05e 0.05f 1.00 0.39 0.18 Pest/Pathogen Low 0.30a 0.30g 0.30c

High 0.75a 0.75g 0.75c Values are proportions of total biomass. CWD = coarse woody debris

aTransferred to litter

bStem transferred off-site, associated branches to go coarse woody debris (CWD)

cTransferred to coarse woody debris

dAll transfers are to the atmosphere

eResidual mortality is transferred to snag

fResidual mortality is transferred to coarse woody debris

gTransferred to snag

Turner et al. Carbon Balance and Management (2015) 10:12 Page 10 of 12

disturbance processes caused by pests/pathogens. Each disturbance event is characterized by year, type, magnitude, and duration. The year is indicated by an inflection point in the trajectory of a spectral vegetation index in each pixel-level Landsat time series. The type of disturbance can be harvest (thinning), fire, or pest/pathogen driven. The distinction for abrupt disturbances between fire and harvest is based on reference to the Monitoring Trends in Burn Severity (MTBS) dataset [88]. Variables in the LandTrendr output are reported in a continuous format but here we binned the outputs into 3 magnitude classes for fire and harvest (see Table 1), and 2 magnitude classes (see Table 1), each with 3 possible durations (9, 17, 25 years), for slow disturbances. Bin midpoints are assigned to represent stand age, disturbance magnitude, and disturbance duration.

The LandTrendr analysis covered the interval from 1985 to 2012, but here we used multiyear buffers at the beginning and end when reporting carbon fluxes to limit artifacts associated with identifying inflection points in slow disturbances.

Disturbances previous to 1986 were prescribed on the basis of mapped stand age class from gradient nearest neighbor analysis (GNN) [89]. GNN integrates data from Landsat and Forest Inventory and Analysis plots [26] to estimate stand age. Conifer stands not having LandTrendr-based disturbance since 1985 were binned into 4 age classes and assigned the age class midpoint (45, 80, 150, 250 years). The binning of disturbance magnitude and duration as well as stand age class was necessary to constrain the number of unique disturbance histories within each 1 km climate grid cell. Model runs were made for the 10 most frequent combinations of cover type and disturbance history in each 1 km grid cell, and fluxes were reported as the weighted mean per 1 km2 cell. Using this approach, we covered >90 % of the forest area under consideration.

The meteorological inputs to Biome-BGC are daily minimum and maximum temperature, precipitation, humidity, and solar radiation. We obtained 25 years (1986-2010) of climate data at 1 km resolution over our study ecoregion through the North American Carbon Program [90]. The distributed climate data (Fig. 1) are based on interpolations of meteorological station observations using a digital elevation model and general meteorological principles [9193]. Uncertainty in the interpolations has previously been evaluated in our region and elsewhere [94, 95].

Soil Data Input. Soil texture and depth are specified from US Geological Survey soil maps [96].

Carbon flux reporting

NEP was calculated as NPP minus Rh. NEE was NEP minus direct fire emissions and an estimate for emissions

from harvested products [97, 68]. This later estimate was made by reference to the change in the stock of previously harvested wood and crop products (based on inputs associated with harvests and product-specific turnover rates) and population distribution [97]. River stream evasion is also a significant term in regional NEE, however much of it derives from CO2 originating in Rh of soil organic matter that is flushed into streams with the flow of the soil solution [98]. Since our Biome-BGC simulation includes Rh of soil organic matter, we did not attempt to additionally account for river/stream evasion here. NECB was calculated as NEP minus fire emissions and harvest removals (from our own simulations). Our NECB did not explicitly account for land use change, but the rate of land use conversion is quite low in the WC ecoregion [72].

AbbreviationsFIA: Forest inventory and analysis; GNN: Gradient nearest neighbor analysis; MTBS: Monitoring trends in burn severity; NEP: Net ecosystem production; NECB: Net ecosystem carbon balance; NEE: Net ecosystem exchange;

Rh: Heterotrophic respiration; NPP: Net primary production; NWFP: Northwest forest plan; U.S.: United States; USDA: United States Department of Agriculture; WC: West Cascades.

Competing interestsThe authors declare that they have no competing interests.

Authors contributionsDT designed the study and contributed to data analysis and writing; WR contributed to assembling the model inputs and implementing model runs; RK provided the stand age map and developed the disturbance maps; AG organized and analyzed the FIA data; and ZY contributed to model coding and implementing the model runs. All authors read and approved the final manuscript.

AcknowledgmentsSupport was provided by the NASA Terrestrial Ecology Program (NNX12AK59G). Climate data from the North American Carbon Program and data from the U.S.D.A. Forest Service Forest Inventory and Analysis Program were essential to this study.

Author details

1Department of Forest Ecosystems and Society, Oregon State University, 97331 Corvallis, OR, USA. 2College of Earth, Ocean, and Atmospheric Sciences, Oregon State University, 97331 Corvallis, OR, USA. 3USDA Forest Service, Pacific Northwest Station, 97331 Corvallis, OR, USA.

Received: 5 March 2015 Accepted: 8 May 2015

References1. Pan YD, Birdsey RA, Fang JY, Houghton R, Kauppi PE, Kurz WA, et al. A large and persistent carbon sink in the worlds forests. Science. 2011;333(6045):98893.

2. Stinson G, Kurz WA, Smyth CE, Neilson ET, Dymond CC, Metsaranta JM, et al. An inventory-based analysis of Canadas managed forest carbon dynamics, 1990 to 2008. Global Change Biol. 2010;17(6):222744.

3. McKinley DC, Ryan MG, Birdsey RA, Giardina CP, Harmon ME, Heath LS, et al. A synthesis of current knowledge on forests and carbon storage in the United States. Ecol Appl. 2011;21(6):190224.

4. Lu XL, Kicklighter DW, Melillo JM, Yang P, Rosenzweig B, Vorosmarty CJ, et al. A contemporary carbon balance for the Northeast Region of the United States. Environ Sci Technol. 2013;47(23):132308.

5. Chen B, Arain MA, Khomik M, Trofymow JA, Grant RF, Kurz WA, et al. Evaluating the impacts of climate variability and disturbance regimes on the historic carbon budget of a forest landscape. Agr Forest Meteorol. 2013;180:26580.

Turner et al. Carbon Balance and Management (2015) 10:12 Page 11 of 12

6. Liu SG, Bond-Lamberty B, Hicke JA, Vargas R, Zhao SQ, Chen J, et al. Simulating the impacts of disturbances on forest carbon cycling in North America: processes, data, models, and challenges. J Geophys Res Biogeosci. 2011;116:22.

7. Schwalm CR, Williams CA, Schaefer K, Baldocchi D, Black TA, Goldstein AH, et al. Reduction in carbon uptake during turn of the century drought in western North America. Nat Geosci. 2012;5(8):5516.

8. Williams CA, Collatz GJ, Masek J, Goward SN. Carbon consequences of forest disturbance and recovery across the conterminous U.S. Global Biogeochem Cy. 2012;26:GB1005.

9. USCCSP. United States Carbon Cycle Science Plan, http://downloads.globalchange.gov/carbon-cycle/us-carbon-cycle-science-plan.pdf

Web End =http:// http://downloads.globalchange.gov/carbon-cycle/us-carbon-cycle-science-plan.pdf

Web End =downloads.globalchange.gov/carbon-cycle/us-carbon-cycle-science- http://downloads.globalchange.gov/carbon-cycle/us-carbon-cycle-science-plan.pdf

Web End =plan.pdf . 2011.

10. Densai AR, Moorcroft PR, Bolstad PV, Davis KJ. Regional carbon fluxes from an observationally constrained dynamic ecosystem model: impacts of disturbance, CO2 fertilization, and heterogeneous land cover. J Geophys Res Biogeosci. 2007;112:G1.

11. Campbell J, Alberti G, Martin J, Law BE. Carbon dynamics of a ponderosa pine plantation following a thinning treatment in the northern Sierra Nevada. Forest Ecol Manag. 2009;257(2):45363.

12. Harmon ME, Bond-Lamberty B, Tang JW, Vargas R. Heterotrophic respiration in disturbed forests: a review with examples from North America. J Geophys Res Biogeosci. 2011;116:17.

13. Acker SA, Kertis J, Bruner H, OConnell K, Sexton J. Dynamics of coarse woody debris following wildfire in a mountain hemlock (Tsuga mertensiana) forest. Forest Ecol Manag. 2013;302:2319.

14. Edburg SL, Hicke JA, Lawrence DM, Thornton PE. Simulating coupled carbon and nitrogen dynamics following mountain pine beetle outbreaks in the western United States. J Geophys Res Biogeosci. 2011;116:15.

15. Kurz WA, Dymond CC, Stinson G, Rampley GJ, Neilson ET, Carroll AL, et al. Mountain pine beetle and forest carbon feedback to climate change. Nature. 2008;452:98790.

16. Cohen WB, Harmon ME, Wallin DO, Fiorella M. Two decades of carbon flux from forests of the Pacific Northwest. Bio Sci. 1996;46:83644.

17. Meigs GW, Turner DP, Ritts WD, Yang ZQ, Law BE. Landscape-scale simulation of heterogeneous fire effects on pyrogenic carbon emissions, tree mortality, and net ecosystem production. Ecosystems. 2011;14(5):75875.

18. Turner DP, Ritts WD, Yang ZQ, Kennedy RE, Cohen WB, Duane MV, et al. Decadal trends in net ecosystem production and net ecosystem carbon balance for a regional socioecological system. Forest Ecol Manag. 2011;262:131825.

19. Cohen WB, Goward SN. Landsats role in ecological applications of remote sensing. Bio Sci. 2004;54(6):53545.

20. Goward SN, Masek JG, Cohen WB, Moisen G, Collatz GJ, Healey SP, et al. Forest disturbance and North American carbon flux. EOS. 2008;89:10516.

21. Turner DP, Ollinger SV, Kimball JS. Integrating remote sensing and ecosystem process models for landscape to regional scale analysis of the carbon cycle. Bio Sci. 2004;54:57384.

22. Turner DP, Guzy M, Lefsky M, Ritts W, VanTuyl S, Law BE. Monitoring forest carbon sequestration with remote sensing and carbon cycle modeling. Environ Manage. 2004;4:45766.

23. Kennedy RE, Yang ZG, Cohen WB. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 1. LandTrendr - Temporal segmentation algorithms. Remote Sens Environ. 2010;114(12):2897910.

24. Kennedy RE, Yang ZQ, Cohen WB, Pfaff E, Braaten J, Nelson P. Spatial and temporal patterns of forest disturbance and regrowth within the area of the Northwest Forest Plan. Remote Sens Environ. 2012;122:11733.

25. Thornton PE, Law BE, Gholz HL, Clark KL, Falge E, Ellsworth DS, et al. Modeling and measuring the effects of disturbance history and climate on carbon and water budgets in evergreen needleleaf forests. Agr Forest Meteorol. 2002;113:185222.

26. Bechtold WA, Patterson PL. The enhanced Forest Inventory and Analysis program - national sampling design and estimation procedures, USDA Forest Service General Technical Report SRS-GTR-80. 2005.

27. Woudenberg SW, Conkling BL, OConnell BM, LaPoint EB, Turner JA, Waddell KL. The Forest Inventory and Analysis database: Description and user manual version 4.0 for Phase 2, USDA Forest Service, General Techncial Report RMRS-GTR-245. 2010.

28. Chapin FS, Woodwell GM, Randerson JT, Rastetter EB, Lovett GM, Baldocchi DD, et al. Reconciling carbon-cycle concepts, terminology, and methods. Ecosystems. 2006;9(7):104150.

29. Peters W, Jacobson AR, Sweeney C, Andrews AE, Conway TJ, Masarie K,et al. An atmospheric perspective on North American carbon dioxide exchange: CarbonTracker. Proc Natl Acad Sci U S A. 2007;104(48):1892530.

30. Hayes DJ, Turner DP. The need for apples-to-apples comparisons of carbon dioxide source and sink estimates. EOS. 2012;93:4045.

31. USDA. Record of Decision for amendments to Forest Service and Bureau of Land Management planning documents within the range of the northern spotted owl. 74 p. (plus Attachment A: Standards and guidelines). 1994.

32. Cohen WB, Spies TA, Alig RJ, Oetter DR, Maiersperger TK, Fiorella M. Characterizing 23 years (1972-95) of stand replacement disturbance in western Oregon forests with Landsat imagery. Ecosystems. 2002;5:12237.

33. Moeur M, Ohmann JL, Kennedy RE, Cohen WB, Gregory MJ, Yang Z, et al. Status and trends of late-successional and old-growth forests, General Technical Report PNW-GTR-853. 2011.

34. Franklin JF, Johnson KN. A restoration framework for federal forests in the Pacific Northwest. J Forest. 2012;110(8):42939.

35. Gray AN, Whittier TR. Carbon stocks and changes on Pacific Northwest national forests and the role of disturbance, management, and growth. Forest Ecol Manag. 2014;328:16778.

36. Weisberg PJ, Swanson FJ. Regional synchroneity in fire regimes of western Oregon and Washington. USA Forest Ecol Manag. 2003;172(1):1728.

37. Weisberg PJ. Historical fire frequency on contrasting slope facets along the McKenzie River, Western Oregon Cascades. West North Am Naturalist. 2009;69(2):20614.

38. Westerling AL, Hidalgo HG, Cayan DR, Swetnam TW. Warming and earlier spring increases western U.S. forest wildfire activity. Science. 2006;313:9403.

39. Mote PW, Salathe EP. Future climate in the Pacific Northwest. Clim Change. 2010;102(1-2):2950.

40. Littell JS, Oneil EE, McKenzie D, Hicke JA, Lutz JA, Norheim RA, et al. Forest ecosystems, disturbance, and climatic change in Washington State, USA. Climatic Change. 2010;102(1-2):12958.

41. Rogers BM, Neilson RP, Drapek R, Lenihan JM, Wells JR, Bachelet D, et al. Impacts of climate change on fire regimes and carbon stocks of the U.S. Pacific Northwest. J Geophys Res Biogeosci. 2011;116:13.

42. NCWFMS. The National Strategy (National Cohesive Wildland Fire Management Strategy). United States Department of Interior and Department of Agriculture http://www.forestsandrangelands.gov/strategy/thestrategy.shtml

Web End =http://www.forestsandrangelands.gov/strategy/thestrategy.shtml , accessed: May 1, 2015. 2014.

43. Meigs GW, Kennedy RE, Gray AN, Gregory MJ. Spatiotemporal dynamics of recent mountain pine beetle and western spruce budworm outbreaks across the Pacific Northwest region, USA. Remote Sens Environ. 2015;339:7186.

44. Meigs GW, Kennedy RE, Cohen WB. A Landsat time series approach to characterize bark beetle and defoliator impacts on tree mortality and surface fuels in conifer forests. Remote Sens Environ. 2011;115(12):370718.

45. Hummel S, Agee JK. Western spruce budworm defoliation effects on forest structure and potential fire behavior. Northwest Sci. 2003;77(2):15969.

46. Simard M, Romme WH, Griffin JM, Turner MG. Do mountain pine beetle outbreaks change the probability of active crown fire in lodgepole pine forests? Ecol Monogr. 2011;81(1):324.

47. Grant RF, Barr AG, Black TA, Margolis HA, McCaughey JH, Trofymow JA. Net ecosystem productivity of temperate and boreal forests after clearcutting-a Fluxnet-Canada measurement and modelling synthesis. Tellus B. 2010;62(5):47596.

48. Kasischke ES, Amiro BD, Barger NN, French NHF, Goetz SJ, Grosse G, et al. Impacts of disturbance on the terrestrial carbon budget of North America. J Geophys Res Biogeosci. 2013;118(1):30316.

49. Turner DP, Guzy M, Lefsky MA, Van Tuyl S, Sun O, Daly C, et al. Effects of land use and fine-scale environmental heterogeneity on net ecosystem production over a temperate coniferous forest landscape. Tellus B. 2003;55:65768.

50. Janisch JE, Harmon ME. Successional changes in live and dead wood carbon stores: implications for net ecosystem productivity. Tree Physiol. 2002;22:7789.

51. Campbell J, Donato D, Azuma DL, Law B. Pyrogenic carbon emission from a large wildfire in Oregon United States. J Geophys Res Biogeosci. 2007;12:G04014.

52. Harmon ME, Moreno A, Domingo JB. Effects of partial harvest on the carbon stores in Douglas-fir/Western Hemlock forests: a simulation study. Ecosystems. 2009;12(5):77791.

53. Schmitz OJ, Raymond PA, Estes JA, Kurz WA, Holtgrieve GW, Ritchie ME, et al. Animating the carbon cycle. Ecosystems. 2014;17(2):34459.

54. Dore S, Montes-Helu M, Hart SC, Hungate BA, Koch GW, Moon JB, et al. Recovery of ponderosa pine ecosystem carbon and water fluxes from thinning and stand-replacing fire. Global Change Biol. 2012;18(10):317185.

Turner et al. Carbon Balance and Management (2015) 10:12 Page 12 of 12

55. Harmon ME, Bible K, Ryan MG, Shaw DC, Chen H, Klopatek J, et al. Production, respiration, and overall carbon balance in an old-growth Pseudotsuga-Tsuga forest ecosystem. Ecosystems. 2004;7(5):498512.

56. Luyssaert S, Schulze ED, Borner A, Knohl A, Hessenmoller D, Law BE, et al. Old-growth forests as global carbon sinks. Nature. 2008;455(7210):2135.

57. Waring RH, Franklin JF. Evergreen forests of the Pacific Northwest. Science. 1979;204:13806.

58. Wharton S, Falk M, Bible K, Schroeder M, Paw KT. Old-growth CO2 flux measurements reveal high sensitivity to climate anomalies across seasonal, annual and decadal time scales. Agr Forest Meteorol. 2012;161:114.

59. Falk M, Wharton S, Schroeder M, Ustin SL, Paw UKT. Flux partitioning in an old-growth forest: seasonal and interannual dynamics. Tree Physiol. 2008;28:50920.

60. Case MJ, Peterson DL. Fine-scale variability in growth-climate relationships of Douglas-fir, North Cascade Range, Washington. Can J For Res-Rev Can Rech For. 2005;35(11):274355.

61. Case MJ, Peterson DL. Growth-climate relations of lodgepole pine in the North Cascades National Park. Washington Northwest Sci. 2007;81(1):6275.

62. Beedlow PA, Lee EH, Tingey DT, Waschmann RS, Burdick CA. The importance of seasonal temperature and moisture patterns on growth of Douglas-fir in western Oregon. USA Agr Forest Meteorol. 2013;169:17485.

63. Wharton S, Chasmer L, Falk M, U KTP. Strong links between teleconnections and ecosystem exchange found at a Pacific Northwest old-growthforest from flux tower and MODIS EVI data. Global Change Biol. 2009;15(9):2187205.

64. Turner DP, Ritts WD, Law BE, Cohen WB, Yang Z, Hudiburg T, et al. Scaling net ecosystem production and net biome production over a heterogeneous region in the western United States. Biogeosci. 2007;4:597612.

65. Keenan TF, Hollinger DY, Bohrer G, Dragoni D, Munger JW, Schmid HP, et al. Increase in forest water-use efficiency as atmospheric carbon dioxide concentrations rise. Nature. 2013;499(7458):324.

66. CT. CarbonTracker. http://www.esrl.noaa.gov/gmd/ccgg/carbontracker/

Web End =http://www.esrl.noaa.gov/gmd/ccgg/carbontracker/ . 2014.

67. Desai AR, Helliker BR, Moorcroft PR, Andrews AE, Berry JA. Climatic controls of interannual variability in regional carbon fluxes from top-down and bottom-up perspectives. J Geophys Res Biogeosci. 2010;115:15.

68. Turner DP, Jacobson AR, Ritts WD, Wang WL, Nemani R. A large proportion of North American net ecosystem production is offset by emissions from harvested products, river/stream evasion, and biomass burning. Global Change Biol. 2013;19(11):351628.

69. Turner DP, Gockede M, Law BE, Ritts WD, Cohen WB, Yang Z, et al. Multiple constraint analysis of regional land-surface carbon flux. Tellus. 2011;63B:20721.

70. Zhu Z, Reed BC, editors. Baseline and projected future carbon storage and greenhouse-gas fluxes in ecosystems of the Western United States. U.S. Geological Survey Professional Paper 1797. 2012.

71. Alig RJ, Krankina ON, Yost A, Kuzminykh J. Forest carbon dynamics in the Pacific Northwest (USA) and the St. Petersburg region of Russia: comparisons and policy implications. Clim Change. 2006;79:33560.

72. Gray AN, Whittier TR, Azuma DL. Estimation of aboveground forest carbon flux in Oregon: adding components of change to stock-difference assessments. Forest Sci. 2014;60(2):31726.

73. Malmsheimer RW, Bowyer JL, Fried JS, Gee E, Izlar RL, Miner RA, et al. Managing forests because carbon matters: integrating energy, products, and land management policy. J Forest. 2011;109(7):S748.

74. Olander LP, Cooley DM, Galik CS. The potential role for management of U.S. public lands in greenhouse gas mitigation and climate policy. Environ Manage. 2012;49(3):52333.

75. Williams CA, Collatz GJ, Masek J, Huang CQ, Goward SN. Impacts of disturbance history on forest carbon stocks and fluxes: Merging satellite disturbance mapping with forest inventory data in a carbon cycle model framework. Remote Sens Environ. 2014;151:5771.

76. Metsaranta JM, Dymond CC, Kurz WA, Spittlehouse DL. Uncertainty of 21st century growing stocks and GHG balance of forests in British Columbia, Canada resulting from potential climate change impacts on ecosystem processes. Forest Ecol Manag. 2011;262(5):82737.

77. Law BE, Turner D, Campbell J, Van Tuyl S, Ritts WD, Cohen WB. Disturbance and climate effects on carbon stocks and fluxes across Western Oregon USA. Global Change Biol. 2004;10:142944.

78. Law BE, Turner DP, Lefsky M, Campbell J, Guzy M, Sun O, et al. Carbon fluxes across regions: observational constraints at multiple scales, Scaling

and uncertainty analysis in ecology: methods and applications. New York: Columbia University Press; 2006.79. Busse MD. Downed bole-wood decomposition in Lodgepole pine forests of central Oregon. Soil Sci Soc Am J. 1994;58(1):2217.

80. Mitchell RG, Priesler HK. Fall rate of lodgepole pine killed by the mountain pine beetle in Central Oregon. West J Appl For. 1998;13:236.

81. Brown M, Black TA, Nesic Z, Foord VN, Spittlehouse DL, Fredeen AL, et al. Impact of mountain pine beetle on the net ecosystem production of lodgepole pine stands in British Columbia. Agr Forest Meteorol. 2011;150(2):25464.

82. Turner DP, Koerper GJ, Harmon ME, Lee JJ. A carbon budget for forests of the conterminous United States. Ecol Appl. 1995;5:42136.

83. Law BE, Sun OJ, Campbell J, Van Tuyl S, Thornton PE. Changes in carbon storage and fluxes in a chronosequence of ponderosa pine. Global Change Biol. 2003;9(4):51024.

84. Sun OJ, Campbell J, Law BE, Wolf V. Dynamics of carbon stocks in soils and detritus across chronosequences of different forest types in the Pacific Northwest. USA Global Change Biol. 2004;10(9):147081.

85. White MA, Thornton PE, Running SW, Nemani RR. Parameterization and sensitivity analysis of the BIOME-BGC terrestrial ecosystem model: net primary production controls. Earth Interact. 2000;4:185.

86. Jin S, Yang L, Danielson P, Homer C, Fry J, Xian G. A comprehensive change detection method for updating the National land cover database to circa 2011. Remote Sens Environ. 2013;132:15975.

87. GAP. US Geological Survey, Gap Analysis Program (GAP). National Land Cover, Version 2. http://gapanalysis.usgs.gov/gaplandcover/data/

Web End =http://gapanalysis.usgs.gov/gaplandcover/data/ . 2014.

88. Eidenshink J, Schwind B, Brewer K, Zhu ZL, Quayle B, Howard SM. A project for monitoring trends in burn severity. Fire Ecol. 2007;3:321.

89. Ohmann JL, Gregory MJ. Predictive mapping of forest composition and structure with direct gradient analysis and nearest-neighbor imputation in coastal Oregon. USA Can J For Res. 2002;32:72541.

90. Daymet: Daily Surface Weather Data on a 1-km Grid for North America, Version2. Data set. Available on-line [http://daac.ornl.gov

Web End =http://daac.ornl.gov] from Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee, USA [database on the Internet] 2014. Available from: http://dx.doi.org/10.3334/ORNLDAAC/1219

Web End =http://dx.doi.org/10.3334/ http://dx.doi.org/10.3334/ORNLDAAC/1219

Web End =ORNLDAAC/1219 . Accessed: January 1, 2015.91. Thornton PE, Running SW, White MA. Generating surfaces of daily meteorological variables over large regions of complex terrain. J Hydrol. 1997;190(190):21451.

92. Thornton PE, Running SW. An improved algorithm for estimating incident daily solar radiation from measurements of temperature, humidity, and precipitation. Agr Forest Meteorol. 1999;93:21128.

93. Thornton PE, Hasenauer H, White MA. Simultaneous estimation of daily solar radiation and humidity from observed temperature and precipitation: an application over complex terrain in Austria. Agr Forest Meteorol. 2000;104:25571.

94. Hasenauer H, Merganicova K, Petritsch R, Pietsch SA, Thornton PE. Validating daily climate interpolations over complex terrain in Austria. Agr Forest Meteorol. 2003;119:87107.

95. Daly C, Halbleib M, Smith JI, Gibson WP, Doggett MK, Taylor GH, et al. Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States. Int J Climatol. 2008;28(15):203164.

96. CONUS. Conterminous United States multi-layer soil characteristics data set for regional climate and hydrology modeling. 2007. http://www.soilinfo.psu.edu/index.cgi?soil_data&conus

Web End =http://www.soilinfo.p http://www.soilinfo.psu.edu/index.cgi?soil_data&conus

Web End =su.edu/index.cgi?soil_data&conus .

97. Hayes DJ, Turner DP, Stinson G, McGuire AD, Wei YX, West TO, et al. Reconciling estimates of the contemporary North American carbon balance among terrestrial biosphere models, atmospheric inversions, and a new approach for estimating net ecosystem exchange from inventory-based data. Global Change Biol. 2012;18(4):128299.

98. Butman D, Raymond PA. Significant efflux of carbon dioxide from streams and rivers in the United States. Nat Geosci. 2011;4(12):83942.