Synthesis, structure and reactivity of novel CuI, CuII and CuIII complexes containing triaza and hexaaza macrocyclic ligands

El treball de tesi s'emmarca dins del camp de la bioinorgànica, disciplina que estudia les propietats estructurals i de reactivitat dels centres actius dels enzims, servint-se de models síntètics de baix pes molecular per tal d'intentar reproduïr la reactivitat presentada per l'enzim i conèixer els mecanismes de reacció a nivell molecular que tenen lloc en els processos biològics.1 Més concretament el treball posa especial èmfasi en els processos d'activació d'oxigen molecular que tenen lloc en les metaloproteïnes de Coure del Tipus 3, com són l'hemocianina i la tirosinasa, ambdues presentant un complex dinuclear de Cu(I)) en el centre actiu de la forma reduïda, capaç d'activar l'O2 cap a espècies de tipus peròxid.2 Un altre camp d'interès ha estat l'estudi dels processos d'activació d'enllaços C-H no activats en hidrocarburs, tant per la seva importàcia a nivell industrial com per comprendre els mecanismes intrínsecs d'aquesta activació a través de metalls de trancisió.3,4 Durant el treball de tesi presentat s'ha desenvolupat la síntesi de nous complexes de Coure(I), Coure(II) y Cu(III) utilitzant lligands macrocíclics de tipus triaza i hexaaza, i s'han estudiat la seves propietats estructurals així com la seva reactivitat. La reacció dels lligands triazacíclics H32m, H2Me33m i H33m amb sals de coure(II) dóna lloc a una reacció de desproporció de Cu(II) per obtenir-se en quantitats equimolars un complex organometàl·lic de Cu(III) i un complex de Cu(I). La caracterizació estructural exhaustiva dels complexes del tipus aryl-Cu(III) evidencia la formació d'un enllaç organometàl·lic entre l'àtom de Cu(III) i el carboni més próxim de l'anell aromàtic del lligand. Aquesta reacció, a més de representar una nova forma de desproporció en la química del Cu, suposa l'activació d'un enllaç C-H aromàtic a temperatura ambient que, mitjançant l'estudi cinètic d'aquesta desproporció per espectroscòpia UV-Vis, dels càlcul de l'efecte cinètic isotòpic utilitzant el lligand deuterat en el C-H de l'anell, juntament amb el recolzament teòrics dels càlculs DFT per a la optimització de geometries d'intermedis de reacció, ens permeten proposar un mecanisme de reacció pel nostre sistema, on l'activació de l'enllaç C-H aromàtic transcorre per la formació d'un enllaç de tipus agòstic C-H ? Cu(II),5 seguit de la desprotonació del C-H aromàtic per acció d'una base i posterior transferència electrònica per obtenir el complex organometàlic de Cu(III) i el complex de de Cu(I). En quant a la reactivitat d'aquests complexes organometàl·lics aryl-Cu(III) s'ha observat que una base en medi aquós causa la inestabilitat d'aquests compostos, evolucionant cap a la inserció d'un àtom d'oxigen sobre la posició activada de l'anell aromàtic, per a donar lloc a un complex dinuclear de Cu(II) amb dos grups fenoxo actuant de pont entre els àtoms metàl·lics. La reacció transcorre per un intermedi colorejat, caracteritzat com el complex ayl-Cu(III) monodesprotonat en una de les seves amines benzíliques, els quals s'observen igualment en la reacció dels correponents complexos de Cu(I) amb oxigen molecular (O2). És en els nostres sistemes en els quals es descriu per primera vegada la participació d'intermedis organometàl·lics Cu(III)-C en processos d'hidroxilació aromàtica, tals com el desenvolupat per l'enzim tirosinasa o per alguns dels seus models químics de síntesi.6,7,8 S'han estudiat les propietats magnètiques dels quatre bis(fenoxo)complexes de Cu(II) descrits, obtenint-se uns acoplaments de tipus antiferromagnètic o ferromagnètic de diversa magnitud, depenent del solapament orbitalari a l'enllaç Cu-O, a través del qual es produeix el superintercanvi. Nous complexos de Cu(I) sintetitzats amb lligands hexaazamacrocíclics han estat estudiats, i posant especial èmfasi a la seva reactivitat respecta a l'activació d'oxigen molecular (O2). S'ha observat una reactivitat diferenciada segons la concentració de complex de Cu(I) utilitzada, de manera que a altes concentracions s'obté un carbonato complex tetranuclear de Cu(II) per fixació de CO2 atmosfèric, mentre que a baixes concentracions s'observa la hidroxilació aromàtica intramolecular d'un dels anells benzílics del lligand, reacció que presumiblement transcorre per atac electrofílic d'un peroxo complex intermedi sobre el sistema ? de l'anell.6 Els resultats obtinguts en aquest treball ens mostren la facilitat per activar enllaços C-H aromàtics per metalls de transició de la primera sèrie (Cu, Ni) quan aquests estan suficientment pròxims a l'enllaç C-H, en unes condicions de reacció molt suaus (1atm., temperatura ambient). Els nous complexos organometàl·lics Aryl-Cu(III) són el producte d'una nova reacció de desproporció de Cu(II), així com un posició aromàtica activada que podria ser el punt de partida per l'estudi de funcionalització selectiva d'aquests grups aromàtics.

