The European storm Kyrill in January 2007: synoptic evolution, meteorological impacts and some considerations with respect to climate change

The synoptic evolution and some meteorological impacts of the European winter storm Kyrill that swept across Western, Central, and Eastern Europe between 17 and 19 January 2007 are investigated. The intensity and large storm damage associated with Kyrill is explained based on synoptic and mesoscale environmental storm features, as well as on comparisons to previous storms. Kyrill appeared on weather maps over the US state of Arkansas about four days before it hit Europe. It underwent an explosive intensification over the Western North Atlantic Ocean while crossing a very intense zonal polar jet stream. A superposition of several favourable meteorological conditions west of the British Isles caused a further deepening of the storm when it started to affect Western Europe. Evidence is provided that a favourable alignment of three polar jet streaks and a dry air intrusion over the occlusion and cold fronts were causal factors in maintaining Kyrill's low pressure very far into Eastern Europe. Kyrill, like many other strong European winter storms, was embedded in a pre-existing, anomalously wide, north-south mean sea-level pressure (MSLP) gradient field. In addition to the range of gusts that might be expected from the synoptic-scale pressure field, mesoscale features associated with convective overturning at the cold front are suggested as the likely causes for the extremely damaging peak gusts observed at many lowland stations during the passage of Kyrill's cold front. Compared to other storms, Kyrill was by far not the most intense system in terms of core pressure and circulation anomaly. However, the system moved into a pre-existing strong MSLP gradient located over Central Europe which extended into Eastern Europe. This fact is considered determinant for the anomalously large area affected by Kyrill. Additionally, considerations of windiness in climate change simulations using two state-of-the-art regional climate models driven by ECHAM5 indicate that not only Central, but also Eastern Central Europe may be affected by higher surface wind speeds at the end of the 21st century. These changes are partially associated with the increased pressure gradient over Europe which is identified in the ECHAM5 simulations. Thus, with respect to the area affected, as well as to the synoptic and mesoscale storm features, it is proposed that Kyrill may serve as an interesting study case to assess future storm impacts.

The integrated WRF/urban modelling system: development, evaluation, and applications to urban environmental problems

To bridge the gaps between traditional mesoscale modelling and microscale modelling, the National Center for Atmospheric Research, in collaboration with other agencies and research groups, has developed an integrated urban modelling system coupled to the weather research and forecasting (WRF) model as a community tool to address urban environmental issues. The core of this WRF/urban modelling system consists of the following: (1) three methods with different degrees of freedom to parameterize urban surface processes, ranging from a simple bulk parameterization to a sophisticated multi-layer urban canopy model with an indoor–outdoor exchange sub-model that directly interacts with the atmospheric boundary layer, (2) coupling to fine-scale computational fluid dynamic Reynolds-averaged Navier–Stokes and Large-Eddy simulation models for transport and dispersion (T&D) applications, (3) procedures to incorporate high-resolution urban land use, building morphology, and anthropogenic heating data using the National Urban Database and Access Portal Tool (NUDAPT), and (4) an urbanized high-resolution land data assimilation system. This paper provides an overview of this modelling system; addresses the daunting challenges of initializing the coupled WRF/urban model and of specifying the potentially vast number of parameters required to execute the WRF/urban model; explores the model sensitivity to these urban parameters; and evaluates the ability of WRF/urban to capture urban heat islands, complex boundary-layer structures aloft, and urban plume T&D for several major metropolitan regions. Recent applications of this modelling system illustrate its promising utility, as a regional climate-modelling tool, to investigate impacts of future urbanization on regional meteorological conditions and on air quality under future climate change scenarios. Copyright © 2010 Royal Meteorological Society

A GCSS model intercomparison for a tropical squall line observed during toga-coare. II: Intercomparison of single-column models and a cloud-resolving model

