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.

Improving flood models using LiDAR: progress and problems

Airborne scanning laser altimetry (LiDAR) is an important new data source for river flood modelling. LiDAR can give dense and accurate DTMs of floodplains for use as model bathymetry. Spatial resolutions of 0.5m or less are possible, with a height accuracy of 0.15m. LiDAR gives a Digital Surface Model (DSM), so vegetation removal software (e.g. TERRASCAN) must be used to obtain a DTM. An example used to illustrate the current state of the art will be the LiDAR data provided by the EA, which has been processed by their in-house software to convert the raw data to a ground DTM and separate vegetation height map. Their method distinguishes trees from buildings on the basis of object size. EA data products include the DTM with or without buildings removed, a vegetation height map, a DTM with bridges removed, etc. Most vegetation removal software ignores short vegetation less than say 1m high. We have attempted to extend vegetation height measurement to short vegetation using local height texture. Typically most of a floodplain may be covered in such vegetation. The idea is to assign friction coefficients depending on local vegetation height, so that friction is spatially varying. This obviates the need to calibrate a global floodplain friction coefficient. It’s not clear at present if the method is useful, but it’s worth testing further. The LiDAR DTM is usually determined by looking for local minima in the raw data, then interpolating between these to form a space-filling height surface. This is a low pass filtering operation, in which objects of high spatial frequency such as buildings, river embankments and walls may be incorrectly classed as vegetation. The problem is particularly acute in urban areas. A solution may be to apply pattern recognition techniques to LiDAR height data fused with other data types such as LiDAR intensity or multispectral CASI data. We are attempting to use digital map data (Mastermap structured topography data) to help to distinguish buildings from trees, and roads from areas of short vegetation. The problems involved in doing this will be discussed. A related problem of how best to merge historic river cross-section data with a LiDAR DTM will also be considered. LiDAR data may also be used to help generate a finite element mesh. In rural area we have decomposed a floodplain mesh according to taller vegetation features such as hedges and trees, so that e.g. hedge elements can be assigned higher friction coefficients than those in adjacent fields. We are attempting to extend this approach to urban area, so that the mesh is decomposed in the vicinity of buildings, roads, etc as well as trees and hedges. A dominant points algorithm is used to identify points of high curvature on a building or road, which act as initial nodes in the meshing process. A difficulty is that the resulting mesh may contain a very large number of nodes. However, the mesh generated may be useful to allow a high resolution FE model to act as a benchmark for a more practical lower resolution model. A further problem discussed will be how best to exploit data redundancy due to the high resolution of the LiDAR compared to that of a typical flood model. Problems occur if features have dimensions smaller than the model cell size e.g. for a 5m-wide embankment within a raster grid model with 15m cell size, the maximum height of the embankment locally could be assigned to each cell covering the embankment. But how could a 5m-wide ditch be represented? Again, this redundancy has been exploited to improve wetting/drying algorithms using the sub-grid-scale LiDAR heights within finite elements at the waterline.

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.

Reactive uptake coefficients for heterogeneous reaction of N2O5 with submicron aerosols of NaCl and natural sea salt

The kinetics of uptake of gaseous N2O5 on submicron aerosols containing NaCl and natural sea salt have been investigated in a flow reactor as a function of relative humidity (RH) in the range 30-80% at 295±2K and a total pressure of 1bar. The measured uptake coefficients, γ, were larger on the aerosols containing sea salt compared to those of pure NaCl, and in both cases increased with increasing RH. These observations are explained in terms of the variation in the size of the salt droplets, which leads to a limitation in the uptake rate into small particles. After correction for this effect the uptake coefficients are independent of relative humidity, and agree with those measured previously on larger droplets. A value of γ=0.025 is recommended for the reactive uptake coefficient for N2O5 on deliquesced sea salt droplets at 298K and RH>40%.

Parameter estimation using a particle method: inferring mixing coefficients from sea level observations

This paper presents a first attempt to estimate mixing parameters from sea level observations using a particle method based on importance sampling. The method is applied to an ensemble of 128 members of model simulations with a global ocean general circulation model of high complexity. Idealized twin experiments demonstrate that the method is able to accurately reconstruct mixing parameters from an observed mean sea level field when mixing is assumed to be spatially homogeneous. An experiment with inhomogeneous eddy coefficients fails because of the limited ensemble size. This is overcome by the introduction of local weighting, which is able to capture spatial variations in mixing qualitatively. As the sensitivity of sea level for variations in mixing is higher for low values of mixing coefficients, the method works relatively well in regions of low eddy activity.