Diversitat de Bacteris Verds del Sofre en llacunes salines i osmoadaptació amb N(epsilon)-acetil-ß-lisina

L'objectiu principal del treball de recerca és aprofundir en el coneixement del grup de Bacteris Verds del Sofre (BVS) des del punt de vista de la relació amb la salinitat. L'estudi s'ha restringit als representants halotolerants o halòfils del grup amb tres línies de treball que corresponen a diferents àmbits de coneixement, com són la descripció del biòtop i la riquesa específica de les comunitats de BVS en els ambients naturals, l'avaluació de la significació taxonòmica de la capacitat de desenvolupament en ambients salins en el marc de discussió de la taxonomia-filogenètica del grup, i l'anàlisi de les estratègies fisiològiques desenvolupades per a l'osmoadaptació en medi salí.

Selective submap joining SLAM for autonomous vehicles

Simultaneous Localization and Mapping (SLAM) do not result in consistent maps of large areas because of gradual increase of the uncertainty for long term missions. In addition, as the size of the map grows the computational cost increases, making SLAM solutions unsuitable for on-line applications. This thesis surveys SLAM approaches paying special attention to those approaches aimed to work on large scenarios. Special focus is given to existing underwater SLAM applications. A technique based on using independent local maps together with a global stochastic map is presented. This technique is called Selective Submap Joining SLAM (SSJS). A global map contains relative transformations between local maps, which are updated once a new loop is detected. Maps sharing several features are fused, maintaining the correlation between landmarks and vehicle's pose. The use of local maps reduces computational costs and improves map consistency as compared to state of the art techniques.

Desenvolupament computacional de la semblança molecular quàntica

