Synthetic RapidEye data used for the detection of area-based spruce tree mortality induced by bark beetles

Tree mortality caused by outbreaks of the bark beetle Ips typographus (L.) plays an important role in the natural dynamics of Norway spruce (Picea abies L.) stands, which could cause far-reaching changes in the occurrence and duration of vegetation phenology. Field-based early detection of tree disturbances is hampered by logistic, terrain, and technical shortcomings, and by the inability to continuously monitor disturbances over large areas. Despite achievements in remote mapping of bark-beetle-induced tree mortalities, early warning has been mostly unsuccessful mainly because of the lack of spectral sensitivity and discrepancies in definitions of field- and image-based disturbance classes. Here we applied a method based on inter-annual phenology of Norway spruce stands derived from synthetic multispectral data to part of the Bavarian Forest National Park in Germany. We fused temporally continuous Moderate Resolution Imaging Spectroradiometer and discrete RapidEye data using a flexible spatiotemporal data fusion method to achieve validated 8-day RapidEye-like composites of normalized difference vegetation index for 2011. We assumed that the dead trees delineated on 2012 aerial photographs were those in which bark beetle infestations were initiated in 2011. Samples were drawn with variable-sized buffering to represent the areas prone to infestations and their surroundings. We applied a conditional inference random forest to select the best image date among the entire 46 synthetic datasets to best discriminate between the core infestation patches and their surroundings from the subsequent year. Of the discrete time points identified, day 281 of the year represented the highest discrepancy between aerial image-based dead trees and their surroundings. Classification results were significantly correlated with beetle count data obtained using pheromone traps. Our method provided valuable information for management purposes and enabled wall-to-wall mapping of stands prone to infestation and its uncertainty. The results offer potential implications for rapid and cost-effective monitoring of bark beetle outbreaks using satellite data, which would be of great benefit for both management and research tasks.     

Lu-Hf garnet geochronology of eclogites from the Balma Unit (Pennine Alps): implications for Alpine paleotectonic reconstructions

Three samples of eclogite from the Balma Unit, an ophiolite sheet on top of the Monte Rosa Nappe in the Pennine Alps, were investigated in terms of their P-T evolution, geochemistry, and Lu-Hf geochronology. The paleogeographic origin of this unit is controversial (North Penninic vs. South Penninic). It has been interpreted as a piece of Late Cretaceous oceanic crust, on the basis of ca. 93 Ma U-Pb SHRIMP ages of synmagmatic zircon cores in an eclogite. Trace element and isotope data suggest a mid ocean ridge (MOR) rather than an intraplate or OIB setting for the protoliths of the eclogites. Electron microprobe analyses of representative garnets show typical prograde zoning profiles. Estimated peak metamorphic temperatures of 550–600 Cº most likely did not exceed the closure temperature of the Lu-Hf system. Hence, Lu-Hf ages most likely reflect garnet growth in the studied samples. To minimize inclusion effects on age determinations, a selective digestion procedure for garnet was applied, in which zircon and rutile inclusions are not dissolved. The ages obtained for three samples, 42.3 ± 0.6 Ma (MSWD: 0.47), 42 ± 1 Ma (MSWD: 3.0) and 45.5 ± 0.3 Ma (MSWD: 0.33), are younger than all Lu-Hf ages reported so far for South Penninic Units. Metamorphic zircon domains of the 42.3 Ma sample (PIS1) were previously dated by U-Pb SHRIMP at 40.4 ± 0.7 Ma, indicating that the growth of metamorphic zircon post-dated the onset of garnet growth.These new data put important constraints on the paleogeographic reconstruction of the Alps. The MORB character of the rocks, together with their previously published protolith age, imply that oceanic spreading was still taking place in the Late Cretaceous. This supports a North Penninic origin for our samples because plate tectonic models predict Cretaceous spreading in the North Penninic but not in the South Penninic Ocean. If the Balma Unit is indeed North Penninic, the new Lu-Hf data, in combination with published geochronological data, require that two independent subduction zones consumed the South and North Penninic oceans.     

