Vertical propagation of information in a middle atmosphere data assimilation system by gravity-wave drag feedbacks

The mesospheric response to the 2002 Antarctic Stratospheric Sudden Warming (SSW) is analysed using the Canadian Middle Atmosphere Model Data Assimilation System (CMAM-DAS), where it represents a vertical propagation of information from the observations into the data-free mesosphere. The CMAM-DAS simulates a cooling in the lowest part of the mesosphere which is accomplished by resolved motions, but which is extended to the mid- to upper mesosphere by the response of the model's non-orographic gravity-wave drag parameterization to the change in zonal winds. The basic mechanism is that elucidated by Holton consisting of a net eastward wave-drag anomaly in the mesosphere during the SSW, although in this case there is a net upwelling in the polar mesosphere. Since the zonal-mean mesospheric response is shown to be predictable, this demonstrates that variations in the mesospheric state can be slaved to the lower atmosphere through gravity-wave drag.

Tropical cyclone climatology in a 10-km global atmospheric GCM: toward weather-resolving climate modeling

Northern Hemisphere tropical cyclone (TC) activity is investigated in multiyear global climate simulations with theECMWFIntegrated Forecast System (IFS) at 10-km resolution forced by the observed records of sea surface temperature and sea ice. The results are compared to analogous simulationswith the 16-, 39-, and 125-km versions of the model as well as observations. In the North Atlantic, mean TC frequency in the 10-km model is comparable to the observed frequency, whereas it is too low in the other versions. While spatial distributions of the genesis and track densities improve systematically with increasing resolution, the 10-km model displays qualitatively more realistic simulation of the track density in the western subtropical North Atlantic. In the North Pacific, the TC count tends to be too high in thewest and too low in the east for all resolutions. These model errors appear to be associated with the errors in the large-scale environmental conditions that are fairly similar in this region for all model versions. The largest benefits of the 10-km simulation are the dramatically more accurate representation of the TC intensity distribution and the structure of the most intense storms. The model can generate a supertyphoon with a maximum surface wind speed of 68.4 m s21. The life cycle of an intense TC comprises intensity fluctuations that occur in apparent connection with the variations of the eyewall/rainband structure. These findings suggest that a hydrostatic model with cumulus parameterization and of high enough resolution could be efficiently used to simulate the TC intensity response (and the associated structural changes) to future climate change.

Ice formation and development in aged, wintertime cumulus over the UK: observations and modelling

In situ high resolution aircraft measurements of cloud microphysical properties were made in coordination with ground based remote sensing observations of a line of small cumulus clouds, using Radar and Lidar, as part of the Aerosol Properties, PRocesses And InfluenceS on the Earth's climate (APPRAISE) project. A narrow but extensive line (~100 km long) of shallow convective clouds over the southern UK was studied. Cloud top temperatures were observed to be higher than −8 °C, but the clouds were seen to consist of supercooled droplets and varying concentrations of ice particles. No ice particles were observed to be falling into the cloud tops from above. Current parameterisations of ice nuclei (IN) numbers predict too few particles will be active as ice nuclei to account for ice particle concentrations at the observed, near cloud top, temperatures (−7.5 °C). The role of mineral dust particles, consistent with concentrations observed near the surface, acting as high temperature IN is considered important in this case. It was found that very high concentrations of ice particles (up to 100 L−1) could be produced by secondary ice particle production providing the observed small amount of primary ice (about 0.01 L−1) was present to initiate it. This emphasises the need to understand primary ice formation in slightly supercooled clouds. It is shown using simple calculations that the Hallett-Mossop process (HM) is the likely source of the secondary ice. Model simulations of the case study were performed with the Aerosol Cloud and Precipitation Interactions Model (ACPIM). These parcel model investigations confirmed the HM process to be a very important mechanism for producing the observed high ice concentrations. A key step in generating the high concentrations was the process of collision and coalescence of rain drops, which once formed fell rapidly through the cloud, collecting ice particles which caused them to freeze and form instant large riming particles. The broadening of the droplet size-distribution by collision-coalescence was, therefore, a vital step in this process as this was required to generate the large number of ice crystals observed in the time available. Simulations were also performed with the WRF (Weather, Research and Forecasting) model. The results showed that while HM does act to increase the mass and number concentration of ice particles in these model simulations it was not found to be critical for the formation of precipitation. However, the WRF simulations produced a cloud top that was too cold and this, combined with the assumption of continual replenishing of ice nuclei removed by ice crystal formation, resulted in too many ice crystals forming by primary nucleation compared to the observations and parcel modelling.

