T. Andrew Black

Observed increase in local cooling effect of deforestation at higher latitudes

Deforestation in mid- to high latitudes is hypothesized to have the potential to cool the Earth’s surface by altering biophysical processes. 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. This feedback is crucial in the model predictions; without it other biophysical processes may overwhelm the albedo effect to generate warming instead. 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. 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 ± 0.44 K (mean ± one standard deviation) northwards of 45° N and 0.21 ± 0.53 K southwards. Below 35° 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.

Remote sensing of photosynthetic-light-use efficiency of boreal forest

Using a helicopter-mounted portable spectroradiometer and continuous eddy covariance data we were able to evaluate the photochemical reflectance index (PRI) as an indicator of canopy photosynthetic light-use efficiency (LUE) in four boreal forest species during the Boreal Ecosystem Atmosphere experiment (BOREAS). PRI was calculated from narrow waveband reflectance data and correlated with LUE calculated from eddy covariance data. Significant linear correlations were found between PRI and LUE when the four species were grouped together and when divided into functional type: coniferous and deciduous. Data from the helicopter-mounted spectroradiometer were then averaged to represent data generated by the Airborne Visible Infrared Imaging Spectrometer (AVIRIS). We calculated PRI from these data and relationships with canopy LUE were investigated. The relationship between PRI and LUE was weakened for deciduous species but strengthened for the coniferous species. The robust nature of this relationship suggests that relative photosynthetic rates may be derived from remotely-sensed reflectance measurements. ©2000 Elsevier Science B.V. All rights reserved.

Atmospheric turbulence within and above a douglas-fir stand. Part II: Eddy fluxes of sensible heat and water vapour

This is the second paper describing a study of the turbulence regimes and exchange processes within and above an extensive Douglas-fir stand. The experiment was conducted on Vancouver Island during a two-week rainless period in July and August 1990. Two eddy correlation units were operated in the daytime to measure the fluxes of sensible heat and water vapour and other turbulence statistics at various heights within and above the stand. Net radiation was measured above the overstory using a stationary net radiometer and beneath the overstory using a tram system. Supplementary measurements included soil heat flux, humidity above and beneath the overstory, profiles of wind speed and air temperature, and the spatial variation of sensible heat flux near the forest floor.The sum of sensible and latent heat fluxes above the stand accounted for, on average, 83% of the available energy flux. On some days, energy budget closure was far better than on others. The average value of the Bowen ratio was 2.1 above the stand and 1.4 beneath the overstory. The mid-morning value of the canopy resistance was 150–450 s/m during the experiment and mid-day value of the Omega factor was about 0.20. The daytime mean canopy resistance showed a strong dependence on the mean saturation deficit during the two-week experimental period.The sum of sensible and latent heat fluxes beneath the overstory accounted for 74% of the available energy flux beneath the overstory. One of the reasons for this energy imbalance was that the small number of soil heat flux plates and the short pathway of the radiometer tram system was unable to account for the large horizontal heterogeneity in the available energy flux beneath the overstory. On the other hand, good agreement was obtained among the measurements of sensible heat flux made near the forest floor at four positions 15 m apart.There was a constant flux layer in the trunk space, a large flux divergence in the canopy layer, and a constant flux layer above the stand. Counter-gradient flux of sensible heat constantly occurred at the base of the canopy.The transfer of sensible heat and water vapour was dominated by intermittent cool downdraft and warm updraft events and dry downdraft and moist updraft events, respectively, at all levels. For sensible heat flux, the ratio of the contribution of cool downdrafts to that of warm updrafts was greater than one in the canopy layer and less than one above the stand and near the forest floor.

Looking deeper into the soil: biophysical controls and seasonal lags of soil CO2 production and efflux.

We seek to understand how biophysical factors such as soil temperature (Ts), soil moisture (theta), and gross primary production (GPP) influence CO2 fluxes across terrestrial ecosystems. Recent advancements in automated measurements and remote-sensing approaches have provided time series in which lags and relationships among variables can be explored. The purpose of this study is to present new applications of continuous measurements of soil CO2 efflux (F0) and soil CO2 concentrations measurements. Here we explore how variation in Ts, theta, and GPP (derived from NASAs moderate-resolution imaging spectroradiometer [MODIS]) influence F0 and soil CO2 production (Ps). We focused on seasonal variation and used continuous measurements at a daily timescale across four vegetation types at 13 study sites to quantify: (1) differences in seasonal lags between soil CO2 fluxes and Ts, theta, and GPP and (2) interactions and relationships between CO2 fluxes with Ts, theta, and GPP. Mean annual Ts did not explain annual F0 and Ps among vegetation types, but GPP explained 73% and 30% of the variation, respectively. We found evidence that lags between soil CO2 fluxes and Ts or GPP provide insights into the role of plant phenology and information relevant about possible timing of controls of autotrophic and heterotrophic processes. The influences of biophysical factors that regulate daily F0 and Ps are different among vegetation types, but GPP is a dominant variable for explaining soil CO2 fluxes. The emergence of long-term automated soil CO2 flux measurement networks provides a unique opportunity for extended investigations into F0 and Ps processes in the near future.

