An assessment of the impact of local processes on dust lifting in martian climate models

Simulation of the lifting of dust from the planetary surface is of substantially greater importance on Mars than on Earth, due to the fundamental role that atmospheric dust plays in the former’s climate, yet the dust emission parameterisations used to date in martian global climate models (MGCMs) lag, understandably, behind their terrestrial counterparts in terms of sophistication. Recent developments in estimating surface roughness length over all martian terrains and in modelling atmospheric circulations at regional to local scales (less than O(100 km)) presents an opportunity to formulate an improved wind stress lifting parameterisation. We have upgraded the conventional scheme by including the spatially varying roughness length in the lifting parameterisation in a fully consistent manner (thereby correcting a possible underestimation of the true threshold level for wind stress lifting), and used a modification to account for deviations from neutral stability in the surface layer. Following these improvements, it is found that wind speeds at typical MGCM resolution never reach the lifting threshold at most gridpoints: winds fall particularly short in the southern midlatitudes, where mean roughness is large. Sub-grid scale variability, manifested in both the near-surface wind field and the surface roughness, is then considered, and is found to be a crucial means of bridging the gap between model winds and thresholds. Both forms of small-scale variability contribute to the formation of dust emission ‘hotspots’: areas within the model gridbox with particularly favourable conditions for lifting, namely a smooth surface combined with strong near-surface gusts. Such small-scale emission could in fact be particularly influential on Mars, due both to the intense positive radiative feedbacks that can drive storm growth and a strong hysteresis effect on saltation. By modelling this variability, dust lifting is predicted at the locations at which dust storms are frequently observed, including the flushing storm sources of Chryse and Utopia, and southern midlatitude areas from which larger storms tend to initiate, such as Hellas and Solis Planum. The seasonal cycle of emission, which includes a double-peaked structure in northern autumn and winter, also appears realistic. Significant increases to lifting rates are produced for any sensible choices of parameters controlling the sub-grid distributions used, but results are sensitive to the smallest scale of variability considered, which high-resolution modelling suggests should be O(1 km) or less. Use of such models in future will permit the use of a diagnosed (rather than prescribed) variable gustiness intensity, which should further enhance dust lifting in the southern hemisphere in particular.

Radiative forcing and climate metrics for ozone precursor emissions: the impact of multi-model averaging

Multi-model ensembles are frequently used to assess understanding of the response of ozone and methane lifetime to changes in emissions of ozone precursors such as NOx, VOCs (volatile organic compounds) and CO. When these ozone changes are used to calculate radiative forcing (RF) (and climate metrics such as the global warming potential (GWP) and global temperature-change potential (GTP)) there is a methodological choice, determined partly by the available computing resources, as to whether the mean ozone (and methane) concentration changes are input to the radiation code, or whether each model's ozone and methane changes are used as input, with the average RF computed from the individual model RFs. We use data from the Task Force on Hemispheric Transport of Air Pollution source–receptor global chemical transport model ensemble to assess the impact of this choice for emission changes in four regions (East Asia, Europe, North America and South Asia). We conclude that using the multi-model mean ozone and methane responses is accurate for calculating the mean RF, with differences up to 0.6% for CO, 0.7% for VOCs and 2% for NOx. Differences of up to 60% for NOx 7% for VOCs and 3% for CO are introduced into the 20 year GWP. The differences for the 20 year GTP are smaller than for the GWP for NOx, and similar for the other species. However, estimates of the standard deviation calculated from the ensemble-mean input fields (where the standard deviation at each point on the model grid is added to or subtracted from the mean field) are almost always substantially larger in RF, GWP and GTP metrics than the true standard deviation, and can be larger than the model range for short-lived ozone RF, and for the 20 and 100 year GWP and 100 year GTP. The order of averaging has most impact on the metrics for NOx, as the net values for these quantities is the residual of the sum of terms of opposing signs. For example, the standard deviation for the 20 year GWP is 2–3 times larger using the ensemble-mean fields than using the individual models to calculate the RF. The source of this effect is largely due to the construction of the input ozone fields, which overestimate the true ensemble spread. Hence, while the average of multi-model fields are normally appropriate for calculating mean RF, GWP and GTP, they are not a reliable method for calculating the uncertainty in these fields, and in general overestimate the uncertainty.

Climate versus carbon dioxide controls on biomass burning: a model analysis of the glacial–interglacial contrast

