Quantifying particle size and turbulent scale dependence of dust uplift in the Sahara using aircraft measurements

The first size-resolved airborne measurements of dust fluxes and the first dust flux measurements from the central Sahara are presented and compared with a parameterization by Kok (2011a). High-frequency measurements of dust size distribution were obtained from 0.16 to 300 µm diameter, and eddy covariance fluxes were derived. This is more than an order of magnitude larger size range than previous flux estimates. Links to surface emission are provided by analysis of particle drift velocities. Number flux is described by a −2 power law between 1 and 144 µm diameter, significantly larger than the 12 µm upper limit suggested by Kok (2011a). For small particles, the deviation from a power law varies with terrain type and the large size cutoff is correlated with atmospheric vertical turbulent kinetic energy, suggesting control by vertical transport rather than emission processes. The measured mass flux mode is in the range 30–100 µm. The turbulent scales important for dust flux are from 0.1 km to 1–10 km. The upper scale increases during the morning as boundary layer depth and eddy size increase. All locations where large dust fluxes were measured had large topographical variations. These features are often linked with highly erodible surface features, such as wadis or dunes. We also hypothesize that upslope flow and flow separation over such features enhance the dust flux by transporting large particles out of the saltation layer. The tendency to locate surface flux measurements in open, flat terrain means these favored dust sources have been neglected in previous studies.

Joint data assimilation of satellite reflectance and net ecosystem exchange data constrains ecosystem carbon fluxes at a high-elevation subalpine forest

We utilized an ecosystem process model (SIPNET, simplified photosynthesis and evapotranspiration model) to estimate carbon fluxes of gross primary productivity and total ecosystem respiration of a high-elevation coniferous forest. The data assimilation routine incorporated aggregated twice-daily measurements of the net ecosystem exchange of CO2 (NEE) and satellite-based reflectance measurements of the fraction of absorbed photosynthetically active radiation (fAPAR) on an eight-day timescale. From these data we conducted a data assimilation experiment with fifteen different combinations of available data using twice-daily NEE, aggregated annual NEE, eight-day f AP AR, and average annual fAPAR. Model parameters were conditioned on three years of NEE and fAPAR data and results were evaluated to determine the information content from the different combinations of data streams. Across the data assimilation experiments conducted, model selection metrics such as the Bayesian Information Criterion and Deviance Information Criterion obtained minimum values when assimilating average annual fAPAR and twice-daily NEE data. Application of wavelet coherence analyses showed higher correlations between measured and modeled fAPAR on longer timescales ranging from 9 to 12 months. There were strong correlations between measured and modeled NEE (R2, coefficient of determination, 0.86), but correlations between measured and modeled eight-day fAPAR were quite poor (R2 = −0.94). We conclude that this inability to determine fAPAR on eight-day timescale would improve with the considerations of the radiative transfer through the plant canopy. Modeled fluxes when assimilating average annual fAPAR and annual NEE were comparable to corresponding results when assimilating twice-daily NEE, albeit at a greater uncertainty. Our results support the conclusion that for this coniferous forest twice-daily NEE data are a critical measurement stream for the data assimilation. The results from this modeling exercise indicate that for this coniferous forest, average annuals for satellite-based fAPAR measurements paired with annual NEE estimates may provide spatial detail to components of ecosystem carbon fluxes in proximity of eddy covariance towers. Inclusion of other independent data streams in the assimilation will also reduce uncertainty on modeled values.

CO2 emission estimation in the urban environment: measurement of the CO2 storage term

Eddy covariance has been used in urban areas to evaluate the net exchange of CO2 between the surface and the atmosphere. Typically, only the vertical flux is measured at a height 2–3 times that of the local roughness elements; however, under conditions of relatively low instability, CO2 may accumulate in the airspace below the measurement height. This can result in inaccurate emissions estimates if the accumulated CO2 drains away or is flushed upwards during thermal expansion of the boundary layer. Some studies apply a single height storage correction; however, this requires the assumption that the response of the CO2 concentration profile to forcing is constant with height. Here a full seasonal cycle (7th June 2012 to 3rd June 2013) of single height CO2 storage data calculated from concentrations measured at 10 Hz by open path gas analyser are compared to a data set calculated from a concurrent switched vertical profile measured (2 Hz, closed path gas analyser) at 10 heights within and above a street canyon in central London. The assumption required for the former storage determination is shown to be invalid. For approximately regular street canyons at least one other measurement is required. Continuous measurements at fewer locations are shown to be preferable to a spatially dense, switched profile, as temporal interpolation is ineffective. The majority of the spectral energy of the CO2 storage time series was found to be between 0.001 and 0.2 Hz (500 and 5 s respectively); however, sampling frequencies of 2 Hz and below still result in significantly lower CO2 storage values. An empirical method of correcting CO2 storage values from under-sampled time series is proposed.

