A high-order arbitrarily unstructured finite-volume model of the global atmosphere: tests solving the shallow-water equations

Simulations of the global atmosphere for weather and climate forecasting require fast and accurate solutions and so operational models use high-order finite differences on regular structured grids. This precludes the use of local refinement; techniques allowing local refinement are either expensive (eg. high-order finite element techniques) or have reduced accuracy at changes in resolution (eg. unstructured finite-volume with linear differencing). We present solutions of the shallow-water equations for westerly flow over a mid-latitude mountain from a finite-volume model written using OpenFOAM. A second/third-order accurate differencing scheme is applied on arbitrarily unstructured meshes made up of various shapes and refinement patterns. The results are as accurate as equivalent resolution spectral methods. Using lower order differencing reduces accuracy at a refinement pattern which allows errors from refinement of the mountain to accumulate and reduces the global accuracy over a 15 day simulation. We have therefore introduced a scheme which fits a 2D cubic polynomial approximately on a stencil around each cell. Using this scheme means that refinement of the mountain improves the accuracy after a 15 day simulation. This is a more severe test of local mesh refinement for global simulations than has been presented but a realistic test if these techniques are to be used operationally. These efficient, high-order schemes may make it possible for local mesh refinement to be used by weather and climate forecast models.

Potential predictability of rapid changes in the Atlantic meridional overturning circulation

We explore the potential predictability of rapid changes in the Atlantic meridional overturning circulation (MOC) using a coupled global climate model (HadCM3). Rapid changes in the temperature and salinity of surface water in the Nordic Seas, and the flow of dense water through Denmark Strait, are found to be precursors to rapid changes in the model's MOC, with a lead time of around 10 years. The mechanism proposed to explain this potential predictability relies on the development of density anomalies in the Nordic Seas which propagate through Denmark Strait and along the deep western boundary current, affecting the overturning. These rapid changes in the MOC have significant, and widespread, climate impacts which are potentially predictable a few years ahead. Whilst the flow through Denmark Strait is too strong in HadCM3, the presence of such potential predictability motivates the monitoring of water properties in the Nordic Seas and Denmark Strait.

An Integrated Assessment of changes in the thermohaline circulation

This paper discusses the risks of a shutdown of the thermohaline circulation (THC) for the climate system, for ecosystems in and around the North Atlantic as well as for fisheries and agriculture by way of an Integrated Assessment. The climate model simulations are based on greenhouse gas scenarios for the 21st century and beyond. A shutdown of the THC, complete by 2150, is triggered if increased freshwater input from inland ice melt or enhanced runoff is assumed. The shutdown retards the greenhouse gas-induced atmospheric warming trend in the Northern Hemisphere, but does not lead to a persistent net cooling. Due to the simulated THC shutdown the sea level at the North Atlantic shores rises by up to 80 cm by 2150, in addition to the global sea level rise. This could potentially be a serious impact that requires expensive coastal protection measures. A reduction of marine net primary productivity is associated with the impacts of warming rather than a THC shutdown. Regional shifts in the currents in the Nordic Seas could strongly deteriorate survival chances for cod larvae and juveniles. This could lead to cod fisheries becoming unprofitable by the end of the 21st century. While regional socioeconomic impacts might be large, damages would be probably small in relation to the respective gross national products. Terrestrial ecosystem productivity is affected much more by the fertilization from the increasing CO2 concentration than by a THC shutdown. In addition, the level of warming in the 22nd to 24th century favours crop production in northern Europe a lot, no matter whether the THC shuts down or not. CO2 emissions corridors aimed at limiting the risk of a THC breakdown to 10% or less are narrow, requiring departure from business-as-usual in the next few decades. The uncertainty about THC risks is still high. This is seen in model analyses as well as in the experts’ views that were elicited. The overview of results presented here is the outcome of the Integrated Assessment project INTEGRATION.

Relative humidity changes in a warmer climate