Under what conditions does embedded convection enhance orographic precipitation?

Idealised convection-permitting simulations are used to quantify the impact of embedded convection on the precipitation generated by moist flow over midlatitude mountain ridges. A broad range of mountain dimensions and moist stabilities are considered to encompass a spectrum of physically plausible flows. The simulations reveal that convection only enhances orographic precipitation in cap clouds that are otherwise unable to efficiently convert cloud condensate into precipitate. For tall and wide mountains (e.g. the Washington Cascades or the southern Andes), precipitate forms efficiently through vapour deposition and collection, even in the absence of embedded convection. When embedded convection develops in such clouds, it produces competing effects (enhanced condensation in updraughts and enhanced evaporation through turbulent mixing and compensating subsidence) that cancel to yield little net change in precipitation. By contrast, convection strongly enhances precipitation over short and narrow mountains (e.g. the UK Pennines or the Oregon Coastal Range) where precipitation formation is otherwise highly inefficient. Although cancellation between increased condensation and evaporation still occurs, the enhanced precipitation formation within the convective updraughts leads to a net increase in precipitation efficiency. The simulations are physically interpreted through non-dimensional diagnostics and relevant time-scales that govern advective, microphysical, and convective processes.

A global perspective on Langmuir turbulence in the ocean surface boundary layer.

The turbulent mixing in thin ocean surface boundary layers (OSBL), which occupy the upper 100 m or so of the ocean, control the exchange of heat and trace gases between the atmosphere and ocean. Here we show that current parameterizations of this turbulent mixing lead to systematic and substantial errors in the depth of the OSBL in global climate models, which then leads to biases in sea surface temperature. One reason, we argue, is that current parameterizations are missing key surface-wave processes that force Langmuir turbulence that deepens the OSBL more rapidly than steady wind forcing. Scaling arguments are presented to identify two dimensionless parameters that measure the importance of wave forcing against wind forcing, and against buoyancy forcing. A global perspective on the occurrence of waveforced turbulence is developed using re-analysis data to compute these parameters globally. The diagnostic study developed here suggests that turbulent energy available for mixing the OSBL is under-estimated without forcing by surface waves. Wave-forcing and hence Langmuir turbulence could be important over wide areas of the ocean and in all seasons in the Southern Ocean. We conclude that surfacewave- forced Langmuir turbulence is an important process in the OSBL that requires parameterization.

Precipitation distributions for explicit versus parameterized convection in a large-domain high-resolution tropical case study

