Estimating urban above ground biomass with multi-scale LiDAR

Background

Urban trees have long been valued for providing ecosystem services (mitigation of the “heat island” effect, suppression of air pollution, etc.); more recently the potential of urban forests to store significant above ground biomass (AGB) has also be recognised. However, urban areas pose particular challenges when assessing AGB due to plasticity of tree form, high species diversity as well as heterogeneous and complex land cover. Remote sensing, in particular light detection and ranging (LiDAR), provide a unique opportunity to assess urban AGB by directly measuring tree structure. In this study, terrestrial LiDAR measurements were used to derive new allometry for the London Borough of Camden, that incorporates the wide range of tree structures typical of an urban setting. Using a wall-to-wall airborne LiDAR dataset, individual trees were then identified across the Borough with a new individual tree detection (ITD) method. The new allometry was subsequently applied to the identified trees, generating a Borough-wide estimate of AGB.

Results

Camden has an estimated median AGB density of 51.6 Mg ha^–1 where maximum AGB density is found in pockets of woodland; terrestrial LiDAR-derived AGB estimates suggest these areas are comparable to temperate and tropical forest. Multiple linear regression of terrestrial LiDAR-derived maximum height and projected crown area explained 93% of variance in tree volume, highlighting the utility of these metrics to characterise diverse tree structure. Locally derived allometry provided accurate estimates of tree volume whereas a Borough-wide allometry tended to overestimate AGB in woodland areas. The new ITD method successfully identified individual trees; however, AGB was underestimated by ≤?25% when compared to terrestrial LiDAR, owing to the inability of ITD to resolve crown overlap. A Monte Carlo uncertainty analysis identified assigning wood density values as the largest source of uncertainty when estimating AGB.

Conclusion

Over the coming century global populations are predicted to become increasingly urbanised, leading to an unprecedented expansion of urban land cover. Urban areas will become more important as carbon sinks and effective tools to assess carbon densities in these areas are therefore required. Using multi-scale LiDAR presents an opportunity to achieve this, providing a spatially explicit map of urban forest structure and AGB.

     

European floods during the winter 1783/1784: scenarios of an extreme event during the ‘Little Ice Age’

The Lakagígar eruption in Iceland during 1783 was followed by the severe winter of 1783/1784, which was characterised by low temperatures, frozen soils, ice-bound watercourses and high rates of snow accumulation across much of Europe. Sudden warming coupled with rainfall led to rapid snowmelt, resulting in a series of flooding phases across much of Europe. The first phase of flooding occurred in late December 1783–early January 1784 in England, France, the Low Countries and historical Hungary. The second phase at the turn of February–March 1784 was of greater extent, generated by the melting of an unusually large accumulation of snow and river ice, affecting catchments across France and Central Europe (where it is still considered as one of the most disastrous known floods), throughout the Danube catchment and in southeast Central Europe. The third and final phase of flooding occurred mainly in historical Hungary during late March and early April 1784. The different impacts and consequences of the above floods on both local and regional scales were reflected in the economic and societal responses, material damage and human losses. The winter of 1783/1784 can be considered as typical, if severe, for the Little Ice Age period across much of Europe.     

Transient Divergent Flow and Transport in an Infinite Anisotropic Porous Formation

When seeking to predict plume geometry resulting from fluid injection through partially penetrating wells, it is common to assume a steady-state spherically diverging flow field. In reality, the flow field is transient. The steady-flow assumption is likely to cause overestimation of injection plume radius since the accommodation of fluid by increases in porosity and fluid density is ignored. In this paper, a transient solution is developed, resulting in a nonlinear ordinary differential equation expressing plume radius as a function of time. It is shown that the problem can be fully described by one type curve. A critical time, t_c, is identified at which the percentage error of the steady-state flow solution compared to the fully dynamic problem is less than 1%. Only for large injection rates and low permeabilities, does t_c become greater than 1 h. Nevertheless, an improved approximate solution is obtained by a simple linearization procedure. The critical time, t_c for the new approximate solution is 0.3% of that required for the steady-state flow solution.     

Monitoring the regional and temporal variability of winter snowfall in the arid Andes using digital NOAA/AVHRR data

Abstract

This paper deals with the spatial distribution and the temporal variability of snowfall in the most arid part of the Andes (18°‐ 28°S) during southern hemisphere winter (May‐September). As the official precipitation data is of poor quality, analyses were carried out by means of digital image processing techniques, using NOAA/AVHRR satellite‐data. Through analysis of 24 different snowfall events from six winters, a previously unknown spatial and temporal precipitation pattern in this remote and unexplored area was revealed. Snowfall is most abundant in the southernmost part of the research area and on the western side of the Andes, indicating the Pacific origin of the snowfall.

Nevertheless, the typical snowfall pattern is modified during different periods of the winter. Three typical time periods could be defined and distinguished from one another. Each of these three periods is characterized by typical weather conditions (cold fronts and “cut‐offs “) leading to a distinct snowfall pattern.

As this study is part of a broader paleoclimatic project, the results will serve as a basis for paleoclimatic reconstruction of past climate. Only by knowing the modern circulation and precipitation patterns is it possible to interpret paleoclimatic signals and archives found in the study area (e.g. paleosol, moraines) correctly.     

Processing and quality control of flux data during LITFASS-2003

Different aspects of the quality assurance and quality control (QA/QC) of micrometeorological measurements were combined to create a comprehensive algorithm which was then applied to experimental data from LITFASS-2003 (Lindenberg Inhomogeneous Terrain—Fluxes between Atmosphere and Surface: a long term Study). Eddy-covariance measurements of the latent heat flux were the main focus of the QA/QC efforts. The results of a turbulence sensor intercomparison experiment showed deviations between the different eddy-covariance systems on the order of 15%, or less than 30 W m⁻², for the latent heat flux and 5%, or less than 10 W m⁻², for the sensible heat flux. In order to avoid uncertainties due to the post-processing of turbulence data, a comprehensive software package was used for the analysis of experimental data from LITFASS-2003, including all necessary procedures for corrections and quality control. An overview of the quality test results shows that for most of the days more than 80% of the available latent heat flux data are of high quality so long as there are no instrumental problems. The representativeness of a flux value for the target land-use type was analysed using a stochastic footprint model. Different methods to calculate soil heat fluxes at the surface are discussed and a sensitivity analysis is conducted to select the most robust method for LITFASS-2003. The lack of energy balance closure, which was found for LITFASS-2003, can probably be attributed to the presence of low-frequency flux contributions that cannot be resolved with an averaging time of 30 min. Though the QA/QC system has been developed for the requirements of LITFASS-2003, it can also be applied to other experiments dealing with similar objectives.     

On the breakdown of zircon upon “dry” thermal annealing

Zircon samples without and with secondary chemical alteration from diverse sources were subjected to heat treatment at 1400 °C for 96 h. Resulting new phases and textures suggest that decomposition of zircon into component oxides occurred in all experiments to various degrees. The crucible material was found to have a strong influence on the extent of breakdown, especially in the case of altered starting materials. In this study the progressive stages of the breakdown of zircon grains are described. The factors that may govern the decomposition are discussed, including radiation damage, secondary alteration and external reaction conditions (sample container, atmosphere). Alumina crucibles should generally be avoided in dry annealing of zircon, to minimise uncontrolled breakdown into oxides.     

Well placement optimization subject to realistic field development constraints

This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.