The impact of global freshwater forcing on the thermohaline circulation: adjustment of North Atlantic convection sites in a CGCM

On the time scale of a century, the Atlantic thermohaline circulation (THC) is sensitive to the global surface salinity distribution. The advection of salinity toward the deep convection sites of the North Atlantic is one of the driving mechanisms for the THC. There is both a northward and a southward contributions. The northward salinity advection (Nsa) is related to the evaporation in the subtropics, and contributes to increased salinity in the convection sites. The southward salinity advection (Ssa) is related to the Arctic freshwater forcing and tends on the contrary to diminish salinity in the convection sites. The THC changes results from a delicate balance between these opposing mechanisms. In this study we evaluate these two effects using the IPSL-CM4 ocean-atmosphere-sea-ice coupled model (used for IPCC AR4). Perturbation experiments have been integrated for 100 years under modern insolation and trace gases. River runoff and evaporation minus precipitation are successively set to zero for the ocean during the coupling procedure. This allows the effect of processes Nsa and Ssa to be estimated with their specific time scales. It is shown that the convection sites in the North Atlantic exhibit various sensitivities to these processes. The Labrador Sea exhibits a dominant sensitivity to local forcing and Ssa with a typical time scale of 10 years, whereas the Irminger Sea is mostly sensitive to Nsa with a 15 year time scale. The GIN Seas respond to both effects with a time scale of 10 years for Ssa and 20 years for Nsa. It is concluded that, in the IPSL-CM4, the global freshwater forcing damps the THC on centennial time scales.

A model study of oceanic mechanisms affecting Equatorial Pacific sea surface temperature during the 1997-98 El Niño

In this study, the processes affecting sea surface temperature variability over the 1992–98 period, encompassing the very strong 1997–98 El Niño event, are analyzed. A tropical Pacific Ocean general circulation model, forced by a combination of weekly ERS1–2 and TAO wind stresses, and climatological heat and freshwater fluxes, is first validated against observations. The model reproduces the main features of the tropical Pacific mean state, despite a weaker than observed thermal stratification, a 0.1 m s−1 too strong (weak) South Equatorial Current (North Equatorial Countercurrent), and a slight underestimate of the Equatorial Undercurrent. Good agreement is found between the model dynamic height and TOPEX/Poseidon sea level variability, with correlation/rms differences of 0.80/4.7 cm on average in the 10°N–10°S band. The model sea surface temperature variability is a bit weak, but reproduces the main features of interannual variability during the 1992–98 period. The model compares well with the TAO current variability at the equator, with correlation/rms differences of 0.81/0.23 m s−1 for surface currents. The model therefore reproduces well the observed interannual variability, with wind stress as the only interannually varying forcing. This good agreement with observations provides confidence in the comprehensive three-dimensional circulation and thermal structure of the model. A close examination of mixed layer heat balance is thus undertaken, contrasting the mean seasonal cycle of the 1993–96 period and the 1997–98 El Niño. In the eastern Pacific, cooling by exchanges with the subsurface (vertical advection, mixing, and entrainment), the atmospheric forcing, and the eddies (mainly the tropical instability waves) are the three main contributors to the heat budget. In the central–western Pacific, the zonal advection by low-frequency currents becomes the main contributor. Westerly wind bursts (in December 1996 and March and June 1997) were found to play a decisive role in the onset of the 1997–98 El Niño. They contributed to the early warming in the eastern Pacific because the downwelling Kelvin waves that they excited diminished subsurface cooling there. But it is mainly through eastward advection of the warm pool that they generated temperature anomalies in the central Pacific. The end of El Niño can be linked to the large-scale easterly anomalies that developed in the western Pacific and spread eastward, from the end of 1997 onward. In the far-western Pacific, because of the shallower than normal thermocline, these easterlies cooled the SST by vertical processes. In the central Pacific, easterlies pushed the warm pool back to the west. In the east, they led to a shallower thermocline, which ultimately allowed subsurface cooling to resume and to quickly cool the surface layer.

The role of lateral ocean physics in the upper ocean thermal balance of a coupled ocean-atmosphere GCM