La present Tesi Doctoral, titulada desenvolupament computacional de la semblança molecular quàntica, tracta, fonamentalment, els aspectes de càlcul de mesures de semblança basades en la comparació de funcions de densitat electrònica.El primer capítol, Semblança quàntica, és introductori. S'hi descriuen les funcions de densitat de probabilitat electrònica i llur significança en el marc de la mecànica quàntica. Se n'expliciten els aspectes essencials i les condicions matemàtiques a satisfer, cara a una millor comprensió dels models de densitat electrònica que es proposen. Hom presenta les densitats electròniques, mencionant els teoremes de Hohenberg i Kohn i esquematitzant la teoria de Bader, com magnituds fonamentals en la descripció de les molècules i en la comprensió de llurs propietats.En el capítol Models de densitats electròniques moleculars es presenten procediments computacionals originals per l'ajust de funcions densitat a models expandits en termes de gaussianes 1s centrades en els nuclis. Les restriccions físico-matemàtiques associades a les distribucions de probabilitat s'introdueixen de manera rigorosa, en el procediment anomenat Atomic Shell Approximation (ASA). Aquest procediment, implementat en el programa ASAC, parteix d'un espai funcional quasi complert, d'on se seleccionen variacionalment les funcions o capes de l'expansió, d'acord als requisits de no negativitat. La qualitat d'aquestes densitats i de les mesures de semblança derivades es verifica abastament. Aquest model ASA s'estén a representacions dinàmiques, físicament més acurades, en quant que afectades per les vibracions nuclears, cara a una exploració de l'efecte de l'esmorteïment dels pics nuclears en les mesures de semblança molecular. La comparació de les densitats dinàmiques respecte les estàtiques evidencia un reordenament en les densitats dinàmiques, d'acord al que constituiria una manifestació del Principi quàntic de Le Chatelier. El procediment ASA, explícitament consistent amb les condicions de N-representabilitat, s'aplica també a la determinació directe de densitats electròniques hidrogenoides, en un context de teoria del funcional de la densitat.El capítol Maximització global de la funció de semblança presenta algorismes originals per la determinació de la màxima sobreposició de les densitats electròniques moleculars. Les mesures de semblança molecular quàntica s'identifiquen amb el màxim solapament, de manera es mesuri la distància entre les molècules, independentment dels sistemes de referència on es defineixen les densitats electròniques. Partint de la solució global en el límit de densitats infinitament compactades en els nuclis, es proposen tres nivells de aproximació per l'exploració sistemàtica, no estocàstica, de la funció de semblança, possibilitant la identificació eficient del màxim global, així com també dels diferents màxims locals. Es proposa també una parametrització original de les integrals de recobriment a través d'ajustos a funcions lorentzianes, en quant que tècnica d'acceleració computacional. En la pràctica de les relacions estructura-activitat, aquests avenços possibiliten la implementació eficient de mesures de semblança quantitatives, i, paral·lelament, proporcionen una metodologia totalment automàtica d'alineació molecular. El capítol Semblances d'àtoms en molècules descriu un algorisme de comparació dels àtoms de Bader, o regions tridimensionals delimitades per superfícies de flux zero de la funció de densitat electrònica. El caràcter quantitatiu d'aquestes semblances possibilita la mesura rigorosa de la noció química de transferibilitat d'àtoms i grups funcionals. Les superfícies de flux zero i els algorismes d'integració usats han estat publicats recentment i constitueixen l'aproximació més acurada pel càlcul de les propietats atòmiques. Finalment, en el capítol Semblances en estructures cristal·lines hom proposa una definició original de semblança, específica per la comparació dels conceptes de suavitat o softness en la distribució de fonons associats a l'estructura cristal·lina. Aquests conceptes apareixen en estudis de superconductivitat a causa de la influència de les interaccions electró-fonó en les temperatures de transició a l'estat superconductor. En aplicar-se aquesta metodologia a l'anàlisi de sals de BEDT-TTF, s'evidencien correlacions estructurals entre sals superconductores i no superconductores, en consonància amb les hipòtesis apuntades a la literatura sobre la rellevància de determinades interaccions.Conclouen aquesta tesi un apèndix que conté el programa ASAC, implementació de l'algorisme ASA, i un capítol final amb referències bibliogràfiques.

Contribución al estudio de la extracción química del Fe3+ contaminante de la madera arqueológica subacuática impregnada con polietilenglicol

Due to the large quantity of contamination by Fe in the waterlogged wood and the composites previously processed with poliethilenglycol, PEG , we have thought the possible extraction and formation of complexes Fe-L (L= 2-Phosphonobuta-1,2,4- tricarboxylic acid, (PBTC)) and their effect on (i) structure of the wood, (ii) structure and physiques characteristics of PEG and (iii) post-treatment contamination by Fe coming from atmospheric agents. This work is a project for to study the formation of these complexes and possible modifications in structures.

Imagem por ressonância magnética da vasculatura da pele humana : viabilidade no contexto clínico

Imagiologia por Ressonância Magnética (IRM) é uma modalidade de imagem médica que está a recuperar o interesse como uma técnica não invasiva no estudo da pele. Tipicamente campos magnéticos de elevada densidade e quipamentos específicos são usados. Este facto limita o usos da técnica a laboratórios e centros de investigação especializados. Neste trabalho estudou-se a viabilidade do uso da IRM no estudo da pele e da sua vasculatura usando equipamento convencional disponível em contexto clínico. Sequências IRM para imagem estrutural e veascular foram optimizadas e testadas para obtenção de imagens da pele do punho de 6 voluntários saudáveis. As sequências observáveis dos vasos, razão sinal-ruído, e razão contraste-ruído. Foi observado que duas sequências volumétricas baseadas em eco de gradiente e com ponderações T1 e T2 forneciam informação complementar em respeito à vasculatura da pele com resoluções espaciais da ordem dos micrómetros, podendo ainda esta informação ser fundida com imagens estruturais das cadamas da pele. Foi igualmente observado que estas sequências fornecem informação útil usando equipamento convencional e perspectiva-se a sua utilização no estudo das vasculatura de tumores cutâneos e na doença vascular periférica.

Benchmarking 2D hydraulic models for urban flooding

