Meteorology and climate inspire secondary school students

As part of its National Science and Engineering Week activities in 2009 and 2010, the University of Reading organised two open days for 60 local key stage 4 pupils. The theme of both open days was ‘How do we predict weather and climate?’ Making use of the students’ familiarity with weather and climate, several concepts of relevance to secondary science were investigated. The open days also provided an opportunity for more than 30 research staff from the university to interact with the students. Feedback from students and teachers was extremely positive. This article shows how meteorological science can be used to illustrate elements of the secondary science and mathematics curricula.

Large-scale and synoptic meteorology in the south-east Pacific during the observations campaign VOCALS-REx in austral Spring 2008

We present a descriptive overview of the meteorology in the south eastern subtropical Pacific (SEP) during the VOCALS-REx intensive observations campaign which was carried out between October and November 2008. Mainly based on data from operational analyses, forecasts, reanalysis, and satellite observations, we focus on spatio-temporal scales from synoptic to planetary. A climatological context is given within which the specific conditions observed during the campaign are placed, with particular reference to the relationships between the large-scale and the regional circulations. The mean circulations associated with the diurnal breeze systems are also discussed. We then provide a summary of the day-to-day synoptic-scale circulation, air-parcel trajectories, and cloud cover in the SEP during VOCALS-REx. Three meteorologically distinct periods of time are identified and the large-scale causes for their different character are discussed. The first period was characterised by significant variability associated with synoptic-scale systems interesting the SEP; while the two subsequent phases were affected by planetary-scale disturbances with a slower evolution. The changes between initial and later periods can be partly explained from the regular march of the annual cycle, but contributions from subseasonal variability and its teleconnections were important. Across the whole of the two months under consideration we find a significant correlation between the depth of the inversion-capped marine boundary layer (MBL) and the amount of low cloud in the area of study. We discuss this correlation and argue that at least as a crude approximation a typical scaling may be applied relating MBL and cloud properties with the large-scale parameters of SSTs and tropospheric temperatures. These results are consistent with previously found empirical relationships involving lower-tropospheric stability.

Resolution of sharp fronts in the presence of model error in variational data assimilation

We show that the four-dimensional variational data assimilation method (4DVar) can be interpreted as a form of Tikhonov regularization, a very familiar method for solving ill-posed inverse problems. It is known from image restoration problems that L1-norm penalty regularization recovers sharp edges in the image more accurately than Tikhonov, or L2-norm, penalty regularization. We apply this idea from stationary inverse problems to 4DVar, a dynamical inverse problem, and give examples for an L1-norm penalty approach and a mixed total variation (TV) L1–L2-norm penalty approach. For problems with model error where sharp fronts are present and the background and observation error covariances are known, the mixed TV L1–L2-norm penalty performs better than either the L1-norm method or the strong constraint 4DVar (L2-norm)method. A strength of the mixed TV L1–L2-norm regularization is that in the case where a simplified form of the background error covariance matrix is used it produces a much more accurate analysis than 4DVar. The method thus has the potential in numerical weather prediction to overcome operational problems with poorly tuned background error covariance matrices.

An extratropical cyclone database: A tool for illustrating cyclone structure and evolution characteristics

Extratropical cyclone lifecycles have been studied extensively with the aim of understanding the dynamical mechanisms involved in their development. Previous work has often been based on subjective analysis of individual case studies. Such case studies have contributed heavily to the generation of conceptual models of extratropical cyclones that provide a framework for understanding the dynamical evolution of cyclones. These conceptual models are widely used in educational meteorology courses throughout the world to illustrate the basic structure and evolution of extratropical cyclones. This article presents a database of extratropical cyclone composites which highlight the average structure and evolution of 20 years of extratropical cyclones, as opposed to individual case studies. The composite fields are achieved by combining a database containing cyclone tracks from the ERA-Interim reanalysis (1989-2009, 6 hourly) with the full 3D ERA-Interim reanalysis fields. Vertical and horizontal composites of cyclone structure for cyclones generated in the Atlantic and Pacific regions identifying features such as the relative positions of cold, warm and occluded fronts and their associated wind and cloud patterns are shown. In addition the evolution of cyclonic flows such as the warm and cold conveyor belts and dry intrusion are illustrated. A webpage containing an archive of the composited data is freely available for educational purposes.

Problem-based learning approaches in meteorology

Problem-Based Learning, despite recent controversies about its effectiveness, is used extensively as a teaching method throughout higher education. In meteorology, there has been little attempt to incorporate Problem-Based Learning techniques into the curriculum. Motivated by a desire to enhance the reflective engagement of students within a current field course module, this project describes the implementation of two test Problem-Based Learning activities and testing and improvement using several different and complementary means of evaluation. By the end of a 2-year program of design, implementation, testing, and reflection and re-evaluation, two robust, engaging activities have been developed that provide an enhanced and diverse learning environment in the field course. The results suggest that Problem-Based Learning techniques would be a useful addition to the meteorology curriculum and suggestions for courses and activities that may benefit from this approach are included in the conclusions.

