Wave activity in the tropical tropopause layer in seven reanalysis and four chemistry climate model data sets

Sub-seasonal variability including equatorial waves significantly influence the dehydration and transport processes in the tropical tropopause layer (TTL). This study investigates the wave activity in the TTL in 7 reanalysis data sets (RAs; NCEP1, NCEP2, ERA40, ERA-Interim, JRA25, MERRA, and CFSR) and 4 chemistry climate models (CCMs; CCSRNIES, CMAM, MRI, and WACCM) using the zonal wave number-frequency spectral analysis method with equatorially symmetric-antisymmetric decomposition. Analyses are made for temperature and horizontal winds at 100 hPa in the RAs and CCMs and for outgoing longwave radiation (OLR), which is a proxy for convective activity that generates tropopause-level disturbances, in satellite data and the CCMs. Particular focus is placed on equatorial Kelvin waves, mixed Rossby-gravity (MRG) waves, and the Madden-Julian Oscillation (MJO). The wave activity is defined as the variance, i.e., the power spectral density integrated in a particular zonal wave number-frequency region. It is found that the TTL wave activities show significant difference among the RAs, ranging from ∼0.7 (for NCEP1 and NCEP2) to ∼1.4 (for ERA-Interim, MERRA, and CFSR) with respect to the averages from the RAs. The TTL activities in the CCMs lie generally within the range of those in the RAs, with a few exceptions. However, the spectral features in OLR for all the CCMs are very different from those in the observations, and the OLR wave activities are too low for CCSRNIES, CMAM, and MRI. It is concluded that the broad range of wave activity found in the different RAs decreases our confidence in their validity and in particular their value for validation of CCM performance in the TTL, thereby limiting our quantitative understanding of the dehydration and transport processes in the TTL.

Observational challenges in evaluating climate models

Activities like the Coupled Model Intercomparison Project (CMIP) have revolutionized climate modelling in terms of our ability to compare models and to process information about climate projections and their uncertainties. The evaluation of models against observations is now considered a key component of multi-model studies. While there are a number of outstanding scientific issues surrounding model evaluation, notably the open question of how to link model performance to future projections, here we highlight a specific but growing problem in model evaluation - that of uncertainties in the observational data that are used to evaluate the models. We highlight the problem using an example obtained from studies of the South Asian Monsoon but we believe the problem is a generic one which arises in many different areas of climate model evaluation and which requires some attention by the community.

An urban parameterization for a global climate model. Part I: Formulation and evaluation for two cities

Urbanization, the expansion of built-up areas, is an important yet less-studied aspect of land use/land cover change in climate science. To date, most global climate models used to evaluate effects of land use/land cover change on climate do not include an urban parameterization. Here, the authors describe the formulation and evaluation of a parameterization of urban areas that is incorporated into the Community Land Model, the land surface component of the Community Climate System Model. The model is designed to be simple enough to be compatible with structural and computational constraints of a land surface model coupled to a global climate model yet complex enough to explore physically based processes known to be important in determining urban climatology. The city representation is based upon the “urban canyon” concept, which consists of roofs, sunlit and shaded walls, and canyon floor. The canyon floor is divided into pervious (e.g., residential lawns, parks) and impervious (e.g., roads, parking lots, sidewalks) fractions. Trapping of longwave radiation by canyon surfaces and solar radiation absorption and reflection is determined by accounting for multiple reflections. Separate energy balances and surface temperatures are determined for each canyon facet. A one-dimensional heat conduction equation is solved numerically for a 10-layer column to determine conduction fluxes into and out of canyon surfaces. Model performance is evaluated against measured fluxes and temperatures from two urban sites. Results indicate the model does a reasonable job of simulating the energy balance of cities.

Climate impacts on river flow: projections for the Medway catchment UK with UKCP09 and CATCHMOD

