Forest Flux Research Articles

Improved assessment of gross and net primary productivity of Canada's landmass

AbstractWe assess Canada's gross primary productivity (GPP) and net primary productivity (NPP) using boreal ecosystem productivity simulator (BEPS) at 250 m spatial resolution with improved input parameter and driver fields and phenology and nutrient release parameterization schemes. BEPS is a process‐based two‐leaf enzyme kinetic terrestrial ecosystem model designed to simulate energy, water, and carbon (C) fluxes using spatial data sets of meteorology, remotely sensed land surface variables, soil properties, and photosynthesis and respiration rate parameters. Two improved key land surface variables, leaf area index (LAI) and land cover type, are derived at 250 m from Moderate Resolution Imaging Spectroradiometer sensor. For diagnostic error assessment, we use nine forest flux tower sites where all measured C flux, meteorology, and ancillary data sets are available. The errors due to input drivers and parameters are then independently corrected for Canada‐wide GPP and NPP simulations. The optimized LAI use, for example, reduced the absolute bias in GPP from 20.7% to 1.1% for hourly BEPS simulations. Following the error diagnostics and corrections, daily GPP and NPP are simulated over Canada at 250 m spatial resolution, the highest resolution simulation yet for the country or any other comparable region. Total NPP (GPP) for Canada's land area was 1.27 (2.68) Pg C for 2008, with forests contributing 1.02 (2.2) Pg C. The annual comparisons between measured and simulated GPP show that the mean differences are not statistically significant (p > 0.05, paired t test). The main BEPS simulation error sources are from the driver fields.

Dry deposition of nitrogen compounds (NO<sub>2</sub>, HNO<sub>3</sub>, NH<sub>3</sub>), sulfur dioxide and ozone in west and central African ecosystems using the inferential method

Abstract. This work is part of the IDAF program (IGAC-DEBITS-AFRICA) and is based on the long-term monitoring of gas concentrations (1998–2007) established at seven remote sites representative of major African ecosystems. Dry deposition fluxes were estimated by the inferential method using on the one hand surface measurements of gas concentrations (NO2, HNO3, NH3, SO2 and O3) and on the other hand modeled exchange rates. Dry deposition velocities (Vd) were calculated using the big-leaf model of Zhang et al. (2003b). The bidirectional approach is used for NH3 surface–atmosphere exchange (Zhang et al., 2010). Surface and meteorological conditions specific to IDAF sites have been used in the models of deposition. The seasonal and annual mean variations of gaseous dry deposition fluxes (NO2, HNO3, NH3, O3 and SO2) are analyzed. Along the latitudinal transect of ecosystems, the annual mean dry deposition fluxes of nitrogen compounds range from −0.4 to −0.8 kg N ha−1 yr−1 for NO2, from −0.7 to −1.0 kg N ha−1 yr−1 for HNO3 and from −0.7 to −8.3 kg N ha−1 yr−1 for NH3 over the study period (1998–2007). The total nitrogen dry deposition flux (NO2+HNO3+NH3) is more important in forests (−10 kg N ha−1 yr−1) than in wet and dry savannas (−1.6 to −3.9 kg N ha−1 yr−1). The annual mean dry deposition fluxes of ozone range between −11 and −19 kg ha−1 yr−1 in dry and wet savannas, and −11 and −13 kg ha−1 yr−1 in forests. Lowest O3 dry deposition fluxes in forests are correlated to low measured O3 concentrations, lower by a factor of 2–3, compared to other ecosystems. Along the ecosystem transect, the annual mean of SO2 dry deposition fluxes presents low values and a small variability (−0.5 to −1 kg S ha−1 yr−1). No specific trend in the interannual variability of these gaseous dry deposition fluxes is observed over the study period.

Greenhouse gas fluxes (CO2, N2O and CH4) from forest soils in the Basque Country: Comparison of different tree species and growth stages

