Summer drought in northern midlatitudes in a time-dependent CO2 climate experiment

A time-dependent climate-change experiment with a coupled ocean–atmosphere general circulation model has been used to study changes in the occurrence of drought in summer in southern Europe and central North America. In both regions, precipitation and soil moisture are reduced in a climate of greater atmospheric carbon dioxide. A detailed investigation of the hydrology of the model shows that the drying of the soil comes about through an increase in evaporation in winter and spring, caused by higher temperatures and reduced snow cover, and a decrease in the net input of water in summer. Evaporation is reduced in summer because of the drier soil, but the reduction in precipitation is larger. Three extreme statistics are used to define drought, namely the frequency of low summer precipitation, the occurrence of long dry spells, and the probability of dry soil. The last of these is arguably of the greatest practical importance, but since it is based on soil moisture, of which there are very few observations, the authors’ simulation of it has the least confidence. Furthermore, long time series for daily observed precipitation are not readily available from a sufficient number of stations to enable a thorough evaluation of the model simulation, especially for the frequency of long dry spells, and this increases the systematic uncertainty of the model predictions. All three drought statistics show marked increases owing to the sensitivity of extreme statistics to changes in their distributions. However, the greater likelihood of long dry spells is caused by a tendency in the character of daily rainfall toward fewer events, rather than by the reduction in mean precipitation. The results should not be taken as firm predictions because extreme statistics for small regions cannot be calculated reliably from the output of the current generation of GCMs, but they point to the possibility of large increases in the severity of drought conditions as a consequence of climate change caused by increased CO2.

Vertical heat transports in the ocean and their effect on time-dependent climate change

In response to increasing atmospheric con- centrations of greenhouse gases, the rate of time- dependent climate change is determined jointly by the strength of climate feedbacks and the e�ciency of pro- cesses which remove heat from the surface into the deep ocean. This work examines the vertical heat transport processes in the ocean of the HADCM2 atmosphere± ocean general circulation model (AOGCM) in experi- ments with CO2 held constant (control) and increasing at 1% per year (anomaly). The control experiment shows that global average heat exchanges between the upper and lower ocean are dominated by the Southern Ocean, where heat is pumped downwards by the wind- driven circulation and di�uses upwards along sloping isopycnals. This is the reverse of the low-latitude balance used in upwelling±di�usion ocean models, the global average upward di�usive transport being against the temperature gradient. In the anomaly experiment, weakened convection at high latitudes leads to reduced diffusive and convective heat loss from the deep ocean, and hence to net heat uptake, since the advective heat input is less a�ected. Reduction of deep water produc- tion at high latitudes results in reduced upwelling of cold water at low latitudes, giving a further contribution to net heat uptake. On the global average, high-latitude processes thus have a controlling in¯uence. The impor- tant role of di�usion highlights the need to ensure that the schemes employed in AOGCMs give an accurate representation of the relevant sub-grid-scale processes.

On the determinants of optimal border enforcement

We extend the current immigration-enforcement literature by incorporating both the practice of people smuggling and a role for non-wage income into a two-country, dynamic general equilibrium model. We use the model economy to examine three questions. First, how does technological progress in the smuggling industry affect the level of migration and capital accumulation for a given level of enforcement? Second, do changes in border enforcement affect the level of migration, capital accumulation, and smuggling activity? Third, is the optimal level of enforcement sensitive to technological progress in the smuggling industry? We show that the government chooses to devote resources to border enforcement only if the deterrent effect on smugglers is large enough. Otherwise, it is not worth taxing host-country natives as the taxes paid will more than offset any income gain resulting from fewer migrants.

Sea-ice anomalies in the Sea of Okhotsk and the relationship with storm tracks in the Northern Hemisphere during winter

Previous studies have shown that sea-ice in the Sea of Okhotsk can be affected by local storms; in turn, the resultant sea-ice changes can affect the downstream development of storm tracks in the Pacific and possibly dampen a pre-existing North Atlantic Oscillation (NAO) signal in late winter. In this paper, a storm tracking algorithm was applied to the six hourly horizontal winds from the National Centers for Environmental Prediction (NCEP) reanalysis data from 1978(9) to 2007 and output from the atmospheric general circulation model (AGCM) ECHAM5 forced by sea-ice anomalies in the Sea of Okhotsk. The life cycle response of storms to sea-ice anomalies is investigated using various aspects of storm activity—cyclone genesis, lysis, intensity and track density. Results show that, for enhanced positive sea-ice concentrations in the Sea of Okhotsk, there is a decrease in secondary cyclogenesis, a westward shift in cyclolysis and changes in the subtropical jet are seen in the North Pacific. In the Atlantic, a pattern resembling the negative phase of the NAO is observed. This pattern is confirmed by the AGCM ECHAM5 experiments driven with above normal sea-ice anomalies in the Sea of Okhotsk