Coherent eddies and temperature structure functions for three contrasting surfaces. Part I: Ramp model with finite microfront time

Air temperature time series within and above canopies reveal ramp patternsassociated with coherent eddies that are responsible for most of thevertical transport of sensible heat. Van Atta used a simple step-changeramp model to analyse the coherent part of air temperature structurefunctions. However, his ocean data, and our own measurements for aDouglas-fir forest, straw mulch, and bare soil, reveal that even withoutlinearization his model cannot account for the observed decrease of thecubic structure function for small time lag. We found that a ramp model inwhich the rapid change at the end of the ramp occurs in a finite microfronttime can describe this decrease very well, and predict at least relativemagnitudes of microfront times between different surfaces. Averagerecurrence time for ramps, determined by analysis of the cubic structurefunction with the new ramp model, agreed well with values determined usingthe Mexican Hat wavelet transform, except at lower levels within theforest. Ramp frequency above the forest and mulch scaled very well withwind speed at the canopy top divided by canopy height. Within the forest,ramp frequency did not vary systematically with height. This is inaccordance with the idea that large-scale canopy turbulence is mostlygenerated by instability of the mean canopy wind profile, similar to aplane mixing layer. The straw mulch and bare soil experiments uniquelyextend measurements of temperature structure functions and ramp frequencyto the smallest scales possible in the field.

Linking foliage spectral responses to canopy-level ecosystem photosynthetic light-use efficiency at a Douglas-fir forest in Canada.

The light-use efficiency (LUE) of a mature Canadian Douglas-fir forest (DF49) was studied using high-resolution in situ temporal, spatial, and spectral measurements in conjunction with fluxes acquired from an instrumented tower. We examined the photochemical reflectance index (PRI), a spectral index responsive to high light conditions that alters reflectance at 531 nm, in combination with several alternative reference bands at 551, 570, and 488 nm. These indices were derived from directional reflectance spectra acquired by a hyperspectral radiometer system mounted on the DF49 tower, viewing the canopy through almost 360° rotations multiple times an hour daily throughout the 2006 growing season. From canopy structure information, three foliage sectors within the canopy were delineated according to instantaneous illumination conditions (sunlit, shaded, and mixed shaded–sunlit). On sunny days, the PRI indices for the sunlit foliage sector captured high light-induced stress responses, expressed as significantly different PRI values for sunlit versus shaded foliage. This difference was not observed on highly diffuse or overcast days. PRIs on sunny days tracked the diurnal photoregulation responses throughout the growing season in concert with illumination intensity. We computed the effective instantaneous LUE for the three foliage groups (LUEfoliage) using modeled and measured information. We provide convincing evidence that LUEfoliage can be well described and strongly related to all variations of the PRI within this coniferous forest under relatively clear skies (0.59 > r2 > 0.80, P < 0.0001). LUEfoliage varied through the growing season between 0.015 and 0.075 µmol C µmol–1 absorbed photosynthetically active radiation (APAR), and the lowest daily values were associated with the sunlit foliage group. The mixed sunlit-shaded foliage was the only group to exhibit monthly averages close to the maximum ecosystem LUE parameter (εmax) used in LUE models for evergreen needle forests (0.0196 µmol C µmol–1 APAR). Implications for remote sensing of carbon uptake dynamics and the interaction of canopy structure and physiology are discussed.

Climatic and Phenological Controls of the Carbon and Energy Balances of Three Contrasting Boreal Forest Ecosystems in Western Canada