CISK or WISHE as the mechanism for tropical cyclone intensification

Examination of conditional instability of the second kind (CISK) and wind-induced surface heat exchange (WISHE), two proposed mechanisms for tropical cyclone and polar low intensification, suggests that the sensitivity of the intensification rate of these disturbances to surface properties, such as surface friction and moisture supply, will be different for the two mechanisms. These sensitivities were examined by perturbing the surface characteristics in a numerical model with explicit convection. The intensification rate was found to have a strong positive dependence on the heat and moisture transfer coefficients, while remaining largely insensitive to the frictional drag coefficient. CISK does not predict the observed dependence of vortex intensification rate on the heat and moisture transfer coefficients, nor the insensitivity to the frictional drag coefficient since it anticipates that intensification rate is controlled by frictional convergence in the boundary layer. Since neither conditional instability nor boundary moisture content showed any significant sensitivity to the transfer coefficients, this is true of CISK using both the convective closures of Ooyama and of Charney and Eliassen. In comparison, the WISHE intensification mechanism does predict the observed increase in intensification rate with heat and moisture transfer coefficients, while not anticipating a direct influence from surface friction.

Time-varying coefficients in the almost ideal demand system: An empirical appraisal

This paper provides a generalisation of the structural time series version of the Almost Ideal Demand System (AIDS) that allows for time-varying coefficients (TVC/AIDS) in the presence of cross-equation constraints. An empirical appraisal of the TVC/AIDS is made using a dynamic AIDS with trending intercept as the baseline model with a data set from the Italian Household Budget Survey (1986-2001). The assessment is based on four criteria: adherence to theoretical constraints, statistical diagnostics on residuals, forecasting performance and economic meaningfulness. No clear evidence is found for superior performance of the TVC/AIDS, apart from improved short-term forecasts.

Octanol-water partition coefficients for predicting the effects of tannins in ruminant nutrition

Tannins can cause beneficial or harmful nutritional effects, but their great diversity has until now prevented a rational distinction between tannin structures and their nutritional responses. An attempt has been made to study this problem by examining the octanol-water solubilities of tannins. A relatively simple HPLC method has been developed for screening mixtures of plant tannins for their octanol-water partition coefficients (K-ow coefficients). Tannins were isolated from the fruits and leaves of different Acacia, Calliandra, Dichrostachys, and Piliostigma species, which are known to produce beneficial or harmful effects. The K-ow coefficients of these tannins ranged from 0.061 to 13.9, average coefficients of variation were 9.2% and recoveries were 107%. Acacia nilotica fruits and leaves had the highest K-ow coefficients, that is, 2.0 and 13.9, respectively. These A. nilotica products also have high concentrations of tannins. The combined effects of high octanol solubilities and high tannin concentrations may explain their negative effects on animal nutrition and health. It is known that compounds with high octanol solubilities are more easily absorbed into tissues, and it is, therefore, proposed that such compounds are more likely to cause toxicity problems especially if consumed in large quantities. According to the literature, tannins in human foods tend to have low K-ow coefficients, and this was confirmed for the tannins in Piliostigma thonningii fruits. Therefore, unconventional feeds or browse products should be screened not only for their tannin concentrations but also for low octanol-water partition coefficients in order to identify nutritionally safe feeds and to avoid potentially toxic feeds.

New magnetic nano-absorbent for the determination of n-octanol/water partition coefficients