Global climate and weather models tend to produce rainfall that is too light and too regular over the tropical ocean. This is likely because of convective parametrizations, but the problem is not well understood. Here, distributions of precipitation rates are analyzed for high-resolution UK Met Office Unified Model simulations of a 10 day case study over a large tropical domain (∼20°S–20°N and 42°E–180°E). Simulations with 12 km grid length and parametrized convection have too many occurrences of light rain and too few of heavier rain when interpolated onto a 1° grid and compared with Tropical Rainfall Measuring Mission (TRMM) data. In fact, this version of the model appears to have a preferred scale of rainfall around 0.4 mm h−1 (10 mm day−1), unlike observations of tropical rainfall. On the other hand, 4 km grid length simulations with explicit convection produce distributions much more similar to TRMM observations. The apparent preferred scale at lighter rain rates seems to be a feature of the convective parametrization rather than the coarse resolution, as demonstrated by results from 12 km simulations with explicit convection and 40 km simulations with parametrized convection. In fact, coarser resolution models with explicit convection tend to have even more heavy rain than observed. Implications for models using convective parametrizations, including interactions of heating and moistening profiles with larger scales, are discussed. One important implication is that the explicit convection 4 km model has temperature and moisture tendencies that favour transitions in the convective regime. Also, the 12 km parametrized convection model produces a more stable temperature profile at its extreme high-precipitation range, which may reduce the chance of very heavy rainfall. Further study is needed to determine whether unrealistic precipitation distributions are due to some fundamental limitation of convective parametrizations or whether parametrizations can be improved, in order to better simulate these distributions.

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.

The effects of explicit versus parameterized convection on the MJO in a large-domain high-resolution tropical case study. Part I: Characterization of large-scale organization and propagation

High-resolution simulations over a large tropical domain (∼20◦S–20◦N and 42◦E–180◦E) using both explicit and parameterized convection are analyzed and compared to observations during a 10-day case study of an active Madden-Julian Oscillation (MJO) event. The parameterized convection model simulations at both 40 km and 12 km grid spacing have a very weak MJO signal and little eastward propagation. A 4 km explicit convection simulation using Smagorinsky subgrid mixing in the vertical and horizontal dimensions exhibits the best MJO strength and propagation speed. 12 km explicit convection simulations also perform much better than the 12 km parameterized convection run, suggesting that the convection scheme, rather than horizontal resolution, is key for these MJO simulations. Interestingly, a 4 km explicit convection simulation using the conventional boundary layer scheme for vertical subgrid mixing (but still using Smagorinsky horizontal mixing) completely loses the large-scale MJO organization, showing that relatively high resolution with explicit convection does not guarantee a good MJO simulation. Models with a good MJO representation have a more realistic relationship between lower-free-tropospheric moisture and precipitation, supporting the idea that moisture-convection feedback is a key process for MJO propagation. There is also increased generation of available potential energy and conversion of that energy into kinetic energy in models with a more realistic MJO, which is related to larger zonal variance in convective heating and vertical velocity, larger zonal temperature variance around 200 hPa, and larger correlations between temperature and ascent (and between temperature and diabatic heating) between 500–400 hPa.

Impact of melt ponds on Arctic sea ice simulations from 1990 to 2007

The extent and thickness of the Arctic sea ice cover has decreased dramatically in the past few decades with minima in sea ice extent in September 2007 and 2011 and climate models did not predict this decline. One of the processes poorly represented in sea ice models is the formation and evolution of melt ponds. Melt ponds form on Arctic sea ice during the melting season and their presence affects the heat and mass balances of the ice cover, mainly by decreasing the value of the surface albedo by up to 20%. We have developed a melt pond model suitable for forecasting the presence of melt ponds based on sea ice conditions. This model has been incorporated into the Los Alamos CICE sea ice model, the sea ice component of several IPCC climate models. Simulations for the period 1990 to 2007 are in good agreement with observed ice concentration. In comparison to simulations without ponds, the September ice volume is nearly 40% lower. Sensitivity studies within the range of uncertainty reveal that, of the parameters pertinent to the present melt pond parameterization and for our prescribed atmospheric and oceanic forcing, variations of optical properties and the amount of snowfall have the strongest impact on sea ice extent and volume. We conclude that melt ponds will play an increasingly important role in the melting of the Arctic ice cover and their incorporation in the sea ice component of Global Circulation Models is essential for accurate future sea ice forecasts.

Climate warming–related strengthening of the tropical hydrological cycle