Forest systems are considered quintessential terrestrial systems for atmospheric CO2 sequestration to mitigate the effect of global warming. Temperate forest soils also present the highest rates of methane uptake among all the natural systems, while may represent a significant source of N2O. Despite of the large area occupied by forest in the Basque Country, no data is yet available regarding greenhouse gas fluxes under these edaphoclimatic conditions. In this manuscript we present a 2-year study which determined the magnitude of CO2, N2O, and CH4 soil gas fluxes in radiata pine, beech and Douglas fir forests using a closed chamber technique. The magnitude of these gas fluxes was additionally compared between different growth stages of radiata pine and beech forest, and the edaphoclimatic parameters that control these gas fluxes in the different forest systems and growth stages were studied. Measured greenhouse gas fluxes were in a low range as already published elsewhere for temperate forest ecosystems. A nitrogen deficit appears to be responsible for these relatively low gas fluxes. Apparently, the different forest species play a key role as controllers responsible for the differences of soil gas-exchange fluxes between the three different forest type systems. The mature pine plantation soil was emitting the most CO2 (1.5 and 2.5 times more than the mature beech and the Douglas fir, respectively), while the Douglas fir forest soil was emitting the most N2O (3 and 17 times more than the mature pine and the mature beech, respectively) and the mature beech forest was the soil type showing the highest CH4 consumption rates (2 and 5.5 times more than the mature pine and the Douglas fir, respectively). The stage of growth and its usual management appear to be important concerning the soil gas-exchange behavior within one forest type. The young beech forest soil emitted 9 times more N2O than the mature, and the new pine and the mature pine plantation soils emitted 2.5 and 2 times more CO2 than the young, respectively. The ground vegetation cover percentage, the organic matter accumulation and the soil porosity seem to be factors which merit a closer look in future studies, as possibly responsible for the differences in gas fluxes among forest types and growth stages.

Land use and management change under climate change adaptation and mitigation strategies: a U.S. case study

We examine the effects of crop management adaptation and climate mitigation strategies on land use and land management, plus on related environmental and economic outcomes. We find that crop management adaptation (e.g. crop mix, new species) increases Greenhouse gas (GHG) emissions by 1.7 % under a more severe climate projection while a carbon price reduces total forest and agriculture GHG annual flux by 15 % and 9 %, respectively. This shows that trade-offs are likely between mitigation and adaptation. Climate change coupled with crop management adaptation has small and mostly negative effects on welfare; mitigation, which is implemented as a carbon price starting at $15 per metric ton carbon dioxide (CO2) equivalent with a 5 % annual increase rate, bolsters welfare carbon payments. When both crop management adaptation and carbon price are implemented the effects of the latter dominates.

Ammonia emissions from deciduous forest after leaf fall

Abstract. The understanding of biochemical feedback mechanisms in the climate system is lacking knowledge in relation to bi-directional ammonia (NH3) exchange between natural ecosystems and the atmosphere. We therefore study the atmospheric NH3 fluxes during a 25-day period during autumn 2010 (21 October to 15 November) for the Danish beech forest Lille Bøgeskov to address the hypothesis that NH3 emissions occur from deciduous forests in relation to leaf fall. This is accomplished by using observations of vegetation status, NH3 fluxes and model calculations. Vegetation status was observed using plant area index (PAI) and leaf area index (LAI). NH3 fluxes were measured using the relaxed eddy accumulation (REA) method. The REA-based NH3 concentrations were compared to NH3 denuder measurements. Model calculations of the atmospheric NH3 concentration were obtained with the Danish Ammonia MOdelling System (DAMOS). The relative contribution from the forest components to the atmospheric NH3 flux was assessed using a simple two-layer bi-directional canopy compensation point model. A total of 57.7% of the fluxes measured showed emission and 19.5% showed deposition. A clear tendency of the flux going from deposition of −0.25 ± 0.30 μg NH3-N m−2 s−1 to emission of up to 0.67 ± 0.28 μg NH3-N m−2 s−1 throughout the measurement period was found. In the leaf fall period (23 October to 8 November), an increase in the atmospheric NH3 concentrations was related to the increasing forest NH3 flux. Following leaf fall, the magnitude and temporal structure of the measured NH3 emission fluxes could be adequately reproduced with the bi-directional resistance model; it suggested the forest ground layer (soil and litter) to be the main contributing component to the NH3 emissions. The modelled concentration from DAMOS fits well the measured concentrations before leaf fall, but during and after leaf fall, the modelled concentrations are too low. The results indicate that the missing contribution to atmospheric NH3 concentration from vegetative surfaces related to leaf fall are of a relatively large magnitude. We therefore conclude that emissions from deciduous forests are important to include in model calculations of atmospheric NH3 for forest ecosystems. Finally, diurnal variations in the measured NH3 concentrations were related to meteorological conditions, forest phenology and the spatial distribution of local anthropogenic NH3 sources. This suggests that an accurate description of ammonia fluxes over forest ecosystems requires a dynamic description of atmospheric and vegetation processes.

