Nitrogen management is essential to prevent tropical oil palm plantations from causing ozone pollution

More than half the world's rainforest has been lost to agriculture since the Industrial Revolution. Among the most widespread tropical crops is oil palm (Elaeis guineensis): global production now exceeds 35 million tonnes per year. In Malaysia, for example, 13% of land area is now oil palm plantation, compared with 1% in 1974. There are enormous pressures to increase palm oil production for food, domestic products, and, especially, biofuels. Greater use of palm oil for biofuel production is predicated on the assumption that palm oil is an “environmentally friendly” fuel feedstock. Here we show, using measurements and models, that oil palm plantations in Malaysia directly emit more oxides of nitrogen and volatile organic compounds than rainforest. These compounds lead to the production of ground-level ozone (O3), an air pollutant that damages human health, plants, and materials, reduces crop productivity, and has effects on the Earth's climate. Our measurements show that, at present, O3 concentrations do not differ significantly over rainforest and adjacent oil palm plantation landscapes. However, our model calculations predict that if concentrations of oxides of nitrogen in Borneo are allowed to reach those currently seen over rural North America and Europe, ground-level O3 concentrations will reach 100 parts per billion (109) volume (ppbv) and exceed levels known to be harmful to human health. Our study provides an early warning of the urgent need to develop policies that manage nitrogen emissions if the detrimental effects of palm oil production on air quality and climate are to be avoided.

Nitrogen management is essential to prevent tropical oil palm plantations from causing ground-level ozone pollution

More than half the world's rainforest has been lost to agriculture since the Industrial Revolution. Among the most widespread tropical crops is oil palm (Elaeis guineensis): global production now exceeds 35 million tonnes per year. In Malaysia, for example, 13% of land area is now oil palm plantation, compared with 1% in 1974. There are enormous pressures to increase palm oil production for food, domestic products, and, especially, biofuels. Greater use of palm oil for biofuel production is predicated on the assumption that palm oil is an "environmentally friendly'' fuel feedstock. Here we show, using measurements and models, that oil palm plantations in Malaysia directly emit more oxides of nitrogen and volatile organic compounds than rainforest. These compounds lead to the production of ground-level ozone (O-3), an air pollutant that damages human health, plants, and materials, reduces crop productivity, and has effects on the Earth's climate. Our measurements show that, at present, O-3 concentrations do not differ significantly over rainforest and adjacent oil palm plantation landscapes. However, our model calculations predict that if concentrations of oxides of nitrogen in Borneo are allowed to reach those currently seen over rural North America and Europe, ground-level O-3 concentrations will reach 100 parts per billion (10(9)) volume (ppbv) and exceed levels known to be harmful to human health. Our study provides an early warning of the urgent need to develop policies that manage nitrogen emissions if the detrimental effects of palm oil production on air quality and climate are to be avoided.

Sensitivity of large-scale atmospheric analyses to humidity observations and its impact on the global water cycle and tropical and extratropical weather systems in ERA40

Reanalysis data obtained from data assimilation are increasingly used for diagnostic studies of the general circulation of the atmosphere, for the validation of modelling experiments and for estimating energy and water fluxes between the Earth surface and the atmosphere. Because fluxes are not specifically observed, but determined by the data assimilation system, they are not only influenced by the utilized observations but also by model physics and dynamics and by the assimilation method. In order to better understand the relative importance of humidity observations for the determination of the hydrological cycle, in this paper we describe an assimilation experiment using the ERA40 reanalysis system where all humidity data have been excluded from the observational data base. The surprising result is that the model, driven by the time evolution of wind, temperature and surface pressure, is able to almost completely reconstitute the large-scale hydrological cycle of the control assimilation without the use of any humidity data. In addition, analysis of the individual weather systems in the extratropics and tropics using an objective feature tracking analysis indicates that the humidity data have very little impact on these systems. We include a discussion of these results and possible consequences for the way moisture information is assimilated, as well as the potential consequences for the design of observing systems for climate monitoring. It is further suggested, with support from a simple assimilation study with another model, that model physics and dynamics play a decisive role for the hydrological cycle, stressing the need to better understand these aspects of model parametrization. .

Feature-point detection using distance transforms: Application to tracking tropical convective complexes

The identification, tracking, and statistical analysis of tropical convective complexes using satellite imagery is explored in the context of identifying feature points suitable for tracking. The feature points are determined based on the shape of complexes using the distance transform technique. This approach has been applied to the determination feature points for tropical convective complexes identified in a time series of global cloud imagery. The feature points are used to track the complexes, and from the tracks statistical diagnostic fields are computed. This approach allows the nature and distribution of organized deep convection in the Tropics to be explored.

A comparison of recent reanalysis datasets using objective feature tracking: Storm tracks and tropical easterly waves