The authors estimate climate warming–related twenty-first-century changes of moisture transports from the descending into the ascending regions in the tropics. Unlike previous studies that employ time and space averaging, here homogeneous high horizontal and vertical resolution data from an Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4) climate model are used. This allows for estimating changes in much greater detail (e.g., the estimation of the distribution of ascending and descending regions, changes in the vertical profile, and separating changes of the inward and outward transports). Low-level inward and midlevel outward moisture transports of the convective regions in the tropics are found to increase in a simulated anthropogenically warmed climate as compared to a simulated twentieth-century atmosphere, indicating an intensification of the hydrological cycle. Since an increase of absolute inward transport exceeds the absolute increase of outward transport, the resulting budget is positive, meaning that more water is projected to converge in the moist tropics. The intensification is found mainly to be due to the higher amount of water in the atmosphere, while the contribution of weakening wind counteracts this response marginally. In addition the changing statistical properties of the vertical profile of the moisture transport are investigated and the importance of the substantial outflow of moisture from the moist tropics at midlevels is demonstrated.

The effect of regional changes in anthropogenic aerosols on rainfall of the East Asian summer monsoon

The response of East Asian Summer Monsoon (EASM) precipitation to long term changes in regional anthropogenic aerosols (sulphate and black carbon) is explored in an atmospheric general circulation model, the atmospheric component of the UK High-Resolution Global Environment Model v1.2 (HiGAM). Separately, sulphur dioxide (SO2) and black carbon (BC) emissions in 1950 and 2000 over East Asia are used to drive model simulations, while emissions are kept constant at year 2000 level outside this region. The response of the EASM is examined by comparing simulations driven by aerosol emissions representative of 1950 and 2000. The aerosol radiative effects are also determined using an off-line radiative transfer model. During June, July and August, the EASM was not significantly changed as either SO2 or BC emissions increased from 1950 to 2000 levels. However, in September, precipitation is significantly decreased by 26.4% for sulphate aerosol and 14.6% for black carbon when emissions are at the 2000 level. Over 80% of the decrease is attributed to changes in convective precipitation. The cooler land surface temperature over China in September (0.8 °C for sulphate and 0.5 °C for black carbon) due to increased aerosols reduces the surface thermal contrast that supports the EASM circulation. However, mechanisms causing the surface temperature decrease in September are different between sulphate and BC experiments. In the sulphate experiment, the sulphate direct and the 1st indirect radiative effects contribute to the surface cooling. In the BC experiment, the BC direct effect is the main driver of the surface cooling, however, a decrease in low cloud cover due to the increased heating by BC absorption partially counteracts the direct effect. This results in a weaker land surface temperature response to BC changes than to sulphate changes. The resulting precipitation response is also weaker, and the responses of the monsoon circulation are different for sulphate and black carbon experiments. This study demonstrates a mechanism that links regional aerosol emission changes to the precipitation changes of the EASM, and it could be applied to help understand the future changes in EASM precipitation in CMIP5 simulations.

Predictability and uncertainty in a regional climate model

The evaluation of the quality and usefulness of climate modeling systems is dependent upon an assessment of both the limited predictability of the climate system and the uncertainties stemming from model formulation. In this study a methodology is presented that is suited to assess the performance of a regional climate model (RCM), based on its ability to represent the natural interannual variability on monthly and seasonal timescales. The methodology involves carrying out multiyear ensemble simulations (to assess the predictability bounds within which the model can be evaluated against observations) and multiyear sensitivity experiments using different model formulations (to assess the model uncertainty). As an example application, experiments driven by assimilated lateral boundary conditions and sea surface temperatures from the ECMWF Reanalysis Project (ERA-15, 1979–1993) were conducted. While the ensemble experiment demonstrates that the predictability of the regional climate varies strongly between different seasons and regions, being weakest during the summer and over continental regions, important sensitivities of the modeling system to parameterization choices are uncovered. In particular, compensating mechanisms related to the long-term representation of the water cycle are revealed, in which summer dry and hot conditions at the surface, resulting from insufficient evaporation, can persist despite insufficient net solar radiation (a result of unrealistic cloud-radiative feedbacks).