Biodegradation of forest floor organic matter bound to minerals via different binding mechanisms

Mineral-associated organic matter (OM) represents a large reservoir of organic carbon (OC) in natural environments. The factors controlling the extent of the mineral-mediated OC stabilization, however, are poorly understood. The protection of OM against biodegradation upon sorption to mineral phases is assumed to result from the formation of strong bonds that limit desorption. To test this, we studied the biodegradation of OM bound to goethite (α-FeOOH), pyrophyllite, and vermiculite via specific mechanisms as estimated from OC uptake in different background electrolytes and operationally defined as ‘ligand exchange’, ‘Ca²⁺ bridging’, and ‘van der Waals forces’. Organic matter extracted from an Oa forest floor horizon under Norway spruce (Picea abies (L.) Karst) was reacted with minerals at dissolved OC concentrations of ∼5-130 mg/L at pH 4. Goethite retained up to 30.1 mg OC/g predominantly by ‘ligand exchange’; pyrophyllite sorbed maximally 12.5 mg OC/g, largely via ‘van der Waals forces’ and ‘Ca²⁺ bridging’, while sorption of OM to vermiculite was 7.3 mg OC/g, mainly due to the formation of ‘Ca²⁺ bridges’. Aromatic OM components were selectively sorbed by all minerals (goethite ? phyllosilicates). The sorption of OM was strongly hysteretic with the desorption into 0.01 M NaCl being larger for OM held by ‘Ca²⁺ bridges’ and ‘van der Waals forces’ than by ‘ligand exchange’. Incubation experiments under aerobic conditions (initial pH 4; 90 days) revealed that OM mainly bound to minerals by ‘ligand exchange’ was more resistant against mineralization than OM held by non-columbic interactions (‘van der Waals forces’). Calcium bridges enhanced the stability of sorbed OM, especially for vermiculite, but less than the binding via ‘ligand exchange’. Combined evidence suggests that the extent and rate of mineralization of mineral-associated OM are governed by desorption. The intrinsic stability of sorbed OM as related to the presence of resistant, lignin-derived aromatic components appears less decisive for the sorptive stabilization of OM than the involved binding mechanisms. In a given environment, the type of minerals present and the solution chemistry determine the operating binding mechanisms, thereby the extent of OM sorption and desorption, and thus ultimately the bioavailability of mineral-associated OM.     

Understanding the time‐varying importance of different uncertainty sources in hydrological modelling using global sensitivity analysis

Simulations from hydrological models are affected by potentially large uncertainties stemming from various sources, including model parameters and observational uncertainty in the input/output data. Understanding the relative importance of such sources of uncertainty is essential to support model calibration, validation and diagnostic evaluation and to prioritize efforts for uncertainty reduction. It can also support the identification of ‘disinformative data’ whose values are the consequence of measurement errors or inadequate observations. Sensitivity analysis (SA) provides the theoretical framework and the numerical tools to quantify the relative contribution of different sources of uncertainty to the variability of the model outputs. In traditional applications of global SA (GSA), model outputs are aggregations of the full set of a simulated variable. For example, many GSA applications use a performance metric (e.g. the root mean squared error) as model output that aggregates the distances of a simulated time series to available observations. This aggregation of propagated uncertainties prior to GSA may lead to a significant loss of information and may cover up local behaviour that could be of great interest. Time‐varying sensitivity analysis (TVSA), where the aggregation and SA are repeated at different time steps, is a viable option to reduce this loss of information. In this work, we use TVSA to address two questions: (1) Can we distinguish between the relative importance of parameter uncertainty versus data uncertainty in time? (2) Do these influences change in catchments with different characteristics? To our knowledge, the results present one of the first quantitative investigations on the relative importance of parameter and data uncertainty across time. We find that the approach is capable of separating influential periods across data and parameter uncertainties, while also highlighting significant differences between the catchments analysed. Copyright © 2016 The Authors. Hydrological Processes. Published by John Wiley & Sons Ltd.     

From a Weak to a Strong Comet – 3d Global Hybrid Simulation Studies