The sensitivity of the upper ocean thermal balance of an ocean-atmosphere coupled GCM to lateral ocean physics is assessed. Three 40-year simulations are performed using horizontal mixing, isopycnal mixing, and isopycnal mixing plus eddy induced advection. The thermal adjustment of the coupled system is quite different between the simulations, confirming the major role of ocean mixing on the heat balance of climate. The initial adjustment phase of the upper ocean (SST) is used to diagnose the physical mechanisms involved in each parametrisation. When the lateral ocean physics is modified, significant changes of SST are seen, mainly in the southern ocean. A heat budget of the annual mixed layer (defined as the “bowl”) shows that these changes are due to a modified heat transfer between the bowl and the ocean interior. This modified heat intake of the ocean interior is directly due to the modified lateral ocean physics. In isopycnal diffusion, this heat exchange, especially marked at mid-latitudes, is both due to an increased effective surface of diffusion and to the sign of the isopycnal gradients of temperature at the base of the bowl. As this gradient is proportional to the isopycnal gradient of salinity, this confirms the strong role of salinity in the thermal balance of the coupled system. The eddy induced advection also leads to increased exchanges between the bowl and the ocean interior. This is both due to the shape of the bowl and again to the existence of a salinity structure. The lateral ocean physics is shown to be a significant contributor to the exchanges between the diabatic and the adiabatic parts of the ocean.

Modelling convective processes during the suppressed phase of the a Madden-Julian Oscillation: Comparing single-column models with cloud resolving models

The role of convective processes in moistening the atmosphere during suppressed periods of the suppressed phase of a Madden-Julian oscillation is investigated in cloud-resolving model (CRM) simulations, and the impact of moistening on the subsequent evolution of convection is assessed as part of a Global Energy and Water Cycle Experiment Cloud System Study (GCSS) intercomparison project. The ability of single-column model (SCM) versions of a number of state-of-the-art climate and numerical weather prediction models to capture these convective processes is also evaluated. During the suppressed periods, the CRMs are found to simulate a maximum moistening around 3 km, which is associated with a predominance of shallow convection. All SCMs produce adequate amounts of shallow convection during the suppressed periods, comparable to that seen in CRMs, but the relatively drier SCMs have higher precipitation rates than the relatively wetter SCMs and CRMs. The relatively drier SCMs dry, rather than moisten, the lower troposphere below the melting level. During the transition periods, convective processes act to moisten the atmosphere above the level at which mean advection changes from moistening to drying, despite an overall drying effect for the column. The SCMs capture some essence of this moistening at upper levels. A gradual transition from shallow to deep convection is simulated by the CRMs and the wetter SCMs during the transition periods, but the onset of deep convection is delayed in the drier SCMs. This results in lower precipitation rates for these SCMs during the active periods, although much better agreement exists between the models at this time.

The propagation mechanism of a vortex on the β plane

The propagation velocity and propagation mechanism for vortices on a β plane are determined for a reduced-gravity model by integrating the momentum equations over the β plane. Isolated vortices, vortices in a background current, and initial vortex propagation from rest are studied. The propagation mechanism for isolated anticyclones as well as cyclones, which has been lacking up to now, is presented. It is shown that, to first order, the vortex moves to generate a Coriolis force on the mass anomaly of the vortex to compensate for the force on the vortex due to the variation of the Coriolis parameter. Only the mass anomaly of the vortex is of importance, because the Coriolis force due to the motion of the bulk of the layer moving with the vortex is almost fully compensated by the Coriolis force on the motion of the exterior flow. Because the mass anomaly of a cyclone is negative the force and acceleration have opposite sign. The role of dipolar structures in steadily moving vortices is discussed, and it is shown that their overall structure is fixed by the steady westward motion of the mass anomaly. Furthermore, it is shown that reduced-gravity vortices are not advected with a background flow. The reason for this behavior is that the background flow changes the ambient vorticity gradient such that the vortex obtains an extra self-propagation term that exactly cancels the advection by the background flow. Last, it is shown that a vortex initially at rest will accelerate equatorward first, after which a westward motion is generated. This result is independent of the sign of the vortex.

Synoptic controls on boundary-layer characteristics

