The role of environmental driving factors in historical and projected carbon dynamics of wetland ecosystems in Alaska.

Wetlands are critical terrestrial ecosystems in Alaska, covering ~177,000 km2 , an area greater than all the wetlands in the remainder of the United States. To assess the relative influence of changing climate, atmospheric carbon dioxide (CO2 ) concentration, and fire regime on carbon balance in wetland ecosystems of Alaska, a modeling framework that incorporates a fire disturbance model and two biogeochemical models was used. Spatially explicit simulations were conducted at 1-km resolution for the historical period (1950-2009) and future projection period (2010-2099). Simulations estimated that wetland ecosystems of Alaska lost 175 Tg carbon (C) in the historical period. Ecosystem C storage in 2009 was 5,556 Tg, with 89% of the C stored in soils. The estimated loss of C as CO2 and biogenic methane (CH4 ) emissions resulted in wetlands of Alaska increasing the greenhouse gas forcing of climate warming. Simulations for the projection period were conducted for six climate change scenarios constructed from two climate models forced under three CO2 emission scenarios. Ecosystem C storage averaged among climate scenarios increased 3.94 Tg C/yr by 2099, with variability among the simulations ranging from 2.02 to 4.42 Tg C/yr. These increases were driven primarily by increases in net primary production (NPP) that were greater than losses from increased decomposition and fire. The NPP increase was driven by CO2 fertilization (~5% per 100 parts per million by volume increase) and by increases in air temperature (~1% per °C increase). Increases in air temperature were estimated to be the primary cause for a projected 47.7% mean increase in biogenic CH4 emissions among the simulations (~15% per °C increase). Ecosystem CO2 sequestration offset the increase in CH4 emissions during the 21st century to decrease the greenhouse gas forcing of climate warming. However, beyond 2100, we expect that this forcing will ultimately increase as wetland ecosystems transition from being a sink to a source of atmospheric CO2 because of (1) decreasing sensitivity of NPP to increasing atmospheric CO2 , (2) increasing availability of soil C for decomposition as permafrost thaws, and (3) continued positive sensitivity of biogenic CH4 emissions to increases in soil temperature.

The role of driving factors in historical and projected carbon dynamics of upland ecosystems in Alaska.

It is important to understand how upland ecosystems of Alaska, which are estimated to occupy 84% of the state (i.e., 1,237,774 km2 ), are influencing and will influence state-wide carbon (C) dynamics in the face of ongoing climate change. We coupled fire disturbance and biogeochemical models to assess the relative effects of changing atmospheric carbon dioxide (CO2 ), climate, logging and fire regimes on the historical and future C balance of upland ecosystems for the four main Landscape Conservation Cooperatives (LCCs) of Alaska. At the end of the historical period (1950-2009) of our analysis, we estimate that upland ecosystems of Alaska store ~50 Pg C (with ~90% of the C in soils), and gained 3.26 Tg C/yr. Three of the LCCs had gains in total ecosystem C storage, while the Northwest Boreal LCC lost C (-6.01 Tg C/yr) because of increases in fire activity. Carbon exports from logging affected only the North Pacific LCC and represented less than 1% of the state's net primary production (NPP). The analysis for the future time period (2010-2099) consisted of six simulations driven by climate outputs from two climate models for three emission scenarios. Across the climate scenarios, total ecosystem C storage increased between 19.5 and 66.3 Tg C/yr, which represents 3.4% to 11.7% increase in Alaska upland's storage. We conducted additional simulations to attribute these responses to environmental changes. This analysis showed that atmospheric CO2 fertilization was the main driver of ecosystem C balance. By comparing future simulations with constant and with increasing atmospheric CO2 , we estimated that the sensitivity of NPP was 4.8% per 100 ppmv, but NPP becomes less sensitive to CO2 increase throughout the 21st century. Overall, our analyses suggest that the decreasing CO2 sensitivity of NPP and the increasing sensitivity of heterotrophic respiration to air temperature, in addition to the increase in C loss from wildfires weakens the C sink from upland ecosystems of Alaska and will ultimately lead to a source of CO2 to the atmosphere beyond 2100. Therefore, we conclude that the increasing regional C sink we estimate for the 21st century will most likely be transitional.