Seasonal surface urban energy balance and wintertime stability simulated using three land-surface models in the high-latitude city Helsinki

The performance of three urban land surface models, run in offline mode, with their default external parameters, is evaluated for two distinctly different sites in Helsinki: Torni and Kumpula. The former is a dense city centre site with 22% vegetation, while the latter is a suburban site with over 50% vegetation. At both locations the models are compared against sensible and latent heat fluxes measured using the eddy covariance technique, along with snow depth observations. The cold climate experienced by the city causes strong seasonal variations that include snow cover and stable atmospheric conditions. Most of the time the three models are able to account for the differences between the study areas as well as the seasonal and diurnal variability of the energy balance components. However, the performances are not systematic across the modelled components, season and surface type. The net all-wave radiation is well simulated, with the greatest uncertainties related to snowmelt timing, when the fraction of snow cover has a key role, particularly in determining the surface albedo. For the turbulent fluxes, more variation between the models is seen which can partly be explained by the different methods in their calculation and partly by surface parameter values. For the sensible heat flux, simulation of wintertime values was the main problem, which also leads to issues in predicting near-surface stabilities particularly at the dense city centre site. All models have the most difficulties in simulating latent heat flux. This study particularly emphasizes that improvements are needed in the parameterization of anthropogenic heat flux and thermal parameters in winter, snow cover in spring and evapotranspiration in order to improve the surface energy balance modelling in cold climate cities.

Surface Urban Energy and Water Balance Scheme (SUEWS): development and evaluation at two UK sites

The Surface Urban Energy and Water Balance Scheme (SUEWS) is evaluated at two locations in the UK: a dense urban site in the centre of London and a residential suburban site in Swindon. Eddy covariance observations of the turbulent fluxes are used to assess model performance over a twoyear period (2011-2013). The distinct characteristics of the sites mean their surface energy exchanges differ considerably. The model suggests the largest differences can be attributed to surface cover (notably the proportion of vegetated versus impervious area) and the additional energy supplied by human activities. SUEWS performs better in summer than winter, and better at the suburban site than the dense urban site. One reason for this is the bias towards suburban summer field campaigns in observational data used to parameterise this (and other) model(s). The suitability of model parameters (such as albedo, energy use and water use) for the UK sites is considered and, where appropriate, alternative values are suggested. An alternative parameterisation for the surface conductance is implemented, which permits greater soil moisture deficits before evaporation is restricted at non-irrigated sites. Accounting for seasonal variation in the estimation of storage heat flux is necessary to obtain realistic wintertime fluxes.

Fluxes of CO(2) above a plantation of Eucalyptus in southeast Brazil

Eddy-covariance measurements of net ecosystem exchange of CO(2) (NEE) and estimates of gross ecosystem productivity (GEP) and ecosystem respiration (R(E)) were obtained in a 2-4 year old Eucalyptus plantation during two years with very different winter rainfall In the first (drier) year the annual NEE GEP and RE were lower than the sums in the second (normal) year and conversely the total respiratory costs of assimilated carbon were higher in the dry year than in the normal year Although the net primary production (NPP) in the first year was 23% lower than that of the second year the decrease in the carbon use efficiency (CUE = NPP/GEP) was 11% and autotrophic respiration utilized more resources in the first dry year than in the second normal year The time variations in NEE were followed by NPP because in these young Eucalyptus plantations NEE is very largely dominated by NPP and heterotrophic respiration plays only a relatively minor role During the dry season a pronounced hysteresis was observed in the relationship between NEE and photosynthetically active radiation and NEE fluxes were inversely proportional to humidity saturation deficit values greater than 0 8 kPa Nighttime fluxes of CO(2) during calm conditions when the friction velocity (u) was below the threshold (0 25 ms(-1)) were estimated based on a Q(10) temperature-dependence relationship adjusted separately for different classes of soil moisture content which regulated the temperature sensitivity of ecosystem respiration (C) 2010 Elsevier B V All rights reserved

The energy and water balance of a Eucalyptus plantation in southeast Brazil