Plasma structures resulting from the solar wind interaction with weakcomets are discussed. Numericalsimulations using a newly developed hybrid code are presented. The simulations are primarily applied to quantitative data for cometWirtanen, which will be the target of the Rosetta mission. It isexpected that Wirtanen is very weak during the first encounter. The main purpose is the discussion of the different features of the plasmaenvironment, such as the structured cycloidal plasma tail and non-linear Mach cones typical for weak comets and their relation tostructures like shocklets, bow shock, diamagnetic cavity and the “classical”magnetotail found at stronger comets. Furthermore, the sensitivity ofthese various features in dependence on the plasma parameters is investigated.     

Effects of a fully submerged boulder within a boulder array on the mean and turbulent flow fields: Implications to bedload transport

The objective of this coupled experimental and numerical study is to provide insight into the mean and turbulent flow fields within an array of fully submerged, isolated, immobile boulders. Our study showed that the velocity defect law performed well for describing the mean flow around the boulder within the array. A prerequisite, however, was to accurately estimate the spatial variability of u* around the boulder, which was achieved via the boundary characteristics method. The u* exhibited considerable spatial variability within the array and form roughness was shown to be up to 2 times larger than the skin roughness in the boulder near-wake region. Because the boulders bear a significant amount of the flow shear, the available bed shear stress for entrainment of the mobile sediment, ?? _ols , near the boulders was roughly 50% lower than the ambient ?? _ols . The ?? _ols variability induced by the boulders could lead to a threefold overestimation of the sediment transport rate.     

Test of statistical means for the extrapolation of soil depth point information using overlays of spatial environmental data and bootstrapping techniques

Hydrological modelling depends highly on the accuracy and uncertainty of model input parameters such as soil properties. Since most of these data are field surveyed, geostatistical techniques such as kriging, classification and regression trees or more sophisticated soil‐landscape models need to be applied to interpolate point information to the area. Most of the existing interpolation techniques require a random or regular distribution of points within the study area but are not adequate to satisfactorily interpolate soil catena or transect data. The soil landscape model presented in this study is predicting soil information from transect or catena point data using a statistical mean (arithmetic, geometric and harmonic mean) to calculate the soil information based on class means of merged spatial explanatory variables. A data set of 226 soil depth measurements covering a range of 0–6·5 m was used to test the model. The point data were sampled along four transects in the Stubbetorp catchment, SE‐Sweden. We overlaid a geomorphology map (8 classes) with digital elevation model‐derived topographic index maps (2–9 classes) to estimate the range of error the model produces with changing sample size and input maps. The accuracy of the soil depth predictions was estimated with the root mean square error (RMSE) based on a testing and training data set. RMSE ranged generally between 0·73 and 0·83 m ± 0·013 m depending on the amount of classes the merged layers had, but were smallest for a map combination with a low number of classes predicted with the harmonic mean (RMSE = 0·46 m). The results show that the prediction accuracy of this method depends on the number of point values in the sample, the value range of the measured attribute and the initial correlations between point values and explanatory variables, but suggests that the model approach is in general scale invariant. Copyright © 2009 John Wiley & Sons, Ltd.     

Garnet‐controlled very low velocities in the lower mantle transition zone at sites of mantle upwelling

Deep mantle plumes and associated increased geotherms are expected to cause an upward deflection of the lower–upper mantle boundary and an overall thinning of the mantle transition zone between about 410 and 660 km depth. We use subsequent forward modelling of mineral assemblages, seismic velocities, and receiver functions to explain the common paucity of such observations in receiver function data. In the lower mantle transition zone, large horizontal differences in seismic velocities may result from temperature‐dependent assemblage variations. At this depth, primitive mantle compositions are dominated by majoritic garnet at high temperatures. Associated seismic velocities are expected to be much lower than for ringwoodite‐rich assemblages at undisturbed thermal conditions. Neglecting this ultralow‐velocity zone at upwelling sites can cause a miscalculation of the lower–upper mantle boundary on the order of 20 km.