Life-cycle simulations of shallow frontal waves and the impact of deformation strain

The life-cycle of shallow frontal waves and the impact of deformation strain on their development is investigated using the idealised version of the Met Office non-hydrostatic Unified Model which includes the same physics and dynamics as the operational forecast model. Frontal wave development occurs in two stages; first, a deformation strain is applied to a front and a positive potential vorticity (PV) strip forms, generated by latent heat release in the frontal updraft; second, as the deformation strain is reduced the PV strip breaks up into individual anomalies. The circulations associated with the PV anomalies cause shallow frontal waves to form. The structure of the simulated frontal waves is consistent with the conceptual model of a frontal cyclone. Deeper frontal waves are simulated if the stability of the atmosphere is reduced. Deformation strain rates of different strengths are applied to the PV strip to determine whether a deformation strain threshold exists above which frontal wave development is suppressed. An objective method of frontal wave activity is defined and frontal wave development was found to be suppressed by deformation strain rates $\ge 0.4\times10^{-5}\mbox{s}^{-1}$. This value compares well with observed deformation strain rate thresholds and the analytical solution for the minimum deformation strain rate needed to suppress barotropic frontal wave development. The deformation strain rate threshold is dependent on the strength of the PV strip with strong PV strips able to overcome stronger deformation strain rates (leading to frontal wave development) than weaker PV strips.

The role of the basic state in the ENSO-monsoon relationship and implications for predictability

The impact of systematic model errors on a coupled simulation of the Asian Summer monsoon and its interannual variability is studied. Although the mean monsoon climate is reasonably well captured, systematic errors in the equatorial Pacific mean that the monsoon-ENSO teleconnection is rather poorly represented in the GCM. A system of ocean-surface heat flux adjustments is implemented in the tropical Pacific and Indian Oceans in order to reduce the systematic biases. In this version of the GCM, the monsoon-ENSO teleconnection is better simulated, particularly the lag-lead relationships in which weak monsoons precede the peak of El Nino. In part this is related to changes in the characteristics of El Nino, which has a more realistic evolution in its developing phase. A stronger ENSO amplitude in the new model version also feeds back to further strengthen the teleconnection. These results have important implications for the use of coupled models for seasonal prediction of systems such as the monsoon, and suggest that some form of flux correction may have significant benefits where model systematic error compromises important teleconnections and modes of interannual variability.

An objective climatology of the dynamical forcing of polar lows in the Nordic seas

A new objective climatology of polar lows in the Nordic (Norwegian and Barents) seas has been derived from a database of diagnostics of objectively identified cyclones spanning the period January 2000 to April 2004. There are two distinct parts to this study: the development of the objective climatology and a characterization of the dynamical forcing of the polar lows identified. Polar lows are an intense subset of polar mesocyclones. Polar mesocyclones are distinguished from other cyclones in the database as those that occur in cold air outbreaks over the open ocean. The difference between the wet-bulb potential temperature at 700 hPa and the sea surface temperature (SST) is found to be an effective discriminator between the atmospheric conditions associated with polar lows and other cyclones in the Nordic seas. A verification study shows that the objective identification method is reliable in the Nordic seas region. After demonstrating success at identifying polar lows using the above method, the dynamical forcing of the polar lows in the Nordic seas is characterized. Diagnostics of the ratio of mid-level vertical motion attributable to quasi-geostrophic forcing from upper and lower levels (U/L ratio) are used to determine the prevalence of a recently proposed category of extratropical cyclogenesis, type C, for which latent heat release is crucial to development. Thirty-one percent of the objectively identified polar low events (36 from 115) exceeded the U/L ratio of 4.0, previously identified as a threshold for type C cyclones. There is a contrast between polar lows to the north and south of the Nordic seas. In the southern Norwegian Sea, the population of polar low events is dominated by type C cyclones. These possess strong convection and weak low-level baroclinicity. Over the Barents and northern Norwegian seas, the well-known cyclogenesis types A and B dominate. These possess stronger low-level baroclinicity and weaker convection.

How many facets are needed to represent the surface energy balance of an urban area?

