Simulating the interannual variability of major dust storms on Mars using variable lifting thresholds

The redistribution of a finite amount of martian surface dust during global dust storms and in the intervening periods has been modelled in a dust lifting version of the UK Mars General Circulation Model. When using a constant, uniform threshold in the model’s wind stress lifting parameterisation and assuming an unlimited supply of surface dust, multiannual simulations displayed some variability in dust lifting activity from year to year, arising from internal variability manifested in surface wind stress, but dust storms were limited in size and formed within a relatively short seasonal window. Lifting thresholds were then allowed to vary at each model gridpoint, dependent on the rates of emission or deposition of dust. This enhanced interannual variability in dust storm magnitude and timing, such that model storms covered most of the observed ranges in size and initiation date within a single multiannual simulation. Peak storm magnitude in a given year was primarily determined by the availability of surface dust at a number of key sites in the southern hemisphere. The observed global dust storm (GDS) frequency of roughly one in every 3 years was approximately reproduced, but the model failed to generate these GDSs spontaneously in the southern hemisphere, where they have typically been observed to initiate. After several years of simulation, the surface threshold field—a proxy for net change in surface dust density—showed good qualitative agreement with the observed pattern of martian surface dust cover. The model produced a net northward cross-equatorial dust mass flux, which necessitated the addition of an artificial threshold decrease rate in order to allow the continued generation of dust storms over the course of a multiannual simulation. At standard model resolution, for the southward mass flux due to cross-equatorial flushing storms to offset the northward flux due to GDSs on a timescale of ∼3 years would require an increase in the former by a factor of 3–4. Results at higher model resolution and uncertainties in dust vertical profiles mean that quasi-periodic redistribution of dust on such a timescale nevertheless appears to be a plausible explanation for the observed GDS frequency.

Heat and moisture budgets from airborne measurements and high-resolution model simulations

High-resolution simulations with a mesoscale model are performed to estimate heat and moisture budgets of a well-mixed boundary layer. The model budgets are validated against energy budgets obtained from airborne measurements over heterogeneous terrain in Western Germany. Time rate of change, vertical divergence, and horizontal advection for an atmospheric column of air are estimated. Results show that the time trend of specific humidity exhibits some deficiencies, while the potential temperature trend is matched accurately. Furthermore, the simulated turbulent surface fluxes of sensible and latent heat are comparable to the measured fluxes, leading to similar values of the vertical divergence. The analysis of different horizontal model resolutions exhibits improved surface fluxes with increased resolution, a fact attributed to a reduced aggregation effect. Scale-interaction effects could be identified: while time trends and advection are strongly influenced by mesoscale forcing, the turbulent surface fluxes are mainly controlled by microscale processes.

The effect of non-uniform radiative damping on the zonal-mean dynamics of the extratropical middle atmosphere

The effect of spatial and temporal variations in the radiative damping rate on the response to an imposed forcing or diabatic heating is examined in a zonal-mean model of the middle atmosphere. Attention is restricted to the extratropics, where a linear approach is viable. It is found that regions with weak radiative damping rates are more sensitive in terms of temperature to the remote influence of the diabatic circulation. The delay in the response in such regions can mean that ‘downward’ control is not achieved on seasonal time-scales. A seasonal variation in the radiative damping rate modulates the evolution of the response and leaves a transient-like signature in the annual mean temperature field. Several idealized examples are considered, motivated by topical questions. It is found that wave drag outside the polar vortex can significantly affect the temperatures in its interior, so that high-latitude, high-altitude gravity-wave drag is not the only mechanism for warming the southern hemisphere polar vortex. Diabatic mass transport through the 100 hPa surface is found to lag the seasonal evolution of the wave drag that drives the transport, and thus cannot be considered to be in the downward control regime. On the other hand, the seasonal variation of the radiative damping rate is found to make only a weak contribution to the annual mean temperature increase that has been observed above the ozone hole. Copyright © 2002 Royal Meteorological Society.

What processes drive the ocean heat transport?

