Multiscale divergence between Landsat- and lidar-based biomass mapping is related to regional variation in canopy cover and composition

Background

Satellite-based aboveground forest biomass maps commonly form the basis of forest biomass and carbon stock mapping and monitoring, but biomass maps likely vary in performance by region and as a function of spatial scale of aggregation. Assessing such variability is not possible with spatially-sparse vegetation plot networks. In the current study, our objective was to determine whether high-resolution lidar-based and moderate-resolution Landsat-base aboveground live forest biomass maps converged on similar predictions at stand- to landscape-levels (10 s to 100 s ha) and whether such differences depended on biophysical setting. Specifically, we examined deviations between lidar- and Landsat-based biomass mapping methods across scales and ecoregions using a measure of error (normalized root mean square deviation), a measure of the unsystematic deviations, or noise (Pearson correlation coefficient), and two measures related to systematic deviations, or biases (intercept and slope of a regression between the two sets of predictions).

Results

Compared to forest inventory data (0.81-ha aggregate-level), lidar and Landsat-based mean biomass predictions exhibited similar performance, though lidar predictions exhibited less normalized root mean square deviation than Landsat when compared with the reference plot data. Across aggregate-levels, the intercepts and slopes of regression equations describing the relationships between lidar- and Landsat-based biomass predictions stabilized (i.e., little additional change with increasing area of aggregates) at aggregate-levels between 10 and 100 ha, suggesting a consistent relationship between the two maps at landscape-scales. Differences between lidar- and Landsat-based biomass maps varied as a function of forest canopy heterogeneity and composition, with systematic deviations (regression intercepts) increasing with mean canopy cover and hardwood proportion within forests and correlations decreasing with hardwood proportion.

Conclusions

Deviations between lidar- and Landsat-based maps indicated that satellite-based approaches may represent general gradients in forest biomass. Ecoregion impacted deviations between lidar and Landsat biomass maps, highlighting the importance of biophysical setting in determining biomass map performance across aggregate scales. Therefore, regardless of the source of remote sensing (e.g., Landsat vs. lidar), factors affecting the measurement and prediction of forest biomass, such as species composition, need to be taken into account whether one is estimating biomass at the plot, stand, or landscape scale.

     

Retreat dynamics of the eastern sector of the British–Irish Ice Sheet during the last glaciation

The findings of BRITICE-CHRONO Transect 2 through the North Sea Basin and eastern England are reported. We define ice-sheet marginal oscillation between ~31 and 16 ka, with seven distinctive former ice-sheet limits (L1–7) constrained by Bayesian statistical analysis. The southernmost limit of the North Sea Lobe is recorded by the Bolders Bank Formation (L1; 25.8–24.6 ka). L2 represents ice-sheet oscillation and early retreat to the northern edge of the Dogger Bank (23.5–22.2 ka), with the Garret Hill Moraine in north Norfolk recording a significant regional readvance to L3 at 21.5–20.8 ka. Ice-marginal oscillations at ~26–21 ka resulted in L1, L2 and L3 being partially to totally overprinted. Ice-dammed lakes related to L1–3, including Lake Humber, are dated at 24.1–22.3 ka. Ice-sheet oscillation and retreat from L4 to L5 occurred between 19.7 and 17.3 ka, with grounding zone wedges marking an important transition from terrestrial to marine tidewater conditions, triggered by the opening of the Dogger Lake spillway between 19.9 and 17.5 ka. L6 relates to ice retreat under glacimarine conditions and final ice retreat into the Firth of Forth by 15.8 ka. L7 (~15 ka) represents an ice retreat from Bosies Bank into the Moray Firth.     

Nanogeochemistry of hydrothermal magnetite

Magnetite from hydrothermal ore deposits can contain up to tens of thousands of parts per million (ppm) of elements such as Ti, Si, V, Al, Ca, Mg, Na, which tend to either structurally incorporate into growth and sector zones or form mineral micro- to nano-sized particles. Here, we report micro- to nano-structural and chemical data of hydrothermal magnetite from the Los Colorados iron oxide–apatite deposit in Chile, where magnetite displays both types of trace element incorporation. Three generations of magnetites (X–Z) were identified with concentrations of minor and trace elements that vary significantly: SiO₂, from below detection limit (bdl) to 3.1 wt%; Al₂O₃, 0.3–2.3 wt%; CaO, bdl–0.9 wt%; MgO, 0.02–2.5 wt%; TiO₂, 0.1–0.4 wt%; MnO, 0.04–0.2 wt%; Na₂O, bdl–0.4 wt%; and K₂O, bdl–0.4 wt%. An exception is V₂O₃, which is remarkably constant, ranging from 0.3 to 0.4 wt%. Six types of crystalline nanoparticles (NPs) were identified by means of transmission electron microscopy in the trace element-rich zones, which are each a few micrometres wide: (1) diopside, (2) clinoenstatite; (3) amphibole, (4) mica, (5) ulvöspinel, and (6) Ti-rich magnetite. In addition, Al-rich nanodomains, which contain 2–3 wt% of Al, occur within a single crystal of magnetite. The accumulation of NPs in the trace element-rich zones suggest that they form owing to supersaturation from a hydrothermal fluid, followed by entrapment during continuous growth of the magnetite surface. It is also concluded that mineral NPs promote exsolution of new phases from the mineral host, otherwise preserved as structurally bound trace elements. The presence of abundant mineral NPs in magnetite points to a complex incorporation of trace elements during growth, and provides a cautionary note on the interpretation of micron-scale chemical data of magnetite.     