Foliar respiration and its temperature sensitivity in trees and lianas: in situ measurements in the upper canopy of a tropical forest

Leaf dark respiration (R) and its temperature sensitivity are essential for efforts to model carbon fluxes in tropical forests under current and future temperature regimes, but insufficient data exist to generalize patterns of R in species-rich tropical forests. Here, we tested the hypothesis that R and its temperature sensitivity (expressed as Q10, the proportional increase in R with a 10 °C rise in temperature) vary in relation to leaf functional traits, and among plant functional types (PFTs). We conducted in situ measurements of R of 461 leaves of 26 species of tree and liana in the upper canopy of a tropical forest in Panama. A construction crane allowed repeated non-destructive access to measure leaves kept in the dark since the previous night and equilibrated to the ambient temperature of 23-31 °C in the morning. R at 25 °C (R25) varied among species (mean 1.11 μmol m(-2) s(-1); range 0.72-1.79 μmol m(-2) s(-1)) but did not differ significantly among PFTs. R25 correlated positively with photosynthetic capacity, leaf mass per unit area, concentrations of nitrogen and phosphorus, and negatively with leaf lifespan. Q10 estimated for each species was on average higher than the 2.0 often assumed in coupled climate-vegetation models (mean 2.19; range 1.24-3.66). Early-successional tree species had higher Q10 values than other functional types, but interspecific variation in Q10 values was not correlated with other leaf traits. Similarity in respiration characteristics across PFTs, and relatively strong correlations of R with other leaf functional traits offer potential for trait-based vegetation modeling in species-rich tropical forests.

Modeling light use efficiency in a subtropical mangrove forest equipped with CO<sub>2</sub> eddy covariance

Abstract. Despite the importance of mangrove ecosystems in the global carbon budget, the relationships between environmental drivers and carbon dynamics in these forests remain poorly understood. This limited understanding is partly a result of the challenges associated with in situ flux studies. Tower-based CO2 eddy covariance (EC) systems are installed in only a few mangrove forests worldwide, and the longest EC record from the Florida Everglades contains less than 9 years of observations. A primary goal of the present study was to develop a methodology to estimate canopy-scale photosynthetic light use efficiency in this forest. These tower-based observations represent a basis for associating CO2 fluxes with canopy light use properties, and thus provide the means for utilizing satellite-based reflectance data for larger scale investigations. We present a model for mangrove canopy light use efficiency utilizing the enhanced green vegetation index (EVI) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) that is capable of predicting changes in mangrove forest CO2 fluxes caused by a hurricane disturbance and changes in regional environmental conditions, including temperature and salinity. Model parameters are solved for in a Bayesian framework. The model structure requires estimates of ecosystem respiration (RE), and we present the first ever tower-based estimates of mangrove forest RE derived from nighttime CO2 fluxes. Our investigation is also the first to show the effects of salinity on mangrove forest CO2 uptake, which declines 5% per each 10 parts per thousand (ppt) increase in salinity. Light use efficiency in this forest declines with increasing daily photosynthetic active radiation, which is an important departure from the assumption of constant light use efficiency typically applied in satellite-driven models. The model developed here provides a framework for estimating CO2 uptake by these forests from reflectance data and information about environmental conditions.

Fine root dynamics along an elevational gradient in tropical Amazonian and Andean forests