This paper describes benchmark testing of six two-dimensional (2D) hydraulic models (DIVAST, DIVASTTVD, TUFLOW, JFLOW, TRENT and LISFLOOD-FP) in terms of their ability to simulate surface flows in a densely urbanised area. The models are applied to a 1·0 km × 0·4 km urban catchment within the city of Glasgow, Scotland, UK, and are used to simulate a flood event that occurred at this site on 30 July 2002. An identical numerical grid describing the underlying topography is constructed for each model, using a combination of airborne laser altimetry (LiDAR) fused with digital map data, and used to run a benchmark simulation. Two numerical experiments were then conducted to test the response of each model to topographic error and uncertainty over friction parameterisation. While all the models tested produce plausible results, subtle differences between particular groups of codes give considerable insight into both the practice and science of urban hydraulic modelling. In particular, the results show that the terrain data available from modern LiDAR systems are sufficiently accurate and resolved for simulating urban flows, but such data need to be fused with digital map data of building topology and land use to gain maximum benefit from the information contained therein. When such terrain data are available, uncertainty in friction parameters becomes a more dominant factor than topographic error for typical problems. The simulations also show that flows in urban environments are characterised by numerous transitions to supercritical flow and numerical shocks. However, the effects of these are localised and they do not appear to affect overall wave propagation. In contrast, inertia terms are shown to be important in this particular case, but the specific characteristics of the test site may mean that this does not hold more generally.

Detection of flood extent in urban areas using high resolution TerraSAR-X data

Flooding is a major hazard in both rural and urban areas worldwide, but it is in urban areas that the impacts are most severe. An investigation of the ability of high resolution TerraSAR-X data to detect flooded regions in urban areas is described. An important application for this would be the calibration and validation of the flood extent predicted by an urban flood inundation model. To date, research on such models has been hampered by lack of suitable distributed validation data. The study uses a 3m resolution TerraSAR-X image of a 1-in-150 year flood near Tewkesbury, UK, in 2007, for which contemporaneous aerial photography exists for validation. The DLR SETES SAR simulator was used in conjunction with airborne LiDAR data to estimate regions of the TerraSAR-X image in which water would not be visible due to radar shadow or layover caused by buildings and taller vegetation, and these regions were masked out in the flood detection process. A semi-automatic algorithm for the detection of floodwater was developed, based on a hybrid approach. Flooding in rural areas adjacent to the urban areas was detected using an active contour model (snake) region-growing algorithm seeded using the un-flooded river channel network, which was applied to the TerraSAR-X image fused with the LiDAR DTM to ensure the smooth variation of heights along the reach. A simpler region-growing approach was used in the urban areas, which was initialized using knowledge of the flood waterline in the rural areas. Seed pixels having low backscatter were identified in the urban areas using supervised classification based on training areas for water taken from the rural flood, and non-water taken from the higher urban areas. Seed pixels were required to have heights less than a spatially-varying height threshold determined from nearby rural waterline heights. Seed pixels were clustered into urban flood regions based on their close proximity, rather than requiring that all pixels in the region should have low backscatter. This approach was taken because it appeared that urban water backscatter values were corrupted in some pixels, perhaps due to contributions from side-lobes of strong reflectors nearby. The TerraSAR-X urban flood extent was validated using the flood extent visible in the aerial photos. It turned out that 76% of the urban water pixels visible to TerraSAR-X were correctly detected, with an associated false positive rate of 25%. If all urban water pixels were considered, including those in shadow and layover regions, these figures fell to 58% and 19% respectively. These findings indicate that TerraSAR-X is capable of providing useful data for the calibration and validation of urban flood inundation models.

Rational modification of the hierarchy of intermolecular interactions in molecular crystal structures by using tunable halogen bonds