The ocean contributes to regulating the Earth’s climate through its ability to transport heat from the equator to the poles. In this study we use long simulations of an ocean model to investigate whether the heat transport is carried primarily by wind-driven gyres or whether it is dominated by deep circulations associated with abyssal mixing and high latitude convection. The heat transport is computed as a function of temperature classes. In the Pacific and Indian ocean, the bulk of the heat transport is associated with wind-driven gyres confined to the thermocline. In the Atlantic, the thermocline gyres account for only 40% of the total heat transport. The remaining 60% is associated with a circulation reaching down to cold waters below the thermocline. Using a series of sensitivity experiments, we show that this deep heat transport is primarily set by the strength and patterns of surface winds and only secondarily by diabatic processes at high latitudes in the North Atlantic. Abyssal mixing below 2000 m has hardly any impact on ocean heat transport. A major implication is that the role of the ocean in regulating Earth’s climate strongly depends on how surface winds change across different climates in both hemispheres at low and high latitudes.

Localization of deep water formation: role of atmospheric moisture transport and geometrical constraints on ocean circulation

A series of coupled atmosphere–ocean–ice aquaplanet experiments is described in which topological constraints on ocean circulation are introduced to study the role of ocean circulation on the mean climate of the coupled system. It is imagined that the earth is completely covered by an ocean of uniform depth except for the presence or absence of narrow barriers that extend from the bottom of the ocean to the sea surface. The following four configurations are described: Aqua (no land), Ridge (one barrier extends from pole to pole), Drake (one barrier extends from the North Pole to 35°S), and DDrake (two such barriers are set 90° apart and join at the North Pole, separating the ocean into a large basin and a small basin, connected to the south). On moving from Aqua to Ridge to Drake to DDrake, the energy transports in the equilibrium solutions become increasingly “realistic,” culminating in DDrake, which has an uncanny resemblance to the present climate. Remarkably, the zonal-average climates of Drake and DDrake are strikingly similar, exhibiting almost identical heat and freshwater transports, and meridional overturning circulations. However, Drake and DDrake differ dramatically in their regional climates. The small and large basins of DDrake exhibit distinctive Atlantic-like and Pacific-like characteristics, respectively: the small basin is warmer, saltier, and denser at the surface than the large basin, and is the main site of deep water formation with a deep overturning circulation and strong northward ocean heat transport. A sensitivity experiment with DDrake demonstrates that the salinity contrast between the two basins, and hence the localization of deep convection, results from a deficit of precipitation, rather than an excess of evaporation, over the small basin. It is argued that the width of the small basin relative to the zonal fetch of atmospheric precipitation is the key to understanding this salinity contrast. Finally, it is argued that many gross features of the present climate are consequences of two topological asymmetries that have profound effects on ocean circulation: a meridional asymmetry (circumpolar flow in the Southern Hemisphere; blocked flow in the Northern Hemisphere) and a zonal asymmetry (a small basin and a large basin).

Ocean heat transport and water vapor greenhouse in a warm equable climate: a new look at the low gradient paradox

The authors study the role of ocean heat transport (OHT) in the maintenance of a warm, equable, ice-free climate. An ensemble of idealized aquaplanet GCM calculations is used to assess the equilibrium sensitivity of global mean surface temperature and its equator-to-pole gradient (ΔT) to variations in OHT, prescribed through a simple analytical formula representing export out of the tropics and poleward convergence. Low-latitude OHT warms the mid- to high latitudes without cooling the tropics; increases by 1°C and ΔT decreases by 2.6°C for every 0.5-PW increase in OHT across 30° latitude. This warming is relatively insensitive to the detailed meridional structure of OHT. It occurs in spite of near-perfect atmospheric compensation of large imposed variations in OHT: the total poleward heat transport is nearly fixed. The warming results from a convective adjustment of the extratropical troposphere. Increased OHT drives a shift from large-scale to convective precipitation in the midlatitude storm tracks. Warming arises primarily from enhanced greenhouse trapping associated with convective moistening of the upper troposphere. Warming extends to the poles by atmospheric processes even in the absence of high-latitude OHT. A new conceptual model for equable climates is proposed, in which OHT plays a key role by driving enhanced deep convection in the midlatitude storm tracks. In this view, the climatic impact of OHT depends on its effects on the greenhouse properties of the atmosphere, rather than its ability to increase the total poleward energy transport.

