Use of fused airborne scanning laser altimetry and digital map data for urban flood modelling

Flood modelling of urban areas is still at an early stage, partly because until recently topographic data of sufficiently high resolution and accuracy have been lacking in urban areas. However, Digital Surface Models (DSMs) generated from airborne scanning laser altimetry (LiDAR) having sub-metre spatial resolution have now become available, and these are able to represent the complexities of urban topography. The paper describes the development of a LiDAR post-processor for urban flood modelling based on the fusion of LiDAR and digital map data. The map data are used in conjunction with LiDAR data to identify different object types in urban areas, though pattern recognition techniques are also employed. Post-processing produces a Digital Terrain Model (DTM) for use as model bathymetry, and also a friction parameter map for use in estimating spatially-distributed friction coefficients. In vegetated areas, friction is estimated from LiDAR-derived vegetation height, and (unlike most vegetation removal software) the method copes with short vegetation less than ~1m high, which may occupy a substantial fraction of even an urban floodplain. The DTM and friction parameter map may also be used to help to generate an unstructured mesh of a vegetated urban floodplain for use by a 2D finite element model. The mesh is decomposed to reflect floodplain features having different frictional properties to their surroundings, including urban features such as buildings and roads as well as taller vegetation features such as trees and hedges. This allows a more accurate estimation of local friction. The method produces a substantial node density due to the small dimensions of many urban features.

A note on boundary layer friction in baroclinic cyclones

The interaction between extratropical cyclones and the underlying boundary layer has been a topic of recent discussion in papers by Adamson et al (2006) and Beare (2007). Their results emphasise different mechanisms through which the boundary layer dynamics may modify the growth of a baroclinic cyclone. By using different sea-surface temperature distributions and comparing the low-level winds, the differences are exposed and both of the proposed mechanisms appear to be acting within a single simulation.

Numerical simulation of baroclinic waves with a parameterized boundary layer

A dry three-dimensional baroclinic life cycle model is used to investigate the role of turbulent fluxes of heat and momentum within the boundary layer on mid-latitude cyclones. Simulations are performed of life cycles for two basic states, both with and without turbulent fluxes. The different basic states produce cyclones with contrasting frontal and mesoscale-flow structures. The analysis focuses on the generation of potential-vorticity (PV) in the boundary layer and its subsequent transport into the free troposphere. The dynamic mechanism through which friction mitigates a barotropic vortex is that of Ekman pumping. This has often been assumed to be also the dominant mechanism for baroclinic developments. The PV framework highlights an additional, baroclinic mechanism. Positive PV is generated baroclinically due to friction to the north-east of a surface low and is transported out of the boundary layer by a cyclonic conveyor belt flow. The result is an anomaly of increased static stability in the lower troposphere which restricts the growth of the baroclinic wave. The reduced coupling between lower and upper levels can be sufficient to change the character of the upper-level evolution of the mature wave. The basic features of the baroclinic damping mechanism are robust for different frontal structures, with and without turbulent heat fluxes, and for the range of surface roughness found over the oceans.

Occurrence of Kelvin-Helmholtz Billows in Sea-breeze Circulations

Centred at the interface between the sea-breeze and the return flow aloft, Kelvin-Helmholtz billows (KHB) are an important feature of the turbulent structure of some sea-breeze circulations (SBCs). In other SBCs, there are no prominent KHBs observed. Factors governing the appearance of billows are determined from a database of 139 sea breezes, constructed from two years of summertime surface observations at a site on the south coast of England. Post-frontal oscillations occur in the surface data for some SBCs and are interpreted as indicating possible KHBs aloft. The SBCs are formed under a wide range of synoptic conditions, enabling various measures of possible billow occurrence to be related to properties of the large-scale, ambient flow. Consistent with laboratory experiments of density currents, KHBs are suppressed for propagation into a head wind and enhanced with a tail wind. They are also found to be enhanced for stronger ambient wind speeds, while large-scale coast-parallel flow is effective in suppressing the billows.

Improved simulation of water vapour and clear-sky radiation using 24 hour forecasts from ERA40

Monthly mean water vapour and clear-sky radiation extracted from the European Centre for Medium Range Weather Forecasts 40-year reanalysis (ERA40) forecasts are assessed using satellite observations and additional reanalysis data. There is a marked improvement in the interannual variability of column-integrated water vapour (CWV) over the oceans when using the 24-hour forecasts compared with the standard 6-hour forecasts products. The spatial distribution of CWV are well simulated by the 6-hour forecasts; using the 24-hour forecasts does not degrade this simulation substantially and in many cases improves on the quality. There is also an improved simulation of clear-sky radiation from the 24-hour forecasts compared with the 6-hour forecasts based on comparison with satellite observations and empirical estimates. Further work is required to assess the quality of water vapour simulation by reanalyses over land regions. Over the oceans, it is recommended that 24-hour forecasts of CWV and clear-sky radiation are used in preference to the standard 6-hour forecast products from ERA40