AbstractThe key role of tropical forest belowground carbon stocks and fluxes is well recognised as one of the main components of the terrestrial ecosystem carbon cycle. This study presents the first detailed investigation of spatial and temporal patterns of fine root stocks and fluxes in tropical forests along an elevational gradient, ranging from the Peruvian Andes (3020 m) to lowland Amazonia (194 m), with mean annual temperatures of 11.8°C to 26.4 °C and annual rainfall values of 1900 to 1560 mm yr‐1, respectively. Specifically, we analyse abiotic parameters controlling fine root dynamics, fine root growth characteristics, and seasonality of net primary productivity along the elevation gradient. Root and soil carbon stocks were measured by means of soil cores, and fine root productivity was recorded using rhizotron chambers and ingrowth cores. We find that mean annual fine root below ground net primary productivity in the montane forests (0–30 cm depth) ranged between 4.27±0.56 Mg C ha‐1 yr‐1 (1855 m) and 1.72±0.87 Mg C ha‐1 yr‐1 (3020 m). These values include a correction for finest roots (<0.6 mm diameter), which we suspect are under sampled, resulting in an underestimation of fine roots by up to 31% in current ingrowth core counting methods. We investigate the spatial and seasonal variation of fine root dynamics using soil depth profiles and an analysis of seasonal amplitude along the elevation gradient. We report a stronger seasonality of NPPFineRoot within the cloud immersion zone, most likely synchronised to seasonality of solar radiation. Finally, we provide the first insights into root growth characteristics along a tropical elevation transect: fine root area and fine root length increase significantly in the montane cloud forest. These insights into belowground carbon dynamics of tropical lowland and montane forests have significant implications for our understanding of the global tropical forest carbon cycle.

A toolbox of permutation tests for structural change

The sup $$LM$$ test for structural change is embedded into a permutation test framework for a simple location model. The resulting conditional permutation distribution is compared to the usual (unconditional) asymptotic distribution, showing that the power of the test can be clearly improved in small samples. Furthermore, the permutation test is embedded into a general framework that encompasses tools for binary and multivariate dependent variables as well as model-based permutation testing for structural change. It is also demonstrated that the methods can not only be employed for analyzing structural changes in time series data but also for recursive partitioning of cross-section data. The procedures suggested are illustrated using both artificial data and empirical applications (number of youth homicides, employment discrimination data, carbon flux in tropical forests, stock returns, and demand for economics journals).

Mineral soil carbon fluxes in forests and implications for carbon balance assessments

AbstractForest carbon cycles play an important role in efforts to understand and mitigate climate change. Large amounts of carbon (C) are stored in deep mineral forest soils, but are often not considered in accounting for global C fluxes because mineral soil C is commonly thought to be relatively stable. We explore C fluxes associated with forest management practices by examining existing data on forest C fluxes in the northeastern US. Our findings demonstrate that mineral soil C can play an important role in C emissions, especially when considering intensive forest management practices. Such practices are known to cause a high aboveground C flux to the atmosphere, but there is evidence that they can also promote comparably high and long‐term belowground C fluxes. If these additional fluxes are widespread in forests, recommendations for increased reliance on forest biomass may need to be reevaluated. Furthermore, existing protocols for the monitoring of forest C often ignore mineral soil C due to lack of data. Forest C analyses will be incomplete until this problem is resolved.

Predicting deciduous forest carbon uptake phenology by upscaling FLUXNET measurements using remote sensing data

Terrestrial ecosystems are highly sensitive to climatic changes in early and late growing seasons. Land surface phenology (LSP), the study of the timing of recurring seasonal pattern of variation in vegetated land surfaces observed from synoptic sensors, has thus received much attention due to its role as a surrogate in detecting the impact of climate change. Although several studies have been conducted on the growing season LSP, studies on the net carbon uptake phenology (CUP) defined as the detrended zero-crossing timing of net ecosystem productivity from a source to a sink in spring and vice versa in autumn, have been scarce. Here we present a CUP determination approach using the commonly available remote sensing data in four selected temperate and boreal deciduous forest CO2 flux tower sites. We test a hypothesis that the mean monthly surface temperature and LSP derived from remote sensing observations explain the CUP both in spring and autumn seasons. Our approach predicts the observed CUP in spring and autumn within 8 day mean errors, equivalent to the temporal resolution of the 8-day composite remote sensing dataset used in this study for the four flux tower sites. The results from this study will have a large implication for global change studies with increasing amount of valuable remote sensing data to be used for monitoring CUP beyond the footprints of CO2 flux towers.

Variation in mass and nutrient concentration of leaf litter across years and sites in a northern hardwood forest