Climate controls fire regimes through its influence on the amount and types of fuel present and their dryness. CO2 concentration constrains primary production by limiting photosynthetic activity in plants. However, although fuel accumulation depends on biomass production, and hence on CO2 concentration, the quantitative relationship between atmospheric CO2 concentration and biomass burning is not well understood. Here a fire-enabled dynamic global vegetation model (the Land surface Processes and eXchanges model, LPX) is used to attribute glacial–interglacial changes in biomass burning to an increase in CO2, which would be expected to increase primary production and therefore fuel loads even in the absence of climate change, vs. climate change effects. Four general circulation models provided last glacial maximum (LGM) climate anomalies – that is, differences from the pre-industrial (PI) control climate – from the Palaeoclimate Modelling Intercomparison Project Phase~2, allowing the construction of four scenarios for LGM climate. Modelled carbon fluxes from biomass burning were corrected for the model's observed prediction biases in contemporary regional average values for biomes. With LGM climate and low CO2 (185 ppm) effects included, the modelled global flux at the LGM was in the range of 1.0–1.4 Pg C year-1, about a third less than that modelled for PI time. LGM climate with pre-industrial CO2 (280 ppm) yielded unrealistic results, with global biomass burning fluxes similar to or even greater than in the pre-industrial climate. It is inferred that a substantial part of the increase in biomass burning after the LGM must be attributed to the effect of increasing CO2 concentration on primary production and fuel load. Today, by analogy, both rising CO2 and global warming must be considered as risk factors for increasing biomass burning. Both effects need to be included in models to project future fire risks.

Impacts of uplift of northern Tibetan Plateau and formation of Asian inland deserts on regional climate and environment

Based on the geological evidence that the northern Tibetan Plateau (NTP) had an uplift of a finite magnitude since the Miocene and the major Asian inland deserts formed in the early Pliocene, a regional climate model (RegCM4.1) with a horizontal resolution of 50 km was used to explore the effects of the NTP uplift and the related aridification of inland Asia on regional climate. We designed three numerical experiments including the control experiment representing the present-day condition, the high-mountain experiment representing the early Pliocene condition with uplifted NTP but absence of the Asian inland deserts, and the low-mountain experiment representing the mid-Miocene condition with reduced topography in the NTP (by as much as 2400 m) and also absence of the deserts. Our simulation results indicated that the NTP uplift caused significant reductions in annual precipitation in a broad region of inland Asia north of the Tibetan Plateau (TP) mainly due to the enhanced rain shadow effect of the mountains and changes in the regional circulations. However, four mountainous regions located in the uplift showed significant increases in precipitation, stretching from the Pamir Plateau in the west to the Qilian Mountains in the east. These mountainous areas also experienced different changes in the rainfall seasonality with the greatest increases occurring during the respective rainy seasons, predominantly resulted from the enhanced orographically forced upwind ascents. The appearance of the major deserts in the inland Asia further reduced precipitation in the region and led to increased dust emission and deposition fluxes, while the spatial patterns of dust deposition were also changed, not only in the regions of uplift-impacted topography, but also in the downwind regions. One major contribution from this study is the comparison of the simulation results with 11 existing geological records representing the moisture conditions from Miocene to Pliocene. The comparisons revealed good matches between the simulation results and the published geological records. Therefore, we conclude that the NTP uplift and the related formation of the major deserts played a controlling role in the evolution of regional climatic conditions in a broad region in inland Asia since the Miocene.

Evidence of resilience to past climate change in Southwest Asia: early farming communities and the 9.2 and 8.2 ka events

Climate change is often cited as a major factor in social change. The so-called 8.2 ka event was one of the most pronounced and abrupt Holocene cold and arid events. The 9.2 ka event was similar, albeit of a smaller magnitude. Both events affected the Northern Hemisphere climate and caused cooling and aridification in Southwest Asia. Yet, the impacts of the 8.2 and 9.2 ka events on early farming communities in this region are not well understood. Current hypotheses for an effect of the 8.2 ka event vary from large-scale site abandonment and migration (including the Neolithisation of Europe) to continuation of occupation and local adaptation, while impacts of the 9.2 ka have not previously been systematically studied. In this paper, we present a thorough assessment of available, quality-checked radiocarbon (14C) dates for sites from Southwest Asia covering the time interval between 9500 and 7500 cal BP, which we interpret in combination with archaeological evidence. In this way, the synchronicity between changes observed in the archaeological record and the rapid climate events is tested. It is shown that there is no evidence for a simultaneous and widespread collapse, large-scale site abandonment, or migration at the time of the events. However, there are indications for local adaptation. We conclude that early farming communities were resilient to the abrupt, severe climate changes at 9250 and 8200 cal BP.