This paper presents single-column model (SCM) simulations of a tropical squall-line case observed during the Coupled Ocean-Atmosphere Response Experiment of the Tropical Ocean/Global Atmosphere Programme. This case-study was part of an international model intercomparison project organized by Working Group 4 ‘Precipitating Convective Cloud Systems’ of the GEWEX (Global Energy and Water-cycle Experiment) Cloud System Study. Eight SCM groups using different deep-convection parametrizations participated in this project. The SCMs were forced by temperature and moisture tendencies that had been computed from a reference cloud-resolving model (CRM) simulation using open boundary conditions. The comparison of the SCM results with the reference CRM simulation provided insight into the ability of current convection and cloud schemes to represent organized convection. The CRM results enabled a detailed evaluation of the SCMs in terms of the thermodynamic structure and the convective mass flux of the system, the latter being closely related to the surface convective precipitation. It is shown that the SCMs could reproduce reasonably well the time evolution of the surface convective and stratiform precipitation, the convective mass flux, and the thermodynamic structure of the squall-line system. The thermodynamic structure simulated by the SCMs depended on how the models partitioned the precipitation between convective and stratiform. However, structural differences persisted in the thermodynamic profiles simulated by the SCMs and the CRM. These differences could be attributed to the fact that the total mass flux used to compute the SCM forcing differed from the convective mass flux. The SCMs could not adequately represent these organized mesoscale circulations and the microphysicallradiative forcing associated with the stratiform region. This issue is generally known as the ‘scale-interaction’ problem that can only be properly addressed in fully three-dimensional simulations. Sensitivity simulations run by several groups showed that the time evolution of the surface convective precipitation was considerably smoothed when the convective closure was based on convective available potential energy instead of moisture convergence. Finally, additional SCM simulations without using a convection parametrization indicated that the impact of a convection parametrization in forced SCM runs was more visible in the moisture profiles than in the temperature profiles because convective transport was particularly important in the moisture budget.

An assessment of Indian monsoon seasonal forecasts and mechanisms underlying monsoon interannual variability in the Met Office GloSea5-GC2 system

We assess Indian summer monsoon seasonal forecasts in GloSea5-GC2, the Met Office fully coupled subseasonal to seasonal ensemble forecasting system. Using several metrics, GloSea5-GC2 shows similar skill to other state-of-the-art forecast systems. The prediction skill of the large-scale South Asian monsoon circulation is higher than that of Indian monsoon rainfall. Using multiple linear regression analysis we evaluate relationships between Indian monsoon rainfall and five possible drivers of monsoon interannual variability. Over the time period studied (1992-2011), the El Nino-Southern Oscillation (ENSO) and the Indian Ocean dipole (IOD) are the most important of these drivers in both observations and GloSea5-GC2. Our analysis indicates that ENSO and its teleconnection with the Indian rainfall are well represented in GloSea5-GC2. However, the relationship between the IOD and Indian rainfall anomalies is too weak in GloSea5-GC2, which may be limiting the prediction skill of the local monsoon circulation and Indian rainfall. We show that this weak relationship likely results from a coupled mean state bias that limits the impact of anomalous wind forcing on SST variability, resulting in erroneous IOD SST anomalies. Known difficulties in representing convective precipitation over India may also play a role. Since Indian rainfall responds weakly to the IOD, it responds more consistently to ENSO than in observations. Our assessment identifies specific coupled biases that are likely limiting GloSea5-GC2 prediction skill, providing targets for model improvement.

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

A new method for the characterisation and verification of local spatial predictability for convective scale ensembles

The use of kilometre-scale ensembles in operational forecasting provides new challenges for forecast interpretation and evaluation to account for uncertainty on the convective scale. A new neighbourhood based method is presented for evaluating and characterising the local predictability variations from convective scale ensembles. Spatial scales over which ensemble forecasts agree (agreement scales, S^A) are calculated at each grid point ij, providing a map of the spatial agreement between forecasts. By comparing the average agreement scale obtained from ensemble member pairs (S^A(mm)_ij), with that between members and radar observations (S^A(mo)_ij), this approach allows the location-dependent spatial spread-skill relationship of the ensemble to be assessed. The properties of the agreement scales are demonstrated using an idealised experiment. To demonstrate the methods in an operational context the S^A(mm)_ij and S^A(mo)_ij are calculated for six convective cases run with the Met Office UK Ensemble Prediction System. The S^A(mm)_ij highlight predictability differences between cases, which can be linked to physical processes. Maps of S^A(mm)_ij are found to summarise the spatial predictability in a compact and physically meaningful manner that is useful for forecasting and for model interpretation. Comparison of S^A(mm)_ij and S^A(mo)_ij demonstrates the case-by-case and temporal variability of the spatial spread-skill, which can again be linked to physical processes.