An ocean modelling and assimilation guide to using GOCE geoid products

We review the procedures and challenges that must be considered when using geoid data derived from the Gravity and steady-state Ocean Circulation Explorer (GOCE) mission in order to constrain the circulation and water mass representation in an ocean 5 general circulation model. It covers the combination of the geoid information with timemean sea level information derived from satellite altimeter data, to construct a mean dynamic topography (MDT), and considers how this complements the time-varying sea level anomaly, also available from the satellite altimeter. We particularly consider the compatibility of these different fields in their spatial scale content, their temporal rep10 resentation, and in their error covariances. These considerations are very important when the resulting data are to be used to estimate ocean circulation and its corresponding errors. We describe the further steps needed for assimilating the resulting dynamic topography information into an ocean circulation model using three different operational fore15 casting and data assimilation systems. We look at methods used for assimilating altimeter anomaly data in the absence of a suitable geoid, and then discuss different approaches which have been tried for assimilating the additional geoid information. We review the problems that have been encountered and the lessons learned in order the help future users. Finally we present some results from the use of GRACE geoid in20 formation in the operational oceanography community and discuss the future potential gains that may be obtained from a new GOCE geoid.

Mountain ranges favour vigorous Atlantic Meridional Overturning

We use a global Ocean-Atmosphere General Circulation Model (OAGCM) to show that the major mountain ranges of the world have a significant role in maintenance of the Atlantic Meridional Overturning Circulation (AMOC). A simulation with mountains has a maximum AMOC of 18 Sv (1 Sv=106 m3 s-1) compared with ~0 Sv for a simulation without mountains. Atlantic heat transport at 25N is 1.1 PW with mountains compared to 0.2 PW without. The difference in AMOC is due to major changes in surface heat and freshwater (FW) fluxes over the Atlantic. In the Pacific changed surface fluxes lead to a meridional overturning circulation of 10 Sv. Our results suggest that the effects of mountains on the large-scale atmospheric circulation is to force the ocean towards a state with a vigorous AMOC and with no overturning in the Pacific.

Accessing extremes of mid-latitudinal wave activity: methodology and application

A statistical methodology is proposed and tested for the analysis of extreme values of atmospheric wave activity at mid-latitudes. The adopted methods are the classical block-maximum and peak over threshold, respectively based on the generalized extreme value (GEV) distribution and the generalized Pareto distribution (GPD). Time-series of the ‘Wave Activity Index’ (WAI) and the ‘Baroclinic Activity Index’ (BAI) are computed from simulations of the General Circulation Model ECHAM4.6, which is run under perpetual January conditions. Both the GEV and the GPD analyses indicate that the extremes ofWAI and BAI areWeibull distributed, this corresponds to distributions with an upper bound. However, a remarkably large variability is found in the tails of such distributions; distinct simulations carried out under the same experimental setup provide sensibly different estimates of the 200-yr WAI return level. The consequences of this phenomenon in applications of the methodology to climate change studies are discussed. The atmospheric configurations characteristic of the maxima and minima of WAI and BAI are also examined.

Changing climate and the Irish landscape

The impacts of current and future changes in climate have been investigated for Irish vegetation. Warming has been observed over the last two decades, with impacts that are also strongly influenced by natural oscillations of the surrounding ocean, seen as fluctuations in the North Atlantic Oscillation and the Atlantic Multidecadal Oscillation. Satellite observations show that vegetation greenness increases in warmer years, a feature mirrored by increases in net ecosystem production observed for a grassland and a plantation forest. An ensemble of general circulation model simulations of future climates indicate temperature rises over the twenty-first century ranging from 1°C to 7°C, depending on future scenarios of greenhouse gas emissions. Net primary production is simulated to increase under all scenarios, due to the positive impacts of rising temperature, a modest rise of precipitation and rising carbon dioxide concentrations. In an optimistic scenario of reducing future emissions, CO2 concentration is simulated to flatten from about 2070, although temperatures continue to increase. Under this scenario Ireland could become a source of carbon, whereas under all other emission scenarios Ireland is a sink for carbon that may increase by up to three-fold over the twenty-first century. A likely and unavoidable impact of changing climate is the arrival of alien plant species, which may disrupt ecosystems and exert negative impacts on native biodiversity. Alien species arrive continually, with about 250 dated arrivals in the twentieth century. A simulation model indicates that this rate of alien arrival may increase by anything between two and ten times, dependent on the future climatic scenario, by 2050. Which alien species may become severely disruptive is, however, not known.