The potential impact of climate change on areas of strategic importance for water resources remains a concern. Here, river flow projections for the River Medway, above Teston in southeast England are presented, which is just such an area of strategic importance. The river flow projections use climate inputs from the Hadley Centre Regional Climate Model (HadRM3) for the time period 1960–2080 (a subset of the early release UKCP09 projections). River flow predictions are calculated using CATCHMOD, the main river flow prediction tool of the Environment Agency (EA) of England and Wales. In order to use this tool in the best way for climate change predictions, model setup and performance are analysed using sensitivity and uncertainty analysis. The model's representation of hydrological processes is discussed and the direct percolation and first linear storage constant parameters are found to strongly affect model results in a complex way, with the former more important for low flows and the latter for high flows. The uncertainty in predictions resulting from the hydrological model parameters is demonstrated and the projections of river flow under future climate are analysed. A clear climate change impact signal is evident in the results with a persistent lowering of mean daily river flows for all months and for all projection time slices. Results indicate that a projection of lower flows under future climate is valid even taking into account the uncertainties considered in this modelling chain exercise. The model parameter uncertainty becomes more significant under future climate as the river flows become lower. This has significant implications for those making policy decisions based on such modelling results. Copyright © 2010 John Wiley & Sons, Ltd.

Toward seamless prediction: calibration of climate change projections using seasonal forecasts

The CWRF is developed as a climate extension of the Weather Research and Forecasting model (WRF) by incorporating numerous improvements in the representation of physical processes and integration of external (top, surface, lateral) forcings that are crucial to climate scales, including interactions between land, atmosphere, and ocean; convection and microphysics; and cloud, aerosol, and radiation; and system consistency throughout all process modules. This extension inherits all WRF functionalities for numerical weather prediction while enhancing the capability for climate modeling. As such, CWRF can be applied seamlessly to weather forecast and climate prediction. The CWRF is built with a comprehensive ensemble of alternative parameterization schemes for each of the key physical processes, including surface (land, ocean), planetary boundary layer, cumulus (deep, shallow), microphysics, cloud, aerosol, and radiation, and their interactions. This facilitates the use of an optimized physics ensemble approach to improve weather or climate prediction along with a reliable uncertainty estimate. The CWRF also emphasizes the societal service capability to provide impactrelevant information by coupling with detailed models of terrestrial hydrology, coastal ocean, crop growth, air quality, and a recently expanded interactive water quality and ecosystem model. This study provides a general CWRF description and basic skill evaluation based on a continuous integration for the period 1979– 2009 as compared with that of WRF, using a 30-km grid spacing over a domain that includes the contiguous United States plus southern Canada and northern Mexico. In addition to advantages of greater application capability, CWRF improves performance in radiation and terrestrial hydrology over WRF and other regional models. Precipitation simulation, however, remains a challenge for all of the tested models.

Optimizing performance of the lightshelf by modifying ceiling geometry in highly luminous climates

The impact of ceiling geometries on the performance of lightshelves was investigated using physical model experiments and radiance simulations. Illuminance level and distribution uniformity were assessed for a working plane in a large space located in sub-tropical climate regions where innovative systems for daylighting and shading are required. It was found that the performance of the lightshelf can be improved by changing the ceiling geometry; the illuminance level increased in the rear of the room and decreased in the front near the window compared to rooms having conventional horizontal ceilings. Moreover, greater uniformity was achieved throughout the room as a result of reducing the difference in the illuminance level between the front and rear of the room. Radiance simulation results were found to be in good agreement with physical model data obtained under a clear sky and high solar radiation. The best ceiling shape was found to be one that is curved in the front and rear of the room.

Climate model benchmarking with glacial and mid-Holocene climates

Past climates provide a test of models’ ability to predict climate change. We present a comprehensive evaluation of state-of-the-art models against Last Glacial Maximum and mid-Holocene climates, using reconstructions of land and ocean climates and simulations from the Palaeoclimate Modelling and Coupled Modelling Intercomparison Projects. Newer models do not perform better than earlier versions despite higher resolution and complexity. Differences in climate sensitivity only weakly account for differences in model performance. In the glacial, models consistently underestimate land cooling (especially in winter) and overestimate ocean surface cooling (especially in the tropics). In the mid-Holocene, models generally underestimate the precipitation increase in the northern monsoon regions, and overestimate summer warming in central Eurasia. Models generally capture large-scale gradients of climate change but have more limited ability to reproduce spatial patterns. Despite these common biases, some models perform better than others.