Data from four recent reanalysis projects [ECMWF, NCEP-NCAR, NCEP - Department of Energy ( DOE), NASA] have been diagnosed at the scale of synoptic weather systems using an objective feature tracking method. The tracking statistics indicate that, overall, the reanalyses correspond very well in the Northern Hemisphere (NH) lower troposphere, although differences for the spatial distribution of mean intensities show that the ECMWF reanalysis is systematically stronger in the main storm track regions but weaker around major orographic features. A direct comparison of the track ensembles indicates a number of systems with a broad range of intensities that compare well among the reanalyses. In addition, a number of small-scale weak systems are found that have no correspondence among the reanalyses or that only correspond upon relaxing the matching criteria, indicating possible differences in location and/or temporal coherence. These are distributed throughout the storm tracks, particularly in the regions known for small-scale activity, such as secondary development regions and the Mediterranean. For the Southern Hemisphere (SH), agreement is found to be generally less consistent in the lower troposphere with significant differences in both track density and mean intensity. The systems that correspond between the various reanalyses are considerably reduced and those that do not match span a broad range of storm intensities. Relaxing the matching criteria indicates that there is a larger degree of uncertainty in both the location of systems and their intensities compared with the NH. At upper-tropospheric levels, significant differences in the level of activity occur between the ECMWF reanalysis and the other reanalyses in both the NH and SH winters. This occurs due to a lack of coherence in the apparent propagation of the systems in ERA15 and appears most acute above 500 hPa. This is probably due to the use of optimal interpolation data assimilation in ERA15. Also shown are results based on using the same techniques to diagnose the tropical easterly wave activity. Results indicate that the wave activity is sensitive not only to the resolution and assimilation methods used but also to the model formulation.

African easterly wave variability and its relationship to Atlantic tropical cyclone activity

Automatic tracking of vorticity centers in European Centre for Medium-Range Weather Forecasts analyses has been used to develop a 20-yr climatology of African easterly wave activity. The tracking statistics at 600 and 850 mb confirm the complicated easterly wave structures present over the African continent. The rainy zone equatorward of 15 degreesN is dominated by 600-mb activity, and the much drier Saharan region poleward of 15 degreesN is more dominated by 850-mb activity. Over the Atlantic Ocean there is just one storm track with the 600- and 850-mb wave activity collocated. Based on growth/decay and genesis statistics, it appears that the 850-mb waves poleward of 15 degreesN over land generally do not get involved with the equatorward storm track over the ocean. Instead, there appears to be significant development of 850-mb activity at the West African coast in the rainy zone around (10 degreesN, 10 degreesW), which, it is proposed, is associated with latent heat release. Based on the tracking statistics, it has been shown that there is marked interannual variability in African easterly wave (AEW) activity. It is especially marked at the 850-mb level at the West African coast between about 10 degrees and 15 degreesN, where the coefficient of variation is 0.29. For the period between 1985 and 1998, a notable positive correlation is seen between this AEW activity and Atlantic tropical cyclone activity. This correlation is particularly strong for the postreanalysis period between 1994 and 1998. This result suggests that Atlantic tropical cyclone activity may be influenced by the number of AEWs leaving the West African coast, which have significant low-level amplitudes, and not simply by the total number of AEWs.

West African storm tracks and their relationship to Atlantic tropical cyclones

The automatic tracking technique used by Thorncroft and Hodges (2001) has been used to identify coherent vorticity structures at 850hPa over West Africa and the tropical Atlantic in the ECMWF 40-year reanalysis. The presence of two dominant source regions, north and south of 15ºN over West Africa, for storm tracks over the Atlantic was confirmed. Results show that the southern storm track provides most of the storms that reach the main development region where most tropical cyclones develop. There exists marked seasonal variability in location and intensity of the storms leaving the West African coast, which may influence the likelihood of downstream intensification and longevity. There exists considerable year-to-year variability in the number of West African storm tracks, both in numbers over the land and continuing out over the tropical Atlantic Ocean. While the low-frequency variability is well correlated with Atlantic tropical cyclone activity, West African rainfall and SSTs, the interannual variability is found to be uncorrelated with these. In contrast, variance of the 2-6-day-filtered meridional wind, which provides a synoptic-scale measure of African Easterly Wave activity, shows a significant, positive correlation with tropical cyclone activity at interannual timescales.

Tropical cyclones in a T159 Resolution Global Climate Model: comparison with observations and re-analyses