Downscaling from GCM precipitation: a benchmark for dynamical and statistical downscaling methods

A precipitation downscaling method is presented using precipitation from a general circulation model (GCM) as predictor. The method extends a previous method from monthly to daily temporal resolution. The simplest form of the method corrects for biases in wet-day frequency and intensity. A more sophisticated variant also takes account of flow-dependent biases in the GCM. The method is flexible and simple to implement. It is proposed here as a correction of GCM output for applications where sophisticated methods are not available, or as a benchmark for the evaluation of other downscaling methods. Applied to output from reanalyses (ECMWF, NCEP) in the region of the European Alps, the method is capable of reducing large biases in the precipitation frequency distribution, even for high quantiles. The two variants exhibit similar performances, but the ideal choice of method can depend on the GCM/reanalysis and it is recommended to test the methods in each case. Limitations of the method are found in small areas with unresolved topographic detail that influence higher-order statistics (e.g. high quantiles). When used as benchmark for three regional climate models (RCMs), the corrected reanalysis and the RCMs perform similarly in many regions, but the added value of the latter is evident for high quantiles in some small regions.

Evaluation of water and energy budgets in regional climate models applied over Europe

This study presents a model intercomparison of four regional climate models (RCMs) and one variable resolution atmospheric general circulation model (AGCM) applied over Europe with special focus on the hydrological cycle and the surface energy budget. The models simulated the 15 years from 1979 to 1993 by using quasi-observed boundary conditions derived from ECMWF re-analyses (ERA). The model intercomparison focuses on two large atchments representing two different climate conditions covering two areas of major research interest within Europe. The first is the Danube catchment which represents a continental climate dominated by advection from the surrounding land areas. It is used to analyse the common model error of a too dry and too warm simulation of the summertime climate of southeastern Europe. This summer warming and drying problem is seen in many RCMs, and to a less extent in GCMs. The second area is the Baltic Sea catchment which represents maritime climate dominated by advection from the ocean and from the Baltic Sea. This catchment is a research area of many studies within Europe and also covered by the BALTEX program. The observed data used are monthly mean surface air temperature, precipitation and river discharge. For all models, these are used to estimate mean monthly biases of all components of the hydrological cycle over land. In addition, the mean monthly deviations of the surface energy fluxes from ERA data are computed. Atmospheric moisture fluxes from ERA are compared with those of one model to provide an independent estimate of the convergence bias derived from the observed data. These help to add weight to some of the inferred estimates and explain some of the discrepancies between them. An evaluation of these biases and deviations suggests possible sources of error in each of the models. For the Danube catchment, systematic errors in the dynamics cause the prominent summer drying problem for three of the RCMs, while for the fourth RCM this is related to deficiencies in the land surface parametrization. The AGCM does not show this drying problem. For the Baltic Sea catchment, all models similarily overestimate the precipitation throughout the year except during the summer. This model deficit is probably caused by the internal model parametrizations, such as the large-scale condensation and the convection schemes.

The Aqua-Planet Experiment (APE): Response to changed meridional SST profile

This paper explores the sensitivity of Atmospheric General Circulation Model (AGCM) simulations to changes in the meridional distribution of sea surface temperature (SST). The simulations are for an aqua-planet, a water covered Earth with no land, orography or sea-ice and with specified zonally symmetric SST. Simulations from 14 AGCMs developed for Numerical Weather Prediction and climate applications are compared. Four experiments are performed to study the sensitivity to the meridional SST profile. These profiles range from one in which the SST gradient continues to the equator to one which is flat approaching the equator, all with the same maximum SST at the equator. The zonal mean circulation of all models shows strong sensitivity to latitudinal distribution of SST. The Hadley circulation weakens and shifts poleward as the SST profile flattens in the tropics. One question of interest is the formation of a double versus a single ITCZ. There is a large variation between models of the strength of the ITCZ and where in the SST experiment sequence they transition from a single to double ITCZ. The SST profiles are defined such that as the equatorial SST gradient flattens, the maximum gradient increases and moves poleward. This leads to a weakening of the mid-latitude jet accompanied by a poleward shift of the jet core. Also considered are tropical wave activity and tropical precipitation frequency distributions. The details of each vary greatly between models, both with a given SST and in the response to the change in SST. One additional experiment is included to examine the sensitivity to an off-equatorial SST maximum. The upward branch of the Hadley circulation follows the SST maximum off the equator. The models that form a single precipitation maximum when the maximum SST is on the equator shift the precipitation maximum off equator and keep it centered over the SST maximum. Those that form a double with minimum on the equatorial maximum SST shift the double structure off the equator, keeping the minimum over the maximum SST. In both situations only modest changes appear in the shifted profile of zonal average precipitation. When the upward branch of the Hadley circulation moves into the hemisphere with SST maximum, the zonal average zonal, meridional and vertical winds all indicate that the Hadley cell in the other hemisphere dominates.

