RACORO extended-term, aircraft observations of boundary-layer clouds

A first-of-a-kind, extended-term cloud aircraft campaign was conducted to obtain an in-situ statistical characterization of continental boundary-layer clouds needed to investigate cloud processes and refine retrieval algorithms. Coordinated by the Atmospheric Radiation Measurement (ARM) Aerial Facility (AAF), the Routine AAF Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations (RACORO) field campaign operated over the ARM Southern Great Plains (SGP) site from 22 January to 30 June 2009, collecting 260 h of data during 59 research flights. A comprehensive payload aboard the Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft measured cloud microphysics, solar and thermal radiation, physical aerosol properties, and atmospheric state parameters. Proximity to the SGP's extensive complement of surface measurements provides ancillary data that supports modeling studies and facilitates evaluation of a variety of surface retrieval algorithms. The five-month duration enabled sampling a range of conditions associated with the seasonal transition from winter to summer. Although about two-thirds of the cloud flights occurred in May and June, boundary-layer cloud fields were sampled under a variety of environmental and aerosol conditions, with about 77% of the flights occurring in cumulus and stratocumulus. Preliminary analyses illustrate use of these data to analyze cloud-aerosol relationships, characterize the horizontal variability of cloud radiative impacts, and evaluate surface-based retrievals. We discuss how an extended-term campaign requires a simplified operating paradigm that is different from that used for typical, short-term, intensive aircraft field programs.

A variational method to retrieve the extinction profile in liquid clouds using multiple field-of-view lidar

Liquid clouds play a profound role in the global radiation budget but it is difficult to remotely retrieve their vertical profile. Ordinary narrow field-of-view (FOV) lidars receive a strong return from such clouds but the information is limited to the first few optical depths. Wideangle multiple-FOV lidars can isolate radiation scattered multiple times before returning to the instrument, often penetrating much deeper into the cloud than the singly-scattered signal. These returns potentially contain information on the vertical profile of extinction coefficient, but are challenging to interpret due to the lack of a fast radiative transfer model for simulating them. This paper describes a variational algorithm that incorporates a fast forward model based on the time-dependent two-stream approximation, and its adjoint. Application of the algorithm to simulated data from a hypothetical airborne three-FOV lidar with a maximum footprint width of 600m suggests that this approach should be able to retrieve the extinction structure down to an optical depth of around 6, and total opticaldepth up to at least 35, depending on the maximum lidar FOV. The convergence behavior of Gauss-Newton and quasi-Newton optimization schemes are compared. We then present results from an application of the algorithm to observations of stratocumulus by the 8-FOV airborne “THOR” lidar. It is demonstrated how the averaging kernel can be used to diagnose the effective vertical resolution of the retrieved profile, and therefore the depth to which information on the vertical structure can be recovered. This work enables exploitation of returns from spaceborne lidar and radar subject to multiple scattering more rigorously than previously possible.

Bells and whistles of convection parameterization

The present workshop constitutes the 5th in the annual series on “Concepts for Convective Parameterizations in Large-Scale Models”. The purpose of the workshop series has been to discuss the fundamental theoretical issues of convection parameterization with a small number of European scientists. The workshop series has been funded by European Cooperation in the Field of Scientific and Technical Research (COST) Action ES0905. The theme of the workshop for the year 2012 was decided from a main conclusion of the previous workshop, which focused on the convective organization problem, seeking a means for implementing such effects into convection parameterizations (Yano et al. 2012). As it turned out, in order to discuss this implementation issue in any concrete manner, we have first to know very well the bells and whistles of convection parameterizations. This was the purpose of the 5th workshop. The title of the workshop is rather metaphorically tagged as “Bulk or Spectrum?”, because this is a typical decision we have to face at the outset of any parameterization development. The following report discusses selected issues of bells and whistles addressed during the meeting.

Controlling the computational modes of the arbitrarily structured C-grid