Global patterns of vegetation response to millennial-scale variability and rapid climate change during the last glacial period

Ninety-four sites worldwide have sufficient resolution and dating to document the impact of millennial-scale climate variability on vegetation and fire regimes during the last glacial period. Although Dansgaard–Oeschger (D–O) cycles all show a basically similar gross structure, they vary in the magnitude and the length of the warm and cool intervals. We illustrate the geographic patterns in the climate-induced changes in vegetation by comparing D–O 6, D–O 8 and D–O 19. There is a strong response to both D–O warming events and subsequent cooling, most marked in the northern extratropics. Pollen records from marine cores from the northern extratropics confirm that there is no lag between the change in climate and the vegetation response, within the limits of the dating resolution (50–100 years). However, the magnitude of the change in vegetation is regionally specific and is not a simple function of either the magnitude or the duration of the change in climate as registered in Greenland ice cores. Fire regimes also show an initial immediate response to climate changes, but during cooling intervals there is a slow recovery of biomass burning after the initial reduction, suggesting a secondary control through the recovery of vegetation productivity. In the extratropics, vegetation changes are largely determined by winter temperatures while in the tropics they are largely determined by changes in plant-available water. Tropical vegetation records show changes corresponding to Heinrich Stadials but the response to D–O warming events is less marked than in the northern extratropics. There are very few high-resolution records from the Southern Hemisphere extratropics, but these records also show both a vegetation and fire response to millennial-scale climate variability. It is not yet possible to determine unequivocally whether terrestrial records reflect the asynchroneity apparent in the ice-core records.

Recovery mechanisms of Arctic summer sea ice

We examine the recovery of Arctic sea ice from prescribed ice-free summer conditions in simulations of 21st century climate in an atmosphere–ocean general circulation model. We find that ice extent recovers typically within two years. The excess oceanic heat that had built up during the ice-free summer is rapidly returned to the atmosphere during the following autumn and winter, and then leaves the Arctic partly through increased longwave emission at the top of the atmosphere and partly through reduced atmospheric heat advection from lower latitudes. Oceanic heat transport does not contribute significantly to the loss of the excess heat. Our results suggest that anomalous loss of Arctic sea ice during a single summer is reversible, as the ice–albedo feedback is alleviated by large-scale recovery mechanisms. Hence, hysteretic threshold behavior (or a “tipping point”) is unlikely to occur during the decline of Arctic summer sea-ice cover in the 21st century.

Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialization of Northern-hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates significantly reduces assimilation error both in identical-twin experiments and when assimilating sea-ice observations, reducing the concentration error by a factor of four to six, and the thickness error by a factor of two. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that the strong dependence of thermodynamic ice growth on ice concentration necessitates an adjustment of mean ice thickness in the analysis update. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that proportional mean-thickness updates are superior to the other two methods considered and enable us to assimilate sea ice in a global climate model using simple Newtonian relaxation.

Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialisation of Northern Hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates leads to good assimilation performance for sea-ice concentration and thickness, both in identical-twin experiments and when assimilating sea-ice observations. The simulation of other Arctic surface fields in the coupled model is, however, not significantly improved by the assimilation. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that an adjustment of mean ice thickness in the analysis update is essential to arrive at plausible state estimates. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that assimilation with proportional mean-thickness updates outperforms the other two methods considered. The method described here is very simple to implement, and gives results that are sufficiently good to be used for initialising sea ice in a global climate model for seasonal to decadal predictions.

The development of a data-driven application benchmarking approach to performance modelling

Performance modelling is a useful tool in the lifeycle of high performance scientific software, such as weather and climate models, especially as a means of ensuring efficient use of available computing resources. In particular, sufficiently accurate performance prediction could reduce the effort and experimental computer time required when porting and optimising a climate model to a new machine. In this paper, traditional techniques are used to predict the computation time of a simple shallow water model which is illustrative of the computation (and communication) involved in climate models. These models are compared with real execution data gathered on AMD Opteron-based systems, including several phases of the U.K. academic community HPC resource, HECToR. Some success is had in relating source code to achieved performance for the K10 series of Opterons, but the method is found to be inadequate for the next-generation Interlagos processor. The experience leads to the investigation of a data-driven application benchmarking approach to performance modelling. Results for an early version of the approach are presented using the shallow model as an example.