The eddy covariance method was used to measure energy and water balance of a plantation of Eucalyptus (grandis x urophylla) hybrids over a 2 year period. The average daily evaporation rates were 5.4 (+/- 2.0) mm day(-1) in summer, but fell to 1.2 (+/- 0.3) mm day(-1) in winter. In contrast, the sensible heat flux was relatively low in summer but dominated the energy balance in winter. Evaporation accounted for 80% and 26% of the available energy, in summer and winter respectively. The annual evaporation was 82% (1124 mm) and 96% (1235 mm) of the annual rainfall recorded during the first and second year, respectively. Daily average canopy and aerodynamic conductance to water vapour were in the summer 51.9 (+/- 38.4) mm s(-1) 84.1 (+/- 25.6) mm s(-1), respectively; and in the winter 6.0 (+/- 10.5) mm s(-1) and 111.6 (+/- 24.6) mm s(-1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

The land-atmosphere water flux in the tropics

Tropical vegetation is a major source of global land surface evapotranspiration, and can thus play a major role in global hydrological cycles and global atmospheric circulation. Accurate prediction of tropical evapotranspiration is critical to our understanding of these processes under changing climate. We examined the controls on evapotranspiration in tropical vegetation at 21 pan-tropical eddy covariance sites, conducted a comprehensive and systematic evaluation of 13 evapotranspiration models at these sites, and assessed the ability to scale up model estimates of evapotranspiration for the test region of Amazonia. Net radiation was the strongest determinant of evapotranspiration (mean evaporative fraction was 0.72) and explained 87% of the variance in monthly evapotranspiration across the sites. Vapor pressure deficit was the strongest residual predictor (14%), followed by normalized difference vegetation index (9%), precipitation (6%) and wind speed (4%). The radiation-based evapotranspiration models performed best overall for three reasons: (1) the vegetation was largely decoupled from atmospheric turbulent transfer (calculated from X decoupling factor), especially at the wetter sites; (2) the resistance-based models were hindered by difficulty in consistently characterizing canopy (and stomatal) resistance in the highly diverse vegetation; (3) the temperature-based models inadequately captured the variability in tropical evapotranspiration. We evaluated the potential to predict regional evapotranspiration for one test region: Amazonia. We estimated an Amazonia-wide evapotranspiration of 1370 mm yr(-1), but this value is dependent on assumptions about energy balance closure for the tropical eddy covariance sites; a lower value (1096 mm yr(-1)) is considered in discussion on the use of flux data to validate and interpolate models.

Atmosphere and hydrological controls of the evapotranspiration over a floodplain forest in the Bananal Island region, Amazonia

This article discusses seasonal and interannual variations of the evapotranspiration (ET) rates in Bananal Island floodplain, Brazil. Measurements included ET and sensible heat flux using the eddy covariance method, atmospheric forcings (net radiation, Rn, vapor pressure deficit, VPD, wind speed and air temperature), soil moisture profiles, groundwater level and flood height, taken from November 2003 to December 2006. For the hydrological years (October-September) of 2003/2004, 2004/2005 and 2005/2006, the accumulated precipitation was 1692, 1471, 1914 mm and the accumulated ET was 1361, 1318 and 1317 mm, respectively. Seasonal analyses indicated that ET decreased in the dry season (average 3.7 mm day(-1)), despite the simultaneous increase in Rn, air temperature and VPD. The increase of ET in the wet season and particularly in the flood period (average 4.1 mm day(-1)) showed that the free water surface evaporation strongly influenced the energy exchange. Soil moisture, which was substantially depleted during the dry season, and adaptative vegetation mechanisms such as leaf senescence contributed to limit the dry season ET. Strong drainage within permeable sandy soils helped to explain the soil moisture depletion. These results suggest that the Bananal flooding area shows a different pattern in relation to the upland Amazon forests, being more similar to the savanna strictu senso areas in central Brazil. For example, seasonal ET variation was not in phase with Rn; the wet season ET was higher than the dry season ET; and the system stored only a tiny memory of the flooding period, being sensitive to extended drought periods.

Análise e estimativa dos componentes do balanço de energia em ecossistema de manguezal amazônico

As medições e estimativas dos componentes do balanço de energia foram feitos acima da copa das árvores no ecossistema de manguezal natural, localizada a 30 km da cidade de Bragança-PA, entre novembro de 2002 e agosto de 2003. Os dados foram utilizados para a análise das variações sazonais e horárias do fluxo de calor sensível e calor latente, bem como a avaliação da partição de energia. Os dados meteorológicos foram coletados pela estação meteorológica automática (EMA) e os fluxos foram calculados utilizando-se a técnica de covariância de vórtices turbulentos. Os modelos de Penman-Monteith e Shuttleworth foram usados para estimar o fluxo de calor sensível e calor latente. O objetivo deste estudo foi analisar o equilíbrio e a partição de energia no manguezal, assim como fazer uma avaliação do comportamento de modelos empíricos para estimar os fluxos de energia. O saldo de radiação apresentou valores mais elevados no período menos chuvoso. A razão de Bowen mostrou valor geralmente baixo, o que indica que uma proporção maior de energia foi utilizada sob a forma de calor latente. O modelo Shuttleworth é mais eficiente na estimativa de fluxos de calor sensível. Para estimar o fluxo de calor latente do modelo de Penman-Monteith foi mais eficiente durante a estação seca e o modelo Shuttleworth durante a estação chuvosa.

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