A novel and generic miniaturization methodology for the determination of partition coefficient values of organic compounds in noctanol/water by using magnetic nanoparticles is, for the first time, described. We have successfully designed, synthesised and characterised new colloidal stable porous silica-encapsulated magnetic nanoparticles of controlled dimensions. These nanoparticles absorbing a tiny amount of n-octanol in their porous silica over-layer are homogeneously dispersed into a bulk aqueous phase (pH 7.40) containing an organic compound prior to magnetic separation. The small size of the particles and the efficient mixing allow a rapid establishment of the partition equilibrium of the organic compound between the solid supported n-octanol nano-droplets and the bulk aqueous phase. UV-vis spectrophotometry is then applied as a quantitative method to determine the concentration of the organic compound in the aqueous phase both before and after partitioning (after magnetic separation). log D values of organic compounds of pharmaceutical interest (0.65-3.50), determined by this novel methodology, were found to be in excellent agreement with the values measured by the shake-flask method in two independent laboratories, which are also consistent with the literature data. It was also found that this new technique gives a number of advantages such as providing an accurate measurement of log D value, a much shorter experimental time and a smaller sample size required. With this approach, the formation of a problematic emulsion, commonly encountered in shake-flask experiments, is eliminated. It is envisaged that this method could be applicable to the high throughput log D screening of drug candidates. (c) 2005 Elsevier B.V. All rights reserved.

Gas-phase rate coefficients for reactions of NO3, OH, O-3 and O(P-3) with unsaturated alcohols and ethers: correlations and structure-activity relations (SARs)

Experimental difficulties sometimes force modellers to use predicted rate coefficients for reactions of oxygenated volatile organic compounds (oVOCs). We examine here methods for making the predictions for reactions of atmospheric initiators of oxidation, NO3, OH, O-3 and O(P-3), with unsaturated alcohols and ethers. Logarithmic correlations are found between measured rate coefficients and calculated orbital energies, and these correlations may be used directly to estimate rate coefficients for compounds where measurements have not been performed. To provide a shortcut that obviates the need to calculate orbital energies, structure-activity relations (SARs) are developed. Our SARs are tested for predictive power against compounds for which experimental rate coefficients exist, and their accuracy is discussed. Estimated atmospheric lifetimes for oVOCs are presented. The SARs for alkenols successfully predict key rate coefficients, and thus can be used to enhance the scope of atmospheric models incorporating detailed chemistry. SARs for the ethers have more limited applicability, but can still be useful in improving tropospheric models. (C) 2008 Elsevier Ltd. All rights reserved.

Gas-phase rate coefficients for the reactions of NO3, OH and O-3 with alpha,beta-unsaturated esters and ketones: structure-activity relations (SARs)

Gas-phase rate coefficients for the atmospherically important reactions of NO3, OH and O-3 are predicted for 55 alpha,beta-unsaturated esters and ketones. The rate coefficients were calculated using a correlation described previously [Pfrang, C., King, M.D., C. E. Canosa-Mas, C.E., Wayne, R.P., 2006. Atmospheric Environment 40, 1170-1179]. These rate coefficients were used to extend structure-activity relations for predicting the rate coefficients for the reactions of NO3, OH or O-3 with alkenes to include alpha,beta-unsaturated esters and ketones. Conjugation of an alkene with an alpha,beta-keto or alpha,beta-ester group will reduce the value of a rate coefficient by a factor of similar to 110, similar to 2.5 and similar to 12 for reaction with NO3, OH or O-3, respectively. The actual identity of the alkyl group, R, in -C(O)R or -C(O)OR has only a small influence. An assessment of the reliability of the SAR is given that demonstrates that it is useful for reactions involving NO3 and OH, but less valuable for those of O-3 or peroxy nitrate esters. (c) 2006 Elsevier Ltd. All rights reserved.

Toughness, the shearing of self-mated junctions and friction

Since ductile fracture (rupture) is the process by which junctions are separated and which prevents ever-increasing plasticity and junction growth, it is argued that models of friction ought to include toughness as well as yield strength. An expression for the coefficient of sliding friction is derived using ductile fracture mechanics. The predictions are quite reasonable.

Extending the friction cone algorithm for arbitrary polygon based haptic objects

Most haptic environments are based on single point interactions whereas in practice, object manipulation requires multiple contact points between the object, fingers, thumb and palm. The Friction Cone Algorithm was developed specifically to work well in a multi-finger haptic environment where object manipulation would occur. However, the Friction Cone Algorithm has two shortcomings when applied to polygon meshes: there is no means of transitioning polygon boundaries or feeling non-convex edges. In order to overcome these deficiencies, Face Directed Connection Graphs have been developed as well as a robust method for applying friction to non-convex edges. Both these extensions are described herein, as well as the implementation issues associated with them.