Deforestation in mid-to high latitudes is hypothesized to have the potential to cool the Earth's surface by altering biophysical processes [1][2][3] . In climate models of continental-scale land clearing, the cooling is triggered by increases in surface albedo and is reinforced by a land albedo-sea ice feedback 4,5 . This feedback is crucial in the model predictions; without it other biophysical processes may overwhelm the albedo effect to generate warming instead 5 . Ongoing land-use activities, such as land management for climate mitigation, are occurring at local scales (hectares) presumably too small to generate the feedback, and it is not known whether the intrinsic biophysical mechanism on its own can change the surface temperature in a consistent manner 6,7 . Nor has the effect of deforestation on climate been demonstrated over large areas from direct observations. Here we show that surface air temperature is lower in open land than in nearby forested land. The effect is 0.85 6 0.44 K (mean 6 one standard deviation) northwards of 456 N and 0.21 6 0.53 K southwards. Below 356 N there is weak evidence that deforestation leads to warming. Results are based on comparisons of temperature at forested eddy covariance towers in the USA and Canada and, as a proxy for small areas of cleared land, nearby surface weather stations. Night-time temperature changes unrelated to changes in surface albedo are an important contributor to the overall cooling effect. The observed latitudinal dependence is consistent with theoretical expectation of changes in energy loss from convection and radiation across latitudes in both the daytime and night-time phase of the diurnal cycle, the latter of which remains uncertain in climate models 8 .The latitudinal gradient of land-use impact is evident in the comparison of the surface air temperature recorded at FLUXNET (www.fluxnet.ornl.gov) forest towers 9 (Supplementary Table 1 and Supplementary Fig. 1) and surface weather stations in North America (Fig. 1a). Here we use the surface stations as proxies for cleared land. In accordance with the requirement of the World Meteorological Organization, these stations are located in open grassy fields that have biophysical characteristics similar to those of open land, such as being covered by snow in northern latitudes in the winter 10 . Latitude accounts for 31% of the variations in the temperature difference DT between the forest sites and the adjacent open lands (number of site pairs n 5 37). The rate of change in DT with latitude is 20.070 6 0.010 K per degree (mean 6 one standard error, s.e., P , 0.005). At these sites, the annual net all-wave radiation R n

The FLUXNET2015 dataset provides ecosystem-scale data on CO 2 , water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

[1] Disturbances are important for renewal of North American forests. Here we summarize more than 180 site years of eddy covariance measurements of carbon dioxide flux made at forest chronosequences in North America. The disturbances included standreplacing fire (Alaska, Arizona, Manitoba, and Saskatchewan) and harvest (British Columbia, Florida, New Brunswick, Oregon, Quebec, Saskatchewan, and Wisconsin) events, insect infestations (gypsy moth, forest tent caterpillar, and mountain pine beetle), Hurricane Wilma, and silvicultural thinning (Arizona, California, and New Brunswick). Net ecosystem production (NEP) showed a carbon loss from all ecosystems following a stand-replacing disturbance, becoming a carbon sink by 20 years for all ecosystems and by 10 years for most. Maximum carbon losses following disturbance (g C m −2 y −1 ) ranged from 1270 in Florida to 200 in boreal ecosystems. Similarly, for forests less than 100 years old, maximum uptake (g C m −2 y −1) was 1180 in Florida mangroves and 210 in boreal ecosystems. More temperate forests had intermediate fluxes. Boreal ecosystems were relatively time invariant after 20 years, whereas western ecosystems tended to increase in carbon gain over time. This was driven mostly by gross photosynthetic production (GPP) because total ecosystem respiration (ER) and heterotrophic respiration were relatively invariant with age. GPP/ER was as low as 0.2 immediately following stand-replacing disturbance reaching a constant value of 1.2 after 20 years. NEP following insect defoliations and silvicultural thinning showed lesser changes than stand-replacing events, with decreases in the year of disturbance followed by rapid recovery. NEP decreased in a mangrove ecosystem following Hurricane Wilma because of a decrease in GPP and an increase in ER.

High-latitude ecosystems store approximately 1700 Pg of soil carbon (C), which is twice as much C as is currently contained in the atmosphere. Permafrost thaw and subsequent microbial decomposition of permafrost organic matter could add large amounts of C to the atmosphere, thereby influencing the global C cycle. The rates at which C is being released from the permafrost zone at different soil depths and across different physiographic regions are poorly understood but crucial in understanding future changes in permafrost C storage with climate change. We assessed the inherent decomposability of C from the permafrost zone by assembling a database of long-term (>1 year) aerobic soil incubations from 121 individual samples from 23 high-latitude ecosystems located across the northern circumpolar permafrost zone. Using a three-pool (i.e., fast, slow and passive) decomposition model, we estimated pool sizes for C fractions with different turnover times and their inherent decomposition rates using a reference temperature of 5 °C. Fast cycling C accounted for less than 5% of all C in both organic and mineral soils whereas the pool size of slow cycling C increased with C : N. Turnover time at 5 °C of fast cycling C typically was below 1 year, between 5 and 15 years for slow turning over C, and more than 500 years for passive C. We project that between 20 and 90% of the organic C could potentially be mineralized to CO2 within 50 incubation years at a constant temperature of 5 °C, with vulnerability to loss increasing in soils with higher C : N. These results demonstrate the variation in the vulnerability of C stored in permafrost soils based on inherent differences in organic matter decomposability, and point toward C : N as an index of decomposability that has the potential to be used to scale permafrost C loss across landscapes.

Results from a systematic investigation of mercury (Hg) concentrations across 14 forest sites in the United States show highest concentrations in litter layers, strongly enriched in Hg compared to aboveground tissues and indicative of substantial postdepositional sorption of Hg. Soil Hg concentrations were lower than in litter, with highest concentrations in surface soils. Aboveground tissues showed no detectable spatial patterns, likely due to 17 different tree species present across sites. Litter and soil Hg concentrations positively correlated with carbon (C), latitude, precipitation, and clay (in soil), which together explained up to 94% of concentration variability. We observed strong latitudinal increases in Hg in soils and litter, in contrast to inverse latitudinal gradients of atmospheric deposition measures. Soil and litter Hg concentrations were closely linked to C contents, consistent with well-known associations between organic matter and Hg, and we propose that C also shapes distribution of Hg in forests at continental scales. The consistent link between C and Hg distribution may reflect a long-term legacy whereby old, C-rich soil and litter layers sequester atmospheric Hg depositions over long time periods. Based on a multiregression model, we present a distribution map of Hg concentrations in surface soils of the United States.