The vertical structure of ocean heat transport

One of the most important contributions the ocean makes to Earth's climate is through its poleward heat transport: about 1.5 PW or more than 30% of that accomplished by the ocean-atmosphere system (Trenberth and Caron, 2001). Recently, concern has arisen over whether global warming could affect this heat transport (Watson et al., 2001), for example, reducing high latitude convection and triggering a collapse of the deep overturning circulation (Rahmstorf, 1995). While the consequences of abrupt changes in oceanic circulation should be of concern, we argue that the attention devoted to deep circulations is disproportionate to their role in heat transport. For this purpose, we introduce a heat function which identifies the contribution to the heat transport by different components of the oceanic circulation. A new view of the ocean emerges in which a shallow surface intensified circulation dominates the poleward heat transport.

The role of anomalous SST and surface fluxes over the Southeastern North Atlantic in the explosive development of windstorm Xynthia

In late February 2010 the extraordinary windstorm Xynthia crossed over Southwestern and Central Europe and caused severe damage, affecting particularly the Spanish and French Atlantic coasts. The storm was embedded in uncommon large-scale atmospheric and boundary conditions prior to and during its development, namely enhanced sea surface temperatures (SST) within the low-level entrainment zone of air masses, an unusual southerly position of the polar jet stream, and a remarkable split jet structure in the upper troposphere. To analyse the processes that led to the rapid intensification of this exceptional storm originating close to the subtropics (30°N), the sensitivity of the cyclone intensification to latent heat release is determined using the regional climate model COSMO-CLM forced with ERA-Interim data. A control simulation with observed SST shows that moist and warm air masses originating from the subtropical North Atlantic were involved in the cyclogenesis process and led to the formation of a vertical tower with high values of potential vorticity (PV). Sensitivity studies with reduced SST or increased laminar boundary roughness for heat led to reduced surface latent heat fluxes. This induced both a weaker and partly retarded development of the cyclone and a weakening of the PV-tower together with reduced diabatic heating rates, particularly at lower and mid levels. We infer that diabatic processes played a crucial role during the phase of rapid deepening of Xynthia and thus to its intensity over the Southeastern North Atlantic. We suggest that windstorms like Xynthia may occur more frequently under future climate conditions due to the warming SSTs and potentially enhanced latent heat release, thus increasing the windstorm risk for Southwestern Europe.

Effects of non-uniform crosswind fields on scintillometry measurements

The effects of a non-uniform wind field along the path of a scintillometer are investigated. Theoretical spectra are calculated for a range of scenarios where the crosswind varies in space or time and compared to the ‘ideal’ spectrum based on a constant uniform crosswind. It is verified that the refractive-index structure parameter relation with the scintillometer signal remains valid and invariant for both spatially and temporally-varying crosswinds. However, the spectral shape may change significantly preventing accurate estimation of the crosswind speed from the peak of the frequency spectrum and retrieval of the structure parameter from the plateau of the power spectrum. On comparison with experimental data, non-uniform crosswind conditions could be responsible for previously unexplained features sometimes seen in observed spectra. By accounting for the distribution of crosswind, theoretical spectra can be generated that closely replicate the observations, leading to a better understanding of the measurements. Spatial variability of wind speeds should be expected for paths other than those that are parallel to the surface and over flat, homogenous areas, whilst fluctuations in time are important for all sites.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

The influence of Antarctica on the momentum budget of the southern extratropics