Key climate feedbacks due to water vapor and clouds rest largely on how relative humidity R changes in a warmer climate, yet this has not been extensively analyzed in models. General circulation models (GCMs) from the CMIP3 archive and several higher resolution atmospheric GCMs examined here generally predict a characteristic pattern of R trend with global temperature that has been reported previously in individual models, including increase around the tropopause, decrease in the tropical upper troposphere, and decrease in midlatitudes. This pattern is very similar to that previously reported for cloud cover in the same GCMs, confirming the role of R in controlling changes in simulated cloud. Comparing different models, the trend in each part of the troposphere is approximately proportional to the upward and/or poleward gradient of R in the present climate. While this suggests that the changes simply reflect a shift of the R pattern upward with the tropopause and poleward with the zonal jets, the drying trend in the subtropics is roughly three times too large to be attributable to shifts of subtropical features, and the subtropical R minima deepen in most models. R trends are correlated with horizontal model resolution, especially outside the tropics, where they show signs of convergence and latitudinal gradients become close to available observations for GCM resolutions near T85 and higher. We argue that much of the systematic change in R can be explained by the local specific humidity having been set (by condensation) in remote regions with different temperature changes, hence the gradients and trends each depend on a model’s ability to resolve moisture transport. Finally, subtropical drying trends predicted from the warming alone fall well short of those observed in recent decades. While this discrepancy supports previous reports of GCMs underestimating Hadley Cell expansion, our results imply that shifts alone are not a sufficient interpretation of changes.

Predicted soil organic carbon stocks and changes in Jordan between 2000 and 2030 made using the GEFSOC modelling system

Estimates of soil organic carbon (SOC) stocks and changes under different land use systems can help determine vulnerability to land degradation. Such information is important for countries in and areas with high susceptibility to desertification. SOC stocks, and predicted changes between 2000 and 2030, were determined at the national scale for Jordan using The Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System. For the purpose of this study, Jordan was divided into three natural regions (The Jordan Valley, the Uplands and the Badia) and three developmental regions (North, Middle and South). Based on this division, Jordan was divided into five zones (based on the dominant land use): the Jordan Valley, the North Uplands, the Middle Uplands, the South Uplands and the Badia. This information was merged using GIS, along with a map of rainfall isohyets, to produce a map with 498 polygons. Each of these was given a unique ID, a land management unit identifier and was characterized in terms of its dominant soil type. Historical land use data, current land use and future land use change scenarios were also assembled, forming major inputs of the modelling system. The GEFSOC Modelling System was then run to produce C stocks in Jordan for the years 1990, 2000 and 2030. The results were compared with conventional methods of estimating carbon stocks, such as the mapping based SOTER method. The results of these comparisons showed that the model runs are acceptable, taking into consideration the limited availability of long-term experimental soil data that can be used to validate them. The main findings of this research show that between 2000 and 2030, SOC may increase in heavily used areas under irrigation and will likely decrease in grazed rangelands that cover most of Jordan giving an overall decrease in total SOC over time if the land is indeed used under the estimated forms of land use. (C) 2007 Elsevier B.V. All rights reserved.

Modelled soil organic carbon stocks and changes in the Indo-Gangetic Plains, India from 1980 to 2030

The Global Environment Facility co-financed Soil Organic Carbon (GEFSOC) Project developed a comprehensive modelling system for predicting soil organic carbon (SOC) stocks and changes over time. This research is an effort to predict SOC stocks and changes for the Indian, Indo-Gangetic Plains (IGP), an area with a predominantly rice (Oryza sativa) - wheat (Triticum aestivum) cropping system, using the GEFSOC Modelling System and to compare output with stocks generated using mapping approaches based on soil survey data. The GEFSOC Modelling System predicts an estimated SOC stock for the IGP, India of 1.27, 1.32 and 1.27 Pg for 1990, 2000 and 2030, respectively, in the top 20 cm of soil. The SOC stock using a mapping approach based on soil survey data was 0.66 and 0.88 Pg for 1980 and 2000, respectively. The SOC stock estimated using the GEFSOC Modelling System is higher than the stock estimated using the mapping approach. This is due to the fact that while the GEFSOC System accounts for variation in crop input data (crop management), the soil mapping approach only considers regional variation in soil texture and wetness. The trend of overall change in the modelled SOC stock estimates shows that the IGP, India may have reached an equilibrium following 30-40 years of the Green Revolution. This can be seen in the SOC stock change rates. Various different estimation methods show SOC stocks of 0.57-1.44 Pg C for the study area. The trend of overall change in C stock assessed from the soil survey data indicates that the soils of the IGP, India may store a projected 1.1 Pg of C in 2030. (C) 2007 Elsevier B.V. All rights reserved.

Simulating SOC changes in 11 land use change chronosequences from the Brazilian Amazon with RothC and Century models