Nonlinearity of ocean heat uptake during warming and cooling in the FAMOUS climate model

Atmospheric CO2 concentration is expected to continue rising in the coming decades, but natural or artificial processes may eventually reduce it. We show that, in the FAMOUS atmosphere-ocean general circulation model, the reduction of ocean heat content as radiative forcing decreases is greater than would be expected from a linear model simulation of the response to the applied forcings. We relate this effect to the behavior of the Atlantic meridional overturning circulation (AMOC): the ocean cools more efficiently with a strong AMOC. The AMOC weakens as CO2 rises, then strengthens as CO2 declines, but temporarily overshoots its original strength. This nonlinearity comes mainly from the accumulated advection of salt into the North Atlantic, which gives the system a longer memory. This implies that changes observed in response to different CO2 scenarios or from different initial states, such as from past changes, may not be a reliable basis for making projections.

The inconstancy of the transient climate response parameter under increasing CO₂

In the Coupled Model Intercomparison Project Phase 5 (CMIP5), the model-mean increase in global mean surface air temperature T under the 1pctCO2 scenario (atmospheric CO2 increasing at 1% yr−1) during the second doubling of CO2 is 40% larger than the transient climate response (TCR), i.e. the increase in T during the first doubling. We identify four possible contributory effects. First, the surface climate system loses heat less readily into the ocean beneath as the latter warms. The model spread in the thermal coupling between the upper and deep ocean largely explains the model spread in ocean heat uptake efficiency. Second, CO2 radiative forcing may rise more rapidly than logarithmically with CO2 concentration. Third, the climate feedback parameter may decline as the CO2 concentration rises. With CMIP5 data, we cannot distinguish the second and third possibilities. Fourth, the climate feedback parameter declines as time passes or T rises; in 1pctCO2, this effect is less important than the others. We find that T projected for the end of the twenty-first century correlates more highly with T at the time of quadrupled CO2 in 1pctCO2 than with the TCR, and we suggest that the TCR may be underestimated from observed climate change.

Sea-breeze dynamics and convection initiation: the influence of convective parameterization in weather and climate model biases

There are some long-established biases in atmospheric models that originate from the representation of tropical convection. Previously, it has been difficult to separate cause and effect because errors are often the result of a number of interacting biases. Recently, researchers have gained the ability to run multiyear global climate model simulations with grid spacings small enough to switch the convective parameterization off, which permits the convection to develop explicitly. There are clear improvements to the initiation of convective storms and the diurnal cycle of rainfall in the convection-permitting simulations, which enables a new process-study approach to model bias identification. In this study, multiyear global atmosphere-only climate simulations with and without convective parameterization are undertaken with the Met Office Unified Model and are analyzed over the Maritime Continent region, where convergence from sea-breeze circulations is key for convection initiation. The analysis shows that, although the simulation with parameterized convection is able to reproduce the key rain-forming sea-breeze circulation, the parameterization is not able to respond realistically to the circulation. A feedback of errors also occurs: the convective parameterization causes rain to fall in the early morning, which cools and wets the boundary layer, reducing the land–sea temperature contrast and weakening the sea breeze. This is, however, an effect of the convective bias, rather than a cause of it. Improvements to how and when convection schemes trigger convection will improve both the timing and location of tropical rainfall and representation of sea-breeze circulations.

Trends in the CERES dataset, 2000–13: the effects of sea ice and jet shifts and comparison to climate models

We review the effects of dynamical variability on clouds and radiation in observations and models and discuss their implications for cloud feedbacks. Jet shifts produce robust meridional dipoles in upper-level clouds and longwave cloud-radiative effect (CRE), but low-level clouds, which do not simply shift with the jet, dominate the shortwave CRE. Because the effect of jet variability on CRE is relatively small, future poleward jet shifts with global warming are only a second-order contribution to the total CRE changes around the midlatitudes, suggesting a dominant role for thermodynamic effects. This implies that constraining the dynamical response is unlikely to reduce the uncertainty in extratropical cloud feedback. However, we argue that uncertainty in the cloud-radiative response does affect the atmospheric circulation response to global warming, by modulating patterns of diabatic forcing. How cloud feedbacks can affect the dynamical response to global warming is an important topic of future research.