Leaf litterfall represents an important nutrient flux in forests, but separating leaves by species and collecting fresh litter annually for nutrient analysis is time-consuming and expensive. To quantify the sources of variation in litterfall nutrient estimates and guide optimal allocation of research effort, we analyzed nutrient concentration (5 years) and mass (6 years) of leaf litter for nine tree species in 13 northern hardwood sites. Coefficients of variation (CVs) in nutrient concentration were higher across sites than over time within sites for most elements; phosphorus was especially variable across sites (56% CV). Thus, to estimate litterfall nutrient fluxes accurately in forests of this type, nutrient analyses should be site-specific as well as species-specific but may not need to be repeated annually (CVs over time averaged 17% for calcium, 21% for magnesium, 28% for potassium, and 32% for phosphorus concentration). Total leaf litterfall mass varied considerably from year to year, ranging from 234 to 370 g·m–2 averaged over 13 sites. We recommend that litter collectors be elevated above the ground to avoid oversampling during extreme wind events. Use of species-specific allometric equations, or even basal area, to estimate the species composition of total litter mass may obviate the need to sort litter by species.

Sustaining forest landscape connectivity under different land cover change scenarios

Managing forest landscapes to sustain functional connectivity is considered one of the key strategies to counteract the negative effects of climate and human-induced changes in forest species pools. With this objective, we evaluated whether a robust network of forest connecting elements can be identified so that it remains efficient when facing different types of potential land cover changes that may affect forest habitat networks and ecological fluxes. For this purpose we considered changes both in the forested areas and in the non-forest intervening landscape matrix. We combined some of the most recent developments in graph theory with models of land cover permeability and least-cost analysis through the forest landscape. We focused on a case of study covering the habitat of a forestdwelling bird (nuthatch, Sitta europaea) in the region of Galicia (NW Spain). Seven land-use change scenarios were analysed for their effects on connecting forest elements (patches and links): one was the simplest case in which the landscape is represented as a binary forest/non-forest pattern (and where matrix heterogeneity is disregarded), four scenarios in which forest lands were converted to other cover types (to scrubland due to wildfires, to extensive and intensive agriculture, and to urban areas), and two scenarios that only involved changes in the non-forested matrix (renaturalization and intensification). Our results show that while the network of connecting elements for the species was very robust to the conversion of the forest habitat patches to different cover types, the different change scenarios in the landscape matrix could more significantly weaken its long-term validity and effectiveness. This is particularly the case when most of the key connectivity providers for the nuthatch are located outside the protected areas or public forests in Galicia, where biodiversity-friendly measures might be more easily implemented. We discuss how the methodology can be applied to a wide range of forest landscape management situations, where both the conservation of the forest critical areas and an adequate management of the landscape matrix between them are of concern to achieve the sustainability of the ecological flows and ecosystem services at the wider forest landscape scale.

The intergalactic medium thermal history at redshiftz= 1.7-3.2 from the Lyα forest: a comparison of measurements using wavelets and the flux distribution

We investigate the thermal history of the intergalactic medium (IGM) in the redshift interval z=1.7--3.2 by studying the small-scale fluctuations in the Lyman alpha forest transmitted flux. We apply a wavelet filtering technique to eighteen high resolution quasar spectra obtained with the Ultraviolet and Visual Echelle Spectrograph (UVES), and compare these data to synthetic spectra drawn from a suite of hydrodynamical simulations in which the IGM thermal state and cosmological parameters are varied. From the wavelet analysis we obtain estimates of the IGM thermal state that are in good agreement with other recent, independent wavelet-based measurements. We also perform a reanalysis of the same data set using the Lyman alpha forest flux probability distribution function (PDF), which has previously been used to measure the IGM temperature-density relation. This provides an important consistency test for measurements of the IGM thermal state, as it enables a direct comparison of the constraints obtained using these two different methodologies. We find the constraints obtained from wavelets and the flux PDF are formally consistent with each other, although in agreement with previous studies, the flux PDF constraints favour an isothermal or inverted IGM temperature-density relation. We also perform a joint analysis by combining our wavelet and flux PDF measurements, constraining the IGM thermal state at z=2.1 to have a temperature at mean density of T0/[10^3 K]=17.3 +/- 1.9 and a power-law temperature-density relation exponent gamma=1.1 +/- 0.1 (1 sigma). Our results are consistent with previous observations that indicate there may be additional sources of heating in the IGM at z<4.

Experimental litterfall manipulation drives large and rapid changes in soil carbon cycling in a wet tropical forest