Land use and land cover changes in the Brazilian Amazon have major implications for regional and global carbon (C) cycling. Cattle pasture represents the largest single use (about 70%) of this once-forested land in most of the region. The main objective of this study was to evaluate the accuracy of the RothC and Century models at estimating soil organic C (SOC) changes under forest-to-pasture conditions in the Brazilian Amazon. We used data from 11 site-specific 'forest to pasture' chronosequences with the Century Ecosystem Model (Century 4.0) and the Rothamsted C Model (RothC 26.3). The models predicted that forest clearance and conversion to well managed pasture would cause an initial decline in soil C stocks (0-20 cm depth), followed in the majority of cases by a slow rise to levels exceeding those under native forest. One exception to this pattern was a chronosequence in Suia-Missu, which is under degraded pasture. In three other chronosequences the recovery of soil C under pasture appeared to be only to about the same level as under the previous forest. Statistical tests were applied to determine levels of agreement between simulated SOC stocks and observed stocks for all the sites within the 11 chronosequences. The models also provided reasonable estimates (coefficient of correlation = 0.8) of the microbial biomass C in the 0-10 cm soil layer for three chronosequences, when compared with available measured data. The Century model adequately predicted the magnitude and the overall trend in delta C-13 for the six chronosequences where measured 813 C data were available. This study gave independent tests of model performance, as no adjustments were made to the models to generate outputs. Our results suggest that modelling techniques can be successfully used for monitoring soil C stocks and changes, allowing both the identification of current patterns in the soil and the projection of future conditions. Results were used and discussed not only to evaluate soil C dynamics but also to indicate soil C sequestration opportunities for the Brazilian Amazon region. Moreover, modelling studies in these 'forest to pasture' systems have important applications, for example, the calculation of CO, emissions from land use change in national greenhouse gas inventories. (0 2007 Elsevier B.V. All rights reserved.

Predicted soil organic carbon stocks and changes in the Brazilian Amazon between 2000 and 2030

Currently we have little understanding of the impacts of land use change on soil C stocks in the Brazilian Amazon. Such information is needed to determine impacts'6n the global C cycle and the sustainability of agricultural systems that are replacing native forest. The aim of this study was to predict soil carbon stocks and changes in the Brazilian Amazon during the period between 2000 and 2030, using the GEFSOC soil carbon (C) modelling system. In order to do so, we devised current and future land use scenarios for the Brazilian Amazon, taking into account: (i) deforestation, rates from the past three decades, (ii) census data on land use from 1940 to 2000, including the expansion and intensification of agriculture in the region, (iii) available information on management practices, primarily related to well managed pasture versus degraded pasture and conventional systems versus no-tillage systems for soybean (Glycine max) and (iv) FAO predictions on agricultural land use and land use changes for the years 2015 and 2030. The land use scenarios were integrated with spatially explicit soils data (SOTER database), climate, potential natural vegetation and land management units using the recently developed GEFSOC soil C modelling system. Results are presented in map, table and graph form for the entire Brazilian Amazon for the current situation (1990 and 2000) and the future (2015 and 2030). Results include soil organic C (SOC) stocks and SOC stock change rates estimated by three methods: (i) the Century ecosystem model, (ii) the Rothamsted C model and (iii) the intergovernmental panel on climate change (IPCC) method for assessing soil C at regional scale. In addition, we show estimated values of above and belowground biomass for native vegetation, pasture and soybean. The results on regional SOC stocks compare reasonably well with those based on mapping approaches. The GEFSOC system provided a means of efficiently handling complex interactions among biotic-edapho-climatic conditions (> 363,000 combinations) in a very large area (similar to 500 Mha) such as the Brazilian Amazon. All of the methods used showed a decline in SOC stock for the period studied; Century and RothC simulated values for 2030 being about 7% lower than those in 1990. Values from Century and RothC (30,430 and 25,000 Tg for the 0-20 cm layer for the Brazilian Amazon region were higher than those obtained from the IPCC system (23,400 Tg in the 0-30 cm layer). Finally; our results can help understand the major biogeochemical cycles that influence soil fertility and help devise management strategies that enhance the sustainability of these areas and thus slow further deforestation. (C) 2007 Elsevier B.V. All rights reserved.

Predicted soil organic carbon stocks and changes in Kenya between 1990 and 2030

Under the United Nations Framework Convention on Climate Change (UNFCCC), Non-Annex 1 countries such as Kenya are obliged to report green house gas (GHG) emissions from all sources where possible, including those from soils as a result of changes in land use or land management. At present, the convention encourages countries to estimate emissions using the most advanced methods possible, given the country circumstances and resources. Estimates of soil organic carbon (SOC) stocks and changes were made for Kenya using the Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System. The tool conducts analysis using three methods: (1) the Century general ecosystem model; (2) the RothC soil C decomposition model; and (3) the Intergovernmental Panel on Climate Change (IPCC) method for assessing soil C at regional scales. The required datasets included: land use history, monthly mean precipitation, monthly mean minimum and maximum temperatures for all the agro-climatic zones of Kenya and historical vegetation cover. Soil C stocks of 1.4-2.0 Pg (0-20 cm), compared well with a Soil and Terrain (SOTER) based approach that estimated similar to .8-2.0 Pg (0-30 cm). In 1990 48% of the country had SOC stocks of < 18 t C ha(-1) and 20% of the country had SOC stocks of 18-30 t C ha(-1), whereas in 2000 56% of the country had SOC stocks of < 18 t C ha(-1) and 31% of the country had SOC stocks of 18-30 t C ha(-1). Conversion of natural vegetation to annual crops led to the greatest soil C losses. Simulations suggest that soil C losses remain substantial throughout the modelling period of 1990-2030. All three methods involved in the GEFSOC System estimated that there would be a net loss of soil C between 2000 and 2030 in Kenya. The decline was more marked with RothC than with Century or the IPCC method. In non-hydric soils the SOC change rates were more pronounced in high sandy soils compared to high clay soils in most land use systems. (C) 2007 Elsevier B.V. All rights reserved.

Changes in the extent of glaciers in the Eastern Terskey Alatoo, the Central Tien-Shan, in response to climatic fluctuations in between the end of the 19th and the beginning of the 21st century

Changes in the extent of glaciers and rates of glacier termini retreat in the eastern Terskey-Alatoo Range, the Tien Shan Mountains, Central Asia have been evaluated using the remote sensing techniques. Changes in the extent of 335 glaciers between the end of the Little Ice Age (LIA; mid-19th century), 1990 and 2003 have been estimated through the delineation of glacier outlines and the LIA moraine positions on the Landsat TM and ASTER imagery for 1990 and 2003 respectively. By 2003, the glacier surface area had decreased by 19% of the LIA value, which constitutes a 76 km(2) reduction in glacier surface area. Mapping of 109 glaciers using the 1965 1:25,000 maps revealed that glacier surface area decreased by 12.6% of the 1965 value between 1965 and 2003. Detailed mapping of 10 glaciers using historical maps and aerial photographs from the 1943-1977 period, has enabled glacier extent variations over the 20th century to be identified with a higher temporal resolution. Glacial retreat was slow in the early 20th century but increased considerably between 1943 and 1956 and then again after 1977. The post-1990 period has been marked by the most rapid glacier retreat since the end of the LIA. The observed changes in the extent of glaciers are in line with the observed climatic warming. The regional weather stations have revealed a strong climatic warming during the ablation season since the 1950s at a rate of 0.02-0.03 degrees Ca-1. At the higher elevations in the study area represented by the Tien Shan meteorological station, the summer warming was accompanied by negative anomalies in annual precipitation in the 1990s enhancing glacier retreat. However, trends in precipitation in the post-1997 period cannot be evaluated due to the change in observational practices at this station. Neither station in the study area exhibits significant long-term trends in precipitation. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

National and sub-national assessments of soil organic carbon stocks and changes: The GEFSOC modelling system

Soil organic carbon (SOC) plays a vital role in ecosystem function, determining soil fertility, water holding capacity and susceptibility to land degradation. In addition, SOC is related to atmospheric CO, levels with soils having the potential for C release or sequestration, depending on land use, land management and climate. The United Nations Convention on Climate Change and its Kyoto Protocol, and other United Nations Conventions to Combat Desertification and on Biodiversity all recognize the importance of SOC and point to the need for quantification of SOC stocks and changes. An understanding of SOC stocks and changes at the national and regional scale is necessary to further our understanding of the global C cycle, to assess the responses of terrestrial ecosystems to climate change and to aid policy makers in making land use/management decisions. Several studies have considered SOC stocks at the plot scale, but these are site specific and of limited value in making inferences about larger areas. Some studies have used empirical methods to estimate SOC stocks and changes at the regional scale, but such studies are limited in their ability to project future changes, and most have been carried out using temperate data sets. The computational method outlined by the Intergovernmental Panel on Climate Change (IPCC) has been used to estimate SOC stock changes at the regional scale in several studies, including a recent study considering five contrasting eco regions. This 'one step' approach fails to account for the dynamic manner in which SOC changes are likely to occur following changes in land use and land management. A dynamic modelling approach allows estimates to be made in a manner that accounts for the underlying processes leading to SOC change. Ecosystem models, designed for site scale applications can be linked to spatial databases, giving spatially explicit results that allow geographic areas of change in SOC stocks to be identified. Some studies have used variations on this approach to estimate SOC stock changes at the sub-national and national scale for areas of the USA and Europe and at the watershed scale for areas of Mexico and Cuba. However, a need remained for a national and regional scale, spatially explicit system that is generically applicable and can be applied to as wide a range of soil types, climates and land uses as possible. The Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System was developed in response to this need. The GEFSOC system allows estimates of SOC stocks and changes to be made for diverse conditions, providing essential information for countries wishing to take part in an emerging C market, and bringing us closer to an understanding of the future role of soils in the global C cycle. (C) 2007 Elsevier B.V. All rights reserved.

Evidence for long-term regional changes in precipitation on the East Coast Mountains in Mauritius

Global climate change and its impacts are being increasingly studied and precipitation trends are one of the measures of quantifying climate change especially in the tropics. This study uses daily rainfall data to determine if there are changes in the long-term trends in rainfall variability in the East Coast Mountains of Mauritius during the last few decades, and to investigate the factors influencing the trends in the inter-annual to inter-decadal rainfall variability. Statistical modelling has been used to investigate the trends in total seasonal rainfall, the number of rain days and the mean amount of rain per rainy days and the local, regional and large-scale factors that affect them on inter-annual to inter-decadal time scales. The strongest inter-decadal trend was found in the number of rain days for both rainfall seasons, and the other variables were found to have weak or insignificant trends. Both local factors, such as the surrounding sea surface temperatures and large-scale phenomena such as Indian Monsoon and the El Niño Southern Oscillation were found to influence rainfall patterns.

Observational constraints on past attributable warming and predictions of future global warming

This paper investigates the impact of aerosol forcing uncertainty on the robustness of estimates of the twentieth-century warming attributable to anthropogenic greenhouse gas emissions. Attribution analyses on three coupled climate models with very different sensitivities and aerosol forcing are carried out. The Third Hadley Centre Coupled Ocean - Atmosphere GCM (HadCM3), Parallel Climate Model (PCM), and GFDL R30 models all provide good simulations of twentieth-century global mean temperature changes when they include both anthropogenic and natural forcings. Such good agreement could result from a fortuitous cancellation of errors, for example, by balancing too much ( or too little) greenhouse warming by too much ( or too little) aerosol cooling. Despite a very large uncertainty for estimates of the possible range of sulfate aerosol forcing obtained from measurement campaigns, results show that the spatial and temporal nature of observed twentieth-century temperature change constrains the component of past warming attributable to anthropogenic greenhouse gases to be significantly greater ( at the 5% level) than the observed warming over the twentieth century. The cooling effects of aerosols are detected in all three models. Both spatial and temporal aspects of observed temperature change are responsible for constraining the relative roles of greenhouse warming and sulfate cooling over the twentieth century. This is because there are distinctive temporal structures in differential warming rates between the hemispheres, between land and ocean, and between mid- and low latitudes. As a result, consistent estimates of warming attributable to greenhouse gas emissions are obtained from all three models, and predictions are relatively robust to the use of more or less sensitive models. The transient climate response following a 1% yr(-1) increase in CO2 is estimated to lie between 2.2 and 4 K century(-1) (5-95 percentiles).

The global response to tropical heating in the Madden-Julian oscillation during the northern winter

A life cycle of the Madden–Julian oscillation (MJO) was constructed, based on 21 years of outgoing long-wave radiation data. Regression maps of NCEP–NCAR reanalysis data for the northern winter show statistically significant upper-tropospheric equatorial wave patterns linked to the tropical convection anomalies, and extratropical wave patterns over the North Pacific, North America, the Atlantic, the Southern Ocean and South America. To assess the cause of the circulation anomalies, a global primitive-equation model was initialized with the observed three-dimensional (3D) winter climatological mean flow and forced with a time-dependent heat source derived from the observed MJO anomalies. A model MJO cycle was constructed from the global response to the heating, and both the tropical and extratropical circulation anomalies generally matched the observations well. The equatorial wave patterns are established in a few days, while it takes approximately two weeks for the extratropical patterns to appear. The model response is robust and insensitive to realistic changes in damping and basic state. The model tropical anomalies are consistent with a forced equatorial Rossby–Kelvin wave response to the tropical MJO heating, although it is shifted westward by approximately 20° longitude relative to observations. This may be due to a lack of damping processes (cumulus friction) in the regions of convective heating. Once this shift is accounted for, the extratropical response is consistent with theories of Rossby wave forcing and dispersion on the climatological flow, and the pattern correlation between the observed and modelled extratropical flow is up to 0.85. The observed tropical and extratropical wave patterns account for a significant fraction of the intraseasonal circulation variance, and this reproducibility as a response to tropical MJO convection has implications for global medium-range weather prediction. Copyright © 2004 Royal Meteorological Society