We investigate the question of how many facets are needed to represent the energy balance of an urban area by developing simplified 3-, 2- and 1-facet versions of a 4-facet energy balance model of two-dimensional streets and buildings. The 3-facet model simplifies the 4-facet model by averaging over the canyon orientation, which results in similar net shortwave and longwave balances for both wall facets, but maintains the asymmetry in the heat fluxes within the street canyon. For the 2-facet model, on the assumption that the wall and road temperatures are equal, the road and wall facets can be combined mathematically into a single street-canyon facet with effective values of the heat transfer coefficient, albedo, emissivity and thermodynamic properties, without further approximation. The 1-facet model requires the additional assumption that the roof temperature is also equal to the road and wall temperatures. Idealised simulations show that the geometry and material properties of the walls and road lead to a large heat capacity of the combined street canyon, whereas the roof behaves like a flat surface with low heat capacity. This means that the magnitude of the diurnal temperature variation of the street-canyon facets are broadly similar and much smaller than the diurnal temperature variation of the roof facets. Consequently, the approximation that the street-canyon facets have similar temperatures is sound, and the road and walls can be combined into a single facet. The roof behaves very differently and a separate roof facet is required. Consequently, the 2-facet model performs similarly to the 4-facet model, while the 1-facet model does not. The models are compared with previously published observations collected in Mexico City. Although the 3- and 2-facet models perform better than the 1-facet model, the present models are unable to represent the phase of the sensible heat flux. This result is consistent with previous model comparisons, and we argue that this feature of the data cannot be produced by a single column model. We conclude that a 2-facet model is necessary, and for numerical weather prediction sufficient, to model an urban surface, and that this conclusion is robust and therefore applicable to more general geometries.

The role of mean ocean salinity in climate

We describe numerical simulations designed to elucidate the role of mean ocean salinity in climate. Using a coupled atmosphere-ocean general circulation model, we study a 100-year sensitivity experiment in which the global-mean salinity is approximately doubled from its present observed value, by adding 35 psu everywhere in the ocean. The salinity increase produces a rapid global-mean sea-surface warming of C within a few years, caused by reduced vertical mixing associated with changes in cabbeling. The warming is followed by a gradual global-mean sea-surface cooling of C within a few decades, caused by an increase in the vertical (downward) component of the isopycnal diffusive heat flux. We find no evidence of impacts on the variability of the thermohaline circulation (THC) or El Niño/Southern Oscillation (ENSO). The mean strength of the Atlantic meridional overturning is reduced by 20% and the North Atlantic Deep Water penetrates less deeply. Nevertheless, our results dispute claims that higher salinities for the world ocean have profound consequences for the thermohaline circulation. In additional experiments with doubled atmospheric carbon dioxide, we find that the amplitude and spatial pattern of the global warming signal are modified in the hypersaline ocean. In particular, the equilibrated global-mean sea-surface temperature increase caused by doubling carbon dioxide is reduced by 10%. We infer the existence of a non-linear interaction between the climate responses to modified carbon dioxide and modified salinity.

On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean-atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10-20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean's thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.

Understanding African easterly waves: a moist singular vector approach

Moist singular vectors (MSV) have been applied successfully to predicting mid-latitude storms growing in association with latent heat of condensation. Tropical cyclone sensitivity has also been assessed. Extending this approach to more general tropical weather systems here, MSVs are evaluated for understanding and predicting African easterly waves, given the importance of moist processes in their development. First results, without initial moisture perturbations, suggest MSVs may be used advantageously. Perturbations bear similar structural and energy profiles to previous idealised non-linear studies and observations. Strong sensitivities prevail in the metrics and trajectories chosen, and benefits of initial moisture perturbations should be appraised. Copyright © 2009 Royal Meteorological Society

Climate entropy budget of the HadCM3 atmosphere-ocean general circulation model and of FAMOUS, its low resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ∼50 mW m−2 K−1, a value intermediate in the range 30–70 mW m−2 K−1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (∼38 mW m−2 K−1). Another 13 mW m−2 K−1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m−2 K−1. The material entropy production within the ocean due to turbulent mixing is ∼1 mW m−2 K−1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m−2 K−1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.