Variability of the Atlantic thermohaline circulation described by three-dimensional empirical orthogonal functions

We describe the use of bivariate 3d empirical orthogonal functions (EOFs) in characterising low frequency variability of the Atlantic thermohaline circulation (THC) in the Hadley Centre global climate model, HadCM3. We find that the leading two modes are well correlated with an index of the meridional overturning circulation (MOC) on decadal timescales, with the leading mode alone accounting for 54% of the decadal variance. Episodes of coherent oscillations in the sub-space of the leading EOFs are identified; these episodes are of great interest for the predictability of the THC, and could indicate the existence of different regimes of natural variability. The mechanism identified for the multi-decadal variability is an internal ocean mode, dominated by changes in convection in the Nordic Seas, which lead the changes in the MOC by a few years. Variations in salinity transports from the Arctic and from the North Atlantic are the main feedbacks which control the oscillation. This mode has a weak feedback onto the atmosphere and hence a surface climatic influence. Interestingly, some of these climate impacts lead the changes in the overturning. There are also similarities to observed multi-decadal climate variability.

Enhancement of ENSO variability by a weakened Atlantic thermohaline circulation in a coupled GCM

A coupled ocean–atmosphere general circulation model is used to investigate the modulation of El Niño–Southern Oscillation (ENSO) variability due to a weakened Atlantic thermohaline circulation (THC). The THC weakening is induced by freshwater perturbations in the North Atlantic, and leads to a well-known sea surface temperature dipole and a southward shift of the intertropical convergence zone (ITCZ) in the tropical Atlantic. Through atmospheric teleconnections and local coupled air–sea feedbacks, a meridionally asymmetric mean state change is generated in the eastern equatorial Pacific, corresponding to a weakened annual cycle, and westerly anomalies develop over the central Pacific. The westerly anomalies are associated with anomalous warming of SST, causing an eastward extension of the west Pacific warm pool particularly in August–February, and enhanced precipitation. These and other changes in the mean state lead in turn to an eastward shift of the zonal wind anomalies associated with El Niño events, and a significant increase in ENSO variability. In response to a 1-Sv (1 Sv ≡ 106 m3 s−1) freshwater input in the North Atlantic, the THC slows down rapidly and it weakens by 86% over years 50–100. The Niño-3 index standard deviation increases by 36% during the first 100-yr simulation relative to the control simulation. Further analysis indicates that the weakened THC not only leads to a stronger ENSO variability, but also leads to a stronger asymmetry between El Niño and La Niña events. This study suggests a role for an atmospheric bridge that rapidly conveys the influence of the Atlantic Ocean to the tropical Pacific and indicates that fluctuations of the THC can mediate not only mean climate globally but also modulate interannual variability. The results may contribute to understanding both the multidecadal variability of ENSO activity during the twentieth century and longer time-scale variability of ENSO, as suggested by some paleoclimate records.

Monitoring present day changes in water vapour and the radiative energy balance using satellite data, reanalyses and models

A combination of satellite data, reanalysis products and climate models are combined to monitor changes in water vapour, clear-sky radiative cooling of the atmosphere and precipitation over the period 1979-2006. Climate models are able to simulate observed increases in column integrated water vapour (CWV) with surface temperature (Ts) over the ocean. Changes in the observing system lead to spurious variability in water vapour and clear-sky longwave radiation in reanalysis products. Nevertheless all products considered exhibit a robust increase in clear-sky longwave radiative cooling from the atmosphere to the surface; clear-sky longwave radiative cooling of the atmosphere is found to increase with Ts at the rate of ~4 Wm-2 K-1 over tropical ocean regions of mean descending vertical motion. Precipitation (P) is tightly coupled to atmospheric radiative cooling rates and this implies an increase in P with warming at a slower rate than the observed increases in CWV. Since convective precipitation depends on moisture convergence, the above implies enhanced precipitation over convective regions and reduced precipitation over convectively suppressed regimes. To quantify this response, observed and simulated changes in precipitation rate are analysed separately over regions of mean ascending and descending vertical motion over the tropics. The observed response is found to be substantially larger than the model simulations and climate change projections. It is currently not clear whether this is due to deficiencies in model parametrizations or errors in satellite retrievals.