Implications of allometric model selection for county-level biomass mapping.

BACKGROUND: Carbon accounting in forests remains a large area of uncertainty in the global carbon cycle. Forest aboveground biomass is therefore an attribute of great interest for the forest management community, but the accuracy of aboveground biomass maps depends on the accuracy of the underlying field estimates used to calibrate models. These field estimates depend on the application of allometric models, which often have unknown and unreported uncertainties outside of the size class or environment in which they were developed. RESULTS: Here, we test three popular allometric approaches to field biomass estimation, and explore the implications of allometric model selection for county-level biomass mapping in Sonoma County, California. We test three allometric models: Jenkins et al. (For Sci 49(1): 12-35, 2003), Chojnacky et al. (Forestry 87(1): 129-151, 2014) and the US Forest Service's Component Ratio Method (CRM). We found that Jenkins and Chojnacky models perform comparably, but that at both a field plot level and a total county level there was a ~ 20% difference between these estimates and the CRM estimates. Further, we show that discrepancies are greater in high biomass areas with high canopy covers and relatively moderate heights (25-45 m). The CRM models, although on average ~ 20% lower than Jenkins and Chojnacky, produce higher estimates in the tallest forests samples (> 60 m), while Jenkins generally produces higher estimates of biomass in forests < 50 m tall. Discrepancies do not continually increase with increasing forest height, suggesting that inclusion of height in allometric models is not primarily driving discrepancies. Models developed using all three allometric models underestimate high biomass and overestimate low biomass, as expected with random forest biomass modeling. However, these deviations were generally larger using the Jenkins and Chojnacky allometries, suggesting that the CRM approach may be more appropriate for biomass mapping with lidar. CONCLUSIONS: These results confirm that allometric model selection considerably impacts biomass maps and estimates, and that allometric model errors remain poorly understood. Our findings that allometric model discrepancies are not explained by lidar heights suggests that allometric model form does not drive these discrepancies. A better understanding of the sources of allometric model errors, particularly in high biomass systems, is essential for improved forest biomass mapping.

Historical and projected trends in landscape drivers affecting carbon dynamics in Alaska.

Modern climate change in Alaska has resulted in widespread thawing of permafrost, increased fire activity, and extensive changes in vegetation characteristics that have significant consequences for socioecological systems. Despite observations of the heightened sensitivity of these systems to change, there has not been a comprehensive assessment of factors that drive ecosystem changes throughout Alaska. Here we present research that improves our understanding of the main drivers of the spatiotemporal patterns of carbon dynamics using in situ observations, remote sensing data, and an array of modeling techniques. In the last 60 yr, Alaska has seen a large increase in mean annual air temperature (1.7°C), with the greatest warming occurring over winter and spring. Warming trends are projected to continue throughout the 21st century and will likely result in landscape-level changes to ecosystem structure and function. Wetlands, mainly bogs and fens, which are currently estimated to cover 12.5% of the landscape, strongly influence exchange of methane between Alaska's ecosystems and the atmosphere and are expected to be affected by thawing permafrost and shifts in hydrology. Simulations suggest the current proportion of near-surface (within 1 m) and deep (within 5 m) permafrost extent will be reduced by 9-74% and 33-55% by the end of the 21st century, respectively. Since 2000, an average of 678 595 ha/yr was burned, more than twice the annual average during 1950-1999. The largest increase in fire activity is projected for the boreal forest, which could result in a reduction in late-successional spruce forest (8-44%) and an increase in early-successional deciduous forest (25-113%) that would mediate future fire activity and weaken permafrost stability in the region. Climate warming will also affect vegetation communities across arctic regions, where the coverage of deciduous forest could increase (223-620%), shrub tundra may increase (4-21%), and graminoid tundra might decrease (10-24%). This study sheds light on the sensitivity of Alaska's ecosystems to change that has the potential to significantly affect local and regional carbon balance, but more research is needed to improve estimates of land-surface and subsurface properties, and to better account for ecosystem dynamics affected by a myriad of biophysical factors and interactions.

Integrating LIDAR and forest inventories to fill the trees outside forests data gap.