Global changes such as variations in plant net primary production are likely to drive shifts in leaf litterfall inputs to forest soils, but the effects of such changes on soil carbon (C) cycling and storage remain largely unknown, especially in C-rich tropical forest ecosystems. We initiated a leaf litterfall manipulation experiment in a tropical rain forest in Costa Rica to test the sensitivity of surface soil C pools and fluxes to different litter inputs. After only 2years of treatment, doubling litterfall inputs increased surface soil C concentrations by 31%, removing litter from the forest floor drove a 26% reduction over the same time period, and these changes in soil C concentrations were associated with variations in dissolved organic matter fluxes, fine root biomass, microbial biomass, soil moisture, and nutrient fluxes. However, the litter manipulations had only small effects on soil organic C (SOC) chemistry, suggesting that changes in C cycling, nutrient cycling, and microbial processes in response to litter manipulation reflect shifts in the quantity rather than quality of SOC. The manipulation also affected soil CO 2 fluxes; the relative decline in CO 2 production was greater in the litter removal plots (-22%) than the increase in the litter addition plots (+15%). Our analysis showed that variations in CO 2 fluxes were strongly correlated with microbial biomass pools, soil C and nitrogen (N) pools, soil inorganic P fluxes, dissolved organic C fluxes, and fine root biomass. Together, our data suggest that shifts in leaf litter inputs in response to localized human disturbances and global environmental change could have rapid and important consequences for belowground C storage and fluxes in tropical rain forests, and highlight differences between tropical and temperate ecosystems, where belowground C cycling responses to changes in litterfall are generally slower and more subtle.

SYSTEMATIC CONTINUUM ERRORS IN THE Lyα FOREST AND THE MEASURED TEMPERATURE–DENSITY RELATION

Continuum fitting uncertainties are a major source of error in estimates of the temperature-density relation (usually parametrized as a power-law, $T \propto \Delta^{\gamma - 1} $) of the inter-galactic medium (IGM) through the flux probability distribution function (PDF) of the Lyman-$\alpha$ forest. Using a simple order-of-magnitude calculation, we show that few percent-level systematic errors in the placement of the quasar continuum due to e.g. a uniform low-absorption Gunn-Peterson component, could lead to errors in $\gamma$ of order unity. This is quantified further using a simple semi-analytic model of the Lya forest flux PDF. We find that under-(over-)estimates in the continuum level can lead to a lower (higher) measured value of $\gamma$. Within current observational uncertainties, continuum biases double the error in $\gamma$ from $\sigma_{\gamma} \approx 0.1$ to $\sigma_{\gamma} \approx 0.2$ within our model. We argue that steps need to be taken to directly estimate the level of continuum bias in order to make recent claims of an inverted \tdr\ more robust.

The Lyman α forest flux probability distribution at z&gt;3★

We present a measurement of the Lyman α flux probability distribution function (PDF) obtained from a set of eight high-resolution quasar spectra with emission redshifts in the range 3.3 ≤z≤ 3.8. We carefully study the effect of metal absorption lines on the shape of the PDF. Metals have a larger impact on the PDF measurements at lower redshift, where there are relatively fewer Lyman α absorption lines. This may be explained by an increase in the number of metal lines that are blended with Lyman α absorption lines towards higher redshift, but may also be due to the presence of fewer metals in the intergalactic medium (IGM) at earlier times. We also provide a new measurement of the redshift evolution of the effective optical depth, τeff, at 2.8 ≤z≤ 3.6, and find no evidence for a deviation from a power-law evolution in the log (τeff)–log (1 +z) plane. The flux PDF measurements are furthermore of interest for studies of the thermal state of the IGM at z≃ 3. By comparing the PDF to state-of-the-art cosmological hydrodynamical simulations, we place constraints on the temperature of the IGM and compare our results with previous measurements of the PDF at lower redshift. At redshift z= 3, our new PDF measurements are consistent with an isothermal temperature–density relation, T=T0Δγ− 1, with a temperature at the mean density of T0= 19 250 ± 4800 K and a slope γ= 0.90 ± 0.21 (1σ uncertainties). In comparison, joint constraints with existing lower redshift PDF measurements at z < 3 favour an inverted temperature–density relation with T0= 17 900 ± 3500 K and γ= 0.70 ± 0.12, in broad agreement with previous analyses.

Data assimilation of photosynthetic light-use efficiency using multi-angular satellite data: I. Model formulation