Lid tectonics-Preface

正The idea that plate tectonics may not have operated deep in Earth's Precambrian past has a long legacy.What predated plate tectonics is unknown,and advances in data-from geochemical,geological and tectonic,to paleomagnetic,as well as modelling ap-     

The role of horizontal resolution in simulating drivers of the global hydrological cycle

The role of atmospheric general circulation model (AGCM) horizontal resolution in representing the global energy budget and hydrological cycle is assessed, with the aim of improving the understanding of model uncertainties in simulating the hydrological cycle. We use two AGCMs from the UK Met Office Hadley Centre: HadGEM1-A at resolutions ranging from 270 to 60 km, and HadGEM3-A ranging from 135 to 25 km. The models exhibit a stable hydrological cycle, although too intense compared to reanalyses and observations. This over-intensity is explained by excess surface shortwave radiation, a common error in general circulation models (GCMs). This result is insensitive to resolution. However, as resolution is increased, precipitation decreases over the ocean and increases over the land. This is associated with an increase in atmospheric moisture transport from ocean to land, which changes the partitioning of moisture fluxes that contribute to precipitation over land from less local to more non-local moisture sources. The results start to converge at 60-km resolution, which underlines the excessive reliance of the mean hydrological cycle on physical parametrization (local unresolved processes) versus model dynamics (large-scale resolved processes) in coarser HadGEM1 and HadGEM3 GCMs. This finding may be valid for other GCMs, showing the necessity to analyze other chains of GCMs that may become available in the future with such a range of horizontal resolutions. Our finding supports the hypothesis that heterogeneity in model parametrization is one of the underlying causes of model disagreement in the Coupled Model Intercomparison Project (CMIP) exercises.     

Generalized Radial Flow in Synthetic Flow Systems

Traditional analysis methods used to determine hydraulic properties from pumping tests work well in many porous media aquifers, but they often do not work in heterogeneous and fractured‐rock aquifers, producing non‐plausible and erroneous results. The generalized radial flow model developed by Barker (1988) can reveal information about heterogeneity characteristics and aquifer geometry from pumping test data by way of a flow dimension parameter. The physical meaning of non‐integer flow dimensions has long been a subject of debate and research. We focus on understanding and interpreting non‐radial flow through high permeability conduits within fractured aquifers. We develop and simulate flow within idealized non‐radial flow conduits and expand on this concept by simulating pumping in non‐fractal random fields with specific properties that mimic persistent sub‐radial flow responses. Our results demonstrate that non‐integer flow dimensions can arise from non‐fractal geometries within aquifers. We expand on these geometric concepts and successfully simulate pumping in random fields that mimic well‐test responses seen in the Culebra Dolomite above the Waste Isolation Pilot Plant.     

The onset of magnetoconvection at large Prandtl number in a rotating layer I. Finite magnetic diffusion

Abstract

This paper develops further a convection model that has been studied several times previously as a very crude idealization of planetary core dynamics. A plane layer of electrically-conducting fluid rotates about the vertical in the presence of a magnetic field. Such a field can be created spontaneously, as in the Childress—Soward dynamo, but here it is uniform, horizontal and externally-applied. The Prandtl number of the fluid is large, but the Ekman, Elsasser and Rayleigh numbers are of order unity, as is the ratio of thermal to magnetic diffusivity. Attention is focused on the onset of convection as the temperature difference applied across the layer is increased, and on the preferred mode, i.e., the planform and time-dependence of small amplitude convection. The case of main interest is the layer confined between electrically-insulating no-slip walls, but the analysis is guided by a parallel study based on illustrative boundary conditions that are mathematically simpler.     

Modeling fractures as interfaces with nonmatching grids

We consider a model for fluid flow in a porous medium with a fracture. In this model, the fracture is treated as an interface between subdomains, on which specific equations have to be solved. In this article, we analyze the discrete problem, assuming that the fracture mesh and the subdomain meshes are completely independent, but that the geometry of the fracture is respected. We show that despite this nonconformity, first-order convergence is preserved with the lowest-order Raviart–Thomas(-Nedelec) mixed finite elements. Numerical simulations confirm this result.     

Buckland,Darwin and the attempted recognition of an Ice Age in Wales, 1837–1842

The concept of a former Ice Age was introduced to Britain by Agassiz, first, through Buckland in 1838 and then by his tour of Britain in 1840. The reception was mixed due to the Iceberg theory, which was held by Darwin, Lyell and Murchison and others. After 1840, Murchison looked for a compromise between Glaciers and Icebergs and this came in the work of Bowman and Buckland in 1841 and Darwin during 1842 in Snowdonia and the Marches. There were three geologists visiting Wales, all familiar with glaciation; Bowman failed to find any glaciation and Buckland and Darwin, who identified both alpine-glacier and “ice-berg” glaciation and reinterpreted their previous work. Thus both a Catastrophist and a Uniformitarian came to similar conclusions, but it was several decades before a consensus was found, which was delayed by Darwin's emphasis on submergence.