Thermopile radiometer signal conditioning for surface atmospheric radiation measurements

A highly stable microvolt amplifier for use with atmospheric broadband thermopile radiometers is described. The amplifier has a nominal gain of 500, for bipolar input signals in the range +/- 10 mV from a floating source. The noise level at the input is less than 5 mu V (at 100 k Omega input impedance), permitting instantaneous diffuse solar radiation measurements to 0.5 W m(-2) resolution with 12 bit analog to digital conversion. The temperature stability of gain is better than 5 ppm/degrees C (-4 to 20 degrees C). Averaged over a decade of use, the long term drift of the amplifier gain is less than similar to 0.02%/yr. As well as radiometers measuring solar and terrestrial radiations, the amplifier has also been successfully used with low level signals from thermocouples and ground heat flux plates.

Moist singular vectors and the predictability of some high impact European cyclones

The ECMWF full-physics and dry singular vector (SV) packages, using a dry energy norm and a 1-day optimization time, are applied to four high impact European cyclones of recent years that were almost universally badly forecast in the short range. It is shown that these full-physics SVs are much more relevant to severe cyclonic development than those based on dry dynamics plus boundary layer alone. The crucial extra ingredient is the representation of large-scale latent heat release. The severe winter storms all have a long, nearly straight region of high baroclinicity stretching across the Atlantic towards Europe, with a tongue of very high moisture content on its equatorward flank. In each case some of the final-time top SV structures pick out the region of the actual storm. The initial structures were generally located in the mid- to low troposphere. Forecasts based on initial conditions perturbed by moist SVs with opposite signs and various amplitudes show the range of possible 1-day outcomes for reasonable magnitudes of forecast error. In each case one of the perturbation structures gave a forecast very much closer to the actual storm than the control forecast. Deductions are made about the predictability of high-impact extratropical cyclone events. Implications are drawn for the short-range forecast problem and suggestions made for one practicable way to approach short-range ensemble forecasting. Copyright © 2005 Royal Meteorological Society.

The influence of physical processes on extratropical singular vectors

An investigation is made of the impact of a full linearized physical (moist) parameterization package on extratropical singular vectors (SVs) using the ECMWF integrated forecasting system (IFS). Comparison is made for one particular period with a dry physical package including only vertical diffusion and surface drag. The crucial extra ingredient in the full package is found to be the large-scale latent heat release. Consistent with basic theory, its inclusion results in a shift to smaller horizontal scales and enhanced growth for the SVs. Whereas, for the dry SVs, T42 resolution is sufficient, the moist SVs require T63 to resolve their structure and growth. A 24-h optimization time appears to be appropriate for the moist SVs because of the larger growth of moist SVs compared with dry SVs. Like dry SVs, moist SVs tend to occur in regions of high baroclinicity, but their location is also influenced by the availability of moisture. The most rapidly growing SVs appear to enhance or reduce large-scale rain in regions ahead of major cold fronts. The enhancement occurs in and ahead of a cyclonic perturbation and the reduction in and ahead of an anticyclonic perturbation. Most of the moist SVs for this situation are slightly modified versions of the dry SVs. However, some occur in new locations and have particularly confined structures. The most rapidly growing SV is shown to exhibit quite linear behavior in the nonlinear model as it grows from 0.5 to 12 hPa in 1 day. For 5 times this amplitude the structure is similar but the growth is about half as the perturbation damps a potential vorticity (PV) trough or produces a cutoff, depending on its sign.

A new perspective on warming of the global oceans

Changes in ocean circulation associated with internal climate variability have a major influence on upper ocean temperatures, particularly in regions such as the North Atlantic, which are relatively well-observed and therefore over-represented in the observational record. As a result, global estimates of upper ocean heat content can give misleading estimates of the roles of natural and anthropogenic factors in causing oceanic warming. We present a method to quantify ocean warming that filters out the natural internal variability from both observations and climate simulations and better isolates externally forced air-sea heat flux changes. We obtain a much clearer picture of the drivers of oceanic temperature changes, being able to detect the effects of both anthropogenic and volcanic influences simultaneously in the observed record. Our results show that climate models are capable of capturing in remarkable detail the externally forced component of ocean temperature evolution over the last five decades.