Freshwater and heat transports from global ocean synthesis

[1] An eddy-permitting ¼° global ocean reanalysis based on the Operational Met Office FOAM data assimilation system has been run for 1989–2010 forced by ERA-Interim meteorology. Freshwater and heat transports are compared with published estimates globally and in each basin, with special focus on the Atlantic. The meridional transports agree with observations within errors at most locations, but where eddies are active the transports by the mean flow are nearly always in better agreement than the total transports. Eddy transports are down gradient and are enhanced relative to a free run. They may oppose or reinforce mean transports and provide 40–50% of the total transport near midlatitude fronts, where eddies with time scales <1 month provide up to 15%. Basin-scale freshwater convergences are calculated with the Arctic/Atlantic, Indian, and Pacific oceans north of 32°S, all implying net evaporation of 0.33 ± 0.04 Sv, 0.65 ± 0.07 Sv, and 0.09 ± 0.04 Sv, respectively, within the uncertainty of observations in the Atlantic and Pacific. The Indian is more evaporative and the Southern Ocean has more precipitation (1.07 Sv). Air-sea fluxes are modified by assimilation influencing turbulent heat fluxes and evaporation. Generally, surface and assimilation fluxes together match the meridional transports, indicating that the reanalysis is close to a steady state. Atlantic overturning and gyre transports are assessed with overturning freshwater transports southward at all latitudes. At 26°N eddy transports are negligible, overturning transport is 0.67 ± 0.19 Sv southward and gyre transport is 0.44 ± 0.17 Sv northward, with divergence between 26°N and the Bering Strait of 0.13 ± 0.23 Sv over 2004–2010.

Meteorology, air quality, and health in London: the ClearfLo project

The ClearfLo project provides integrated measurements of the meteorology, composition and particulate loading of London's urban atmosphere to improve predictive capability for air quality. Air quality and heat are strong health drivers and their accurate assessment and forecast are important in densely populated urban areas. However, the sources and processes leading to high concentrations of main pollutants such as ozone, nitrogen dioxide, and fine and coarse particulate matter in complex urban areas are not fully understood, limiting our ability to forecast air quality accurately. This paper introduces the ClearfLo project's interdisciplinary approach to investigate the processes leading to poor air quality and elevated temperatures. Within ClearfLo (www.clearflo.ac.uk), a large multi-institutional project funded by the UK Natural Environment Research Council (NERC), integrated measurements of meteorology, gaseous and particulate composition/loading within London's atmosphere were undertaken to understand the processes underlying poor air quality. Long-term measurement infrastructure installed at multiple levels (street and elevated), and at urban background, kerbside and rural locations were complemented with high-resolution numerical atmospheric simulations . Combining these (measurement/modeling) enhances understanding of seasonal variations in meteorology and composition together with the controlling processes. Two intensive observation periods (winter 2012 and summer Olympics 2012) focus upon the vertical structure and evolution of the urban boundary layer, chemical controls on nitrogen dioxide and ozone production, in particular the role of volatile organic compounds, and processes controlling the evolution, size, distribution and composition of particulate matter. The paper shows that mixing heights are deeper over London than in the rural surroundings and the seasonality of the urban boundary layer evolution controls when concentrations peak. The composition also reflects the seasonality of sources such as domestic burning and biogenic emissions.

Differences between the CME fronts tracked by an expert, an automated algorithm, and the Solar Stormwatch project

Observations from the Heliospheric Imager (HI) instruments aboard the twin STEREO spacecraft have enabled the compilation of several catalogues of coronal mass ejections (CMEs), each characterizing the propagation of CMEs through the inner heliosphere. Three such catalogues are the Rutherford Appleton Laboratory (RAL)-HI event list, the Solar Stormwatch CME catalogue, and, presented here, the J-tracker catalogue. Each catalogue uses a different method to characterize the location of CME fronts in the HI images: manual identification by an expert, the statistical reduction of the manual identifications of many citizen scientists, and an automated algorithm. We provide a quantitative comparison of the differences between these catalogues and techniques, using 51 CMEs common to each catalogue. The time-elongation profiles of these CME fronts are compared, as are the estimates of the CME kinematics derived from application of three widely used single-spacecraft-fitting techniques. The J-tracker and RAL-HI profiles are most similar, while the Solar Stormwatch profiles display a small systematic offset. Evidence is presented that these differences arise because the RAL-HI and J-tracker profiles follow the sunward edge of CME density enhancements, while Solar Stormwatch profiles track closer to the antisunward (leading) edge. We demonstrate that the method used to produce the time-elongation profile typically introduces more variability into the kinematic estimates than differences between the various single-spacecraft-fitting techniques. This has implications for the repeatability and robustness of these types of analyses, arguably especially so in the context of space weather forecasting, where it could make the results strongly dependent on the methods used by the forecaster.