The arbitrarily structured C-grid, TRiSK (Thuburn, Ringler, Skamarock and Klemp, 2009, 2010) is being used in the ``Model for Prediction Across Scales'' (MPAS) and is being considered by the UK Met Office for their next dynamical core. However the hexagonal C-grid supports a branch of spurious Rossby modes which lead to erroneous grid-scale oscillations of potential vorticity (PV). It is shown how these modes can be harmlessly controlled by using upwind-biased interpolation schemes for PV. A number of existing advection schemes for PV are tested, including that used in MPAS, and none are found to give adequate results for all grids and all cases. Therefore a new scheme is proposed; continuous, linear-upwind stabilised transport (CLUST), a blend between centred and linear-upwind with the blend dependent on the flow direction with respect to the cell edge. A diagnostic of grid-scale oscillations is proposed which gives further discrimination between schemes than using potential enstrophy alone and indeed some schemes are found to destroy potential enstrophy while grid-scale oscillations grow. CLUST performs well on hexagonal-icosahedral grids and unrotated skipped latitude-longitude grids of the sphere for various shallow water test cases. Despite the computational modes, the hexagonal icosahedral grid performs well since these modes are easy and harmless to filter. As a result TRiSK appears to perform better than a spectral shallow water model.

Separating the dynamical effects of climate change and ozone depletion. Part II: Southern Hemisphere troposphere

The separate effects of ozone depleting substances (ODSs) and greenhouse gases (GHGs) on forcing circulation changes in the Southern Hemisphere extratropical troposphere are investigated using a version of the Canadian Middle Atmosphere Model (CMAM) that is coupled to an ocean. Circulation-related diagnostics include zonal wind, tropopause pressure, Hadley cell width, jet location, annular mode index, precipitation, wave drag, and eddy fluxes of momentum and heat. As expected, the tropospheric response to the ODS forcing occurs primarily in austral summer, with past (1960-99) and future (2000-99) trends of opposite sign, while the GHG forcing produces more seasonally uniform trends with the same sign in the past and future. In summer the ODS forcing dominates past trends in all diagnostics, while the two forcings contribute nearly equally but oppositely to future trends. The ODS forcing produces a past surface temperature response consisting of cooling over eastern Antarctica, and is the dominant driver of past summertime surface temperature changes when the model is constrained by observed sea surface temperatures. For all diagnostics, the response to the ODS and GHG forcings is additive: that is, the linear trend computed from the simulations using the combined forcings equals (within statistical uncertainty) the sum of the linear trends from the simulations using the two separate forcings. Space time spectra of eddy fluxes and the spatial distribution of transient wave drag are examined to assess the viability of several recently proposed mechanisms for the observed poleward shift in the tropospheric jet.

Investigating global tropical cyclone activity with a hierarchy of AGCMs: the role of model resolution

The ability to run General Circulation Models (GCMs) at ever-higher horizontal resolutions has meant that tropical cyclone simulations are increasingly credible. A hierarchy of atmosphere-only GCMs, based on the Hadley Centre Global Environmental Model (HadGEM1), with horizontal resolution increasing from approximately 270km to 60km (at 50N), is used to systematically investigate the impact of spatial resolution on the simulation of global tropical cyclone activity, independent of model formulation. Tropical cyclones are extracted from ensemble simulations and reanalyses of comparable resolutions using a feature-tracking algorithm. Resolution is critical for simulating storm intensity and convergence to observed storm intensities is not achieved with the model hierarchy. Resolution is less critical for simulating the annual number of tropical cyclones and their geographical distribution, which are well captured at resolutions of 135km or higher, particularly for Northern Hemisphere basins. Simulating the interannual variability of storm occurrence requires resolutions of 100km or higher; however, the level of skill is basin dependent. Higher resolution GCMs are increasingly able to capture the interannual variability of the large-scale environmental conditions that contribute to tropical cyclogenesis. Different environmental factors contribute to the interannual variability of tropical cyclones in the different basins: in the North Atlantic basin the vertical wind shear, potential intensity and low-level absolute vorticity are dominant, while in the North Pacific basins mid-level relative humidity and low-level absolute vorticity are dominant. Model resolution is crucial for a realistic simulation of tropical cyclone behaviour, and high-resolution GCMs are found to be valuable tools for investigating the global location and frequency of tropical cyclones.