The antarctic plateau acts as a strong heat sink for the global climate, cooling the atmosphere and radiating energy to space. A cold dense atmospheric boundary layer is formed. Strong surface winds are formed as the boundary layer drains off the plateau. These drainage winds and the eddy fluxes necessary to maintain them are analysed in a general circulation model (GCM). The drainage flow is well represented in the GCM. The associated mean meridional circulation is analysed in isentropic coordinates. The momentum budget over Antarctica reveals a balance between the Eliassen-Palm flux convergence and the Coriolis torque exerted by the mean meridional mass flux. Both vertical and horizontal components of the Eliassen-Palm flux contribute, the vertical component being the greater.

Results of PMIP2 coupled simulations of the mid-Holocene and Last Glacial Maximum, Part 2: feedbacks with emphasis on the location of the ITCZ and mid- and high latitudes heat budget

A set of coupled ocean-atmosphere(-vegetation) simulations using state of the art climate models is now available for the Last Glacial Maximum (LGM) and the Mid-Holocene (MH) through the second phase of the Paleoclimate Modeling Intercomparison Project (PMIP2). Here we quantify the latitudinal shift of the location of the Intertropical Convergence Zone (ITCZ) in the tropical regions during boreal summer and the change in precipitation in the northern part of the ITCZ. For both periods the shift is more pronounced over the continents and East Asia. The maritime continent is the region where the largest spread is found between models. We also clearly establish that the larger the increase in the meridional temperature gradient in the tropical Atlantic during summer at the MH, the larger the change in precipitation over West Africa. The vegetation feedback is however not as large as found in previous studies, probably due to model differences in the control simulation. Finally, we show that the feedback from snow and sea-ice at mid and high latitudes contributes for half of the cooling in the Northern Hemisphere for the LGM, with the remaining being achieved by the reduced CO2 and water vapour in the atmosphere. For the MH the snow and albedo feedbacks strengthen the spring cooling and enhance the boreal summer warming, whereas water vapour reinforces the late summer warming. These feedbacks are modest in the Southern Hemisphere. For the LGM most of the surface cooling is due to CO2 and water vapour.

Multi-scale sensible heat fluxes in the urban environment from large aperture scintillometry and eddy covariance

Sensible heat fluxes (QH) are determined using scintillometry and eddy covariance over a suburban area. Two large aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 km and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2-4 km2, whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5-10 km2. These large-scale heat fluxes are compared with local-scale eddy covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly QH is related to the meteorological conditions. At shorter time scales the response of QH to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation (Q*) coincide with greater differences between heat fluxes. For clear days QH lags Q*, thus the ratio of QH to Q* increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

Effect of cattle urine addition on the surface emissions and subsurface concentrations of greenhouse gases in a UK peat grassland

Grazing systems represent a substantial percentage of the global anthropogenic flux of nitrous oxide (N2O) as a result of nitrogen addition to the soil. The pool of available carbon that is added to the soil from livestock excreta also provides substrate for the production of carbon dioxide (CO2) and methane (CH4) by soil microorganisms. A study into the production and emission of CO2, CH4 and N2O from cattle urine amended pasture was carried out on the Somerset Levels and Moors, UK over a three-month period. Urine-amended plots (50 g N m−2) were compared to control plots to which only water (12 mg N m−2) was applied. CO2 emission peaked at 5200 mg CO2 m−2 d−1 directly after application. CH4 flux decreased to −2000 μg CH4 m−2 d−1 two days after application; however, net CH4 flux was positive from urine treated plots and negative from control plots. N2O emission peaked at 88 mg N2O m−2 d−1 12 days after application. Subsurface CH4 and N2O concentrations were higher in the urine treated plots than the controls. There was no effect of treatment on subsurface CO2 concentrations. Subsurface N2O peaked at 500 ppm 12 days after and 1200 ppm 56 days after application. Subsurface NO3− concentration peaked at approximately 300 mg N kg dry soil−1 12 days after application. Results indicate that denitrification is the key driver for N2O release in peatlands and that this production is strongly related to rainfall events and water-table movement. N2O production at depth continued long after emissions were detected at the surface. Further understanding of the interaction between subsurface gas concentrations, surface emissions and soil hydrological conditions is required to successfully predict greenhouse gas production and emission.