Decadal predictability of regional-scale peak winds over Europe using the Earth System Model of the Max-Planck-Institute for Meteorology

The predictability of high impact weather events on multiple time scales is a crucial issue both in scientific and socio-economic terms. In this study, a statistical-dynamical downscaling (SDD) approach is applied to an ensemble of decadal hindcasts obtained with the Max-Planck-Institute Earth System Model (MPI-ESM) to estimate the decadal predictability of peak wind speeds (as a proxy for gusts) over Europe. Yearly initialized decadal ensemble simulations with ten members are investigated for the period 1979–2005. The SDD approach is trained with COSMO-CLM regional climate model simulations and ERA-Interim reanalysis data and applied to the MPI-ESM hindcasts. The simulations for the period 1990–1993, which was characterized by several windstorm clusters, are analyzed in detail. The anomalies of the 95 % peak wind quantile of the MPI-ESM hindcasts are in line with the positive anomalies in reanalysis data for this period. To evaluate both the skill of the decadal predictability system and the added value of the downscaling approach, quantile verification skill scores are calculated for both the MPI-ESM large-scale wind speeds and the SDD simulated regional peak winds. Skill scores are predominantly positive for the decadal predictability system, with the highest values for short lead times and for (peak) wind speeds equal or above the 75 % quantile. This provides evidence that the analyzed hindcasts and the downscaling technique are suitable for estimating wind and peak wind speeds over Central Europe on decadal time scales. The skill scores for SDD simulated peak winds are slightly lower than those for large-scale wind speeds. This behavior can be largely attributed to the fact that peak winds are a proxy for gusts, and thus have a higher variability than wind speeds. The introduced cost-efficient downscaling technique has the advantage of estimating not only wind speeds but also estimates peak winds (a proxy for gusts) and can be easily applied to large ensemble datasets like operational decadal prediction systems.

Rossby wave propagation on potential vorticity fronts with finite width

The horizontal gradient of potential vorticity (PV) across the tropopause typically declines with lead time in global numerical weather forecasts and tends towards a steady value dependent on model resolution. This paper examines how spreading the tropopause PV contrast over a broader frontal zone affects the propagation of Rossby waves. The approach taken is to analyse Rossby waves on a PV front of finite width in a simple single-layer model. The dispersion relation for linear Rossby waves on a PV front of infinitesimal width is well known; here an approximate correction is derived for the case of a finite width front, valid in the limit that the front is narrow compared to the zonal wavelength. Broadening the front causes a decrease in both the jet speed and the ability of waves to propagate upstream. The contribution of these changes to Rossby wave phase speeds cancel at leading order. At second order the decrease in jet speed dominates, meaning phase speeds are slower on broader PV fronts. This asymptotic phase speed result is shown to hold for a wide class of single-layer dynamics with a varying range of PV inversion operators. The phase speed dependence on frontal width is verified by numerical simulations and also shown to be robust at finite wave amplitude, and estimates are made for the error in Rossby wave propagation speeds due to the PV gradient error present in numerical weather forecast models.

A numerical investigation of the atmosphere-ocean thermal contrast over the coastal upwelling region of Cabo Frio, Brazil

An one-dimensional atmospheric second order closure model, coupled to an oceanic mixed layer model, is used to investigate the short term variation of the atmospheric and oceanic boundary layers in the coastal upwelling area of Cabo Frio, Brazil (23 degrees S, 42 degrees 08`W). The numerical simulations were carried out to evaluate the impact caused by the thermal contrast between atmosphere and ocean on the vertical extent and other properties of both atmospheric and oceanic boundary layers. The numerical simulations were designed taking as reference the observations carried out during the passage of a cold front that disrupted the upwelling regime in Cabo Frio in July of 1992. The simulations indicated that in 10 hours the mechanical mixing, sustained by a constant background flow of 10 in s(-1), increases the atmospheric boundary layer in 214 in when the atmosphere is initially 2 K warmer than the ocean (positive thermal contrast observed during upwelling regime). For an atmosphere initially -2 K colder than the ocean (negative thermal contrast observed during passage of the cold front), the incipient thermal convection intensifies the mechanical mixing increasing the vertical extent of the atmospheric boundary layer in 360 in. The vertical evolution of the atmospheric boundary layer is consistent with the observations carried out in Cabo Frio during upwelling condition. When the upwelling is disrupted, the discrepancy between the simulated and observed atmospheric boundary layer heights in Cabo Frio during July of 1992 increases considerably. During the period of 10 hours, the simulated oceanic mixed layer deepens 2 in and 5.4 in for positive and negative thermal contrasts of 2 K and -2 K, respectively. In the latter case, the larger vertical extent of the oceanic mixed layer is due to the presence of thermal convection in the atmospheric boundary layer, which in turn is associated to the absence of upwelling caused by the passage of cold fronts in Cabo Frio.