We report the characteristics of the three-dimensional, time evolving, atmospheric boundary layer that develops beneath an idealised, dry, baroclinic weather system. The boundary-layer structure is forced by thermal advection associated with the weather system. Large positive heat fluxes behind the cold front drive a vigorous convective boundary layer, whereas moderate negative heat fluxes in the warm sector between the cold and warm fronts generate shallow, stably stratified or neutral boundary layers. The forcing of the boundary-layer structure is quantified by forming an Eulerian mass budget integrated over the depth of the boundary layer. The mass budget indicates that tropospheric air is entrained into the boundary layer both in the vicinity of the high-pressure centre, and behind the cold front. It is then transported horizontally within the boundary layer and converges towards the cyclone's warm sector, whence it is ventilated out into the troposphere. This cycling of air is likely to be important for the ventilation of pollution out of the boundary layer, and for the transformation of the properties of large-scale air masses.

Categorisation of synoptic environments associated with mesoscale convective systems over the UK

Mesoscale convective systems (MCSs) are relatively rare events in the UK but, when they do occur, can be associated with weather that is considered extreme with respect to climatology (as indicated by the number of such events that have been analysed as case studies). These case studies usually associate UK MCSs with a synoptic environment known as the Spanish plume. Here a previously published 17 year climatology of UK MCS events is extended to the present day (from 1998 to 2008) and these events classified according to the synoptic environment in which they form. Three distinct synoptic environments have been identified, here termed the classical Spanish plume, modified Spanish plume, and European easterly plume. Detailed case studies of the two latter, newly defined, environments are presented. Composites produced for each environment further reveal the differences between them. The classical Spanish plume is associated with an eastward propagating baroclinic cyclone that evolves according to idealised life cycle 1. Conditional instability is released from a warm moist plume of air advected northeastwards from Iberia that is capped by warmer, but very dry air, from the Spanish plateau. The modified Spanish plume is associated with a slowly moving mature frontal system associated with a forward tilting trough (and possibly cut-off low) at 500 hPa that evolves according to idealised life cycle 2. As in the classical Spanish plume, conditional instability is released from a warm plume of air advected northwards from Iberia. The less frequent European easterly plume is associated with an omega block centred over Scandinavia at upper levels. Conditional instability is released from a warm plume of air advected westwards across northern continental Europe. Unlike the Spanish plume environments, the European easterly plume is not a warm sector phenomena associated with a baroclinic cyclone. However, in all environments the organisation of convection is associated with the interaction of an upper-level disturbance with a low-level region of warm advection.

Controls on boundary layer ventilation: Boundary layer processes and large-scale dynamics

Midlatitude cyclones are important contributors to boundary layer ventilation. However, it is uncertain how efficient such systems are at transporting pollutants out of the boundary layer, and variations between cyclones are unexplained. In this study 15 idealized baroclinic life cycles, with a passive tracer included, are simulated to identify the relative importance of two transport processes: horizontal divergence and convergence within the boundary layer and large-scale advection by the warm conveyor belt. Results show that the amount of ventilation is insensitive to surface drag over a realistic range of values. This indicates that although boundary layer processes are necessary for ventilation they do not control the magnitude of ventilation. A diagnostic for the mass flux out of the boundary layer has been developed to identify the synoptic-scale variables controlling the strength of ascent in the warm conveyor belt. A very high level of correlation (R-2 values exceeding 0.98) is found between the diagnostic and the actual mass flux computed from the simulations. This demonstrates that the large-scale dynamics control the amount of ventilation, and the efficiency of midlatitude cyclones to ventilate the boundary layer can be estimated using the new mass flux diagnostic. We conclude that meteorological analyses, such as ERA-40, are sufficient to quantify boundary layer ventilation by the large-scale dynamics.

NOx and O3 above a tropical rainforest: an analysis with a global and box model

A cross-platform field campaign, OP3, was conducted in the state of Sabah in Malaysian Borneo between April and July of 2008. Among the suite of observations recorded, the campaign included measurements of NOx and O3 – crucial outputs of any model chemistry mechanism. We describe the measurements of these species made from both the ground site and aircraft. We then use the output from two resolutions of the chemistry transport model p-TOMCAT to illustrate the ability of a global model chemical mechanism to capture the chemistry at the rainforest site. The basic model performance is good for NOx and poor for ozone. A box model containing the same chemical mechanism is used to explore the results of the global model in more depth and make comparisons between the two. Without some parameterization of the nighttime boundary layer – free troposphere mixing (i.e. the use of a dilution parameter), the box model does not reproduce the observations, pointing to the importance of adequately representing physical processes for comparisons with surface measurements. We conclude with a discussion of box model budget calculations of chemical reaction fluxes, deposition and mixing, and compare these results to output from p-TOMCAT. These show the same chemical mechanism behaves similarly in both models, but that emissions and advection play particularly strong roles in influencing the comparison to surface measurements.