Integrated analysis of climate, soil, topography and vegetative growth in Iberian viticultural regions

The Iberian viticultural regions are convened according to the Denomination of Origin (DO) and present different climates, soils, topography and management practices. All these elements influence the vegetative growth of different varieties throughout the peninsula, and are tied to grape quality and wine type. In the current study, an integrated analysis of climate, soil, topography and vegetative growth was performed for the Iberian DO regions, using state-of-the-art datasets. For climatic assessment, a categorized index, accounting for phenological/thermal development, water availability and grape ripening conditions was computed. Soil textural classes were established to distinguish soil types. Elevation and aspect (orientation) were also taken into account, as the leading topographic elements. A spectral vegetation index was used to assess grapevine vegetative growth and an integrated analysis of all variables was performed. The results showed that the integrated climate-soil-topography influence on vine performance is evident. Most Iberian vineyards are grown in temperate dry climates with loamy soils, presenting low vegetative growth. Vineyards in temperate humid conditions tend to show higher vegetative growth. Conversely, in cooler/warmer climates, lower vigour vineyards prevail and other factors, such as soil type and precipitation acquire more important roles in driving vigour. Vines in prevailing loamy soils are grown over a wide climatic diversity, suggesting that precipitation is the primary factor influencing vigour. The present assessment of terroir characteristics allows direct comparison among wine regions and may have great value to viticulturists, particularly under a changing climate.

Methane in the atmosphere of the transiting hot Neptune GJ436b?

We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5, and 8 μm obtained with the Infrared Array Camera on the Spitzer Space Telescope. After correcting for systematic effects, we fitted the light curves using the Markov Chain Monte Carlo technique. Combining these new data with the EPOXI, Hubble Space Telescope, and ground-based V, I, H, and Ks published observations, the range 0.5-10 μm can be covered. Due to the low level of activity of GJ436, the effect of starspots on the combination of transits at different epochs is negligible at the accuracy of the data set. Representative climate models were calculated by using a three-dimensional, pseudospectral general circulation model with idealized thermal forcing. Simulated transit spectra of GJ436b were generated using line-by-line radiative transfer models including the opacities of the molecular species expected to be present in such a planetary atmosphere. A new, ab-initio-calculated, line list for hot ammonia has been used for the first time. The photometric data observed at multiple wavelengths can be interpreted with methane being the dominant absorption after molecular hydrogen, possibly with minor contributions from ammonia, water, and other molecules. No clear evidence of carbon monoxide and carbon dioxide is found from transit photometry. We discuss this result in the light of a recent paper where photochemical disequilibrium is hypothesized to interpret secondary transit photometric data. We show that the emission photometric data are not incompatible with the presence of abundant methane, but further spectroscopic data are desirable to confirm this scenario.

Observed climate variations and change in Vietnam

In this study, observed changes of temperature, rainfall, and some extreme climate indices in Vietnam were investigated by using daily observations during the period 1961-2012. The observed data were collected from 80 meteorological stations for temperature, and from 170 stations for rainfall over the seven climatological sub-regions of Vietnam. Results show that there were insignificant differences between the trends of changes obtained from the 1961-2011 and 1979-2012 periods. Near-surface temperature, including mean (T2m), maximum (Tx) and minimum temperature (Tm), increased consistently at almost all stations. Tm increased faster than Tx. Temperature also increased faster in winter than in summer. Consequently, the number of hot days and warm nights increased whereas the number of cold days, cold nights and cool days decreased. In the northern regions, temperature tended to slightly decrease in May but significantly increased in June. Annual rainfall decreased in the northern area of Vietnam, while it increased at almost all stations in the central regions, and had insignificant trends in the southern sub-region. Changes in some extreme rainfall indices were likely consistent with changes in annual rainfall. Monthly rainfall in the central regions significantly increased from August to December. Rainfall generally increased in May and decreased in June over almost all country.