The effect of climate change on the variability of the Northern Hemisphere stratospheric polar vortex

With extreme variability of the Arctic polar vortex being a key link for stratosphere–troposphere influences, its evolution into the twenty-first century is important for projections of changing surface climate in response to greenhouse gases. Variability of the stratospheric vortex is examined using a state-of-the-art climate model and a suite of specifically developed vortex diagnostics. The model has a fully coupled ocean and a fully resolved stratosphere. Analysis of the standard stratospheric zonal mean wind diagnostic shows no significant increase over the twenty-first century in the number of major sudden stratospheric warmings (SSWs) from its historical value of 0.7 events per decade, although the monthly distribution of SSWs does vary, with events becoming more evenly dispersed throughout the winter. However, further analyses using geometric-based vortex diagnostics show that the vortex mean state becomes weaker, and the vortex centroid is climatologically more equatorward by up to 2.5°, especially during early winter. The results using these diagnostics not only characterize the vortex structure and evolution but also emphasize the need for vortex-centric diagnostics over zonally averaged measures. Finally, vortex variability is subdivided into wave-1 (displaced) and -2 (split) components, and it is implied that vortex displacement events increase in frequency under climate change, whereas little change is observed in splitting events.

The role of atmosphere feedbacks during ENSO in the CMIP3 models. Part III: the shortwave flux feedback

Previous studies using coupled general circulation models (GCMs) suggest that the atmosphere model plays a dominant role in the modeled El Nin ̃ o–Southern Oscillation (ENSO), and that intermodel differences in the thermodynamical damping of sea surface temperatures (SSTs) are a dominant contributor to the ENSO amplitude diversity. This study presents a detailed analysis of the shortwave flux feedback (aSW) in 12 Coupled Model Intercomparison Project phase 3 (CMIP3) simulations, motivated by findings that aSW is the primary contributor to model thermodynamical damping errors. A ‘‘feedback decomposition method,’’ developed to elucidate the aSW biases, shows that all models un- derestimate the dynamical atmospheric response to SSTs in the eastern equatorial Pacific, leading to un- derestimated aSW values. Biases in the cloud response to dynamics and the shortwave interception by clouds also contribute to errors in aSW. Changes in the aSW feedback between the coupled and corresponding atmosphere-only simulations are related to changes in the mean dynamics. A large nonlinearity is found in the observed and modeled SW flux feedback, hidden when linearly cal- culating aSW. In the observations, two physical mechanisms are proposed to explain this nonlinearity: 1) a weaker subsidence response to cold SST anomalies than the ascent response to warm SST anomalies and 2) a nonlinear high-level cloud cover response to SST. The shortwave flux feedback nonlinearity tends to be underestimated by the models, linked to an underestimated nonlinearity in the dynamical response to SST. The process-based methodology presented in this study may help to correct model ENSO atmospheric biases, ultimately leading to an improved simulation of ENSO in GCMs.

The effect of wind shear and curvature on the gravity wave drag produced by a ridge

The analytical model proposed by Teixeira, Miranda, and Valente is modified to calculate the gravity wave drag exerted by a stratified flow over a 2D mountain ridge. The drag is found to be more strongly affected by the vertical variation of the background velocity than for an axisymmetric mountain. In the hydrostatic approximation, the corrections to the drag due to this effect do not depend on the detailed shape of the ridge as long as this is exactly 2D. Besides the drag, all the perturbed quantities of the flow at the surface, including the pressure, may be calculated analytically.

An analytical model of mountain wave drag for wind profiles with shear and curvature