Tropical cyclones have been investigated in a T159 version of the MPI ECHAM5 climate model using a novel technique to diagnose the evolution of the 3-dimensional vorticity structure of tropical cyclones, including their full life cycle from weak initial vortex to their possible extra-tropical transition. Results have been compared with reanalyses (ERA40 and JRA25) and observed tropical storms during the period 1978-1999 for the Northern Hemisphere. There is no indication of any trend in the number or intensity of tropical storms during this period in ECHAM5 or in re-analyses but there are distinct inter-annual variations. The storms simulated by ECHAM5 are realistic both in space and time, but the model and even more so the re-analyses, underestimate the intensities of the most intense storms (in terms of their maximum wind speeds). There is an indication of a response to ENSO with a smaller number of Atlantic storms during El Niño in agreement with previous studies. The global divergence circulation responds to El Niño by setting up a large-scale convergence flow, with the center over the central Pacific with enhanced subsidence over the tropical Atlantic. At the same time there is an increase in the vertical wind shear in the region of the tropical Atlantic where tropical storms normally develop. There is a good correspondence between the model and ERA40 except that the divergence circulation is somewhat stronger in the model. The model underestimates storms in the Atlantic but tends to overestimate them in the Western Pacific and in the North Indian Ocean. It is suggested that the overestimation of storms in the Pacific by the model is related to an overly strong response to the tropical Pacific SST anomalies. The overestimation in 2 the North Indian Ocean is likely to be due to an over prediction in the intensity of monsoon depressions, which are then classified as intense tropical storms. Nevertheless, overall results are encouraging and will further contribute to increased confidence in simulating intense tropical storms with high-resolution climate models.

How may tropical cyclones change in a warmer climate?

Tropical Cyclones (TC) under different climate conditions in the Northern Hemisphere have been investigated with the Max Planck Institute (MPI) coupled (ECHAM5/MPIOM) and atmosphere (ECHAM5) climate models. The intensity and size of the TC depend crucially on resolution with higher wind speed and smaller scales at the higher resolutions. The typical size of the TC is reduced by a factor of 2.3 from T63 to T319 using the distance of the maximum wind speed from the centre of the storm as a measure. The full three dimensional structure of the storms becomes increasingly more realistic as the resolution is increased. For the T63 resolution, three ensemble runs are explored for the period 1860 until 2100 using the IPCC SRES scenario A1B and evaluated for three 30 year periods at the end of the 19th, 20th and 21st century, respectively. While there is no significant change between the 19th and the 20th century, there is a considerable reduction in the number of the TC by some 20% in the 21st century, but no change in the number of the more intense storms. Reduction in the number of storms occurs in all regions. A single additional experiment at T213 resolution was run for the two latter 30-year periods. The T213 is an atmospheric only experiment using the transient Sea Surface Temperatures (SST) of the T63 resolution experiment. Also in this case, there is a reduction by some 10% in the number of simulated TC in the 21st century compared to the 20th century but a marked increase in the number of intense storms. The number of storms with maximum wind speeds greater than 50ms-1 increases by a third. Most of the intensification takes place in 2 the Eastern Pacific and in the Atlantic where also the number of storms more or less stays the same. We identify two competing processes effecting TC in a warmer climate. First, the increase in the static stability and the reduced vertical circulation is suggested to contribute to the reduction in the number of storms. Second, the increase in temperature and water vapor provide more energy for the storms so that when favorable conditions occur, the higher SST and higher specific humidity will contribute to more intense storms. As the maximum intensity depends crucially on resolution, this will require higher resolution to have its full effect. The distribution of storms between different regions does not, at first approximation, depend on the temperature itself but on the distribution of the SST anomalies and their influence on the atmospheric circulation. Two additional transient experiments at T319 resolution where run for 20 years at the end of the 20th and 21st century, respectively using the same conditions as in the T213 experiments. The results are consistent with the T213 study. The total number of tropical cyclones were similar to the T213 experiment but were generally more intense. The change from the 20th to the 21st century was also similar with fewer TC in total but with more intense cyclones.

Deep Convective System Evolution over Africa and the Tropical Atlantic

In the tropical African and neighboring Atlantic region there is a strong contrast in the properties of deep convection between land and ocean. Here, satellite radar observations are used to produce a composite picture of the life cycle of convection in these two regions. Estimates of the broadband thermal flux from the geostationary Meteosat-8 satellite are used to identify and track organized convective systems over their life cycle. The evolution of the system size and vertical extent are used to define five life cycle stages (warm and cold developing, mature, cold and warm dissipating), providing the basis for the composite analysis of the system evolution. The tracked systems are matched to overpasses of the Tropical Rainfall Measuring Mission satellite, and a composite picture of the evolution of various radar and lightning characteristics is built up. The results suggest a fundamental difference in the convective life cycle between land and ocean. African storms evolve from convectively active systems with frequent lightning in their developing stages to more stratiform conditions as they dissipate. Over the Atlantic, the convective fraction remains essentially constant into the dissipating stages, and lightning occurrence peaks late in the life cycle. This behavior is consistent with differences in convective sustainability in land and ocean regions as proposed in previous studies. The area expansion rate during the developing stages of convection is used to provide an estimate of the intensity of convection. Reasonable correlations are found between this index and the convective system lifetime, size, and depth.