Forest inventories are commonly used to estimate total tree biomass of forest land even though they are not traditionally designed to measure biomass of trees outside forests (TOF). The consequence may be an inaccurate representation of all of the aboveground biomass, which propagates error to the outputs of spatial and process models that rely on the inventory data. An ideal approach to fill this data gap would be to integrate TOF measurements within a traditional forest inventory for a parsimonious estimate of total tree biomass. In this study, Light Detection and Ranging (LIDAR) data were used to predict biomass of TOF in all "nonforest" Forest Inventory and Analysis (FIA) plots in the state of Maryland. To validate the LIDAR-based biomass predictions, a field crew was sent to measure TOF on nonforest plots in three Maryland counties, revealing close agreement at both the plot and county scales between the two estimates. Total tree biomass in Maryland increased by 25.5 Tg, or 15.6%, when biomass of TOF were included. In two counties (Carroll and Howard), there was a 47% increase. In contrast, counties located further away from the interstate highway corridor showed only a modest increase in biomass when TOF were added because nonforest conditions were less common in those areas. The advantage of this approach for estimating biomass of TOF is that it is compatible with, and explicitly separates TOF biomass from, forest biomass already measured by FIA crews. By predicting biomass of TOF at actual FIA plots, this approach is directly compatible with traditionally reported FIA forest biomass, providing a framework for other states to follow, and should improve carbon reporting and modeling activities in Maryland.

Local discrepancies in continental scale biomass maps: a case study over forested and non-forested landscapes in Maryland, USA.

BACKGROUND: Continental-scale aboveground biomass maps are increasingly available, but their estimates vary widely, particularly at high resolution. A comprehensive understanding of map discrepancies is required to improve their effectiveness in carbon accounting and local decision-making. To this end, we compare four continental-scale maps with a recent high-resolution lidar-derived biomass map over Maryland, USA. We conduct detailed comparisons at pixel-, county-, and state-level. RESULTS: Spatial patterns of biomass are broadly consistent in all maps, but there are large differences at fine scales (RMSD 48.5-92.7 Mg ha-1). Discrepancies reduce with aggregation and the agreement among products improves at the county level. However, continental scale maps exhibit residual negative biases in mean (33.0-54.6 Mg ha-1) and total biomass (3.5-5.8 Tg) when compared to the high-resolution lidar biomass map. Three of the four continental scale maps reach near-perfect agreement at ~4 km and onward but do not converge with the high-resolution biomass map even at county scale. At the State level, these maps underestimate biomass by 30-80 Tg in forested and 40-50 Tg in non-forested areas. CONCLUSIONS: Local discrepancies in continental scale biomass maps are caused by factors including data inputs, modeling approaches, forest/non-forest definitions and time lags. There is a net underestimation over high biomass forests and non-forested areas that could impact carbon accounting at all levels. Local, high-resolution lidar-derived biomass maps provide a valuable bottom-up reference to improve the analysis and interpretation of large-scale maps produced in carbon monitoring systems.

Integrating forest inventory and analysis data into a LIDAR-based carbon monitoring system.

BACKGROUND: Forest Inventory and Analysis (FIA) data may be a valuable component of a LIDAR-based carbon monitoring system, but integration of the two observation systems is not without challenges. To explore integration methods, two wall-to-wall LIDAR-derived biomass maps were compared to FIA data at both the plot and county levels in Anne Arundel and Howard Counties in Maryland. Allometric model-related errors were also considered. RESULTS: In areas of medium to dense biomass, the FIA data were valuable for evaluating map accuracy by comparing plot biomass to pixel values. However, at plots that were defined as "nonforest", FIA plots had limited value because tree data was not collected even though trees may be present. When the FIA data were combined with a previous inventory that included sampling of nonforest plots, 21 to 27% of the total biomass of all trees was accounted for in nonforest conditions, resulting in a more accurate benchmark for comparing to total biomass derived from the LIDAR maps. Allometric model error was relatively small, but there was as much as 31% difference in mean biomass based on local diameter-based equations compared to regional volume-based equations, suggesting that the choice of allometric model is important. CONCLUSIONS: To be successfully integrated with LIDAR, FIA sampling would need to be enhanced to include measurements of all trees in a landscape, not just those on land defined as "forest". Improved GPS accuracy of plot locations, intensifying data collection in small areas with few FIA plots, and other enhancements are also recommended.