The dynamics of a polar low assessed using potential vorticity inversion

The dynamics of a polar low are examined using a piecewise potential vorticity (PV) inversion method. In previous studies of this and other polar lows, structural evolution has been described in terms of regions of anomalous PV. In this study the relative importance of different PV anomalies and the interactions between them have been quantified using PV diagnostics. The intensification of the polar low occurred in three stages (in contrast to previous studies of polar lows that have only identified two stages). The dynamical characteristics of stages one and two are consistent with the proposed type C cyclogenesis mechanism. A diabatically-generated lower-tropospheric PV anomaly dominated intensification after initial triggering by a positive upper-level PV anomaly. A phase tilt between the upper and lower levels was maintained through retardation of the positive upper-level anomaly by the effects of latent heat release. Stage three was a period of growth dominated by wind-induced surface heat exchange (WISHE), which contributed at least 18% to the amplitude of the mature surface polar low.

Influence of vegetation on the local climate and hydrology in the tropics: sensitivity to soil parameters

Land use change with accompanying major modifications to the vegetation cover is widespread in the tropics, due to increasing demands for agricultural land, and may have significant impacts on the climate. This study investigates (1) the influence of vegetation on the local climate in the tropics; (2) how that influence varies from region to region; and (3) how the sensitivity of the local climate to vegetation, and hence land use change, depends on the hydraulic characteristics of the soil. A series of idealised experiments with the Hadley Centre atmospheric model, HadAM3, are described in which the influence of vegetation in the tropics is assessed by comparing the results of integrations with and without tropical vegetation. The sensitivity of the results to the soil characteristics is then explored by repeating the experiments with a differing, but equally valid, description of soil hydraulic parameters. The results have shown that vegetation has a significant moderating effect on the climate throughout the tropics by cooling the surface through enhanced latent heat fluxes. The influence of vegetation is, however, seasonally dependent, with much greater impacts during the dry season when the availability of surface moisture is limited. Furthermore, there are significant regional variations both in terms of the magnitude of the cooling and in the response of the precipitation. Not all regions show a feedback of vegetation on the local precipitation; this result has been related both to vegetation type and to the prevailing meteorological conditions. An important finding has been the sensitivity of the results to the specification of the soil hydraulic parameters. The introduction of more freely draining soils has changed the soil-moisture contents of the control, vegetated system and has reduced, significantly, the climate sensitivity to vegetation and by implication, land use change. Changes to the soil parameters have also had an impact on the soil hydrology and its interaction with vegetation, by altering the partitioning between fast and slow runoff processes. These results raise important questions about the representation of highly heterogeneous soil characteristics in climate models, as well as the potential influence of land use change on the soil characteristics themselves.

Modelling the partitioning of evapotranspiration in a maize-sunflower intercrop

The primary purpose of this study was to model the partitioning of evapotranspiration in a maize-sunflower intercrop at various canopy covers. The Shuttleworth-Wallace (SW) model was extended for intercropping systems to include both crop transpiration and soil evaporation and allowing interaction between the two. To test the accuracy of the extended SW model, two field experiments of maize-sunflower intercrop were conducted in 1998 and 1999. Plant transpiration and soil evaporation were measured using sap flow gauges and lysimeters, respectively. The mean prediction error (simulated minus measured values) for transpiration was zero (which indicated no overall bias in estimation error), and its accuracy was not affected by the plant growth stages, but simulated transpiration during high measured transpiration rates tended to be slightly underestimated. Overall, the predictions for daily soil evaporation were also accurate. Model estimation errors were probably due to the simplified modelling of soil water content, stomatal resistances and soil heat flux as well as due to the uncertainties in characterising the 2 micrometeorological conditions. The SW’s prediction of transpiration was most sensitive to parameters most directly related to the canopy characteristics such as the partitioning of captured solar radiation, canopy resistance, and bulk boundary layer resistance.