ENSO representation in climate models: from CMIP3 to CMIP5

We analyse the ability of CMIP3 and CMIP5 coupled ocean–atmosphere general circulation models (CGCMs) to simulate the tropical Pacific mean state and El Niño-Southern Oscillation (ENSO). The CMIP5 multi-model ensemble displays an encouraging 30 % reduction of the pervasive cold bias in the western Pacific, but no quantum leap in ENSO performance compared to CMIP3. CMIP3 and CMIP5 can thus be considered as one large ensemble (CMIP3 + CMIP5) for multi-model ENSO analysis. The too large diversity in CMIP3 ENSO amplitude is however reduced by a factor of two in CMIP5 and the ENSO life cycle (location of surface temperature anomalies, seasonal phase locking) is modestly improved. Other fundamental ENSO characteristics such as central Pacific precipitation anomalies however remain poorly represented. The sea surface temperature (SST)-latent heat flux feedback is slightly improved in the CMIP5 ensemble but the wind-SST feedback is still underestimated by 20–50 % and the shortwave-SST feedbacks remain underestimated by a factor of two. The improvement in ENSO amplitudes might therefore result from error compensations. The ability of CMIP models to simulate the SST-shortwave feedback, a major source of erroneous ENSO in CGCMs, is further detailed. In observations, this feedback is strongly nonlinear because the real atmosphere switches from subsident (positive feedback) to convective (negative feedback) regimes under the effect of seasonal and interannual variations. Only one-third of CMIP3 + CMIP5 models reproduce this regime shift, with the other models remaining locked in one of the two regimes. The modelled shortwave feedback nonlinearity increases with ENSO amplitude and the amplitude of this feedback in the spring strongly relates with the models ability to simulate ENSO phase locking. In a final stage, a subset of metrics is proposed in order to synthesize the ability of each CMIP3 and CMIP5 models to simulate ENSO main characteristics and key atmospheric feedbacks.

The influence of sea ice dynamics on the climate sensitivity and memory to increased Antarctic sea ice

The study analyzes the sensitivity and memory of the Southern Hemisphere coupled climate system to increased Antarctic sea ice (ASI), taking into account the persistence of the sea ice maxima in the current climate. The mechanisms involved in restoring the climate balance under two sets of experiments, which differ in regard to their sea ice models, are discussed. The experiments are perturbed with extremes of ASI and integrated for 10 yr in a large 30-member ensemble. The results show that an ASI maximum is able to persist for ; 4 yr in the current climate, followed by a negative sea ice phase. The sea ice insulating effect during the positive phase reduces heat fluxes south of 60 8 S, while at the same time these are intensified at the sea ice edge. The increased air stability over the sea ice field strengthens the polar cell while the baroclinicity increases at midlatitudes. The mean sea level pressure is reduced (increased) over high latitudes (midlatitudes), typical of the southern annular mode (SAM) positive phase. The Southern Ocean (SO) becomes colder and fresher as the sea ice melts mainly through sea ice lateral melting, the consequence of which is an increase in the ocean stability by buoyancy and mixing changes. The climate sensitivity is triggered by the sea ice insulating process and the resulting freshwater pulse (fast response), while the climate equilibrium is restored by the heat stored in the SO subsurface layers (long response). It is concluded that the time needed for the ASI anomaly to be dissipated and/or melted is shortened by the sea ice dynamical processes.

The impact of European legislative and technology measures to reduce air pollutants on air quality, human health and climate

European air quality legislation has reduced emissions of air pollutants across Europe since the 1970s, affecting air quality, human health and regional climate. We used a coupled composition-climate model to simulate the impacts of European air quality legislation and technology measures implemented between 1970 and 2010. We contrast simulations using two emission scenarios; one with actual emissions in 2010 and the other with emissions that would have occurred in 2010 in the absence of technological improvements and end-of-pipe treatment measures in the energy, industrial and road transport sectors. European emissions of sulphur dioxide, black carbon (BC) and organic carbon in 2010 are 53%, 59% and 32% lower respectively compared to emissions that would have occurred in 2010 in the absence of legislative and technology measures. These emission reductions decreased simulated European annual mean concentrations of fine particulate matter(PM2.5) by 35%, sulphate by 44%, BC by 56% and particulate organic matter by 23%. The reduction in PM2.5 concentrations is calculated to have prevented 80 000 (37 000–116 000, at 95% confidence intervals) premature deaths annually across the European Union, resulting in a perceived financial benefit to society of US$232 billion annually (1.4% of 2010 EU GDP). The reduction in aerosol concentrations due to legislative and technology measures caused a positive change in the aerosol radiative effect at the top of atmosphere, reduced atmospheric absorption and also increased the amount of solar radiation incident at the surface over Europe. We used an energy budget approximation to estimate that these changes in the radiative balance have increased European annual mean surface temperatures and precipitation by 0.45 ± 0.11 °C and by 13 ± 0.8 mm yr−1 respectively. Our results show that the implementation of European legislation and technological improvements to reduce the emission of air pollutants has improved air quality and human health over Europe, as well as having an unintended impact on the regional radiative balance and climate.

