What drives the seasonality of photosynthesis across the Amazon basin? A cross-site analysis of eddy flux tower measurements from the Brasil flux network

We investigated the seasonal patterns of Amazonian forest photosynthetic activity, and the effects thereon of variations in climate and land-use, by integrating data from a network of ground-based eddy flux towers in Brazil established as part of the ‘Large-Scale Biosphere Atmosphere Experiment in Amazonia’ project. We found that degree of water limitation, as indicated by the seasonality of the ratio of sensible to latent heat flux (Bowen ratio) predicts seasonal patterns of photosynthesis. In equatorial Amazonian forests (5◦ N–5◦ S), water limitation is absent, and photosynthetic fluxes (or gross ecosystem productivity, GEP) exhibit high or increasing levels of photosynthetic activity as the dry season progresses, likely a consequence of allocation to growth of new leaves. In contrast, forests along the southern flank of the Amazon, pastures converted from forest, and mixed forest-grass savanna, exhibit dry-season declines in GEP, consistent with increasing degrees of water limitation. Although previous work showed tropical ecosystem evapotranspiration (ET) is driven by incoming radiation, GEP observations reported here surprisingly show no or negative relationships with photosynthetically active radiation (PAR). Instead, GEP fluxes largely followed the phenology of canopy photosynthetic capacity (Pc), with only deviations from this primary pattern driven by variations in PAR. Estimates of leaf flush at three

Amazon Rainforest Exchange of Carbon and Subcanopy Air Flow: Manaus LBA Site-A Complex Terrain Condition

On the moderately complex terrain covered by dense tropical Amazon Rainforest (Reserva Biologica do Cuieiras-ZF2-02 degrees 36'17.1 '' S, 60 degrees 12'24.4 '' W), subcanopy horizontal and vertical gradients of the air temperature, CO2 concentration and wind field were measured for the dry and wet periods in 2006. We tested the hypothesis that horizontal drainage flow over this study area is significant and can affect the interpretation of the high carbon uptake rates reported by previous works at this site. A similar experimental design as the one by Tota et al. (2008) was used with a network of wind, air temperature, and CO2 sensors above and below the forest canopy. A persistent and systematic subcanopy nighttime upslope (positive buoyancy) and daytime downslope (negative buoyancy) flow pattern on a moderately inclined slope (12%) was observed. The microcirculations observed above the canopy (38 m) over the sloping area during nighttime presents a downward motion indicating vertical convergence and correspondent horizontal divergence toward the valley area. During the daytime an inverse pattern was observed. The microcirculations above the canopy were driven mainly by buoyancy balancing the pressure gradient forces. In the subcanopy space the microcirculations were also driven by the same physical mechanisms but probably with the stress forcing contribution. The results also indicated that the horizontal and vertical scalar gradients (e. g., CO2) were modulated by these micro-circulations above and below the canopy, suggesting that estimates of advection using previous experimental approaches are not appropriate due to the tridimensional nature of the vertical and horizontal transport locally. This work also indicates that carbon budget from tower-based measurement is not enough to close the system, and one needs to include horizontal and vertical advection transport of CO2 into those estimates.

A MODIS-Based Energy Balance to Estimate Evapotranspiration for Clear-Sky Days in Brazilian Tropical Savannas

Evapotranspiration (ET) plays an important role in global climate dynamics and in primary production of terrestrial ecosystems; it represents the mass and energy transfer from the land to atmosphere. Limitations to measuring ET at large scales using ground-based methods have motivated the development of satellite remote sensing techniques. The purpose of this work is to evaluate the accuracy of the SEBAL algorithm for estimating surface turbulent heat fluxes at regional scale, using 28 images from MODIS. SEBAL estimates are compared with eddy-covariance (EC) measurements and results from the hydrological model MGB-IPH. SEBAL instantaneous estimates of latent heat flux (LE) yielded r(2) = 0.64 and r(2) = 0.62 over sugarcane croplands and savannas when compared against in situ EC estimates. At the same sites, daily aggregated estimates of LE were r(2) = 0.76 and r(2) = 0.66, respectively. Energy balance closure showed that turbulent fluxes over sugarcane croplands were underestimated by 7% and 9% over savannas. Average daily ET from SEBAL is in close agreement with estimates from the hydrological model for an overlay of 38,100 km(2) (r(2) = 0.88). Inputs to which the algorithm is most sensitive are vegetation index (NDVI), gradient of temperature (dT) to compute sensible heat flux (H) and net radiation (Re). It was verified that SEBAL has a tendency to overestimate results both at local and regional scales probably because of low sensitivity to soil moisture and water stress. Nevertheless the results confirm the potential of the SEBAL algorithm, when used with MODIS images for estimating instantaneous LE and daily ET from large areas.