Forest photosynthetic exchange rates at landscape scales have proven difficult to either accurately measure or estimate. Recent developments (Hall et al., 2011, 2008; Hilker et al., 2011a, 2010a) permit us to infer photosynthetic forest light use efficiency (ε) using multi-angle measurements of photochemical reflectance index (PRI) from the CHRIS/PROBA satellite imaging spectrometer, thus completing a long sought-after capability to remotely sense the major inputs driving gross primary production GPP i.e., ε and absorbed photosynthetically active radiation (APAR). In this first of two companion papers we introduce the theoretical underpinnings of an innovative approach that utilizes our recent developments to produce remotely sensed and spatially explicit maps of ε and GPP from space, and a data assimilation approach to extend the spatially explicit maps to diurnal, daily and annual time scales. We quantify GPP using the traditional radiation-limited approach of Monteith (1972); however we apply it in an innovative way. [I] Using CHRIS/PROBA we quantify ε at each satellite overpass for a 625km2 area at 30m resolution. [II] We use a novel physiologically-based multivariate function of APAR, temperature and water vapor pressure deficit model (described herein) and use it to down-regulate ε at 30 minute intervals. [III] We use the CHRIS/PROBA images of spatial variation in ε, and NDVI to quantify APAR, hence produce snapshots of GPP. We use a data assimilation approach to extend ε and GPP to temporally continuous and spatially contiguous maps of vegetation carbon uptake.In the second part of this study (Hilker et al., 2011b) we demonstrate and validate our approach over eight different forest flux tower sites in North America.

Evaluating the usefulness of hemispherical photographs as a means to estimate photosynthetic photon flux density during a growing season in the understorey of Nothofagus pumilio forests

Old-growth Nothofagus pumilio forests in Chile are managed employing a shelterwood system. A wide range of canopy openings can be found in old-growth and managed forests. Plant survival and growth in the understorey are influenced by the light available. There are limitations (practical and economic) to monitoring the light in the understorey. The aim of this study was to assess the options to estimate the forest understorey photosynthetic photon flux density (PPFD) measured during the growing season (GS) using canopy openness (CO) estimated by means of hemispherical photographs (HP). PPFD was measured using 31 sensors (Li-190SA quantum sensor) over the course of three GSs (October to March). The sensors were installed in an old-growth stand and another subjected to a regeneration felling under a shelterwood system. One HP was taken above each sensor (during the final GS) and the CO estimated. A comparison of the three seasons revealed that the sum of the PPFD during the GSs did not differ significantly. The CO could be used to effectively predict the sum of the PPFD during a GS (R 2 = 0.959). These results demonstrate the usefulness of HPs as a means to estimate the sum of the PPFD during a GS.

Nitrogen availability links forest productivity, soil nitrous oxide and nitric oxide fluxes of a tropical montane forest in southern Ecuador

[1] Tropical forests are important sources of the greenhouse gas nitrous oxide (N2O) and of nitric oxide (NO), a precursor of ozone. In tropical montane forests nitrogen limitation is common which affects both soil N2O and NO fluxes and forest productivity. Here we present evidence that forest productivity and N-oxide (N2O + NO) fluxes are linked through N availability along elevation and topographic gradients in tropical montane forests. We measured N-oxide fluxes, several indices of N availability, and forest productivity along an elevation gradient from 1000 m to 3000 m and along topographic gradients. Organic layer thickness of the soils increased and N availability decreased with increasing elevation and along the topographic gradient from the lower slope position to the ridges. Annual N2O fluxes ranged from −0.53 μg(N)m−2h−1 to 14.54 μg(N)m−2h−1 while NO fluxes ranged from −0.02 μg(N)m−2h−1 to 1.13 μg(N)m−2h−1. Both N-oxide fluxes and forest productivity increased with increasing N availability and showed close positive correlations with indices of N availability (C/N ratio andδ 15N signature of litterfall). We interpret the close correlations of N-oxide fluxes with total litterfall and tree basal area increment as evidence that N availability links N-oxide fluxes and forest productivity. This opens the possibility to include forest productivity as co-variable in predictions of N-oxide fluxes in nitrogen limited tropical montane forests. Especially increment of tree basal area was a promising proxy to predict soil N-oxide fluxes in these N limited ecosystems, possibly because it better reflects long-term forest productivity than total litterfall.