Simulated variations in atmospheric CO2 over a Wisconsin forest using a coupled ecosystem-atmosphere model

Ecosystem fluxes of energy, water, and CO2 result in spatial and temporal variations in atmospheric properties. In principle, these variations can be used to quantify the fluxes through inverse modelling of atmospheric transport, and can improve the understanding of processes and falsifiability of models. We investigated the influence of ecosystem fluxes on atmospheric CO2 in the vicinity of the WLEF-TV tower in Wisconsin using an ecophysiological model (Simple Biosphere, SiB2) coupled to an atmospheric model (Regional Atmospheric Modelling System). Model parameters were specified from satellite imagery and soil texture data. In a companion paper, simulated fluxes in the immediate tower vicinity have been compared to eddy covariance fluxes measured at the tower, with meteorology specified from tower sensors. Results were encouraging with respect to the ability of the model to capture observed diurnal cycles of fluxes. Here, the effects of fluxes in the tower footprint were also investigated by coupling SiB2 to a high-resolution atmospheric simulation, so that the model physiology could affect the meteorological environment. These experiments were successful in reproducing observed fluxes and concentration gradients during the day and at night, but revealed problems during transitions at sunrise and sunset that appear to be related to the canopy radiation parameterization in SiB2.

The surprising role of orography in the initiation of an isolated thunderstorm in southern England

Many factors, both mesoscale and larger scale, often come together in order for a particular convective initiation to take place. The authors describe a modeling study of a case from the Convective Storms Initiation Project (CSIP) in which a single thunderstorm formed behind a front in the southern United Kingdom. The key features of the case were a tongue of low-level high θw air associated with a forward-sloping split front (overrunning lower θw air above), a convergence line, and a “lid” of high static stability air, which the shower was initially constrained below but later broke through. In this paper, the authors analyze the initiation of the storm, which can be traced back to a region of high ground (Dartmoor) at around 0700 UTC, in more detail using model sensitivity studies with the Met Office Unified Model (MetUM). It is established that the convergence line was initially caused by roughness effects but had a significant thermal component later. Dartmoor had a key role in the development of the thunderstorm. A period of asymmetric flow over the high ground, with stronger low-level descent in the lee, led to a hole in a layer of low-level clouds downstream. The surface solar heating through this hole, in combination with the tongue of low-level high θw air associated with the front, caused the shower to initiate with sufficient lifting to enable it later to break through the lid.

Characteristics of high-resolution versions of the Met Office unified model for forecasting convection over the United Kingdom

With many operational centers moving toward order 1-km-gridlength models for routine weather forecasting, this paper presents a systematic investigation of the properties of high-resolution versions of the Met Office Unified Model for short-range forecasting of convective rainfall events. The authors describe a suite of configurations of the Met Office Unified Model running with grid lengths of 12, 4, and 1 km and analyze results from these models for a number of convective cases from the summers of 2003, 2004, and 2005. The analysis includes subjective evaluation of the rainfall fields and comparisons of rainfall amounts, initiation, cell statistics, and a scale-selective verification technique. It is shown that the 4- and 1-km-gridlength models often give more realistic-looking precipitation fields because convection is represented explicitly rather than parameterized. However, the 4-km model representation suffers from large convective cells and delayed initiation because the grid length is too long to correctly reproduce the convection explicitly. These problems are not as evident in the 1-km model, although it does suffer from too numerous small cells in some situations. Both the 4- and 1-km models suffer from poor representation at the start of the forecast in the period when the high-resolution detail is spinning up from the lower-resolution (12 km) starting data used. A scale-selective precipitation verification technique implies that for later times in the forecasts (after the spinup period) the 1-km model performs better than the 12- and 4-km models for lower rainfall thresholds. For higher thresholds the 4-km model scores almost as well as the 1-km model, and both do better than the 12-km model.