Evaluating the ability of a numerical weather prediction model to forecast tracer concentrations during ETEX 2

In this paper the meteorological processes responsible for transporting tracer during the second ETEX (European Tracer EXperiment) release are determined using the UK Met Office Unified Model (UM). The UM predicted distribution of tracer is also compared with observations from the ETEX campaign. The dominant meteorological process is a warm conveyor belt which transports large amounts of tracer away from the surface up to a height of 4 km over a 36 h period. Convection is also an important process, transporting tracer to heights of up to 8 km. Potential sources of error when using an operational numerical weather prediction model to forecast air quality are also investigated. These potential sources of error include model dynamics, model resolution and model physics. In the UM a semi-Lagrangian monotonic advection scheme is used with cubic polynomial interpolation. This can predict unrealistic negative values of tracer which are subsequently set to zero, and hence results in an overprediction of tracer concentrations. In order to conserve mass in the UM tracer simulations it was necessary to include a flux corrected transport method. Model resolution can also affect the accuracy of predicted tracer distributions. Low resolution simulations (50 km grid length) were unable to resolve a change in wind direction observed during ETEX 2, this led to an error in the transport direction and hence an error in tracer distribution. High resolution simulations (12 km grid length) captured the change in wind direction and hence produced a tracer distribution that compared better with the observations. The representation of convective mixing was found to have a large effect on the vertical transport of tracer. Turning off the convective mixing parameterisation in the UM significantly reduced the vertical transport of tracer. Finally, air quality forecasts were found to be sensitive to the timing of synoptic scale features. Errors in the position of the cold front relative to the tracer release location of only 1 h resulted in changes in the predicted tracer concentrations that were of the same order of magnitude as the absolute tracer concentrations.

A statistical mechanical approach for the computation of the climatic response to general forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. In this paper we show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied with great success to analyze the climatic response to general forcings. The crucial value of the Ruelle theory lies in the fact that it allows to compute the response of the system in terms of expectation values of explicit and computable functions of the phase space averaged over the invariant measure of the unperturbed state. We choose as test bed a classical version of the Lorenz 96 model, which, in spite of its simplicity, has a well-recognized prototypical value as it is a spatially extended one-dimensional model and presents the basic ingredients, such as dissipation, advection and the presence of an external forcing, of the actual atmosphere. We recapitulate the main aspects of the general response theory and propose some new general results. We then analyze the frequency dependence of the response of both local and global observables to perturbations having localized as well as global spatial patterns. We derive analytically several properties of the corresponding susceptibilities, such as asymptotic behavior, validity of Kramers-Kronig relations, and sum rules, whose main ingredient is the causality principle. We show that all the coefficients of the leading asymptotic expansions as well as the integral constraints can be written as linear function of parameters that describe the unperturbed properties of the system, such as its average energy. Some newly obtained empirical closure equations for such parameters allow to define such properties as an explicit function of the unperturbed forcing parameter alone for a general class of chaotic Lorenz 96 models. We then verify the theoretical predictions from the outputs of the simulations up to a high degree of precision. The theory is used to explain differences in the response of local and global observables, to define the intensive properties of the system, which do not depend on the spatial resolution of the Lorenz 96 model, and to generalize the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable for finite as well as infinite time. Finally, we propose a simple yet general methodology to study general Climate Change problems on virtually any time scale by resorting to only well selected simulations, and by taking full advantage of ensemble methods. The specific case of globally averaged surface temperature response to a general pattern of change of the CO2 concentration is discussed. We believe that the proposed approach may constitute a mathematically rigorous and practically very effective way to approach the problem of climate sensitivity, climate prediction, and climate change from a radically new perspective.

Tracer concentration profiles measured in central London as part of the REPARTEE campaign

