Impact of data model and point density on aboveground forest biomass estimation from airborne LiDAR

Background

Accurate estimation of aboveground forest biomass (AGB) and its dynamics is of paramount importance in understanding the role of forest in the carbon cycle and the effective implementation of climate change mitigation policies. LiDAR is currently the most accurate technology for AGB estimation. LiDAR metrics can be derived from the 3D point cloud (echo-based) or from the canopy height model (CHM). Different sensors and survey configurations can affect the metrics derived from the LiDAR data. We evaluate the ability of the metrics derived from the echo-based and CHM data models to estimate AGB in three different biomes, as well as the impact of point density on the metrics derived from them.

Results

Our results show that differences among metrics derived at different point densities were significantly different from zero, with a larger impact on CHM-based than echo-based metrics, particularly when the point density was reduced to 1 point m^?2. Both data models-echo-based and CHM-performed similarly well in estimating AGB at the three study sites. For the temperate forest in the Sierra Nevada Mountains, California, USA, R² ranged from 0.79 to 0.8 and RMSE (relRMSE) from 69.69 (35.59%) to 70.71 (36.12%) Mg ha^?1 for the echo-based model and from 0.76 to 0.78 and 73.84 (37.72%) to 128.20 (65.49%) Mg ha^?1 for the CHM-based model. For the moist tropical forest on Barro Colorado Island, Panama, the models gave R² ranging between 0.70 and 0.71 and RMSE between 30.08 (12.36%) and 30.32 (12.46) Mg ha^?1 [between 0.69–0.70 and 30.42 (12.50%) and 61.30 (25.19%) Mg ha^?1] for the echo-based [CHM-based] models. Finally, for the Atlantic forest in the Sierra do Mar, Brazil, R² was between 0.58–0.69 and RMSE between 37.73 (8.67%) and 39.77 (9.14%) Mg ha^?1 for the echo-based model, whereas for the CHM R² was between 0.37–0.45 and RMSE between 45.43 (10.44%) and 67.23 (15.45%) Mg ha^?1.

Conclusions

Metrics derived from the CHM show a higher dependence on point density than metrics derived from the echo-based data model. Despite the median of the differences between metrics derived at different point densities differing significantly from zero, the mean change was close to zero and smaller than the standard deviation except for very low point densities (1 point m^?2). The application of calibrated models to estimate AGB on metrics derived from thinned datasets resulted in less than 5% error when metrics were derived from the echo-based model. For CHM-based metrics, the same level of error was obtained for point densities higher than 5 points m^?2. The fact that reducing point density does not introduce significant errors in AGB estimates is important for biomass monitoring and for an effective implementation of climate change mitigation policies such as REDD + due to its implications for the costs of data acquisition. Both data models showed similar capability to estimate AGB when point density was greater than or equal to 5 point m^?2.

     

The large earthquake of 8 August 1303 in Crete: seismic scenario and tsunami in the Mediterranean area

By conducting a historical review of this large seismic event in the Mediterranean, it has been possible to identify both the epicentral area and the area in which its effects were principally felt. Ever since the nineteenth century, the seismological tradition has offered a variety of partial interpretations of the earthquake, depending on whether the main sources used were Arabic, Greek or Latin texts. Our systematic research has involved the analysis not only of Arab, Byzantine and Italian chronicle sources, but also and in particular of a large number of never previously used official and public authority documents, preserved in Venice in the State Archive, in the Marciana National Library and in the Library of the Museo Civico Correr. As a result, it has been possible to establish not only chronological parameters for the earthquake (they were previously uncertain) but also its overall effects (epicentral area in Crete, Imax XI MCS). Sources containing information in 41 affected localities and areas were identified. The earthquake also gave rise to a large tsunami, which scholars have seen as having certain interesting elements in common with that of 21 July 365, whose epicentre was also in Crete. As regards methodology, this research made it clear that knowledge of large historical earthquakes in the Mediterranean is dependent upon developing specialised research and going beyond the territorial limits of current national catalogues.     

Determination of Lithium, Beryllium and Boron at Trace Levels by Laser Ablation-Inductively Coupled Plasma-Sector Field Mass Spectrometry

The analytical capabilities of laser ablation (LA)-ICP-MS in determining Li, Be and B at trace levels in geological samples have been tested on a series of glass reference materials and natural samples. The LA-ICP-MS instrument used consisted of a sector-field ICP-MS coupled with a laser ablation microprobe operating at either 266 or 213 nm wavelength. Reference glasses from NIST (SRM 612, 614 and 616) and MPI-DING (KL2-G, ML3B-G, StHs6/80-G, GOR128-G, GOR132-G, T1-G and ATHO-G) were selected to develop the analytical method and to assess the best instrumental configuration. A series of calcic amphiboles with different Li, Be and B concentrations were also analysed using both LA-ICP-MS and SIMS to test the applicability of the method to natural minerals. Results indicated that with a spot size of 40 μm the agreement between measured and reference values of Li, Be and B is generally better than 10% for NIST SRM 612 and 20% for NIST SRM 614. Average reproducibility at the 2s level was 10% for Li, 20% for Be and 15% for B. Limits of detection were approximately 100 ng g^-1 for Be and B and 200 ng g^-1 for Li. These results were confirmed by analyses carried out on natural amphiboles and compared well in terms of precision and accuracy with those commonly achieved by SIMS.     