The medieval climate anomaly and Byzantium: a review of the evidence on climatic fluctuations, economic performance and societal change

At the beginning of the Medieval Climate Anomaly, in the ninth and tenth century, the medieval eastern Roman empire, more usually known as Byzantium, was recovering from its early medieval crisis and experiencing favourable climatic conditions for the agricultural and demographic growth. Although in the Balkans and Anatolia such favourable climate conditions were prevalent during the eleventh century, parts of the imperial territories were facing significant challenges as a result of external political/military pressure. The apogee of medieval Byzantine socio-economic development, around AD 1150, coincides with a period of adverse climatic conditions for its economy, so it becomes obvious that the winter dryness and high climate variability at this time did not hinder Byzantine society and economy from achieving that level of expansion. Soon after this peak, towards the end of the twelfth century, the populations of the Byzantine world were experiencing unusual climatic conditions with marked dryness and cooler phases. The weakened Byzantine socio-political system must have contributed to the events leading to the fall of Constantinople in AD 1204 and the sack of the city. The final collapse of the Byzantine political control over western Anatolia took place half century later, thus contemporaneous with the strong cooling effect after a tropical volcanic eruption in AD 1257. We suggest that, regardless of a range of other influential factors, climate change was also an important contributing factor to the socio-economic changes that took place in Byzantium during the Medieval Climate Anomaly. Crucially, therefore, while the relatively sophisticated and complex Byzantine society was certainly influenced by climatic conditions, and while it nevertheless displayed a significant degree of resilience, external pressures as well as tensions within the Byzantine society more broadly contributed to an increasing vulnerability in respect of climate impacts. Our interdisciplinary analysis is based on all available sources of information on the climate and society of Byzantium, that is textual (documentary), archaeological, environmental, climate and climate model-based evidence about the nature and extent of climate variability in the eastern Mediterranean. The key challenge was, therefore, to assess the relative influence to be ascribed to climate variability and change on the one hand, and on the other to the anthropogenic factors in the evolution of Byzantine state and society (such as invasions, changes in international or regional market demand and patterns of production and consumption, etc.). The focus of this interdisciplinary

The impact of two coupled cirrus microphysics-radiation parameterizations on the temperature and specific humidity biases in the tropical tropopause layer in a climate model

The impact of two different coupled cirrus microphysics-radiation parameterizations on the zonally averaged temperature and humidity biases in the tropical tropopause layer (TTL) of a Met Office climate model configuration is assessed. One parameterization is based on a linear coupling between a model prognostic variable, the ice mass mixing ratio, qi, and the integral optical properties. The second is based on the integral optical properties being parameterized as functions of qi and temperature, Tc, where the mass coefficients (i.e. scattering and extinction) are parameterized as nonlinear functions of the ratio between qi and Tc. The cirrus microphysics parameterization is based on a moment estimation parameterization of the particle size distribution (PSD), which relates the mass moment (i.e. second moment if mass is proportional to size raised to the power of 2 ) of the PSD to all other PSD moments through the magnitude of the second moment and Tc. This same microphysics PSD parameterization is applied to calculate the integral optical properties used in both radiation parameterizations and, thus, ensures PSD and mass consistency between the cirrus microphysics and radiation schemes. In this paper, the temperature-non-dependent and temperature-dependent parameterizations are shown to increase and decrease the zonally averaged temperature biases in the TTL by about 1 K, respectively. The temperature-dependent radiation parameterization is further demonstrated to have a positive impact on the specific humidity biases in the TTL, as well as decreasing the shortwave and longwave biases in the cloudy radiative effect. The temperature-dependent radiation parameterization is shown to be more consistent with TTL and global radiation observations.