There have been relatively few tracer experiments carried out that have looked at vertical plume spread in urban areas. In this paper we present results from two tracer (cyclic perfluorocarbon) experiments carried out in 2006 and 2007 in central London centred on the BT Tower as part of the REPARTEE (Regent’s Park and Tower Environmental Experiment) campaign. The height of the tower gives a unique opportunity to study vertical dispersion profiles and transport times in central London. Vertical gradients are contrasted with the relevant Pasquill stability classes. Estimation of lateral advection and vertical mixing times are made and compared with previous measurements. Data are then compared with a simple operational dispersion model and contrasted with data taken in central London as part of the DAPPLE campaign. This correlates dosage with non-dimensionalised distance from source. Such analyses illustrate the feasibility of the use of these empirical correlations over these prescribed distances in central London.

Moisture transport in mid-latitude cyclones

We discuss how synoptic-scale variability controls the transport of atmospheric water vapour by mid-latitude cyclones. Idealised simulations are used to investigate quantitatively what factors determine the magnitude of cyclone moisture transport. It is demonstrated that large-scale ascent on the warmconveyor belt and shallow cumulus convection are equally important for ventilating moisture from the boundary layer into the free troposphere, and that ventilated moisture can be transported large distances eastwards and polewards by the cyclone, before being returned to the surface as precipitation. The initial relative humidity is shown to have little affect on the ability of the cyclone to transport moisture, whilst the absolute temperature and meridional temperature gradient provide much stronger controls. Scaling arguments are presented to quantify the dependence of moisture transport on large-scale and boundary-layer parameters. It is shown that ventilation by shallow convection and warm-conveyor belt advection vary in the same way with changes to large-scale parameters. However, shallow convective ventilation has a much stronger dependence on boundary-layer parameters than warm-conveyor belt ventilation.

Simulations of the London urban heat island

We present simulations of London's meteorology using the Met Office Unified Model with a new, sophisticated surface energy-balance scheme to represent the urban surfaces, called MORUSES. Simulations are performed with the urban surfaces represented and with the urban surfaces replaced with grass in order to calculate the urban increment on the local meteorology. The local urban effects were moderated to some extent by the passage of an onshore flow that propagated up the Thames estuary and across the city, cooling London slightly in the afternoon. Validations of screen-level temperature show encouraging agreement to within 1–2 K, when the urban increment is up to 5 K. The model results are then used to examine factors shaping the spatial and temporal structure of London's atmospheric boundary layer. The simulations reconcile the differences in the temporal evolution of the urban heat island (UHI) shown in various studies and demonstrate that the variation of UHI with time depends strongly on the urban fetch. The UHI at a location downwind of the city centre shows a decrease in UHI during the night, while the UHI at the city centre stays constant. Finally, the UHI at a location upwind of the city centre increases continuously. The magnitude of the UHI by the time of the evening transition increases with urban fetch. The urban increments are largest at night, when the boundary layer is shallow. The boundary layer experiences continued warming after sunset, as the heat from the urban fabric is released, and a weakly convective boundary layer develops across the city. The urban land-use fraction is the dominant control on the spatial structure in the sensible heat flux and the resulting urban increment, although even the weak advection present in this case study is sufficient to advect the peak temperature increments downwind of the most built-up areas. Copyright © 2011 Royal Meteorological Society and British Crown Copyright, the Met Office

Tephra fallout hazard assessment for a Plinian eruption scenario at Volcán de Colima (Mexico)

Volcanic ash fallout associated with renewal of explosive activity at Colima, represents a serious threat to the surrounding urbanized area. Here we assess the tephra fallout hazard associated with a Plinian eruption scenario. The eruptive history of Volcán de Colima shows that Plinian eruptions occur approximately every 100 years and the last eruption, the 1913, represents the largest historic eruption of this volcano. We used the last eruption as a reference to discuss volcanic hazard and risk scenarios connected with ash fallout. Tephra fallout deposits are modeled using HAZMAP, a model based on a semi-analytical solution of the advection– diffusion–sedimentation equation for volcanic particles. Based on a statistical study of wind profiles at Colima region, we first reconstructed ash loading maps and then computed ground load probability maps for different seasons. The obtained results show that a Plinian eruptive scenario at Volcán de Colima, could seriously damage more than 10 small towns and ranches, and potentially affect big cities located at tens of kilometers from the eruptive center. The probability maps obtained are aimed to give support to the risk mitigation strategies