An analytical model is developed to predict the surface drag exerted by internal gravity waves on an isolated axisymmetric mountain over which there is a stratified flow with a velocity profile that varies relatively slowly with height. The model is linear with respect to the perturbations induced by the mountain, and solves the Taylor–Goldstein equation with variable coefficients using a Wentzel–Kramers–Brillouin (WKB) approximation, formally valid for high Richardson numbers, Ri. The WKB solution is extended to a higher order than in previous studies, enabling a rigorous treatment of the effects of shear and curvature of the wind profile on the surface drag. In the hydrostatic approximation, closed formulas for the drag are derived for generic wind profiles, where the relative magnitude of the corrections to the leading-order drag (valid for a constant wind profile) does not depend on the detailed shape of the orography. The drag is found to vary proportionally to Ri21, decreasing as Ri decreases for a wind that varies linearly with height, and increasing as Ri decreases for a wind that rotates with height maintaining its magnitude. In these two cases the surface drag is predicted to be aligned with the surface wind. When one of the wind components varies linearly with height and the other is constant, the surface drag is misaligned with the surface wind, especially for relatively small Ri. All these results are shown to be in fairly good agreement with numerical simulations of mesoscale nonhydrostatic models, for high and even moderate values of Ri.

The influence of the QBO on the propagation of equatorial waves into the stratosphere

The variation of stratospheric equatorial wave characteristics with the phase of the quasi-biennial oscillation (QBO) is investigated using ECMWF Re-Analysis and NOAA outgoing longwave radiation (OLR) data. The impact of the QBO phases on the upward propagation of equatorial waves is found to be consistent and significant. In the easterly phase, there is larger Kelvin wave amplitude but smaller westward-moving mixed Rossby–gravity (WMRG) and n = 1 Rossby (R1) wave amplitude due to reduced propagation from the upper troposphere into the lower stratosphere, compared with the westerly phase. Differences in the wave amplitude exist in a deeper layer in summer than in winter, consistent with the seasonality of ambient zonal winds. There is a strong evidence of Kelvin wave amplitude peaking just below the descending westerly phase, suggesting that Kelvin waves act to bring the westerly phase downward. However, the corresponding evidence for WMRG and R1 waves is less clear. In the lower stratosphere there is zonal variation in equatorial waves. This reflects the zonal asymmetry of wave amplitudes in the upper troposphere, the source for the lower-stratospheric waves. In easterly winters the upper-tropospheric WMRG and R1 waves over the eastern Pacific region appear to be somewhat stronger compared to climatology, perhaps because of the accumulation of waves that are unable to propagate upward into the lower stratosphere. Vertical propagation features of these waves are generally consistent with theory and suggest a mixture of Doppler shifting by ambient flows and filtering. Some lower-stratosphere equatorial waves have a connection with preceding tropical convection, especially for Kelvin and R1 waves in winter.

Impact of duration thresholds on Atlantic tropical cyclone counts*

Records of Atlantic basin tropical cyclones (TCs) since the late nineteenth century indicate a very large upward trend in storm frequency. This increase in documented TCs has been previously interpreted as resulting from anthropogenic climate change. However, improvements in observing and recording practices provide an alternative interpretation for these changes: recent studies suggest that the number of potentially missed TCs is sufficient to explain a large part of the recorded increase in TC counts. This study explores the influence of another factor—TC duration—on observed changes in TC frequency, using a widely used Atlantic hurricane database (HURDAT). It is found that the occurrence of short-lived storms (duration of 2 days or less) in the database has increased dramatically, from less than one per year in the late nineteenth–early twentieth century to about five per year since about 2000, while medium- to long-lived storms have increased little, if at all. Thus, the previously documented increase in total TC frequency since the late nineteenth century in the database is primarily due to an increase in very short-lived TCs. The authors also undertake a sampling study based upon the distribution of ship observations, which provides quantitative estimates of the frequency of missed TCs, focusing just on the moderate to long-lived systems with durations exceeding 2 days in the raw HURDAT. Upon adding the estimated numbers of missed TCs, the time series of moderate to long-lived Atlantic TCs show substantial multidecadal variability, but neither time series exhibits a significant trend since the late nineteenth century, with a nominal decrease in the adjusted time series. Thus, to understand the source of the century-scale increase in Atlantic TC counts in HURDAT, one must explain the relatively monotonic increase in very short-duration storms since the late nineteenth century. While it is possible that the recorded increase in short-duration TCs represents a real climate signal, the authors consider that it is more plausible that the increase arises primarily from improvements in the quantity and quality of observations, along with enhanced interpretation techniques. These have allowed National Hurricane Center forecasters to better monitor and detect initial TC formation, and thus incorporate increasing numbers of very short-lived systems into the TC database.