A family of 16 isomolecular salts (3-XpyH)(2)[MX'(4)] (3-XpyH=3-halopyridinium; M=Co, Zn; X=(F), Cl, Br, (I); X'=Cl, Br, I) each containing rigid organic cations and tetrahedral halometallate anions has been prepared and characterized by X-ray single crystal and/or powder diffraction. Their crystal structures reflect the competition and cooperation between non-covalent interactions: N-H center dot center dot center dot X'-M hydrogen bonds, C-X center dot center dot center dot X'-M halogen bonds and pi-pi stacking. The latter are essentially unchanged in strength across the series, but both halogen bonds and hydrogen bonds are modified in strength upon changing the halogens involved. Changing the organic halogen (X) from F to I strengthens the C-X center dot center dot center dot X'-M halogen bonds, whereas an analogous change of the inorganic halogen (X') weakens both halogen bonds and N-H center dot center dot center dot X'-M hydrogen bonds. By so tuning the strength of the putative halogen bonds from repulsive to weak to moderately strong attractive interactions, the hierarchy of the interactions has been modified rationally leading to systematic changes in crystal packing. Three classes of crystal structure are obtained. In type A (C-F center dot center dot center dot X'-M) halogen bonds are absent. The structure is directed by N-H center dot center dot center dot X'-M hydrogen bonds and pi-stacking interactions. In type B structures, involving small organic halogens (X) and large inorganic halogens (X'), long (weak) C-X center dot center dot center dot X'-M interactions are observed with type I halogen-halogen interaction geometries (C-X center dot center dot center dot X' approximate to X center dot center dot center dot X'-M approximate to 155 degrees), but hydrogen bonds still dominate. Thus, minor but quite significant perturbations from the type A structure arise. In type C, involving larger organic halogens (X) and smaller inorganic halogens (X'), stronger halogen bonds are formed with a type II halogen-halogen interaction geometry (C-X center dot center dot center dot X' approximate to 180 degrees; X center dot center dot center dot X'-M approximate to 110 degrees) that is electrostatically attractive. The halogen bonds play a major role alongside hydrogen bonds in directing the type C structures, which as a result are quite different from type A and B.

Fe-sulphate-rich evaporative mineral precipitates from the Rio Tinto, southwest Spain

The soluble metal sulphate salts melanterite, rozenite, rhomboclase, szornolnokite, copiapite, coquimbite, hexahydrite and halotrichite, together with gypsum, have been identified, some for the first time oil the banks of the Rio Tinto, SW Spain. Secondary Fe-sulphate minerals call form directly from evaporating acid, sulphate-rich Solutions as a result of pyrite oxidation. Chemical analyses of mixtures of these salt minerals indicate concentrations of Fe (up to 31 wt.%), Mg (up to 4 wt.%), Cu (up to 2 wt.%) and Zn (up to wt.%). These minerals are shown to act as transient storage Cor metals and can store on average up to 10% (9.5 - 11%) and 22% (20-23%) Zn and Cu respectively, of the total discharge of the Rio Tinto during the summer period. Melanterite and rozenite precipitates at Rio Tinto are only found in association with very acidic drainage waters (pH <1.0) draining directly from pyritic waste piles. Copiapite precipitates abundantly oil the banks of the Rio Tinto by (1) direct evaporation of the river water; or (2) as part of a paragenetic sequence with the inclusion of minor halotrichite, indicating natural dehydration and decomposition. The natural occurrences are comparable with the process of paragenesis from the evaporation of Rio Tinto river water under laboratory experiments resulting in the formation of aluminocopiapite, halotrichite, coquimbite, voltaite and gypsum.

Impact of sewage sludge applications on the biogeochemistry of soils

This report describes an investigation into the bioavailability and fate of trace metals and their subsequent impact on important soil microbiological functions such as nitrification, denitrification and methane oxidation in low and high Cu containing soils in the presence and absence of residual organic matter from sewage sludge additions made 10 years earlier. The soils being studied are part of long term sewage sludge trials and include a low Cu soil ( 13.3mg Cu/ kg soil, 4.18 LOI %), left un- amended to serve as a control soil, soil amended with a high Cu sewage sludge ( 278.3mg Cu/ kg soil, 6.52 LOI %) and soil amended with a low Cu sewage sludge ( 46.3mg Cu/ kg soil, 6.18 LOI %). Soil was also amended with inorganic metal salts ( 273.4mg Cu/ kg soil, 4.52 LOI %) to further investigate the impact of Cu in the absence of additional organic matter contained in applied sewage sludge. Data from the first two years of a project are presented which has included field- based studies at long term sewage sludge trials based in Watlington, Oxford, UK and laboratory based studies at the Institute of Grassland & Environmental Research, North Wyke, Devon, UK.

Influence of aluminum hydroxide on potassium release from two Colombian soils