A thermo-plastic constitutive law for brittle-plastic behavior of rocks at high temperatures

Mechanical properties of rocks change under the influence of, temperature. Stress at the onset of yielding, ultimate strength, dilatancy, strain hardening and softening, and the confining pressure at brittle-ductile transition are all reduced by the increasing temperature. This study presents a framework of constitutive modeling of thermo-brittle-plastic behavior of rocks which encompasses these changes. The constitutive law is based on a thermo-plasticity theory first proposed for metals byPrager (1958). Two phenomenological mechanisms have been identified as central for the modeling: temperature dependence of the yield locus (thermal softening), and temperature dependence of the strain-hardening function (thermally enhanced ductility). Material parameters for two rocks, Carrara marble and Westerly granite, were determined on the basis of additional hypotheses. These parameters are used in numerical simulations of triaxial tests at different temperatures. The obtained stress-strain curves compare well to the experimental results. The changes with temperature in the stress at the onset of yielding are more accurately reproduced that the evolution of hardening or softening. Suggestions for possible improvements and future research directions are indicated.     

Assessment of calanchi and rill–interrill erosion susceptibilities using terrain analysis and geostochastics: A case study in the Oltrepo Pavese,Northern Apennines,Italy

Soil erosion is one of the most important environmental problems distributed worldwide. In the last decades, numerous studies have been published on the assessment of soil erosion and the related processes and forms using empirical, conceptual and physically based models. For the prediction of the spatial distribution, more and more sophisticated stochastic modelling approaches have been proposed – especially on smaller spatial scales such as river basins. In this work, we apply a maximum entropy model (MaxEnt) to evaluate badlands (calanchi) and rill–interrill (sheet erosion) areas in the Oltrepo Pavese (Northern Apennines, Italy). The aim of the work is to assess the important environmental predictors that influence calanchi and rill–interrill erosion at the regional scale. We used 13 topographic parameters derived from a 12 m digital elevation model (TanDEM-X) and data on the lithology and land use. Additional information about the vegetation is introduced through the normalized difference vegetation index based on remotely sensed data (ASTER images). The results are presented in the form of susceptibility maps showing the spatial distribution of the occurrence probability for calanchi and rill–interrill erosion. For the validation of the MaxEnt model results, a support vector machine approach was applied. The models show reliable results and highlight several locations of the study area that are potentially prone to future soil erosion. Thus, coping and mitigation strategies may be developed to prevent or fight the soil erosion phenomenon under consideration. © 2020 John Wiley & Sons, Ltd.     

Uncertainty quantification of overpressure buildup through inverse modeling of compaction processes in sedimentary basins

This study illustrates a procedure conducive to a preliminary risk analysis of overpressure development in sedimentary basins characterized by alternating depositional events of sandstone and shale layers. The approach rests on two key elements: (1) forward modeling of fluid flow and compaction, and (2) application of a model-complexity reduction technique based on a generalized polynomial chaos expansion (gPCE). The forward model considers a one-dimensional vertical compaction processes. The gPCE model is then used in an inverse modeling context to obtain efficient model parameter estimation and uncertainty quantification. The methodology is applied to two field settings considered in previous literature works, i.e. the Venture Field (Scotian Shelf, Canada) and the Navarin Basin (Bering Sea, Alaska, USA), relying on available porosity and pressure information for model calibration. It is found that the best result is obtained when porosity and pressure data are considered jointly in the model calibration procedure. Uncertainty propagation from unknown input parameters to model outputs, such as pore pressure vertical distribution, is investigated and quantified. This modeling strategy enables one to quantify the relative importance of key phenomena governing the feedback between sediment compaction and fluid flow processes and driving the buildup of fluid overpressure in stratified sedimentary basins characterized by the presence of low-permeability layers. The results here illustrated (1) allow for diagnosis of the critical role played by the parameters of quantitative formulations linking porosity and permeability in compacted shales and (2) provide an explicit and detailed quantification of the effects of their uncertainty in field settings.     

Origin of peraluminous granites and granodiorites, Iberian massif, Spain: an experimental test of granite petrogenesis