Water use in a sugarcane plantation

The evapotranspiration (E) from a sugarcane plantation in the southeast Brazil was measured by the eddy-covariance method during two consecutive cycles. These represented the second (393 similar to days) and third year (374 similar to days) re-growth (ratoon). The total E in the first cycle was 829 similar to mm, accounting for 69% of rainfall, whereas in the second cycle, it was 690 similar to mm, despite the total rainfall (1353 similar to mm) being 13% greater. The ratio of E to available energy, the evaporative fraction, exhibited a smaller variation between the first and second cycles: 0.58 and 0.51, respectively. The estimated interception losses were 88 and 90 similar to mm, respectively, accounting for approximately 7% of the total rainfall. The sugarcane yield in the second cycle (61.5 similar to +/-similar to 4.0 similar to t similar to ha-1) was 26% lower than the first cycle, as well as lower than the regional average for the third ratoon (76 similar to t similar to ha-1). The below average yield was associated with less available soil water at the beginning of the cycle, with the amount of rainfall recorded during the first 120 similar to days of re-growth in the second cycle being 16% of that recorded in the first (203 similar to mm).

Validation and comparison of two soil-vegetation-atmosphere transfer models for tropical Africa

This study aims to compare and validate two soil-vegetation-atmosphere-transfer (SVAT) schemes: TERRA-ML and the Community Land Model (CLM). Both SVAT schemes are run in standalone mode (decoupled from an atmospheric model) and forced with meteorological in-situ measurements obtained at several tropical African sites. Model performance is quantified by comparing simulated sensible and latent heat fluxes with eddy-covariance measurements. Our analysis indicates that the Community Land Model corresponds more closely to the micrometeorological observations, reflecting the advantages of the higher model complexity and physical realism. Deficiencies in TERRA-ML are addressed and its performance is improved: (1) adjusting input data (root depth) to region-specific values (tropical evergreen forest) resolves dry-season underestimation of evapotranspiration; (2) adjusting the leaf area index and albedo (depending on hard-coded model constants) resolves overestimations of both latent and sensible heat fluxes; and (3) an unrealistic flux partitioning caused by overestimated superficial water contents is reduced by adjusting the hydraulic conductivity parameterization. CLM is by default more versatile in its global application on different vegetation types and climates. On the other hand, with its lower degree of complexity, TERRA-ML is much less computationally demanding, which leads to faster calculation times in a coupled climate simulation.

Chemically Resolved Particle Fluxes Over Tropical and Temperate Forests

Chemically resolved submicron (PM1) particlemass fluxes were measured by eddy covariance with a high resolution time-of-flight aerosolmass spectrometer over temperate and tropical forests during the BEARPEX-07 and AMAZE-08 campaigns. Fluxes during AMAZE-08 were small and close to the detection limit (<1 ng m−2 s−1) due to low particle mass concentrations (<1 μg m−3). During BEARPEX-07, concentrations were five times larger, with mean mid-day deposition fluxes of −4.8 ng m−2 s−1 for total nonrefractory PM1 (Vex,PM1 = −1 mm s−1) and emission fluxes of +2.6 ng m−2 s−1 for organic PM1 (Vex,org = +1 mm s−1). Biosphere–atmosphere fluxes of different chemical components are affected by in-canopy chemistry, vertical gradients in gas-particle partitioning due to canopy temperature gradients, emission of primary biological aerosol particles, and wet and dry deposition. As a result of these competing processes, individual chemical components had fluxes of varying magnitude and direction during both campaigns. Oxygenated organic components representing regionally aged aerosol deposited, while components of fresh secondary organic aerosol (SOA) emitted. During BEARPEX-07, rapid incanopy oxidation caused rapid SOA growth on the timescale of biosphere-atmosphere exchange. In-canopy SOA mass yields were 0.5–4%. During AMAZE-08, the net organic aerosol flux was influenced by deposition, in-canopy SOA formation, and thermal shifts in gas-particle partitioning.Wet deposition was estimated to be an order ofmagnitude larger than dry deposition during AMAZE-08. Small shifts in organic aerosol concentrations from anthropogenic sources such as urban pollution or biomass burning alters the balance between flux terms. The semivolatile nature of the Amazonian organic aerosol suggests a feedback in which warmer temperatures will partition SOA to the gas-phase, reducing their light scattering and thus potential to cool the region.