The effect of sesquioxides on the mechanisms of chemical reactions that govern the transformation between exchangeable potassium (Kex) and non-exchangeable K (Knex) was studied on acid tropical soils from Colombia: Caribia with predominantly 2 : 1 clay minerals and High Terrace with predominantly 1 : 1 clay minerals and sesquioxides. Illite and vermiculite are the main clay minerals in Caribia followed by kaolinite, gibbsite, and plagioclase, and kaolinite is the major clay mineral in High Terrace followed by hydroxyl-Al interlayered vermiculite, quartz, and pyrophyllite. The soils have 1.8 and 0.5% of K2O, respectively. They were used either untreated or prepared by adding AlCl3 and NaOH, which produced aluminum hydroxide. The soils were percolated continuously with 10mM NH4OAc at pH 7.0 and 10 mM CaCl2 at pH 5.8 for 120 h at 6 mL h(-1) to examine the release of Kex and Knex. In the untreated soils, NH4+ and Ca-2(+) released the same amounts of Kex from Caribia, whereas NH4+ released about twice as much Kex as Ca2+ from High Terrace. This study proposes that the small ionic size of NH4+ (0.54nm) enables it to enter more easily into the K sites at the broken edges of the kaolinite where Ca2+ (0.96 nm) cannot have access. As expected for a soil dominated by 2 : 1 clay minerals, Ca2+ caused Knex to be released from Caribia with no release by NH4+. No Knex was released by either ion from High Terrace. After treatment with aluminum hydroxide, K release from the exchangeable fraction was reduced in Caribia due to the blocking of the exchange sites but release of Knex was not affected. The treatment increased the amount of Kex released from the High Terrace soil and the release of Knex remained negligible although with Ca2+ the distinction between Kex and Knex was unclear. The increase in Kex was attributed to the initially acidic conditions produced by adding AlCl3 which may have dissolved interlayered aluminum hydroxide from the vermiculite present, thus exposing trapped K as exchangeable K. The subsequent precipitation of aluminum hydroxide when NaOH was added did not interfere with the release of this K, and so was probably formed mostly on the surface of the dominant kaolinite. Measurement of availability of K by standard methods using NH4 salts could result in overestimates in High Terrace and this may be a more general shortcoming of the methods in kaolinitic soils.

Extraction of tidal channel networks from aerial photographs alone and combined with laser altimetry

Tidal channel networks play an important role in the intertidal zone, exerting substantial control over the hydrodynamics and sediment transport of the region and hence over the evolution of the salt marshes and tidal flats. The study of the morphodynamics of tidal channels is currently an active area of research, and a number of theories have been proposed which require for their validation measurement of channels over extensive areas. Remotely sensed data provide a suitable means for such channel mapping. The paper describes a technique that may be adapted to extract tidal channels from either aerial photographs or LiDAR data separately, or from both types of data used together in a fusion approach. Application of the technique to channel extraction from LiDAR data has been described previously. However, aerial photographs of intertidal zones are much more commonly available than LiDAR data, and most LiDAR flights now involve acquisition of multispectral images to complement the LiDAR data. In view of this, the paper investigates the use of multispectral data for semiautomatic identification of tidal channels, firstly from only aerial photographs or linescanner data, and secondly from fused linescanner and LiDAR data sets. A multi-level, knowledge-based approach is employed. The algorithm based on aerial photography can achieve a useful channel extraction, though may fail to detect some of the smaller channels, partly because the spectral response of parts of the non-channel areas may be similar to that of the channels. The algorithm for channel extraction from fused LiDAR and spectral data gives an increased accuracy, though only slightly higher than that obtained using LiDAR data alone. The results illustrate the difficulty of developing a fully automated method, and justify the semi-automatic approach adopted.

LiDAR mapping of tidal marshes for ecogeomorphological modelling in the TIDE project

The European research project TIDE (Tidal Inlets Dynamics and Environment) is developing and validating coupled models describing the morphological, biological and ecological evolution of tidal environments. The interactions between the physical and biological processes occurring in these regions requires that the system be studied as a whole rather than as separate parts. Extensive use of remote sensing including LiDAR is being made to provide validation data for the modelling. This paper describes the different uses of LiDAR within the project and their relevance to the TIDE science objectives. LiDAR data have been acquired from three different environments, the Venice Lagoon in Italy, Morecambe Bay in England, and the Eden estuary in Scotland. LiDAR accuracy at each site has been evaluated using ground reference data acquired with differential GPS. A semi-automatic technique has been developed to extract tidal channel networks from LiDAR data either used alone or fused with aerial photography. While the resulting networks may require some correction, the procedure does allow network extraction over large areas using objective criteria and reduces fieldwork requirements. The networks extracted may subsequently be used in geomorphological analyses, for example to describe the drainage patterns induced by networks and to examine the rate of change of networks. Estimation of the heights of the low and sparse vegetation on marshes is being investigated by analysis of the statistical distribution of the measured LiDAR heights. Species having different mean heights may be separated using the first-order moments of the height distribution.

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.