The discrimination between potential source materials involved in the genesis of Iberian granites and granodiorites, as well as the role of mantle-crust interactions, are examined using constraints imposed by melting experiments, melting-assimilation experiments and Sr-Nd isotope systematics. The Sr-Nd isotope relationships indicate the existence of different genetic trends in which juvenile mantle materials are involved by different mechanisms: (1) a source trend, traced by a particular evolution of the pre-Hercynian basement and indicating mantle participation at the time of sedimentation; (2) a set of magmatic trends traced by gabbro-tonalite-enclave-granodiorite associations, implying the incorporation of new mantle material at the time of granite generation. These relationships strongly support a pure crustal origin for the peraluminous leucogranites, derived from partial melting of crustal protoliths, and a hybrid origin for the peraluminous granodiorites. These granodiorites are the most abundant granitic rocks of the Central Iberian zone (CIZ) of the Iberian massif, implying that processes of hybridisation by assimilation and/or magma mixing played an important role in granitoid production during the Hercynian orogeny. These hypotheses have been tested by means of melting and assimilation experiments. Melting experiments in the range 800–900 °C and at pressures of 3, 6, 10 and 15 kbar indicate that: (1) several potential source materials such as Bt-Ms gneisses and metagreywackes are suitable for the production of peraluminous leucogranite melts; (2) the melt compositions are always leucogranitic, regardless of pressure; (3) pressure exerts a strong influence on the fertility of the source: experiments at 3 kbar produce more than 20 vol% of melt, compared with less than 5 vol% of melt produced at 10 and 15 kbar and at the same temperature. The melting-assimilation experiments carried out at 1000 °C and 4, 7 and 10 kbar and using a proportion of 50% gabbro and 50% gneiss give high melt proportions (more than 50 vol.%) and noritic residues. These melts have the composition of leucogranodiorites, and overlap with part of the compositional range of peraluminous granodiorites of the Iberian massif. The generation of more mafic granodiorites may be explained by the incorporation of some residual orthopyroxene to the granodiorite magmas. The low solubility of Fe + Mg prevents the generation of granodiorite melts with more than 3 wt% of MgO + FeO at all crustal pressures. The large volumes of peraluminous, hybrid granodiorites, produced by assimilation of crustal rocks by mantle magmas, imply that an important episode of crustal growth took place during the Late-Palaeozoic Hercynian orogeny in the Iberian massif. Received: 30 June 1998 / Accepted: 27 November 1998     

The amphipod Talitrus saltator as a bioindicator of human trampling on sandy beaches

The present paper assesses the use of the supralittoral amphipod Talitrus saltator as a bioindicator of the effects of human trampling on the supralittoral sandy band. Samplings in delimited areas were carried out at sites subjected to different human impact. The results showed a strong negative correlation between the number of swimmers and the sandhopper population density, while there was no clear relationship between sandhopper abundance and the other factors considered: granulometry, compactness and organic carbon content of the sand, and trace metal contents in the sand and sandhoppers. A field test of trampling conducted in a confined space showed its direct negative effect on sandhopper survival. However, trace metal analysis confirmed the ability of T. saltator to bioaccumulate some elements (Hg, Zn, Cu, Cd). Our study demonstrates that T. saltator is a good bioindicator of human impact in the supralittoral zone of sandy shores.     

Global change and relative sea level rise at Venice: what impact in term of flooding

Relative sea level rise (RSLR) due to climate change and geodynamics represents the main threat for the survival of Venice, emerging today only 90 cm above the Northern Adriatic mean sea level (msl). The 25 cm RSLR occurred over the 20th century, consisting of about 12 cm of land subsidence and 13 cm of sea level rise, has increased the flood frequency by more than seven times with severe damages to the urban heritage. Reasonable forecasts of the RSLR expected to the century end must be investigated to assess the suitability of the Mo.S.E. project planned for the city safeguarding, i.e., the closure of the lagoon inlets by mobile barriers. Here we consider three RSLR scenarios as resulting from the past sea level rise recorded in the Northern Adriatic Sea, the IPCC mid-range A1B scenario, and the expected land subsidence. Available sea level measurements show that more than 5 decades are required to compute a meaningful eustatic trend, due to pseudo-cyclic 7–8 year long fluctuations. The period from 1890 to 2007 is characterized by an average rate of 0.12 ± 0.01 cm/year. We demonstrate that linear regression is the most suitable model to represent the eustatic process over these 117 year. Concerning subsidence, at present Venice is sinking due to natural causes at 0.05 cm/year. The RSLR is expected to range between 17 and 53 cm by 2100, and its repercussions in terms of flooding frequency are associated here to each scenario. In particular, the frequency of tides higher than 110 cm, i.e., the value above which the gates would close the lagoon to the sea, will increase from the nowadays 4 times per year to a range between 20 and 250. These projections provide a large spread of possible conditions concerning the survival of Venice, from a moderate nuisance to an intolerable aggression. Hence, complementary solutions to Mo.S.E. may well be investigated.