Evaluation of the Plant–Craig stochastic convection scheme (v2.0) in the ensemble forecasting system MOGREPS-R (24 km) based on the Unified Model (v7.3)

The Plant–Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic scheme only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant–Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant–Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.

New perspectives on the synoptic and mesoscale structure of Hurricane Catarina

This work explores in detail synoptic and mesoscale features of Hurricane Catarina during its life cycle from a decaying baroclinic wave to a tropical depression that underwent tropical transition (TT) and finally to a Category 2 hurricane at landfall over Santa Catarina State coast, southern Brazil. This unique system caused 11 deaths mostly off the Brazilian coast and an estimated half billion dollars in damage in a matter of a few hours on 28 March 2004. Although the closest meteorological station available was tens of kilometres away from the eye, in situ meteorological measurements provided by a work-team sent to the area where the eye made landfall unequivocally reproduces the tropical signature with category 2 strength, adding to previous analysis where this data was not available. Further analyses are based mostly on remote sensing data available at the time of the event. A classic dipole blocking set synoptic conditions for Hurricane Catarina to develop, dynamically contributing to the low wind shear observed. On the other hand, on its westward transit, large scale subsidence limited its strength and vertical development. Catarina had relatively cool SST conditions, but this was mitigated by favourable air-sea fluxes leading to latent heat release-driven processes during the mature phase. The ocean`s dynamic topography also suggested the presence of nearby warm core rings which may have facilitated the transition and post-transition intensification. Since there were no records of such a system at least in the past 30 years and given that SSTs were generally below 26 degrees C and vertical shear was usually strong, despite all satellite data available, the system was initially classified as an extratropical cyclone. Here we hypothesise that this categorization was based oil inadequate regional scale model outputs which did not account for the importance of the latent heat fluxes over the ocean. Hurricane Catarina represents a dramatic event on weather systems in South America. It has attracted attention worldwide and poses questions as whether or not it is a symptom of global warming. (C) 2009 Elsevier B.V. All rights reserved.

The South American Monsoon System and the 1970s climate transition

The South American Monsoon System (SAMS) is characterised by intense convective activity and precipitation during austral summer. This study investigates changes in the onset, demise and duration of SAMS during 1948-2008. The results show a significant change in these characteristics in the early 1970s. Onset becomes steadily earlier from 1948 to early 1970s and has occurred earlier than 23-27 October after 1972-1973. Demise dates have remained later than 21-25 April after the mid-to-late 1970s. SAMS duration shows a statistical changepoint in the summer of 1971-1972 such that the mean duration was similar to 170 days (1948-1972) and 195 days (1972-1982). Vertically integrated moisture flux is used to diagnose changes in mean state and reveal statistically significant increases over South America after 1971-1972. Copyright. (C) 2010 Royal Meteorological Society

Moisture transport and intraseasonal variability in the South America monsoon system

This paper examines moisture transport on intraseasonal timescales over the continent and over the South Atlantic convergence zone (SACZ) during the South America (SA) summer monsoon. Combined Empirical Orthogonal Function analysis (EOFc) of Global Precipitation Climatology Project pentad precipitation, specific humidity, air temperature, zonal and meridional winds at 850 hPa (NCEP/NCAR reanalysis) are performed to identify the large-scale variability of the South America monsoon system and the SACZ. The first EOFc was used as a large-scale index for the South American monsoon (LISAM), whereas the second EOFc characterized the SACZ. LISAM (SACZ) index showed spectral variance on 30-90 (15-20) days and were both band filtered (10-100 days). Intraseasonal wet anomalies were defined when LISAM and SACZ anomalies were above the 75th percentile of their respective distribution. LISAM and SACZ wet events were examined independently of each other and when they occur simultaneously. LISAM wet events were observed with the amplification of wave activity in the Northern Hemisphere and the enhancement of northwesterly cross-equatorial moisture transport over tropical continental SA. Enhanced SACZ was observed with moisture transport from the extratropics of the Southern Hemisphere. Simultaneous LISAM and SACZ wet events are associated with cross-equatorial moisture transport along with moisture transport from Subtropical Southwestern Atlantic.

NONLINEAR DIFFUSION PROCESS IN A BENARD SYSTEM AT THE CRITICAL-POINT FOR THE ONSET OF CONVECTION

The evolution equation governing surface perturbations of a shallow fluid heated from below at the critical Rayleigh number for the onset of convective motion, and with boundary conditions leading to zero critical wave number, is obtained. A solution for negative or cooling perturbations is explicitly exhibited, which shows that the system presents sharp propagating fronts.

A case study of convective organization into precipitating lines in the Southwest Amazon during the WETAMC and TRMM-LBA

A case study of convective development in the Southwest Amazon region during the Wet Season Atmospheric Mesoscale Campaign (WETAMC) and Tropical Rainfall Measuring Mission (TRMM)/Large-Scale Biosphere-Atmosphere (LBA) Experiment in Amazonia is presented. The convective development during 7 February 1999 is shown to occur during a period of very weak large-scale forcing in the presence of topography and deforestation. The available data include dual Doppler radar analysis, radiosonde launches, and surface and boundary layer observations. The observational analysis is complemented with a series of model simulations using the RAMS with 2-km resolution over a 300 km 300 km area forced by a morning radiosonde profile. A comparison of the observed and simulated thermodynamic transformation of the boundary layer and of the formation of convective lines, and of their kinematic and microphysical properties is presented. It is shown that only a few very deep and intense convective cells are necessary to explain the overall precipitating line formation and that discrete propagation and coupling with upper atmosphere circulations may explain the appearance of several lines. The numerical simulation indicates that topography may be the cause of initial convective development, although later on the convective line is parallel to the midlevel shear. There are indications that small-scale deforestation may have an effect on increasing rainfall in the wet season when the large-scale forcing is very weak.

Implementation of convective heating in Companhia Siderúrgica Nacional Blast Furnace runners

This paper describes the main characteristics and advantages of convective heating system for refractory lining, compared with conventional heating systems. In addition the main results obtained are presented with its implementation in CSN Blast Furnace #2 and 3 Runners, in terms of cost and equipment availability, as well as the need for ceramic coating to protect the lining against oxidation, arising from excessive air combustion. © 2012 Elsevier Ltd. All rights reserved.

A two-layer approximation to the Brazil Current-Intermediate Western Boundary Current System between 20 degrees S and 28 degrees S

It has been shown that the vertical structure of the Brazil Current (BC)-Intermediate Western Boundary Current (IWBC) System is dominated by the first baroclinic mode at 22 degrees S-23 degrees S. In this work, we employed the Miami Isopycnic Coordinate Ocean Model to investigate whether the rich mesoscale activity of this current system, between 20 degrees S and 28 degrees S, is reproduced by a two-layer approximation of its vertical structure. The model results showed cyclonic and anticyclonic meanders propagating southwestward along the current axis, resembling the dynamical pattern of Rossby waves superposed on a mean flow. Analysis of the upper layer zonal velocity component, using a space-time diagram, revealed a dominant wavelength of about 450 km and phase velocity of about 0.20 ms(-1) southwestward. The results also showed that the eddy-like structures slowly grew in amplitude as they moved downstream. Despite the simplified design of the numerical experiments conducted here, these results compared favorably with observations and seem to indicate that weakly unstable long baroclinic waves are responsible for most of the variability observed in the BC-IWBC system. (C) 2009 Elsevier Ltd. All rights reserved.