Seasonal and interannual variability in the carbon and energy cycles of boreal forests are controlled by the interaction of climate, ecophysiology and plant phenology. This study analyses eddy-covariance data from mature trembling aspen, black spruce and jack pine stands in western Canada. The seasonal cycles of the surface carbon and energy balances were tightly coupled to the seasonal cycle of soil temperature. The contiguous carbon-uptake period was ∼50 days longer for the black spruce and jack pine stands than the trembling aspen stand, with 30 days difference in spring and 20 days difference in autumn. The black spruce and jack pine carbon-uptake period spanned the warm season, with gross ecosystem photosynthesis beginning during spring thaw and continuing until air temperature dropped to below freezing in autumn. In contrast, the trembling aspen carbon-uptake period was determined by the timing of leaf emergence and senescence, which occurred well after spring thaw and before autumn freeze. Regression analysis identified spring temperature as the primary factor controlling annual net ecosystem production at all three sites, through its influence on the onset of the growing season. Precipitation and soil water content had significant but secondary influences on the annual carbon fluxes. The impact of spring warming on annual net ecosystem production was 2–3 times greater at the deciduous-broadleaf than the evergreen-coniferous sites, confirming the high sensitivity of boreal deciduous-broadleaf forests to spring warming. The analysis confirmed the pivotal role of phenology in the response of northern ecosystems to climate variability and change.

Coherent eddies and temperature structure functions for three contrasting surfaces. Part II: Renewal model for sensible heat flux

Sensible heat, latent heat, and other scalar fluxes cannot be measuredwithin short dense canopies, e.g., straw mulches, with standard approachessuch as eddy correlation, Bowen ratio-energy balance, aerodynamic, andvariance methods. However, recently developed surface renewal models, thatare based on the fact that most of the turbulent transfer within and abovecanopies is associated with large-scale coherent eddies, which are evidentas ramp patterns in scalar time series, offer a feasible solution. Wepresent a new air renewal model that calculates sensible heat flux atdifferent heights within and above a canopy from the average cubictemperature structure function, sampled at a moderate rate, and measuredaverage friction velocity. The model is calibrated and tested with datameasured above and within a Douglas-fir forest and above a straw mulch andbare soil. We show that the model describes half-hour variations ofsensible heat flux very well, both within the canopy and roughnesssublayers and in the inertial sublayer, for stable and unstable atmosphericconditions. The combined empirical coefficient that appears in the modelhas an apparently universal value of about 0.4 for all surfaces andheights, which makes application of the model particularly simple. Themodel is used to predict daytime and nighttime sensible heat flux profileswithin the straw mulch and within a small bare opening in the mulch.

Carbon Dioxide Exchange and Nocturnal Processes Over a Mixed Deciduous Forest

This paper reports the results of the analysis of CO2 exchange from a one-month experiment conducted at a mixed deciduous forest, Camp Borden (80°65′W, 44°19′ N), Canada, in the summer of 1993. The mid-day CO2 flux from the forest under clear sky conditions was around −1.0 mg m−2 s−1, the average light and water use efficiencies 13 mmol CO2(mol photon)−1 and 7.95 mg CO2(gH2O)−1 , and the average nocturnal respiration rate 0.21 mg CO2m−2s−1. We observed different flow features at heights of 34.5 (14.5 m above the canopy) and 22.4 m at night. Wavelike structures were frequently encountered at z = 34.5 m. Depending on the phase angle between the vertical velocity and CO2 concentration time series, they could act to enhance the co-gradient (upward) flux or to create counter-gradient (downward) flux of CO2. We speculate that the wave events were limited to isolated regions in the upwind direction. Near the tree-tops (z = 22.4 m), the strong wind shear was able to maintain turbulence. Inverse temperature ramp structures were very common and flux of sensible heat was well behaved (directed downward).

Thermal optimality of net ecosystem exchange of carbon dioxide and underlying mechanisms.

• It is well established that individual organisms can acclimate and adapt to temperature to optimize their functioning. However, thermal optimization of ecosystems, as an assemblage of organisms, has not been examined at broad spatial and temporal scales. • Here, we compiled data from 169 globally distributed sites of eddy covariance and quantified the temperature response functions of net ecosystem exchange (NEE), an ecosystem-level property, to determine whether NEE shows thermal optimality and to explore the underlying mechanisms. • We found that the temperature response of NEE followed a peak curve, with the optimum temperature (corresponding to the maximum magnitude of NEE) being positively correlated with annual mean temperature over years and across sites. Shifts of the optimum temperature of NEE were mostly a result of temperature acclimation of gross primary productivity (upward shift of optimum temperature) rather than changes in the temperature sensitivity of ecosystem respiration. • Ecosystem-level thermal optimality is a newly revealed ecosystem property, presumably reflecting associated evolutionary adaptation of organisms within ecosystems, and has the potential to significantly regulate ecosystem-climate change feedbacks. The thermal optimality of NEE has implications for understanding fundamental properties of ecosystems in changing environments and benchmarking global models.