Serving two masters: a reflective narrative of reconciling the tensions faced in designing doctoral research

This paper presents a reflective narrative of the process of designing a PhD project. Using the analogy of the play 'One Man, Two Guvnors' , this paper discusses the tensions a beginning researcher faces in reconciling her own vision for a project with the academic demands of doctoral-level study. Focusing on an ethnographic study of a reading group for visually-impaired people, the paper explores how the researcher's developing understanding of the considerations necessary when working with disabled people impacted on the research design. In particular, it focuses on the conflict faced by doctoral students when working in a paradigm that requires actively involving research participants, thereby relinquishing some control over the project. The aim of the paper is to provide an honest narrative that will resonate with other beginning researchers.

Rank histograms of stratified Monte-Carlo ensembles

The application of forecast ensembles to probabilistic weather prediction has spurred considerable interest in their evaluation. Such ensembles are commonly interpreted as Monte Carlo ensembles meaning that the ensemble members are perceived as random draws from a distribution. Under this interpretation, a reasonable property to ask for is statistical consistency, which demands that the ensemble members and the verification behave like draws from the same distribution. A widely used technique to assess statistical consistency of a historical dataset is the rank histogram, which uses as a criterion the number of times that the verification falls between pairs of members of the ordered ensemble. Ensemble evaluation is rendered more specific by stratification, which means that ensembles that satisfy a certain condition (e.g., a certain meteorological regime) are evaluated separately. Fundamental relationships between Monte Carlo ensembles, their rank histograms, and random sampling from the probability simplex according to the Dirichlet distribution are pointed out. Furthermore, the possible benefits and complications of ensemble stratification are discussed. The main conclusion is that a stratified Monte Carlo ensemble might appear inconsistent with the verification even though the original (unstratified) ensemble is consistent. The apparent inconsistency is merely a result of stratification. Stratified rank histograms are thus not necessarily flat. This result is demonstrated by perfect ensemble simulations and supplemented by mathematical arguments. Possible methods to avoid or remove artifacts that stratification induces in the rank histogram are suggested.

Tropical-extratropical interactions over Southern Africa: three cases of heavy summer season rainfall

The synoptic evolution of three tropical–extratropical (TE) interactions, each responsible for extreme rainfall events over southern Africa, is discussed in detail. Along with the consideration of previously studied events, common features of these heavy rainfall producing tropical temperate troughs (TTTs) over southern Africa are discussed. It is found that 2 days prior to an event, northeasterly moisture transports across Botswana, set up by the Angola low, are diverted farther south into the semiarid region of subtropical southern Africa. The TTTs reach full maturity as a TE cloud band, rooted in the central subcontinent, which is triggered by upper-level divergence along the leading edge of an upper-tropospheric westerly wave trough. Convection and rainfall within the cloud band is supported by poleward moisture transports with subtropical air rising as it leaves the continent and joins the midlatitude westerly flow. It is shown that these systems fit within a theoretical framework describing similar TE interactions found globally. Uplift forcing for the extreme rainfall of each event is investigated. Unsurprisingly, quasigeostrophic uplift is found to dominate in the midlatitudes with convective processes strongest in the subtropics. Rainfall in the semiarid interior of South Africa appears to be a result of quasigeostrophically triggered convection. Investigation of TTT formation in the context of planetary waves shows that early development is sometimes associated with previous anticyclonic wave breaking south of the subcontinent, with full maturity of TTTs occurring as a potential vorticity trough approaches the continent from the west. Sensitivity to upstream wave perturbations and effects on anticyclonic wave breaking in the South Indian Ocean are also observed.