Fluxes of CO2 above a sugarcane plantation in Brazil

Fluxes of CO2 were measured above a sugarcane plantation using the eddy-covariance method covering two growth cycles, representing the second and third re-growth (ratoons) harvested with stubble burning. The total net ecosystem exchange (NEE) in the first cycle (second ratoon, 393 days long) was −1964 ± 44 g C m−2; the gross ecosystem productivity (GEP) was 3612 ± 46 g C m−2 and the ecosystem respiration (RE) was 1648 ± 14 g C m−2. The NEE and GEP totals in the second cycle (third ratoon, 374 days long) decreased 51% and 25%, respectively and RE increased 7%. Accounting for the carbon emitted during biomass burning and the removal of stalks at harvest, net ecosystem carbon balance (NECB) totals were 102 ± 130 g C m−2 and 403 ± 84 g C m−2 in each cycle respectively. Thus the sugarcane agrosystem was approximately carbon neutral in the second ratoon. Yield in stalks fresh weight (SFW) attained the regional average (8.3 kg SFW m−2). Although it was a carbon source to the atmosphere, observed productivity (6.2 kg SFW m−2) of the third ratoon was 19% lower than the regional average due to the lower water availability observed during the initial 120 days of re-growth. However, the overall water use efficiency (WUE) achieved in the first cycle (4.3 g C kg−1 H2O) decreased only 5% in the second cycle. © 2013 Elsevier B.V. All rights reserved

Eddy-induced upwelling off Cape Sao Tome (22 degrees S, Brazil)

The regional ocean off southeast Brazil (20 degrees S-28 degrees S) is known as a current-eddy-upwelling region. The proximity of the Brazil Current to the coast in the Cape Sao Tome vicinities, as well as of its quasi-stationary unstable meanders, suggests the possibility of background eddy-induced upwelling. Such phenomenon can intensify the prevalent coastal upwelling due to wind and topographic effects. In this paper, with the help of a numerical simulation, we provide evidence that eddy-induced upwelling in the absence of wind is possible in this region. The simulation was conducted with a regional configuration of the 3-D Princeton Ocean Model initialized by a feature-based implementation of the Brazil Current and Cape Frio eddy, blended with climatology. (C) 2010 Elsevier Ltd. All rights reserved.

Is the Brazil Current eddy-dominated to the north of 20 degrees S?

The Brazil Current (BC) originates with the arrival and bifurcation of the southernmost branch of the South Equatorial Current (sSEC) between 10-20 degrees S. Previous climatological studies showed a stratified sSEC bifurcation and that the resulting southern branch formed a shallow BC - a weak western boundary current. The analysis of three recent synoptic surveys and global model outputs challenge the description of a continuous BC. The sSEC bifurcation signal near the continental margin was unclear in the analyses, and the velocity fields were dominated by mesoscale eddies. Recurrent anticyclones that seemed to be related to the meandering BC led us to construct a picture of a flow strongly influenced by topography and probably very unstable. Given this new emerging scenario, we hypothesize that the Brazil Current is eddy-dominated to the north of 20 degrees S. Citation: Soutelino, R. G., I. C. A. da Silveira, A. Gangopadhyay, and J. A. Miranda (2011), Is the Brazil Current eddy-dominated to the north of 20 S?, Geophys. Res. Lett., 38, L03607, doi:10.1029/2010GL046276.

A Single Statistic for Monitoring the Covariance Matrix of Bivariate Processes

In this article, we present a new control chart for monitoring the covariance matrix in a bivariate process. In this method, n observations of the two variables were considered as if they came from a single variable (as a sample of 2n observations), and a sample variance was calculated. This statistic was used to build a new control chart specifically as a VMIX chart. The performance of the new control chart was compared with its main competitors: the generalized sampled variance chart, the likelihood ratio test, Nagao's test, probability integral transformation (v(t)), and the recently proposed VMAX chart. Among these statistics, only the VMAX chart was competitive with the VMIX chart. For shifts in both variances, the VMIX chart outperformed VMAX; however, VMAX showed better performance for large shifts (higher than 10%) in one variance.