Radiative forcing due to aviation water vapour emissions

Three emissions inventories have been used with a fully Lagrangian trajectory model to calculate the stratospheric accumulation of water vapour emissions from aircraft, and the resulting radiative forcing. The annual and global mean radiative forcing due to present-day aviation water vapour emissions has been found to be 0.9 [0.3 to 1.4] mW m^2. This is around a factor of three smaller than the value given in recent assessments, and the upper bound is much lower than a recently suggested 20 mW m^2 upper bound. This forcing is sensitive to the vertical distribution of emissions, and, to a lesser extent, interannual variability in meteorology. Large differences in the vertical distribution of emissions within the inventories have been identified, which result in the choice of inventory being the largest source of differences in the calculation of the radiative forcing due to the emissions. Analysis of Northern Hemisphere trajectories demonstrates that the assumption of an e-folding time is not always appropriate for stratospheric emissions. A linear model is more representative for emissions that enter the stratosphere far above the tropopause.

Estimation of the friction velocity in stably stratified boundary-layer flows over hills

A method is suggested for the calculation of the friction velocity for stable turbulent boundary-layer flow over hills. The method is tested using a continuous upstream mean velocity profile compatible with the propagation of gravity waves, and is incorporated into the linear model of Hunt, Leibovich and Richards with the modification proposed by Hunt, Richards and Brighton to include the effects of stability, and the reformulated solution of Weng for the near-surface region. Those theoretical results are compared with results from simulations using a non-hydrostatic microscale-mesoscale two-dimensional numerical model, and with field observations for different values of stability. These comparisons show a considerable improvement in the behaviour of the theoretical model when the friction velocity is calculated using the method proposed here, leading to a consistent variation of the boundary-layer structure with stability, and better agreement with observational and numerical data.

Revisiting Asian monsoon formation and change associated with Tibetan Plateau forcing: II. Change

Data analysis based on station observations reveals that many meteorological variables averaged over the Tibetan Plateau (TP) are closely correlated, and their trends during the past decades are well correlated with the rainfall trend of the Asian summer monsoon. However, such correlation does not necessarily imply causality. Further diagnosis confirms the existence of a weakening trend in TP thermal forcing, characterized by weakened surface sensible heat flux in spring and summer during the past decades. This weakening trend is associated with decreasing summer precipitation over northern South Asia and North China and increasing precipitation over northwestern China, South China, and Korea. An atmospheric general circulation model, the HadAM3, is employed to elucidate the causality between the weakening TP forcing and the change in the Asian summer monsoon rainfall. Results demonstrate that a weakening in surface sensible heating over the TP results in reduced summer precipitation in the plateau region and a reduction in the associated latent heat release in summer. These changes in turn result in the weakening of the near-surface cyclonic circulation surrounding the plateau and the subtropical anticyclone over the subtropical western North Pacific, similar to the results obtained from the idealized TP experiment in Part I of this study. The southerly that normally dominates East Asia, ranging from the South China Sea to North China, weakens, resulting in a weaker equilibrated Sverdrup balance between positive vorticity generation and latent heat release. Consequently, the convergence of water vapor transport is confined to South China, forming a unique anomaly pattern in monsoon rainfall, the so-called “south wet and north dry.” Because the weakening trend in TP thermal forcing is associated with global warming, the present results provide an effective means for assessing projections of regional climate over Asia in the context of global warming.

How difficult is it to recover from dangerous levels of global warming?

Climate models provide compelling evidence that if greenhouse gas emissions continue at present rates, then key global temperature thresholds (such as the European Union limit of two degrees of warming since pre-industrial times) are very likely to be crossed in the next few decades. However, there is relatively little attention paid to whether, should a dangerous temperature level be exceeded, it is feasible for the global temperature to then return to safer levels in a usefully short time. We focus on the timescales needed to reduce atmospheric greenhouse gases and associated temperatures back below potentially dangerous thresholds, using a state-of-the-art general circulation model. This analysis is extended with a simple climate model to provide uncertainty bounds. We find that even for very large reductions in emissions, temperature reduction is likely to occur at a low rate. Policy-makers need to consider such very long recovery timescales implicit in the Earth system when formulating future emission pathways that have the potential to 'overshoot' particular atmospheric concentrations of greenhouse gases and, more importantly, related temperature levels that might be considered dangerous.