ALBERT

Description: Abstract Low‐frequency passive integrated transponders (PIT tags), are commonly used for monitoring pebble mobility in gravel‐bed rivers. Although early studies reported high recovery rates for PIT tags used in small streams, recovery rates in larger systems remain low, substantially limiting the possibilities for their use in such rivers. These low recovery rates are potentially due to missed detections caused by tag signal collision, burial in the sediment layer deeper than the maximum detection range and insufficient (but still exhausting) field effort to cover the concerned areas. A potential solution for addressing these problems is to use active ultra‐high frequency (a‐UHF) transponders as these have a greater detection range and anti‐collision protocols. In order to assess the potential of such transponders for pebble tracking in rivers, we used 433.92 MHz COIN‐ID and COIN‐HC models (ELA Innovation Company, Montpellier, France). We completed several tests to (i) characterize transponder detection ranges in the water and in saturated sediment and (ii) develop field protocols for locating tags by combining global positioning systems (GPS) sites and transponder received signal strength indication (RSSI) levels. The results showed that (i) the maximum detection ranges are about 2.4 m in the water column and more than 2.6 m in a column of saturated gravelly‐sandy sediment, (ii) RSSI spatial interpolation can be used to determine transponder position with good accuracy (〈 1 m), (iii) the desired minimal level of accuracy can be adjusted depending on in‐field effort and signal impulse interval, (iv) the RSSI maximal value observed cannot yet be used to determine transponder burial depth because of the multipath propagation of radio frequencies and the semi‐directional emission of the tag signal. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient kv was found to be vary around kv ~ 105 for hard materials (tensile strength ft 〉 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract When fine sediments are present in gravel streambeds (gravel‐framework beds), the gravel can be more easily removed from its original position, compared with gravel in a streambed without fine sediment but otherwise under the same hydraulic conditions. In this study, the effect of the presence of sand on the initiation of gravel motion in gravel riverbeds was investigated using flume experiments. The relationship between the critical Shields stress for gravel motion initiation and the fraction of sand in the bed was determined experimentally. The results can be summarized as follows. (1) When the fraction of sand in the bed is smaller than about 0.4, the critical Shields stress for the initiation of gravel motion decreases with increasing fraction of sand. The critical Shields stress increases, however, with increasing fraction of sand when it is larger than about 0.4. (2) The difference between the value of the critical Shields stress predicted by the Egiazaroff equation and the value obtained from the experimental data becomes maximum at about 0.4 of the fraction of sand. Here an empirical relation between the critical Shields stress and the fraction of sand is proposed so as to consider the effects of the ratio of the characteristic gravel size to the mean size of the bed material on the critical Shields stress. (3) Gravel in armored beds can be more easily mobilized by supplying sand as part of a sediment augmentation scheme. The sand fraction in the subsurface layer of the bed appears to reduce the friction angle of exposed particles. Sediment augmentation using sand has been recently demonstrated to be a viable alternative for mobilizing gravel for the restoration of gravel‐bed rivers downstream of dams. The quantitative evaluation obtained through the experiments reported here may be useful for the design of augmentation schemes. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract The assumption of the equilibrium state of gravel surfaces in flume experiments under feeding or recirculating conditions is generally justified by three equilibrium criteria based on sediment transport, slopes, and bed features. When these parameters become stable, an experiment is expected to reach equilibrium. This equilibrium state, however, is based on a one‐dimensional model, the Exner equation, which may not truly reflect the equilibrium state of the system considering the complex flow and sediment processes. In this paper, the evolutionary process of a gravel surface is investigated based on a large‐scale recirculating flume experiment. The performances of the three equilibrium criteria are evaluated first, and then the evolution of the bed morphology is studied. The key findings include the following: (1) the sediment transport rate, slopes of water and bed surfaces, and one‐dimensional morphological features reach equilibrium roughly simultaneously; (2) two‐dimensional morphology continually evolves after these characteristics reach equilibrium, which is confirmed by the characteristics of the sediment transport process; and (3) the results from a numerical simulation suggest that a much longer time is required to reach an equilibrium state. Our results suggest that sufficient experimental time is required to investigate the equilibrium morphological characteristics of gravel surfaces, which is much longer than the equilibrium time reflected by the one‐dimensional equilibrium criteria. Copyright © 2017 John Wiley & Sons, Ltd.

Description: ABSTRACT The Late Pleistocene Middle Stone Age (MSA) record of Africa provides early examples of standardised stone tool production and complex manufacturing sequences, superficially implying a long-term trend towards greater complexity in MSA technology at a continental scale. However, at this scale, spatial and temporal expressions of technological complexity are uneven. New lithic and chronometric data from the East–Central African record add further regional perspective to these patterns. Stone artefact assemblages from Karonga, northern, Malawi (92–22 ka), persistently lack the complexity demonstrated elsewhere in Africa at the same times, despite similar lithic raw materials. These new data provide an essential avenue for exploring hypotheses about the roles of environmental risk and demography in shaping the expression of MSA technology across the continent, not just at a local scale. When set within this framework, the simplicity of the Karonga MSA is best explained by its position in an environment that was persistently low in relative extrinsic subsistence risk. These results reinforce that motivations to invest in complex tools were variable through space and time, and that this variation, more than factors relating to behavioural capacity, may explain the patchy evidence for lithic complexity in the later MSA.

Description: Abstract In meandering rivers cut into bedrock, erosion across a channel cross‐section can be strongly asymmetric. At a meander apex, deep undercutting of the outer bank can result in the formation of a hanging cliff (which may drive hillslope failure), whereas the inner bank adjoins a slip‐off slope that connects to the hillslope itself. Here we propose a physically‐based model for predicting channel planform migration and incision, point bar and slip‐off slope formation, bedrock abrasion, the spatial distribution of alluvial cover, and adaptation of channel width in a mixed bedrock‐alluvial channel. We simplify the analysis by considering a numerical model of steady, uniform bend flow satisfying cyclic boundary conditions. Thus in our analysis, ‘sediment supply’, i.e. the total volume of alluvium in the system, is conserved. In our numerical simulations, the migration rate of the outer bank is a specified parameter. Our simulations demonstrate the existence of an approximate state of dynamic equilibrium corresponding to a near‐solution of permanent form in which a bend of constant curvature, width, cross‐sectional shape and alluvial cover distribution migrates diagonally downward at constant speed, leaving a bedrock equivalent of a point bar on the inside of the bend. Channel width is set internally by the processes of migration and incision. We find that equilibrium width increases with increasing sediment supply, but is insensitive to outer bank migration rate. The slope of the bedrock point bar varies inversely with both outer bank migration rate and sediment supply. Although the migration rate of the outer bank is externally imposed here, we discuss a model modification that would allow lateral side‐wall abrasion to be treated in a manner similar to the process of bedrock incision. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Abstract Particle dynamics are investigated experimentally in supercritical high‐speed open channel flow over a fixed planar bed of low relative roughness height simulating flows in high‐gradient non‐alluvial mountain streams and hydraulic structures. Non‐dimensional equations were developed for transport mode, particle velocity, hop length and hop height accounting for a wide range of literature data encompassing sub‐ and supercritical flow conditions as well as planar and alluvial bed configurations. Particles were dominantly transported in saltation and particle trajectories on planar beds were rather flat and long compared with alluvial bed data due to (1) increased lift forces by spinning motion, (2) strongly downward directed secondary currents, and (3) a planar flume bed where variation in particle reflection and damping effects were minor. The analysis of particle saltation trajectories revealed that the rising and falling limbs were almost symmetrical contradicting alluvial bed data. Furthermore, no or negligible effect of particle size and shape on particle dynamics were found. Implications of experimental findings for mechanistic saltation‐abrasion models are briefly discussed. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract Large wood along rivers influences entrainment, transport, and storage of mineral sediment and particulate organic matter. We review how wood alters sediment dynamics and explore patterns among volumes of in‐stream wood, sediment storage, and residual pools for dispersed pieces of wood, logjams, and beaver dams. We hypothesized that: volume of sediment per unit area of channel stored in association with wood is inversely proportional to drainage area; the form of sediment storage changes downstream; sediment storage correlates with wood load; the residual volume of pools created in association with wood correlates inversely with drainage area; and volume of sediment stored behind beaver dams correlates with pond area. Lack of data from larger drainage areas limits tests of these hypotheses, but the analyses suggest that sediment volume correlates positively with drainage area and wood volume. The form of sediment storage in relation to wood appears to change downstream, with wedges of sediment upstream from jammed steps most prevalent in small, steep channels and more dispersed sediment storage in lower gradient channels. Pool volume correlates positively with wood volume and negatively with channel gradient. Sediment volume correlates well with beaver pond area. More abundant in‐stream wood and beaver populations present historically equated to greater sediment storage within river corridors and greater residual pool volume. One implication of these changes is that protecting and re‐introducing wood and beavers can be used to restore rivers. This review of the existing literature on wood and sediment dynamics highlights the lack of studies on larger rivers. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Abstract For more than 150 years, humans have tried to limit the geomorphic activity of mountain streams, and the related damage, using torrent control works. Check dams are likely the most emblematic civil engineering structures used in soil conservation programs. Modern mountain societies have inherited thousands of these structures built in upland gullies and streams. To help define their effectiveness and decisions concerning their maintenance or new project designs, a clear understanding of potential effects of check dams on river systems, i.e. their functions, is first needed. The next steps concern quantitative assessments of each function on the flood features and combination of all effects. The present understanding of these sometimes old structures' functions can be complicated because the societal and environmental contexts in which the original structures were built may have changed. To bridge this gap, this paper traces the purposes for which check dams were built, through a detailed analysis of French archives. We first analyze chronologically how each function was theorized and applied in the field. In the nineteenth century, engineers developed a thorough empirical and conceptual knowledge of mountain soil erosion, torrential geomorphology, and sediment transport processes as well as check dam interactions with these natural processes. The second part of this paper synthesizes conceptual descriptions of the check dams' functions, in the light of more than 150 years of experience, with their implication on the features of the structures. The French experience is compared to other countries' pioneering works. Finally, the next steps and remaining research challenges toward a comprehensive analysis of check dams' efficiency in torrent hazard mitigation are presented. This analysis is proposed to remind how, conceptually, check dams may influence geomorphic systems, bearing in mind the knowledge represented in pioneer guidelines and recent works on the subject. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Abstract Structure‐from‐Motion (SfM) photogrammetry is now used widely to study a range of earth surface processes and landforms, and is fast becoming a core tool in fluvial geomorphology. SfM photogrammetry allows extraction of topographic information and orthophotos from aerial imagery. However, one field where it is not yet widely used is that of river restoration. The characterisation of physical habitat conditions pre‐ and post‐restoration is critical for assessing project success, and SfM can be used easily and effectively for this purpose. In this paper we outline a workflow model for the application of SfM photogrammetry to collect topographic data, develop surface models and assess geomorphic change resulting from river restoration actions. We illustrate the application of the model to a river restoration project in the NW of England, to show how SfM techniques have been used to assess whether the project is achieving its geomorphic objectives. We outline the details of each stage of the workflow, which extend from preliminary decision‐making related to the establishment of a ground control network, through fish‐eye lens camera testing and calibration, to final image analysis for the creation of facies maps, the extraction of point clouds, and the development of digital elevation models (DEMs) and channel roughness maps. The workflow enabled us to confidently identify geomorphic changes occurring in the river channel over time, as well as assess spatial variation in erosion and aggradation. Critical to the assessment of change was the high number of ground control points and the application of a minimum level of detection threshold used to assess uncertainties in the topographic models. We suggest that these two things are especially important for river restoration applications. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Abstract Experiments were undertaken to study the nature of granular interaction in running water by examining the influence of fine grain inputs to a coarser sediment bed with a mobile surface. Video recordings of grain sorting by both kinetic sieving and spontaneous percolation are used to diagnose the critical processes controlling the overall bed response. Kinetic sieving takes place in the mobile bed surface, with the finer sediment moving to the bottom of the bedload transport layer at the interface with the underlying quasi‐static coarse bed. We show that the behavior at this interface dictates how a channel responds to a fine sediment input. If, by spontaneous percolation, the fine sediment is able to infiltrate into the underlying quasi‐static bed, the total transport increases and the channel degrades. However, if the fine sediment input rate exceeds the transport capacity or is geometrically unable to infiltrate into the underlying bed, it forms a quasi‐static layer underneath the transport layer that inhibits entrainment from the underlying bed, resulting in aggradation and an increase in bed slope. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Topographic monitoring of a gravel augmentation during a dam‐controlled flood in the Trinity River, CA, indicates that all of the gravel entering the study reach was captured in one of two bed‐material storage reservoirs. Such interactions with storage reservoirs can cause large bed‐material pulses to disperse by fragmenting into multiple smaller pulses. As a refinement to the conceptual model of sediment pulse evolution by dispersion and translation, the concept of pulse fragmentation has practical implications for gravel management. Abstract The geomorphic effect of introducing a gravel augmentation totaling 520 m3 into a gravel‐bed stream during a dam‐controlled flood in May of 2015 was monitored with bedload transport measurements, an array of seismometers, and repeated topographic surveys. Half of the augmented gravel was injected into the flow with front‐end loaders on the rising limb of the flood and the other half was injected on the first day of the peak. Virtually all of the gravel transported past the injection point was deposited within about 7 to 10 channel widths of the injection point. Most of the injected gravel deposited along the left bank of the river whereas the right half of the channel bed was dominated by scour. The downstream third of the depositional area consisted of a small dune field that developed prior to the second gravel injection and subsequently migrated about one channel width downstream. A second depositional front was observed upstream from the gravel injection point, where a delta‐like wedge of bed material developed in the first hours of the flow release and changed little over the remainder of the release. These two depositional areas represent small‐scale bed‐material storage reservoirs with the potential to accumulate and periodically release packets of bed material. Interactions with such storage reservoirs are hypothesized to cause large bed‐material pulses to disperse by fragmenting into multiple smaller pulses. As a refinement to the conceptual model that views sediment pulse evolution in terms of dispersion and translation, the concept of pulse fragmentation has practical implications for gravel management. It implies that gravel augmentations can produce morphologic changes at locations that are separated from the augmentation point by arbitrarily long reaches, and it highlights the dependence of pulse propagation rates on the nature and distribution of the bed‐material storage reservoirs in the channel system. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract The artificial gravel augmentation of river channels is increasingly being used to mitigate the adverse effects of river regulation and sediment starvation. A systematic framework for designing and assessing such gravel augmentations is still lacking, notably on large rivers. Monitoring is required to quantify the movement of augmented gravel, measure bedform changes, assess potential habitat enhancement, and reduce the uncertainty in sediment management. Here we present the results of an experiment conducted in the Rhine River (French and German border). In 2010, 23 000 m3 of sediments (approximately the mean annual bedload transport capacity) were supplied in a by‐passed reach downstream of the Kembs dam to test the feasibility of enhancing sediment transport and bedform changes. A 620‐m‐long and 12‐m‐wide gravel deposit was created 8 km downstream from the dam. Monitoring included topo‐bathymetric surveys, radio‐frequency particle tracking using passive integrated transponder (PIT) tags, bed grain size measurement, and airborne imagery. Six surveys performed since 2009 have been described (before and after gravel augmentation, and after Q2 and Q15 floods). The key findings are that (i) the augmented gravel was partially dispersed by the first flood event of December 2010 (Q1); (ii) PIT tags were found up to 3200 m downstream of the gravel augmentation site after four years, but the effects of gravel augmentation could not be clearly distinguished from the effects of floods and internal remobilization on more than 3500 m downstream; (iii) linear and log‐linear relationships linking bedload transport, particle mobility, and grain size were established; and (iv) combined bathymetry and PIT tag surveys were useful for evaluating potential environmental risks and the first morpho‐ecological responses. This confirmed the complementary nature of such techniques in the monitoring of gravel augmentation in large rivers. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Mountain gravel‐bed rivers are subjected to punctuated floods and sediment delivery. We model this numerically in terms of a cycled hydrograph upstream and a cycled “sedimentograph” midreach. Both the hydrograph and the sedimentograph generate boundary layers, within which the effect of the fluctuations are felt. Outside the boundary layers, steady, uniform conditions representing the overall effects of the hydrograph and sedimentograph are achieved. Our result provides a useful tool for simplifying the modeling of long‐term morphodynamics in mountain rivers. Abstract Sediment often enters rivers in the form of sediment pulses associated with landslides and debris flows. This is particularly so in gravel‐bed rivers in earthquake‐prone mountain regions, such as Southwest China. Under such circumstances, sediment pulses can rapidly change river topography and leave the river in repeated states of gradual recovery. In this paper, we implement a one‐dimensional morphodynamic model of river response to pulsed sediment supply. The model is validated using data from flume experiments, so demonstrating that it can successfully reproduce the overall morphodynamics of experimental pulses. The model is then used to explore the evolution of a gravel‐bed river subject to cycled hydrographs and repeated sediment pulses. These pulses are fed into the channel in a fixed region centered at a point halfway down the calculational domain. The pulsed sediment supply is in addition to a constant sediment supply at the upstream end. Results indicate that the river can reach a mobile‐bed equilibrium in which two regions exist within which bed elevation and surface grain size distribution vary periodically in time. One of these is at the upstream end, where a periodic discharge hydrograph and constant sediment supply are imposed, and the other is in a region about halfway down the channel where periodic sediment pulses are introduced. Outside these two regions, bed elevation and surface grain size distribution reach a mobile‐bed equilibrium that is invariant in time. The zone of fluctuation‐free mobile‐bed equilibrium upstream of the pulse region is not affected by repeated sediment pulses under the scenarios tested, but downstream of the pulse region, the channel reaches different fluctuation‐free mobile‐bed equilibriums under different sediment pulse scenarios. The vertical bed structure predicted by the simulations indicates that the cyclic variation associated with the hydrograph and sediment pulses can affect the substrate stratigraphy to some depth. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract Large wood (LW) is a ubiquitous feature in rivers of forested watersheds worldwide, and its importance for river diversity has been recognized for several decades. Although the role of LW in fluvial dynamics has been extensively documented, there is a need to better quantify the most significant components of LW budgets at the river scale. The purpose of our study was to quantify each component (input, accumulation, and output) of a LW budget at the reach and watershed scales for different time periods (i.e. a 50‐year period, decadal cycle, and interannual cycle). The LW budget was quantified by measuring the volumes of LW inputs, accumulations, and outputs within river sections that were finally evacuated from the watershed. The study site included three unusually large but natural wood rafts in the delta of the Saint‐Jean River (SJR; Québec, Canada) that have accumulated all LW exported from the watershed for the last 50 years. We observed an increase in fluvial dynamics since 2004, which led to larger LW recruitment and a greater LW volume trapped in the river corridor, suggesting that the system is not in equilibrium in terms of the wood budget but is rather recovering from previous human pressures as well as adjusting to hydroclimatic changes. The results reveal the large variability in the LW budget dynamics during the 50‐year period and allow us to examine the eco‐hydromorphological trajectory that highlights key variables (discharge, erosion rates, bar surface area, sinuosity, wood mobility, and wood retention). Knowledge on the dynamics of these variables improves our understanding of the historical and future trajectories of LW dynamics and fluvial dynamics in gravel‐bed rivers. Extreme events (flood and ice‐melt) significantly contribute to LW dynamics in the SJR river system. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract The question: ‘how does a streambed change over a minor flood?’ does not have a clear answer due to lack of measurement methods during high flows. We investigate bedload transport and disentrainment during a 1.5‐year flood by linking field measurements using fiber optic distributed temperature sensing (DTS) cable with sediment transport theory and an existing explicit analytical solution to predict depth of sediment deposition from amplitude and phase changes of the diurnal near‐bed pore‐water temperature. The method facilitates the study of gravel transport by using near‐bed temperature time series to estimate rates of sediment deposition continuously over the duration of a high flow event coinciding with bar formation. The observations indicate that all gravel and cobble particles present were transported along the riffle at a relatively low Shields Number for the median particle size, and were re‐deposited on the lee side of the bar at rates that varied over time during a constant flow. Approximately 1–6% of the bed was predicted to be mobile during the 1.5‐year flood, indicating that large inactive regions of the bed, particularly between riffles, persist between years despite field observations of narrow zones of local transport and bar growth on the order ~3–5 times the median particle size. In contrast, during a seven‐year flood approximately 8–55% of the bed was predicted to become mobile, indicating that the continuous along‐stream mobility required to mobilize coarse gravel through long pools and downstream to the next riffle is infrequent. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract Piping has been recognized as an important geomorphic, soil erosion and hydrologic process. It seems that it is far more widespread than it has often been supposed. However, our knowledge about piping dynamics and its quantification currently relies on a limited number of data for mainly loess‐derived areas and marl badlands. Therefore, this research aimed to recognize piping dynamics in mid‐altitude mountains under a temperate climate, where piping occurs in Cambisols, not previously considered as piping‐prone soils. It has been expressed by the estimation of erosion rates due to piping and elongation of pipes in the Bereźnica Wyżna catchment in the Bieszczady Mountains, eastern Carpathians (305 ha, 188 collapsed pipes). The research was based on the monitoring of selected piping systems (1971–1974, 2013–2016). Changes in soil loss vary significantly between different years (up to 27.36 t ha−1 yr−1), as well as between the mean short‐term erosion rate (up to 13.10 t ha−1 yr−1), and the long‐term (45 years) mean of 1.34 t ha−1 yr−1. The elongation of pipes also differs, from no changes to 36 m during one year. The mean total soil loss is 48.8 t ha−1 in plots, whereas in the whole studied catchment it is 2.0 t ha−1. Hence, piping is both spatially and temporally dependent. The magnitude of piping in the study area is at least three orders of magnitude higher than surface erosion rates (i.e. sheet and rill erosion) under similar land use (grasslands), and it is comparable to the magnitude of surface soil erosion on arable lands. It means that piping constitutes a significant environmental problem and, wherever it occurs, it is an important, or even the main, sediment source. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Abstract This paper reports on a first attempt of using the virtual velocity approach to assess sediment mobility and transport in two wide and complex gravel‐bed rivers of northern Italy. Displacement length and virtual velocity of spray‐painted tracers were measured in the field. Also, the thickness of the sediment active layer during floods was measured using scour chains and post‐flood morphological changes as documented by repeated survey of channel cross‐sections. The effects of eight and seven floods were studied on the Tagliamento and Brenta Rivers, where 259 and 277 spray‐painted areas were surveyed, respectively. In the Tagliamento River 36% of the spray‐painted areas experienced partial transport, whereas in the Brenta River this accounted for 20%. Whereas, full removal/gravel deposition was observed on 37% and 26% of these areas on the Tagliamento and Brenta Rivers, respectively. The mean displacement length of particles, the thickness of the active layer and the extent of partial transport are well correlated with the dimensionless shear stress. The virtual velocity approach allowed calculation of bed material transport over a wide range of flood magnitudes. Annual coarse sediment transport was calculated up to 150 for the Tagliamento, and 30 × 103  m3 yr−1 for the Brenta. The outcomes of this work highlight the relevance of partial transport condition, as it could represent more than 70% of the total bed material transported during low‐magnitude floods, and up to 40% for near‐bankfull events. Results confirm that bed material load tends to be overestimated by traditional formulas. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Abstract During floods, large quantities of wood can be mobilized and transported downstream. At critical sections, such as bridges, the transported wood might be entrapped and a quick succession of backwater effects can occur as a result of the reduction of the cross‐sectional area. The aim of this work is to explore large wood‐related hazards during floods in the gravel‐bed river Czarny Dunajec (Polish Carpathians), where the river flows through the village of Długopole. This work is based on the numerical modelling of large wood transport together with flow dynamics in which inlet and boundary conditions were designed based on field observations. The exploratory approach developed in this study uses multiple scenarios (193) to analyse the factors controlling bridge clogging: wood size, wood supply, flow conditions, morphology and obstacles in the riverbed. Results highlighted the strong control of log length (stronger than that of log diameter) on potential blockage probability; however, according to our results the main factor controlling bridge clogging was the flood discharge. River morphology and wood supply play an important role as well. The river morphology may reduce bridge blockage, as it influences flow velocity and depth, and creates natural retention zones for wood. In addition, the impacts of bridge blockage were analysed in terms of afflux depth and length, and flooded area. Results showed that bridge blockage may result in a significant increase in water depth (up to 0.7 m) and flooded area (up to 33% more), therefore increasing flood risk in the village. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Abstract Regulated rivers generally incise below dams that cut off sediment supply, but how that happens and what the consequences are at different spatial scales is poorly understood. Modern topographic mapping at meter‐scale resolution now enables investigation of the details of spatial processes. In this study, spatial segregation was applied to a meter‐scale raster map of topographic change from 1999 to 2008 on the gravel‐cobble, regulated lower Yuba River in California to answer specific scientific questions about how a decadal hydrograph that included a flood peak of 22 times bankfull discharge affected the river at segment, reach, and morphological unit scales. The results show that the river preferentially eroded sediment from floodplains compared to the channel, and this not only promoted valley‐wide sediment evacuation, but also facilitated the renewal and differentiation of morphological units, especially in the channel. At the reach scale, area of fill and mean net rate of elevational change were directly correlated with better connectivity between the channel and floodplain, while the mean rate of scour in scour areas was influenced by the ratio of slope to bankfull Froude number, a ratio indicative of lateral migration versus vertical downcutting. Hierarchical segregation of topographic change rasters proved useful for understanding multi‐scalar geomorphic dynamics. Copyright © 2016 John Wiley & Sons, Ltd.

Description: Comprehensive research on glacier changes in the Tian Shan is available for the current decade; however, there is limited information about glacier investigations of previous decades and especially before the mid 1970s. The earliest stereo images from the Corona missions were acquired in the 1960s but existing studies dealing with these images focus on single glaciers or small areas only. We developed a workflow to generate digital terrain models (DTMs) and orthophotos from 1964 Corona KH-4 for an entire mountain range (Ak-Shirak) located in the Central Tian Shan. From these DTMs and orthoimages, we calculated geodetic mass balances and length changes in comparison to 1973 and 1980 Hexagon KH-9 data. We found mass budgets between −0.4 ± 0.1 m·w.e.a−1 (1964–1980) and −0.9 ± 0.4 m·w.e.a−1 (1973–1980) for the whole region and individual glaciers. The length changes, on the other hand, vary heterogeneously between +624 ± 18 m (+39.0 ± 1.1 m·a−1) and −923 ± 18 m (−57.7 ± 1.1 m·a−1) for 1964–1980. An automation of the processing line can successively lead to region-wide Corona data processing allowing the analysis and interpretation of glacier changes on a larger scale and supporting a refinement of glacier modelling.

Description: The potato cyst nematodes Globodera pallida and G. rostochiensis are economically important plant pathogens causing losses to UK potato harvests estimated at £50m/ year. Implications of climate change on their future pest status have not been fully considered. Here, we report growth of female G. pallida and G. rostochiensis over the range 15 to 25 °C. Females per plant and their fecundity declined progressively with temperatures above 17.5 °C for G. pallida , whilst females per plant were optimal between 17.5 and 22.5 °C for G. rostochiensis . Relative reproductive success with temperature was confirmed on two potato cultivars infected with either species at 15, 22.5 and 25 °C. The reduced reproductive success of G. pallida at 22.5 °C relative to 15 °C was also recorded for a further seven host cultivars studied. The differences in optimal temperatures for reproductive success may relate to known differences in the altitude of their regions of origin in the Andes. Exposure of G. pallida to a diurnal temperature stress for one week during female growth significantly suppressed subsequent growth for one week at 17.5 °C but had no effect on G. rostochiensis . However, after two weeks of recovery female size was not significantly different from that for the control treatment. Future soil temperatures were simulated for medium and high emissions scenarios and combined with nematode growth data to project future implications of climate change for the two species. Increased soil temperatures associated with climate change may reduce the pest status of G. pallida but benefit G. rostochiensis especially in southern UK. We conclude that plant breeders may be able to exploit the thermal limits of G. pallida by developing potato cultivars able to grow under future warm summer conditions. Existing widely deployed resistance to G. rostochiensis is an important characteristic to retain for new potato cultivars. This article is protected by copyright. All rights reserved.

Description: ABSTRACT Hydrology has advanced considerably as a scientific discipline since its recognized inception in the mid-20th century. Modern water resource related questions have forced adaptation from exclusively physical or engineering science viewpoints toward a deliberate interdisciplinary context. Over the past few decades, many of the eventual manifestations of this evolution were foreseen by prominent expert hydrologists. However, their narrative descriptions have lacked substantial quantification. This study addressed that gap by measuring the prevalence of and analyzing the relationships between the terms most frequently used by hydrologists to define and describe their research. We analyzed 16,591 journal article titles from 1965-2015 in Water Resources Research , through which the scientific dialogue and its time-sensitive progression emerged. Our word frequency and term co-occurrence network results revealed the dynamic timing of the lateral movement of hydrology across multiple disciplines as well as the deepening of scientific discourse with respect to traditional hydrologic questions. The conversation among water resource scientists surrounding the hydrologic sub-disciplines of catchment-hydrology, hydro-meteorology, socio-hydrology, hydro-climatology and eco-hydrology all gained statistically significant momentum in the analyzed time period, while hydro-geology and contaminant-hydrology experienced periods of increase followed by significant decline. This study concludes that formerly exotic disciplines can potentially modify hydrology, prompting new insights and inspiring unconventional perspectives on old questions that may have otherwise become obsolete. This article is protected by copyright. All rights reserved.

Description: Threshold estimation in the Peaks Over Threshold (POT) method, and the impact of the estimation method on the calculation of high return period quantiles and their uncertainty (or confidence intervals) are issues that are still unresolved. In the past, methods based on goodness-of-fit tests and EDF-statistics have yielded satisfactory results, but their use has not yet been systematized. This paper proposes a methodology for automatic threshold estimation, based on the Anderson-Darling EDF-statistic and goodness-of-fit test. When combined with bootstrapping techniques, this methodology can be used to quantify both the uncertainty of threshold estimation and its impact on the uncertainty of high return period quantiles. This methodology was applied to several simulated series and to four precipitation/riverflow data series. The results obtained confirmed its robustness. For the measured series, the estimated thresholds corresponded to those obtained by non-automatic methods. Moreover, even though the uncertainty of the threshold estimation was high, this did not have a significant effect on the width of the confidence intervals of high return period quantiles. This article is protected by copyright. All rights reserved.

Description: Snow availability in Alpine catchments plays an important role in water resources management. In this paper we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical-based methodologies. Different techniques (multiple regression, multi-criteria analysis and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful based on a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) based on multi-criteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent (SWE) resources.

Description: Geographic Object-Based Image Analysis (GEOBIA) mostly uses proprietary software,but the interest in Free and Open-Source Software (FOSS) for GEOBIA is growing. This interest stems not only from cost savings, but also from benefits concerning reproducibility and collaboration. Technical challenges hamper practical reproducibility, especially when multiple software packages are required to conduct an analysis. In this study, we use containerization to package a GEOBIA workflow in a well-defined FOSS environment. We explore the approach using two software stacks to perform an exemplary analysis detecting destruction of buildings in bi-temporal images of a conflict area. The analysis combines feature extraction techniques with segmentation and object-based analysis to detect changes using automatically-defined local reference values and to distinguish disappeared buildings from non-target structures. The resulting workflow is published as FOSS comprising both the model and data in a ready to use Docker image and a user interface for interaction with the containerized workflow. The presented solution advances GEOBIA in the following aspects: higher transparency of methodology; easier reuse and adaption of workflows; better transferability between operating systems; complete description of the software environment; and easy application of workflows by image analysis experts and non-experts. As a result, it promotes not only the reproducibility of GEOBIA, but also its practical adoption.

Description: Digital surface models (DSMs) derived from spaceborne and airborne sensors enable the monitoring of the vertical structures for forests in large areas. Nevertheless, due to the lack of an objective performance assessment for this task, it is difficult to select the most appropriate data source for DSM generation. In order to fill this gap, this paper performs change detection analysis including forest decrease and tree growth. The accuracy of the DSMs is evaluated by comparison with measured tree heights from inventory plots (field data). In addition, the DSMs are compared with LiDAR data to perform a pixel-wise quality assessment. DSMs from four different satellite stereo sensors (ALOS/PRISM, Cartosat-1, RapidEye and WorldView-2), one satellite InSAR sensor (TanDEM-X), two aerial stereo camera systems (HRSC and UltraCam) and two airborne laser scanning datasets with different point densities are adopted for the comparison. The case study is a complex central European temperate forest close to Traunstein in Bavaria, Germany. As a major experimental result, the quality of the DSM is found to be robust to variations in image resolution, especially when the forest density is high. The forest decrease results confirm that besides aerial photogrammetry data, very high resolution satellite data, such as WorldView-2, can deliver results with comparable quality as the ones derived from LiDAR, followed by TanDEM-X and Cartosat DSMs. The quality of the DSMs derived from ALOS and Rapid-Eye data is lower, but the main changes are still correctly highlighted. Moreover, the vertical tree growth and their relationship with tree height are analyzed. The major tree height in the study site is between 15 and 30 m and the periodic annual increments (PAIs) are in the range of 0.30–0.50 m.

Description: This paper presents an automated and effective framework for classifying airborne laser scanning (ALS) point clouds. The framework is composed of four stages: (i) step-wise point cloud segmentation, (ii) feature extraction, (iii) Random Forests (RF) based feature selection and classification, and (iv) post-processing. First, a step-wise point cloud segmentation method is proposed to extract three kinds of segments, including planar, smooth and rough surfaces. Second, a segment, rather than an individual point, is taken as the basic processing unit to extract features. Third, RF is employed to select features and classify these segments. Finally, semantic rules are employed to optimize the classification result. Three datasets provided by Open Topography are utilized to test the proposed method. Experiments show that our method achieves a superior classification result with an overall classification accuracy larger than 91.17%, and kappa coefficient larger than 83.79%.

Description: The response of hillslope processes to changes in precipitation may drive the observed changes in the solute geochemistry of rivers with discharge. This conjecture is most robust when variations in the key environmental factors that affect hillslope processes (e.g., lithology, erosion rate, and climate) are minimal across a river's catchment area. For rivers with heterogenous catchments, temporal variations in the relative contributions of different tributary sub-catchments may modulate variations in solute geochemistry with runoff. In the absence of a dense network of hydrologic gauging stations, alternative approaches are required to distinguish between the different drivers of temporal variability in river solute concentrations. In this contribution, we apportion the water and solute fluxes of a reach of the Madre de Dios River (Peru) between its four major tributary sub-catchments during two sampling campaigns (wet and dry seasons) using spatial variations in conservative tracers. Guided by the results of a mixing model, we identify temporal variations in solute concentrations of the mainstem Madre de Dios that are due to changes in the relative contributions of each tributary. Our results suggest that variations in tributary mixing are, in part, responsible for the observed concentration-discharge (C-Q) relationships. The implications of these results are further explored by re-analyzing previously published C-Q data from this region, developing a theoretical model of tributary mixing, and, in a companion paper, comparing the C-Q behavior of a suite of major and trace elements in the Madre de Dios River system. This article is protected by copyright. All rights reserved.

Description: Solving inverse problems in a complex, geologically realistic, and discrete model space and from a sparse set of observations is a very challenging task. Extensive exploration by Markov chain Monte Carlo (McMC) methods often results in considerable computational efforts. Most optimization methods, on the other hand, are limited to linear (continuous) model spaces and the minimization of an objective function, what often proves to be insufficient. To overcome these problems, we propose a new ensemble based exploration scheme for geostatistical prior models generated by a multiple-point statistics (MPS) tool. The principle of our method is to expand an existing set of models by using posterior facies information for conditioning new MPS realizations. The algorithm is independent of the physical parametrization. It is tested on a simple synthetic inverse problem. When compared to two existing McMC methods (Iterative Spatial Resampling (ISR) and Interrupted Markov chain Monte Carlo (IMcMC)) the required number of forward model runs was divided by a factor of 8-12. This article is protected by copyright. All rights reserved.

Description: Although vegetation is present in many rivers, the bulk of past work concerned with modeling the influence of vegetation on flow has considered vegetation to be morphologically simple, and has generally neglected the complexity of natural plants. Here we report on a combined flume and numerical model experiment which incorporates time-averaged plant posture, collected through Terrestrial Laser Scanning, into a Computational Fluid Dynamics model to predict flow around a submerged riparian plant. For three depth-limited flow conditions (Reynolds number = 65 000 – 110 000), plant dynamics were recorded through high-definition video imagery, and the numerical model was validated against flow velocities collected with an acoustic Doppler velocimeter. The plant morphology shows an 18% reduction in plant height and a 14% increase in plant length, compressing and reducing the volumetric canopy morphology as the Reynolds number increases. Plant shear layer turbulence is dominated by Kelvin–Helmholtz type vortices generated through shear instability, the frequency of which is estimated to be between 0.20 and 0.30 Hz, increasing with Reynolds number. These results demonstrate the significant effect that the complex morphology of natural plants has on in-stream drag, and allows a physically determined, species-dependent drag coefficient to be calculated. Given the importance of vegetation in river corridor management, the approach developed here demonstrates the necessity to account for plant motion when calculating vegetative resistance. This article is protected by copyright. All rights reserved.

Description: The vast majority of continental sediment delivered to the world's oceans moves by suspension in rivers. Depth- or point-integrated bottle sampling are the traditional methods used to determine the mean concentration of suspended sediment in rivers. While there has been some investigation of the error associated with depth-integrated sampling, the representativeness of a point-integrated bottle sample has not been addressed in the literature. Here, we analyze continuous hour-long measurements of suspended sediment and grain size fractions collected using a LISST-SL in the sand-bed portion of the Fraser River, British Columbia to determine an appropriate sampling time. The 2σ uncertainty range of individual 30 s samples varied from ±3% to ±33% about the observed mean, with a systematic increase toward the streambed. Mean concentrations for suspended sediment and grain size fractions were computed over increasing time periods and compared with a long duration mean concentration to determine when a sample becomes representative. A cumulative probability distribution was generated from multiple iterations of this process. All suspended sediment load and grain size fractions bear a low probability of representing the actual mean concentration over standard bottle sample durations. A probability 〉90% of representing the mean concentration and grain-size of various fractions requires ∼570 seconds (9.5 minutes) of sampling. Sampling for a shorter period of 264 seconds (4.4 minutes) can yield a sample with 73% probability of representing the mean concentration. This article is protected by copyright. All rights reserved.

Description: Variations in riverine solute chemistry with changing runoff are used to interrogate catchment hydrology and to investigate chemical reactions in Earth's critical zone. This approach requires some understanding of how spatial and temporal averaging of solute-generating reactions affect the dissolved load of rivers and streams. In this study, we investigate the concentration-runoff (C-Q) dynamics of a suite of major (Na, Mg, Ca, Si, K, and SO 4 ) and trace (Al, Ba, Cd, Co, Cr, Cu, Fe, Ge, Li, Mn, Mo, Nd, Ni, Rb, Sr, U, V, and Zn) elements in nested catchments of variable size, spanning the geomorphic gradient from the Andes mountains to the Amazon foreland-floodplain. The major elements exhibit various degrees of dilution with increasing runoff at all sites, whereas the concentrations of most trace elements either increase or show no relationship with increasing runoff in the three larger catchments (160 to 28 000 km 2 area). We show that the observed mainstem C-Q dynamics are influenced by variable mixing of tributaries with distinct C-Q relationships. Trace element C-Q relationships are more variable among tributaries relative to major elements, which could be the result of variations in geomorphology, lithology, and hydrology of the sub-catchments. Certain trace metals are also lost from solution during in-channel processes (possibly related to colloidal size-partitioning), which may exert an additional control on C-Q dynamics. Overall, we suggest that aggregation effects should be assessed in heterogeneous catchments before C-Q or ratio-Q relationships can be interpreted as reflecting catchment-wide solute generation processes and their relationship to hydrology. This article is protected by copyright. All rights reserved.

Description: Nitrogen is an essential nutrient in many terrestrial ecosystems because it affects vegetation’s primary production. Due to the variety of nitrogen-containing substances and the differences in their composition across species, statistical approaches are now dominant in remote sensing retrieval of leaf nitrogen content. Many studies remove spectral regions characterized by strong water absorptions before retrieving nitrogen content, because water is believed to mask the absorption features of nitrogen. The objectives of this study are to discuss the necessity of this practice and to explore how water absorption affects leaf nitrogen estimation. Spectral measurements and chemical analyses for Maize, Sawtooth Oak, and Sweetgum leaves were carried out in 2014. The leaf optical properties model PROSPECT5 was used to eliminate the influences of water on the measured reflectance spectra. The inversion accuracy of PROPECT5 for chlorophyll, carotenoid, water, and dry matter of Maize was also discussed. Measured, simulated, and water-removed spectra were used to: (1) find the optimal nitrogen-related spectral index; and (2) regress with the area-based leaf nitrogen concentration (LNC) using the partial least square regression technique (PLSR). Two types of spectral indices were selected in this study: Normalized Difference Spectral Index (NDSI) and Ratio Spectral Index (RSI). Additionally, first-order derivative forms of measured, simulated, and water-removed spectra were devised to search for the optimal spectral indices. Finally, species-specific optimal indices and cross-species optimal indices, as well as their root mean square errors (RMSE) and coefficients of determination (R2), were obtained. The Ending Top Percentile (ETP), an indicator of the performance of cross-species optimal indices, was also calculated. PLSR was combined with leave-one-out cross validation (LOOCV) for each species. The predicted root mean square errors (RMSEP) and predicted R2 were finally calculated. The results showed that chlorophyll, carotenoid, and water contents could be estimated with R2 of 0.75, 0.59, and 0.69, respectively, which were acceptable for fresh leaves. The dry matter was retrieved with a relatively lower accuracy because of the fixed absorption coefficients adopted by PROSPECT5. The performances of species-specific optimal indices using water-free spectra were comparable to or worse than the corresponding indices derived with measured or simulated spectra. Compared with measured spectra, ETP did not change much after the effects of water were removed, and the R2 between cross-species optimal spectral indices and area-based LNC for Sawtooth Oak and Sweetgum decreased while it remained almost the same for Maize, suggesting that the water-removed cross-species optimal indices were inferior to the corresponding optimal indices found without water removal. ETP was larger than 30% for all spectra, demonstrating the non-existence of common optimal NDSI or RSI for the three species. After water removal, the accuracy of PLSR for Sawtooth Oak and Sweetgum decreased and increased negligibly for Maize. The results suggest that water absorption has limited effects on reducing the accuracy of leaf nitrogen estimation. On the contrary, the accuracy may decrease due to the loss of spectral information caused by the removal of water-sensitive spectral regions.

Description: No agreement on what constitutes a safe and reproducible anticontamination protocol exists for ancient starch research. Protocols applied to laboratory work may represent ‘symptomatic treatment’ only, as contamination of archaeological materials in the field may be more extensive than realized. This paper is the first systematic study on the impact that modern starches from surface and buried soils, windborne dispersal, human motion, excavation techniques and toolkits, and field attire has on archaeological sample quality. The study area is Olduvai Gorge, Tanzania. We identify seven starch types (discrete granules, n = 788) that embody the starch contamination landscape for the region. This study also demonstrates the various diagenetic changes that buried starch granules undergo in a short time, such as cavitation, fissuring, disruption and gelatinization. There are significant differences in morphotype class representation between the topsoil starches and those collected deeper below ground at excavated sites. Diagenetically transformed granules from underground storage organs dominate in soils, while native starches from cereal endosperm (Panicoideae and Triticeae) abound above ground in airborne samples. Furthermore, we illustrate how lithic samples excavated under standard field conditions can be contaminated, and that when a sample is compromised during excavation, it may be impossible to distinguish between target and introduced starches, especially when granules are identical or morphologically similar. The paper provides field recommendations to control false positives.

Description: A total of three binary tropical cyclone (TC) cases over the west North Pacific are selected to investigate the effects of satellite radiance data assimilation on analyses and forecasts of binary TCs. Two parallel cycling experiments with a 6-h interval are performed for each binary TC case, and the difference between the two experiments is whether satellite radiance observations are assimilated. Satellite radiance observations are assimilated using the Weather Research and Forecasting Data Assimilation (WRFDA)'s three dimensional variational (3D-Var) system, which includes the observation operator, quality control procedures, and bias correction algorithm for radiance observations. On average, radiance assimilation results in slight improvements of environmental fields and track forecasts of binary TC cases, but the detailed effects vary with the case. When there is no direct interaction between binary TCs, radiance assimilation leads to better depictions of environmental fields, and finally it results in improved track forecasts. However, positive effects of radiance assimilation on track forecasts can be reduced when there exists a direct interaction between binary TCs and intensities/structures of binary TCs are not represented well. An initialization method (e.g. dynamic initialization) combined with radiance assimilation and/or more advanced DA techniques (e.g. hybrid method) can be considered to overcome these limitations. This article is protected by copyright. All rights reserved.

Description: The phenology of diameter-growth cessation in trees will likely play a key role in mediating species and ecosystem responses to climate change. A common expectation is that warming will delay cessation, but the environmental and genetic influences on this process are poorly understood. We modeled the effects of temperature, photoperiod and seed-source climate on diameter-growth cessation timing in coast Douglas-fir (an ecologically and economically vital tree) using high-frequency growth measurements across broad environmental gradients for a range of genotypes from different seed sources. Our model suggests that cool temperatures or short photoperiods can induce cessation in autumn. At cool locations (high latitude and elevation), cessation seems to be induced primarily by low temperatures in early autumn (under relatively long photoperiods), so warming will likely delay cessation and extend the growing season. But at warm locations (low latitude or elevation), cessation seems to be induced primarily by short photoperiods later in autumn, so warming will likely lead to only slight extensions of the growing season, reflecting photoperiod limitations on phenological shifts. Trees from seed sources experiencing frequent frosts in autumn or early winter tended to cease growth earlier in the autumn, potentially as an adaptation to avoid frost. Thus, gene flow into populations in warm locations with little frost will likely have limited potential to delay mean cessation dates because these populations already cease growth relatively late. In addition, data from an abnormal heat wave suggested that very high temperatures during long photoperiods in early summer might also induce cessation. Climate change could make these conditions more common in warm locations, leading to much earlier cessation. Thus, photoperiod cues, patterns of genetic variation and summer heat waves could limit the capacity of coast Douglas-fir to extend its growing season in response to climate change in the warm parts of its range. This article is protected by copyright. All rights reserved.

Description: The urban environment modifies the hydrologic cycle resulting in increased runoff rates, volumes and peak flows. Green Infrastructure, which uses best management practices (BMPs), is a natural system approach used to mitigate the impacts of urbanization onto stormwater runoff. Patterns of stormwater runoff from urban environments are complex and it is unclear how efficiently green infrastructure will improve the urban water cycle. These challenges arise from issues of scale, the merits of BMPs depend on changes to small scale hydrologic processes aggregated up from the neighborhood to the urban watershed. Here, we use a hyper-resolution (1 m), physically-based hydrologic model of the urban hydrologic cycle with explicit inclusion of the built environment. This model represents the changes to hydrology at the BMP scale (~1 m), represents each individual BMP explicitly to represent response over the urban watershed. Our study varies both the percentage of BMP emplacement and their spatial location for storm events of increasing intensity in an urban watershed. We develop a metric of effectiveness that indicates a nonlinear relationship is seen between percent BMP emplacement and storm intensity. Results indicate that BMP effectiveness varies with spatial location and that type and emplacement within the urban watershed may be more important than overall percent.

Description: A number of extensive droughts and destructive floods have occurred in Poland in the last 25 years, hence projections of low and high river flows are of considerable interest and importance. In the first part of this paper, projections of low and high flows in the rivers of the Vistula and the Odra basins (VOB region), for two future time horizons, are presented. Projections are based on the SWAT hydrological model simulations driven by results of the EURO-CORDEX experiment under Representative Concentration Pathways 4.5 and 8.5. The VOB region covers most of Poland and parts of five neighboring countries, giving this study an international relevance. In the second part of the paper a review of projections of low and high flows in rivers in Central and Eastern Europe is presented. Despite a substantial spread of flow projections, the main message of the modelling part is that increases of both low and high flows are dominating. The magnitude of increase of low flow is considerably higher than that of high flow. In other words, future streamflow droughts are projected to be less severe, while, in contrast, river floods are projected to increase, which is a challenge for flood risk reduction, water management and climate change adaptation. There is an overall agreement of our findings for the VOB region with projections of hydrological extremes from large-scale models forced by EURO-CORDEX results in the European-scale studies.

Description: Because the capability of terrestrial ecosystems to fix carbon is constrained by nutrient availability, understanding how nutrients limit plant growth is a key contemporary question. However, what drives nutrient limitations at global scale remains to be clarified. Using global data on plant growth, plant nutritive status, and soil fertility, we investigated to which extent soil parent materials explain nutrient limitations. We found that N limitation was not linked to soil parent materials, but was best explained by climate: ecosystems under harsh (i.e. cold and or dry) climates were more N-limited than ecosystems under more favourable climates. Contrary to N limitation, P limitation was not driven by climate, but by soil parent materials. The influence of soil parent materials was the result of the tight link between actual P pools of soils and physical-chemical properties (acidity; P richness) of soil parent materials. Some other ground-related factors (i.e. soil weathering stage, landform) had a noticeable influence on P limitation, but their role appeared to be relatively smaller than that of geology. The relative importance of N limitation versus P limitation was explained by a combination of climate and soil parent material: at global scale, N limitation became prominent with increasing climatic constraints, but this global trend was modulated at lower scales by the effect of parent materials on P limitation, particularly under climates favourable to biological activity. As compared with soil parent materials, atmospheric deposition had only a weak influence on the global distribution of actual nutrient limitation. Our work advances our understanding of the distribution of nutrient limitation at global scale. In particular, it stresses the need to take soil parent materials into account when investigating plant growth response to environment changes. This article is protected by copyright. All rights reserved.

Description: Transformations of precipitation into groundwater and streamflow are fundamental hydrological processes, critical to irrigated agriculture, hydroelectric power generation and ecosystem health. Our understanding of the timing of groundwater recharge and streamflow generation remains incomplete, limiting our ability to predict fresh water, nutrient, and contaminant fluxes, especially in large basins. Here we analyze thousands of rain, snow, groundwater and streamflow δ 18 O and δ 2 H values in the Nelson River basin, which covers 1.2 million km 2 of central Canada. We show that the fraction of precipitation that recharges aquifers is ~1.3-5 times higher for precipitation falling during cold months with subzero mean monthly temperatures than for precipitation falling during warmer months. The near-ubiquity of cold-season-biased groundwater recharge implies that changes to winter water balances may have disproportionate impacts on annual groundwater recharge rates. We also show that young streamflow—defined as precipitation that enters a river in less than ~2.3 months—comprises ~27 % of annual streamflow, but varies widely among tributaries in the Nelson River basin (1-59 %). Young streamflow fractions are lower in steep catchments, and higher in flatter catchments such as the transboundary Red River basin. Our findings imply that flat, lower-permeability, heavily-tiled landscapes favor more rapid transmission of precipitation into rivers, possibly mobilizing excess soluble fertilizers and exacerbating eutrophication events in Lake Winnipeg.

Description: An operational system for producing a global diurnally varying analysis of skin Sea Surface Temperature (SST) has been developed at the Met Office. Skin SST is formulated as the sum of a foundation temperature, a warm layer temperature difference, and a thermal skin temperature difference. Foundation temperature is taken from the Operational Sea surface Temperature and sea Ice Analysis (OSTIA) system, while numerical models are used for the warm layer and thermal skin layer. Both the thermal skin layer and warm layer models are forced using outputs from the Met Office's numerical weather prediction system. Data assimilation is used to improve estimates of the warm layer, with observations coming from the geostationary SEVIRI and GOES-W instruments, as well as the polar orbiting NOAA-AVHRR sensors. The SST observations of these instruments are converted to observations of the warm layer temperature difference by subtracting an instrument specific foundation SST calculated using just night time data from that sensor. A quality control procedure removes warm layer observations where the satellite specific foundation estimates are less robust due to lack of data. Validation of the analysis system has been performed by comparing the model to the assimilated satellite data, and to independent near-surface observations from Argo floats. The mean state of the system was assessed via a comparison to a climatology generated from drifting buoys. Results from the validation show that the system does a good job of replicating the climatology and that assimilation improves the analysis when assessed against Argo. However, the system has been found to underestimate the diurnal range of skin SST by approximately 0.1-0.3 ° C on average. Data from the system is available free of charge from the Copernicus Marine Environment Monitoring Service.

Description: Ensemble methods are an attractive option for data assimilation (DA) in convection-permitting models, providing flow-dependent covariances which respect the complex balances that apply at these scales. The Met Office deterministic UK forecast currently uses variational assimilation based on predominantly climatological covariances, whilst the MOGREPS-UK ensemble interpolates its initial perturbations from the corresponding global members. An experimental convective-scale ensemble DA system has now been built based on a serial ensemble filter, where observations are assimilated sequentially. This has been tested in a series of one-month trials assimilating sonde, surface and aircraft observations, exploring issues such as cycle length, filter type, localisation and inflation. The tests use 44 perturbed ensemble members, plus an unperturbed control forecast that is updated using the same Kalman gain as the ensemble mean. The best configurations provide control forecasts which are competitive with interpolating global initial states like MOGREPS-UK did at the time, even though the local assimilation has seen fewer observation types than the global DA. An hourly cycle performs much better than a six-hourly one, as suggested by theoretical arguments, even though both see the same boundary conditions and essentially the same total set of observations. The serial filter calculates the difference between each observation and its model equivalent as updated by all previous observations. We call this the Innovation After assimilating Prior Observations (IAPO). The mean square IAPO measures the error of the developing analysis using observations which are statistically independent of it, without the need to involve the forecast model. This diagnostic confirms the beneficial impact of assimilation, can be used to cheaply tune parameters such as localisation radii, and contributes to evaluating the correctness of the ensemble spread.

Description: The ability of a newly developed non-hydrostatic regional climate model (RCM) based on the Unified Model of the UK Met Office (MetUM) at a resolution of 12 km is examined for the simulation of tropical cyclone (TC) activity affecting the South China Sea and compared with the current released RCM version of MetUM at the resolution of 25 km. The results show that both the 25 km and 12 km models can reasonably simulate the TC-associated large-scale environments, while the 12 km model has a better ability to simulate the South Asian monsoon. Compared with the 25 km model, the 12 km model generally improves the simulation of track density and the radial wind structure of TCs. However, the annual cycle of simulated TCs show that both models tend to over-estimate the TC frequency in May and November–January while underestimating the frequency in June-September. Compared with the 25 km model, the 12 km model produces fewer intense TCs with 10 m maximum wind speeds 〉 30 m s −1 . It is also found that both the 12 km and 25 km models reproduce the observed modulation of TC activity associated with different phases of the El Niño/Southern Oscillation (ENSO) such as the reduced track density and accumulated cyclonic energy during El Niño events, while the 12 km model better captures the TC-ENSO response including the track density and the large-scale environments than the 25 km model.

Description: The Soil Moisture Active Passive (SMAP) satellite makes coincident global measurements of soil moisture using an L-band radar instrument and an L-band radiometer. It is crucial to evaluate the errors in the newest L-band SMAP satellite-derived soil moisture products, before they are routinely used in scientific research and applications. This study represents the first evaluation of the SMAP radiometer soil moisture product over China. In this paper, a preliminary evaluation was performed using sparse in situ measurements from 655 China Meteorological Administration (CMA) monitoring stations between 1 April 2015 and 31 August 2016. The SMAP radiometer-derived soil moisture product was evaluated against two schemes of original soil moisture and the soil moisture anomaly in different geographical zones and land cover types. Four performance metrics, i.e., bias, root mean square error (RMSE), unbiased root mean square error (ubRMSE), and the correlation coefficient (R), were used in the accuracy evaluation. The results indicated that the SMAP radiometer-derived soil moisture product agreed relatively well with the in situ measurements, with ubRMSE values of 0.058 cm3·cm−3 and 0.039 cm3·cm−3 based on original data and anomaly data, respectively. The values of the SMAP radiometer-based soil moisture product were overestimated in wet areas, especially in the Southwest China, South China, Southeast China, East China, and Central China zones. The accuracies over croplands and in Northeast China were the worst. Soil moisture, surface roughness, and vegetation are crucial factors contributing to the error in the soil moisture product. Moreover, radio frequency interference contributes to the overestimation over the northern portion of the East China zone. This study provides guidelines for the application of the SMAP-derived soil moisture product in China and acts as a reference for improving the retrieval algorithm.

Description: Agricultural land use results in multiple stressors affecting stream ecosystems. Flow reduction due to water abstraction, elevated levels of nutrients and chemical contaminants are common agricultural stressors worldwide. Concurrently, stream ecosystems are also increasingly affected by climate change. Interactions among multiple co-occurring stressors result in biological responses that cannot be predicted from single-stressor effects (i.e. synergisms and antagonisms). At the ecosystem level, multiple-stressor effects can be further modified by biotic interactions (e.g. trophic interactions). We conducted a field experiment using 128 flow-through stream mesocosms to examine the individual and combined effects of water abstraction, nutrient enrichment and elevated levels of the nitrification inhibitor dicyandiamide (DCD) on survival, condition and gut content of juvenile brown trout, and on benthic abundance of their invertebrate prey. Flow velocity reduction decreased fish survival (-12% compared to controls) and condition (-8% compared to initial condition), whereas effects of nutrient and DCD additions and interactions among these stressors were not significant. Negative effects of flow velocity reduction on fish survival and condition were consistent with effects on fish gut content (-25% compared to controls) and abundance of dominant invertebrate prey (-30%), suggesting a negative metabolic balance driving fish mortality and condition decline, which was confirmed by structural equation modelling. Fish mortality under reduced flow velocity increased as maximal daily water temperatures approached the upper limit of their tolerance range, reflecting synergistic interactions between these stressors. Our study highlights the importance of indirect stressor effects such as those transferred through trophic interactions, which need to be considered when assessing and managing fish populations and stream food webs in multiple-stressor situations. However, in real streams compensatory mechanisms, behavioural responses, as well as seasonal and spatial variation may alter the intensity of stressor effects and the sensitivity of trout populations. This article is protected by copyright. All rights reserved.

Description: Saturation-excess runoff is the major runoff mechanism in humid well-vegetated areas where infiltration rates often exceed rainfall intensity. While the Soil and Water Assessment Tool (SWAT) is one of the most widely used models, it predicts runoff based mainly on soil and land use characteristics, and is implicitly an infiltration-excess runoff type of model. Previous attempts to incorporate the saturation-excess runoff mechanism in SWAT fell short due to the inability to distribute water from one Hydrological Response Unit (HRU) to another. This paper introduces a modified version of SWAT, referred to as SWAT-Hillslope (SWAT-HS). This modification improves the simulation of saturation-excess runoff by redefining HRUs based on wetness classes, and by introducing a surface aquifer with the ability to route interflow from “drier” to “wetter” wetness classes. Mathematically, the surface aquifer is a non-linear reservoir that generates rapid subsurface stormflow as the water table in the surface aquifer rises. The SWAT-HS model was tested in the Town Brook watershed in the upper reaches of the West Branch Delaware River in the Catskill region of New York, USA. SWAT-HS predicted discharge well with a Nash-Sutcliffe Efficiency of 0.68 and 0.87 for daily and monthly time steps. Compared to the original SWAT model, SWAT-HS predicted less surface runoff and groundwater flow and more lateral flow. The saturated areas predicted by SWAT-HS were concentrated in locations with a high topographic index and were in agreement with field observations. With the incorporation of topographic characteristics and the addition of the surface aquifer, SWAT-HS improved streamflow simulation and gave a good representation of saturated areas on the date that measurements were available. SWAT-HS is expected to improve water quality model predictions where the location of the surface runoff matters.

Description: Nano-silver and nano-titanium oxide films can be coated over brackets in order to reduce bacterial aggregation and friction. However, their antimicrobial efficacy, surface roughness, and frictional resistance are not assessed before. Fifty-five stainless-steel brackets were divided into 5 groups of 11 brackets each: uncoated brackets, brackets coated with 60 µm silver, 100 µm silver, 60 µm titanium, and 100 µm titanium. Coating was performed using physical vapor deposition method. For friction test, three brackets from each group were randomly selected and tested. For scanning electron microscopy and atomic-force microscopy assessments, one and one brackets were selected from each group. For antibacterial assessment, six brackets were selected from each group. Of them, three were immediately subjected to direct contact with S. mutans . Colonies were counted 3, 6, 24, and 48 h of contact. The other three were stored in water for 3 months. Then were subjected to a similar direct contact test. Results pertaining to both subgroups were combined. Groups were compared statistically. Mean (SD) friction values of the groups 'control, silver-60, silver-100, titanium-60, and titanium-100' were 0.55 ± 0.14, 0.77 ± 0.08, 0.82 ± 0.11, 1.52 ± 0.24, and 1.57 ± 0.41 N, respectively ( p  =   .0004, Kruskal–Wallis). Titanium frictions were significantly greater than control ( p  〈   .05), but silver groups were not ( p  〉   .05, Dunn). In the uncoated group, colony count increased exponentially within 48 h. The coated groups showed significant reductions in colony count ( p  〈   .05, two-way-repeated-measures ANOVA). In conclusions, all four explained coatings reduce surface roughness and bacterial growth. Nano-titanium films are not suitable for friction reduction. Nano-silver results were not conclusive and need future larger studies.

Description: Many papers have claimed the attainment of super-resolution, i.e. resolution beyond that achieved classically, by measurement of the profile of a feature in the image. We argue that measurement of the contrast of the image of a dark bar on a bright background does not give a measure of resolution, but of detection sensitivity. The width of a bar that gives an intensity at the center of the bar of 0.735 that in the bright region (the same ratio as in the Rayleigh resolution criterion) is for the coherent case with central illumination. This figure, which compares with for the Abbe resolution limit with central illumination, holds for the classical case, and so is no indication of super-resolution. Theoretical images for two points, two lines, arrays of lines, arrays of bars, and grating objects are compared. These results can be used a reference for experimental results, to determine if super-resolution has indeed been attained. The history of the development of the theory of microscope resolution is outlined. This shows the image of a dark bar on a coherent bright background, for different bar widths. There is no sharp cutoff. There is always a dip in intensity at the center, even for a bar width of, corresponding to 0.05 Airy units.

Description: As modeling capabilities at regional and global scales improve, questions remain regarding the appropriate process representation required to accurately simulate multichannel river hydraulics. This study uses the hydrodynamic model LISFLOOD-FP to simulate patterns of water surface elevation (WSE), depth, and inundation extent across a ∼90 km, anabranching reach of the Tanana River, Alaska. To provide boundary conditions, we collected field observations of bathymetry and WSE during a two-week field campaign in summer 2013. For the first time at this scale, we test a simple, raster-based model's capabilities to simulate 2D, in-channel patterns of WSE and inundation extent. Additionally, we compare finer resolution (≤ 25 m) 2D models to four other models of lower dimensionality and coarser resolution (100–500 m) to determine the effects of simplifying process representation. Results indicate that simple, raster-based models can accurately simulate 2D, in-channel hydraulics in the Tanana. Also, the fine-resolution, 2D models produce lower errors in spatiotemporal outputs of WSE and inundation extent compared to coarse-resolution, 1D models: 22.6 cm vs. 56.4 cm RMSE for WSE, and 90% vs. 41% Critical Success Index values for simulating inundation extent. Incorporating the anabranching channel network using subgrid representations for smaller channels is important for simulating accurate hydraulics and lowers RMSE in spatially distributed WSE by at least 16%. As a result, better representation of the converging and diverging multichannel network by using subgrid solvers or downscaling techniques in multichannel rivers is needed to improve errors in regional to global scale models. This article is protected by copyright. All rights reserved.

Description: Attempts to estimate the influence of erosion on the carbon (C) cycle are limited by difficulties in accounting for the fate of mobilized organic material and for the uncertainty associated with land management practices. This study proposes a method to quantify the uncertainty introduced by the influence of land management on soil organic C (SOC) generation and decomposition at eroding soils. The framework is implemented in tRIBS-ECO (Triangulated Irregular Network-based Real-time Integrated Basin Simulator-Erosion and Carbon Oxidation). tRIBS-ECO is a spatially- and depth-explicit model of C dynamics coupled with a process-based hydro-geomorphic model. We assess the impact of soil erosion on the net soil-atmosphere CO 2 exchange at the Calhoun Critical Zone Observatory, one of the most severely agriculturally eroded regions in the U.S. Measurements of SOC storage are used from different catena positions. We demonstrate that the spatio-temporal variations of land management practices introduce significant uncertainty in estimates of the erosion-induced CO 2 exchange with the atmosphere. Observations and simulations suggest that a substantial portion of eroded organic material is buried in alluvial sediments at the study site. According to results, recent reforestation led to a partial decline in soil and SOC erosion rates. It is suggested that the representation of the fine spatio-temporal variability of the dynamics of eroded C is important in the computation of C budgets in regional and global scales. This article is protected by copyright. All rights reserved.

Description: Understanding how historical processes modulate the response of ecosystems to perturbations is becoming increasingly important. In contrast to the growing interest in projecting biodiversity and ecosystem functioning under future climate scenarios, how legacy effects originating from historical conditions drive change in ecosystems remains largely unexplored. Using experiments in combination with Stochastic Antecedent Modelling, we evaluated how extreme warming, sediment deposition and grazing events modulated the ecological memory of rocky intertidal epilithic microphytobenthos (EMPB). We found memory effects in the non-clustered scenario of disturbance (60 days apart), where EMPB biomass fluctuated in time, but not under clustered disturbances (15 days apart), where EMPB biomass was consistently low. A massive grazing event impacted on EMPB biomass in a second run of the experiment, also muting ecological memory. Our results provide empirical support to the theoretical expectation that stochastic fluctuations promote ecological memory, but also show that contingencies may lead to memory loss. This article is protected by copyright. All rights reserved.

Description: Vegetation is often represented by the leaf area index (LAI) in many ecological, hydrological and meteorological land surface models. However, the spatio-temporal dynamics of the vegetation are important to represent in these models. While the widely applied methods, such as the Canopy Structure Dynamic Model (CSDM) and the Double Logistic Model (DLM) are solely based on cumulative daily mean temperature data as input, a new spatio-temporal LAI prediction model referred to as the Temperature Precipitation Vegetation Model (TPVM) is developed that also considers cumulative precipitation data as input into the modelling process. TPVM as well as CDSM and DLM model performances are compared and evaluated against filtered LAI data from the Moderate Resolution Imaging Spectroradiometer (MODIS). The calibration/validation results of a cross-validation performed in the meso-scale Attert catchment in Luxembourg indicated that the DLM and TPVM generally provided more realistic and accurate LAI data. The TPVM performed superiorly for the agricultural land cover types compared to the other two models, which only used the temperature data. The Pearson's correlation coefficient (CC) between TPVM and the field measurement is 0.78, compared to 0.73 and 0.69 for the DLM and CSDM, respectively. The phenological metrics were derived from the TPVM model to investigate the interaction between the climate variables and the LAI variations. These interactions illustrated the dominant control of temperature on the LAI dynamics for deciduous forest cover, and a combined influence of temperature with precipitation for the agricultural land use areas.

Description: Water vapour (H2O) is the dominant species in volcanic gas plumes. Therefore, measurements of H2O fluxes could provide valuable constraints on subsurface degassing and magmatic processes. However, due to the large and variable concentration of this species in the background atmosphere, little attention has been devoted to monitoring the emission rates of this species from volcanoes. Instead, the focus has been placed on remote measurements of SO2, which is present in far lower abundances in plumes, and therefore provides poorer single flux proxies for overall degassing conditions. Here, we present a new technique for the measurement of H2O emissions at degassing volcanoes at high temporal resolution (≈1 Hz), via remote sensing with low cost digital cameras. This approach is analogous to the use of dual band ultraviolet (UV) cameras for measurements of volcanic SO2 release, but is focused on near infrared absorption by H2O. We report on the field deployment of these devices on La Fossa crater, Vulcano Island, and the North East Crater of Mt. Etna, during which in-plume calibration was performed using a humidity sensor, resulting in estimated mean H2O fluxes of ≈15 kg·s−1 and ≈34 kg·s−1, respectively, in accordance with previously reported literature values. By combining the Etna data with parallel UV camera and Multi-GAS observations, we also derived, for the first time, a combined record of 1 Hz gas fluxes for the three most abundant volcanic gas species: H2O, CO2, and SO2. Spectral analysis of the Etna data revealed oscillations in the passive emissions of all three species, with periods spanning ≈40–175 s, and a strong degree of correlation between the periodicity manifested in the SO2 and H2O data, potentially related to the similar exsolution depths of these two gases. In contrast, there was a poorer linkage between oscillations in these species and those of CO2, possibly due to the deeper exsolution of carbon dioxide, giving rise to distinct periodic degassing behaviour.

Description: Identifying individual trees and delineating their canopy structures from the forest point clouddataacquiredbyanairborneLiDAR(LightDetectionAndRanging)hassignificantimplications in forestry inventory. Once accurately identified, tree structural attributes such as tree height, crown diameter, canopy based height and diameter at breast height can be derived. This paper focuses on a novel computationally efficient method to adaptively calibrate the kernel bandwidth of a computational scheme based on mean shift—a non-parametric probability density-based clustering technique—to segment the 3D (three-dimensional) forest point clouds and identify individual tree crowns. The basic concept of this method is to partition the 3D space over each test plot into small vertical units (irregular columns containing 3D spatial features from one or more trees) first, by using a fixed bandwidth mean shift procedure and a small square grouping technique, and then rough estimation of crown sizes for distinct trees within a unit, based on an original 2D (two-dimensional) incremental grid projection technique, is applied to provide a basis for dynamical calibration of the kernel bandwidth for an adaptive mean shift procedure performed in each partition. The adaptive mean shift-based scheme, which incorporates our proposed bandwidth calibration method, is validated on 10 test plots of a dense, multi-layered evergreen broad-leaved forest located in South China. Experimental results reveal that this approach can work effectively and when compared to the conventional point-based approaches (e.g., region growing, k-means clustering, fixed bandwidth or multi-scale mean shift), its accuracies are relatively high: it detects 86 percent of the trees (“recall”) and 92 percent of the identified trees are correct (“precision”), showing good potential for use in the area of forest inventory.

Description: A 3D refractive-index matching Lagrangian particle tracking (3D-RIM-LPT) system was developed to study the filtration and clogging process inside a homogeneous porous medium. A small subset of particles flowing through the porous medium was dyed and tracked. As this subset was randomly chosen, its dynamics is representative of all the rest. The statistics of particle locations, number, and velocity vectors were obtained as functions of different volumetric concentrations. It is found that, in our system, the clogging time decays with particle concentration following a power law relationship. As the concentration increases, there is a transition from depth filtration to cake filtration. At high concentration, more clogged pores lead to frequent flow redirections and more transverse migrations of particles. In addition, the velocity distribution of the transverse direction is symmetrical around zero, and it is slightly more intermittent than the random Gaussian curve due to particle-particle interactions and particle-grain interactions. In contrast, as clogging develops, the longitudinal velocity of particles along the main flow direction has peak near zero because of those trapped particles. But at the same time, the remaining open pores will experience larger pressure and, as a result, particles through those pores will have a larger longitudinal velocity. This article is protected by copyright. All rights reserved.

Description: Precipitation patterns are changing across the globe causing more severe and frequent drought for many forest ecosystems. Although research has focused on the resistance of tree populations and communities to these novel precipitation regimes, resilience of forests is also contingent on recovery following drought, which remains poorly understood, especially in aseasonal tropical forests. We used rainfall exclusion shelters to manipulate the inter-annual frequency of drought for diverse seedling communities in a tropical forest and assessed resistance, recovery and resilience of seedling growth and mortality relative to everwet conditions. We found seedlings exposed to recurrent periods of drought altered their growth rates throughout the year relative to seedlings in everwet conditions. During drought periods seedlings grew slower than seedlings in everwet conditions (i.e. resistance phase) while compensating with faster growth after drought (i.e. recovery phase). However, the response to frequent drought was species dependent as some species grew significantly slower with frequent drought relative to everwet conditions while others grew faster with frequent drought due to overcompensating growth during the recovery phase. In contrast, mortality was unrelated to rainfall conditions and instead correlated with differences in light. Intra-annual plasticity of growth and increased annual growth of some species led to an overall maintenance of growth rates of tropical seedling communities in response to more frequent drought. These results suggest these communities can potentially adapt to predicted climate change scenarios and that plasticity in the growth of species, and not solely changes in mortality rates among species, may contribute to shifts in community composition under drought. This article is protected by copyright. All rights reserved.

Description: The Advanced Along-Track Scanning Radiometer (AATSR) land surface temperature (LST) product has a long-term time series of data from 20 May 2002 to 8 April 2012 and is a crucial dataset for global change studies. Accuracy and uncertainty assessment of satellite derived LST is important for its use in studying land–surface–atmosphere interactions. However, the validation of AATSR-derived LST products is scarce in China, especially in arid and semi-arid areas. In this study, we evaluated the accuracy of the AATSR LST product using ground-based measurements from 2007 to 2011 in the Heihe River Basin (HRB), China. The AATSR-derived LST results over Yingke site are closer to ground measurements than those over A’rou site for both daytime and nighttime temperatures. For nighttime, the averaged bias, STD, RMSE and R2 over both sites are 0.67 K, 3.03 K, 3.13 K and 0.93 K, respectively. Based on the accuracy assessment, we analyzed the AATSR-derived annual LST variations both in the HRB region and the two validation sites for the period of 2003 to 2011. The results at the A’rou site show an obvious increasing trend for daytime from 2003 to 2011. For the whole HRB region, the warming trend is clearly shown in the downstream of HRB.

Description: Recently, we reported on the development of low-cost ultraviolet (UV) cameras, based on the modification of sensors designed for the smartphone market. These units are built around modified Raspberry Pi cameras (PiCams; ≈USD 25), and usable system sensitivity was demonstrated in the UVA and UVB spectral regions, of relevance to a number of application areas. Here, we report on the first deployment of PiCam devices in one such field: UV remote sensing of sulphur dioxide emissions from volcanoes; such data provide important insights into magmatic processes and are applied in hazard assessments. In particular, we report on field trials on Mt. Etna, where the utility of these devices in quantifying volcanic sulphur dioxide (SO2) emissions was validated. We furthermore performed side-by-side trials of these units against scientific grade cameras, which are currently used in this application, finding that the two systems gave virtually identical flux time series outputs, and that signal-to-noise characteristics of the PiCam units appeared to be more than adequate for volcanological applications. Given the low cost of these sensors, allowing two-filter SO2 camera systems to be assembled for ≈USD 500, they could be suitable for widespread dissemination in volcanic SO2 monitoring internationally.

Description: The assessment of the soil redistribution and real long-term soil degradation due to erosion on agriculture land is still insufficient in spite of being essential for soil conservation policy. Imaging spectroscopy has been recognized as a suitable tool for soil erosion assessment in recent years. In our study, we bring an approach for assessment of soil degradation by erosion by means of determining soil erosion classes representing soils differently influenced by erosion impact. The adopted methods include extensive field sampling, laboratory analysis, predictive modelling of selected soil surface properties using aerial hyperspectral data and the digital elevation model and fuzzy classification. Different multivariate regression techniques (Partial Least Square, Support Vector Machine, Random forest and Artificial neural network) were applied in the predictive modelling of soil properties. The properties with satisfying performance (R2 > 0.5) were used as input data in erosion classes determination by fuzzy C-means classification method. The study was performed at four study sites about 1 km2 large representing the most extensive soil units of the agricultural land in the Czech Republic (Chernozems and Luvisols on loess and Cambisols and Stagnosols on crystalline rocks). The influence of site-specific conditions on prediction of soil properties and classification of erosion classes was assessed. The prediction accuracy (R2) of the best performing models predicting the soil properties varies in range 0.8–0.91 for soil organic carbon content, 0.21–0.67 for sand content, 0.4–0.92 for silt content, 0.38–0.89 for clay content, 0.73–089 for Feox, 0.59–0.78 for Fed and 0.82 for CaCO3. The performance and suitability of different properties for erosion classes’ classification are highly variable at the study sites. Soil organic carbon was the most frequently used as the erosion classes’ predictor, while the textural classes showed lower applicability. The presented approach was successfully applied in Chernozem and Luvisol loess regions where the erosion classes were assessed with a good overall accuracy (82% and 67%, respectively). The model performance in two Cambisol/Stagnosol regions was rather poor (51%–52%). The results showed that the presented method can be directly and with a good performance applied in pedologically and geologically homogeneous areas. The sites with heterogeneous structure of the soil cover and parent material will require more precise local-fitted models and use of further auxiliary information such as terrain or geological data. The future application of presented approach at a regional scale promises to produce valuable data on actual soil degradation by erosion usable for soil conservation policy purposes.

Description: The methods used for image contrast enhancement in the wavelet domain have been previously documented. The essence of these methods lies in the manipulation of the image during the reconstruction process, by changing the relationship between the components that require transformation. This paper proposes a new variant based on using undecimated wavelet transform and adapting the Gaussian function for scaling the coefficients of detail wavelet components, so that the role of low coefficients in the reconstructed image is greater. The enhanced image is then created by combining the new components. Applying the Haar wavelet minimises the effects of the relationship disturbance between components, and creates a small buffer around the edge. The proposed method was tested using six images at different scales, collected with handheld photo cameras, and aerial and satellite optical sensors. The results of the tests indicate that the method can achieve comparable, or even better enhancement effects for weak edges, than the well-known unsharp masking and Retinex methods. The proposed method can be applied in order to improve the visual interpretation of remote sensing images taken by various sensors at different scales.

Description: Remote Sensing, Vol. 9, Pages 820: Soil Moisture Data Assimilation in a Hydrological Model: A Case Study in Belgium Using Large-Scale Satellite Data Remote Sensing doi: 10.3390/rs9080820 Authors: Pierre Baguis Emmanuel Roulin In the present study, we focus on the assimilation of satellite observations for Surface Soil Moisture (SSM) in a hydrological model. The satellite data are produced in the framework of the EUMETSAT project H-SAF and are based on measurements with the Advanced radar Scatterometer (ASCAT), embarked on the Meteorological Operational satellites (MetOp). The product generated with these measurements has a horizontal resolution of 25 km and represents the upper few centimeters of soil. Our approach is based on the Ensemble Kalman Filter technique (EnKF), where observation and model uncertainties are taken into account, implemented in a conceptual hydrological model. The analysis is carried out in the Demer catchment of the Scheldt River Basin in Belgium, for the period from June 2013–May 2016. In this context, two methodological advances are being proposed. First, the generation of stochastic terms, necessary for the EnKF, of bounded variables like SSM is addressed with the aid of specially-designed probability distributions, so that the bounds are never exceeded. Second, bias due to the assimilation procedure itself is removed using a post-processing technique. Subsequently, the impact of SSM assimilation on the simulated streamflow is estimated using a series of statistical measures based on the ensemble average. The differences from the control simulation are then assessed using a two-dimensional bootstrap sampling on the ensemble generated by the assimilation procedure. Our analysis shows that data assimilation combined with bias correction can improve the streamflow estimations or, at a minimum, produce results statistically indistinguishable from the control run of the hydrological model.

Description: Remote Sensing, Vol. 9, Pages 823: Low-Altitude Aerial Methane Concentration Mapping Remote Sensing doi: 10.3390/rs9080823 Authors: Bara Emran Dwayne Tannant Homayoun Najjaran Detection of leaks of fugitive greenhouse gases (GHGs) from landfills and natural gas infrastructure is critical for not only their safe operation but also for protecting the environment. Current inspection practices involve moving a methane detector within the target area by a person or vehicle. This procedure is dangerous, time consuming, labor intensive and above all unavailable when access to the desired area is limited. Remote sensing by an unmanned aerial vehicle (UAV) equipped with a methane detector is a cost-effective and fast method for methane detection and monitoring, especially for vast and remote areas. This paper describes the integration of an off-the-shelf laser-based methane detector into a multi-rotor UAV and demonstrates its efficacy in generating an aerial methane concentration map of a landfill. The UAV flies a preset flight path measuring methane concentrations in a vertical air column between the UAV and the ground surface. Measurements were taken at 10 Hz giving a typical distance between measurements of 0.2 m when flying at 2 m/s. The UAV was set to fly at 25 to 30 m above the ground. We conclude that besides its utility in landfill monitoring, the proposed method is ready for other environmental applications as well as the inspection of natural gas infrastructure that can release methane with much higher concentrations.

Description: Globally, long-term research is critical to monitor the responses of tropical species to climate and land cover change at the range scale. Citizen science surveys can reveal the long-term persistence of poorly known nomadic tropical birds occupying fragmented forest patches. We applied dynamic occupancy models to 13 years (2002-2014) of citizen science driven presence/absence data on Cape parrot ( Poicephalus robustus ), a food nomadic bird endemic to South Africa. We modelled its underlying range dynamics as a function of resource distribution, and change in climate and land cover through the estimation of colonization and extinction patterns. The range occupancy of Cape parrot changed little over time (ψ = 0.75-0.83) because extinction was balanced by recolonization. Yet, there was considerable regional variability in occupancy and detection probability increased over the years. Colonizations increased with warmer temperature and area of orchards, thus explaining their range shifts southeastwards in recent years. Although colonizations were higher in the presence of nests and yellowwood trees ( Afrocarpus and Podocarpus spp .), the extinctions in small forest patches (≤ 227 ha) and during low precipitation (≤ 41 mm) are attributed to resource constraints and unsuitable climatic conditions. Loss of indigenous forest cover, and artificial lake/water bodies increased extinction probabilities of Cape parrot. The land use matrix (fruit farms, gardens and cultivations) surrounding forest patches provides alternative food sources thereby facilitating spatiotemporal colonization and extinction in the human-modified matrix. Our models show that Cape parrots are vulnerable to extreme climatic conditions such as drought which is predicted to increase under climate change. Therefore, management of optimum sized high quality forest patches is essential for long-term survival of Cape parrot populations. Our novel application of dynamic occupancy models to long-term citizen science monitoring data unfolds the complex relationships between the environmental dynamics and range fluctuations of this food nomadic species. This article is protected by copyright. All rights reserved.

Description: A method is presented for estimating the error covariance of the errors in the model equations in observation space. Estimating model errors in this systematic way opens up the possibility to use data assimilation for systematic model improvement at the level of the model equations, which would be a huge step forward. This model error covariance is perhaps the hardest covariance matrix to estimate. It represents how the missing physics and errors in parameterisations manifest themselves at the scales the model can resolve. A new element is that we use an efficient particle filter to avoid the need to estimate the error covariance of the state as well, which most other data-assimilation methods do require. Starting from a reasonable first estimate, the method generates new estimates iteratively during the data assimilation run, and the method is shown to converge to the correct model error matrix. We also investigate the influence of the accuracy of the observation error covariance on the estimation of the model error covariance and show that when the observation errors are known the model error covariance can be estimated well, but, as expected and perhaps unavoidable, the diagonal elements are estimated too low when the observation errors are estimated too high, and vice versa.

Description: ABSTRACT The Met Office developed a prototype high resolution, hourly-cycling NWP forecasting system using four-dimensional variational data assimilation, covering the southern half of England. This system was known as the Met Office Nowcasting Demonstration Project (NDP), and run in real time from June 2012 to March 2013 (covering the period of the London 2012 Olympics and Paralympic Games). It was principally developed to assess the abilities of a frequent update NWP-based nowcasting system to forecast precipitation at short range such as convective storms for flood forecasting. This article provides a comprehensive assessment of the precipitation forecast skill of the NDP and its comparison against the Met-Office operational nowcasting system. It was found that the NWP-based nowcasting system becomes more skillful than an advanced extrapolation-based nowcasting system from T+1.5 to T+2 depending on weather type.

Description: Remote Sensing, Vol. 9, Pages 815: Validation of Automatically Generated Global and Regional Cropland Data Sets: The Case of Tanzania Remote Sensing doi: 10.3390/rs9080815 Authors: Juan Laso Bayas Linda See Christoph Perger Christina Justice Catherine Nakalembe Jan Dempewolf Steffen Fritz There is a need to validate existing global cropland maps since they are used for different purposes including agricultural monitoring and assessment. In this paper we validate three recent global products (ESA-CCI, GlobeLand30, FROM-GC) and one regional product (Tanzania Land Cover 2010 Scheme II) using a validation data set that was collected by students through the Geo-Wiki tool. The ultimate aim was to understand the usefulness of these products for agricultural monitoring. Data were collected wall-to-wall for Kilosa district and for a sample across Tanzania. The results show that the amount of and spatial extent of cropland in the different products differs considerably from 8% to 42% for Tanzania, with similar values for Kilosa district. The agreement of the validation data with the four different products varied between 36% and 54% and highlighted that cropland is overestimated by the ESA-CCI and underestimated by FROM-GC. The validation data were also analyzed for consistency between the student interpreters and also compared with a sample interpreted by five experts for quality assurance. Regarding consistency between the students, there was more than 80% agreement if one difference in cropland category was considered (e.g., between low and medium cropland) while most of the confusion with the experts was also within one category difference. In addition to the validation of current cropland products, the data set collected by the students also has potential value as a training set for improving future cropland products.

Description: Remote Sensing, Vol. 9, Pages 816: Mapping Aboveground Carbon in Oil Palm Plantations Using LiDAR: A Comparison of Tree-Centric versus Area-Based Approaches Remote Sensing doi: 10.3390/rs9080816 Authors: Matheus Nunes Robert Ewers Edgar Turner David Coomes Southeast Asia is the epicentre of world palm oil production. Plantations in Malaysia have increased 150% in area within the last decade, mostly at the expense of tropical forests. Maps of the aboveground carbon density (ACD) of vegetation generated by remote sensing technologies, such as airborne LiDAR, are vital for quantifying the effects of land use change for greenhouse gas emissions, and many papers have developed methods for mapping forests. However, nobody has yet mapped oil palm ACD from LiDAR. The development of carbon prediction models would open doors to remote monitoring of plantations as part of efforts to make the industry more environmentally sustainable. This paper compares the performance of tree-centric and area-based approaches to mapping ACD in oil palm plantations. We find that an area-based approach gave more accurate estimates of carbon density than tree-centric methods and that the most accurate estimation model includes LiDAR measurements of top-of-canopy height and canopy cover. We show that tree crown segmentation is sensitive to crown density, resulting in less accurate tree density and ACD predictions, but argue that tree-centric approach can nevertheless be useful for monitoring purposes, providing a method to detect, extract and count oil palm trees automatically from images.

Description: Remote Sensing, Vol. 9, Pages 818: Water Optics and Water Colour Remote Sensing Remote Sensing doi: 10.3390/rs9080818 Authors: Yunlin Zhang Claudia Giardino Linhai Li The editorial paper aims to highlight the main topics investigated in the Special Issue (SI) “Water Optics and Water Colour Remote Sensing”. The outcomes of the 21 papers published in the SI are presented, along with a bibliometric analysis in the same field, namely, water optics and water colour remote sensing. This editorial summarises how the research articles of the SI approach the study of bio-optical properties of aquatic systems, the development of remote sensing algorithms, and the application of time-series satellite data for assessing long-term and temporal-spatial dynamics in inland, coastal, and oceanic waters. The SI shows the progress with a focus on: (1) bio-optical properties (three papers); (2) atmospheric correction and data uncertainties (five papers); (3) remote sensing estimation of chlorophyll-a (Chl-a) (eight papers); (4) remote sensing estimation of suspended matter and chromophoric dissolved organic matter (CDOM) (four papers); and (5) water quality and water ecology remote sensing (four papers). Overall, the SI presents a variety of applications at the global scale (with case studies in Europe, Asia, South and North America, and the Antarctic), achieved with different remote sensing instruments, such as hyperspectral field and airborne sensors, ocean colour radiometry, geostationary platforms, and the multispectral Landsat and Sentinel-2 satellites. The bibliometric analysis, carried out to include research articles published from 1900 to 2016, indicates that “chlorophyll-a”, “ocean colour”, “phytoplankton”, “SeaWiFS” (Sea-Viewing Wide Field-of-View Sensor), and “chromophoric dissolved organic matter” were the five most frequently used keywords in the field. The SI contents, along with the bibliometric analysis, clearly suggest that remote sensing of Chl-a is one of the topmost investigated subjects in the field.

Description: Remote Sensing, Vol. 9, Pages 811: Radiometric Cross-Calibration of GF-4 PMS Sensor Based on Assimilation of Landsat-8 OLI Images Remote Sensing doi: 10.3390/rs9080811 Authors: Yepei Chen Kaimin Sun Deren Li Ting Bai Chengquan Huang Earth observation data obtained from remote sensors must undergo radiometric calibration before use in quantitative applications. However, the large view angles of the panchromatic multispectral sensor (PMS) aboard the GF-4 satellite pose challenges for cross-calibration due to the effects of atmospheric radiation transfer and the bidirectional reflectance distribution function (BRDF). To address this problem, this paper introduces a novel cross-calibration method based on data assimilation considering cross-calibration as an optimal approximation problem. The GF-4 PMS was cross-calibrated with the well-calibrated Landsat-8 Operational Land Imager (OLI) as the reference sensor. In order to correct unequal bidirectional reflection effects, an adjustment factor for the BRDF was established, making complex models unnecessary. The proposed method employed the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm to find the optimal calibration coefficients and BRDF adjustment factor through an iterative process. The validation results revealed a surface reflectance error of <5% for the new cross-calibration coefficients. The accuracy of calibration coefficients were significantly improved when compared to the officially published coefficients as well as those derived using conventional methods. The uncertainty produced by the proposed method was less than 7%, meeting the demands for future quantitative applications and research. This method is also applicable to other sensors with large view angles.

Description: Remote Sensing, Vol. 9, Pages 817: Estimation of Satellite-Based SO42− and NH4+ Composition of Ambient Fine Particulate Matter over China Using Chemical Transport Model Remote Sensing doi: 10.3390/rs9080817 Authors: Yidan Si Shenshen Li Liangfu Chen Chao Yu Wende Zhu Epidemiologic and health impact studies have examined the chemical composition of ambient PM2.5 in China but have been constrained by the paucity of long-term ground measurements. Using the GEOS-Chem chemical transport model and satellite-derived PM2.5 data, sulfate and ammonium levels were estimated over China from 2004 to 2014. A comparison of the satellite-estimated dataset with model simulations based on ground measurements obtained from the literature indicated our results are more accurate. Using satellite-derived PM2.5 data with a spatial resolution of 0.1 × 0.1°, we further presented finer satellite-estimated sulfate and ammonium concentrations in anthropogenic polluted regions, including the NCP (the North China Plain), the SCB (the Sichuan Basin) and the PRD (the Pearl River Delta). Linear regression results obtained on a national scale yielded an r value of 0.62, NMB of −35.9%, NME of 48.2%, ARB_50% of 53.68% for sulfate and an r value of 0.63, slope of 0.67, and intercept of 5.14 for ammonium. In typical regions, the satellite-derived dataset was significantly robust. Based on the satellite-derived dataset, the spatial-temporal variation of 11-year annual average satellite-derived SO42− and NH4+ concentrations and time series of monthly average concentrations were also investigated. On a national scale, both exhibited a downward trend each year between 2004 and 2014 (SO42−: −0.61%; NH4+: −0.21%), large values were mainly concentrated in the NCP and SCB. For regions captured at a finer resolution, the inter-annual variation trends presented a positive trend over the periods 2004–2007 and 2008–2011, followed by a negative trend over the period 2012–2014, and sulfate concentrations varied appreciably. Moreover, the seasonal distributions of the 11-year satellite-derived dataset over China were presented. The distribution of both sulfate and ammonium concentrations exhibited seasonal characteristics, with the seasonal concentrations ranking as follows: winter > summer > autumn > spring. High concentrations of these species were concentrated in the NCP and SCB, originating from coal-fired power plants and agricultural activities, respectively. Efforts to reduce sulfur dioxide (SO2) emissions have yielded remarkable results since the government has adopted stricter control measures in recent years. Moreover, ammonia emissions should be controlled while reducing the concentration of sulfur, nitrogen and particulate matter. This study provides an assessment of the population’s exposure to certain chemical components.

Description: Oceanic dimethyl sulfide (DMS) is of interest due to its critical influence on atmospheric sulfur compounds in the marine atmosphere and its hypothesized significant role in global climate. High resolution shipboard underway measurements of surface seawater DMS and the partial pressure of carbon dioxide ( p CO 2 ) were conducted in the Atlantic Ocean and Indian Ocean sectors of the Southern Ocean (SO), the southeast Indian Ocean and the northwest Pacific Ocean from February – April 2014 during the 30th Chinese Antarctic Research Expedition. The SO, particularly in the region south of 58 °S, had the highest mean surface seawater DMS concentration of 4.1 ± 8.3 nM (ranged from 0.1 to 73.2 nM) and lowest mean seawater p CO 2 level of 337 ± 50 μatm (ranged from 221 to 411 μatm) over the entire cruise. Significant variations of surface seawater DMS and p CO 2 in the seasonal ice zone (SIZ) of SO were observed, which are mainly controlled by biological process and sea ice activity. We found a significant negative relationship between DMS and p CO 2 in the SO SIZ using 0.1 degree resolution, [DMS] seawater = - 0.160 [ p CO 2 ] seawater + 61.8 (r 2 = 0.594, n = 924, p 〈 0.001). We anticipate that the relationship may possibly be utilized to reconstruct the surface seawater DMS climatology in the SO SIZ. Further studies are necessary to improve the universality of this approach.

Description: Recent observations confirm the rising temperatures of Atlantic waters transported into the Arctic Ocean via the West Spitsbergen Current (WSC). We studied the overall abundance and population structure of the North Atlantic keystone zooplankton copepod Calanus finmarchicus , which is the main prey for pelagic fish and some seabirds, in relation to selected environmental variables in this area between 2001 and 2011, when warming in the Arctic and Subarctic was particularly pronounced. Sampling within a three-week time window each summer demonstrated that trends in the overall abundance of C. finmarchicus varied between years, with the highest values in “extreme” years, due to high numbers of nauplii and early copepodite stages in colder years (2001, 2004, 2010), and contrary to that, the fifth copepodite stage (C5) peaking in warm years (2006, 2007, 2009). The most influential environmental variable, driving C. finmarchicus life cycle was temperature, which promoted an increased C5 abundance when the temperature was above 6°C, indicating earlier spawning and/or accelerated development, and possibly leading to their development to adults later in the summer and spawning for the second time, given adequate food supply. Based on the presented high interannual and spatial variability, we hypothesize that under a warmer climate, C. finmarchicus may annually produce two generations in the southern part of the WSC, what in turn could lead to food web reorganization of important top predators, such as little auks, and induce northward migrations of fish, especially the Norwegian herring. This article is protected by copyright. All rights reserved.

Description: Remote Sensing, Vol. 9, Pages 837: National BDS Augmentation Service System (NBASS) of China: Progress and Assessment Remote Sensing doi: 10.3390/rs9080837 Authors: Chuang Shi Fu Zheng Yidong Lou Shengfeng Gu Weixing Zhang Xiaolei Dai Xianjie Li Hailin Guo Xiaopeng Gong Abstract: In this contribution, the processing strategies of real-time BeiDou System (BDS) precise orbits, clocks, and ionospheric corrections in the National BDS Augmentation Service System (NBASS) are briefly introduced. The Root Mean Square (RMS) of BDS predicted orbits are better than 10 cm in radial and cross-track components, and the accuracy of the BDS real-time clock is better than 0.5 ns for Inclined Geosynchronous Orbit (IGSO) and Mid Earth Orbit (MEO) satellites. The accuracy of BDS Geostationary Earth Orbit (GEO) orbits and clocks are worse than the IGSO and MEO satellites due to its poor geometry conditions. The real-time ionospheric correction is evaluated by cross-validation, and the average accuracy in the vertical direction is about 4 TECU. With these real-time corrections, the overall single and dual-frequency kinematic precise point positioning (PPP) performance in China are evaluated in terms of positioning accuracy at the 95% confidence level and convergence time. The BDS PPP positioning accuracy shows significant regional characteristics due to the geometry distribution of BDS satellites and the accuracy of ionospheric model in different regions. The BDS dual-frequency PPP positioning accuracy in high-latitude and western fringe region is about 0.5 m and 1.0 m in the horizontal and vertical component, respectively, while the horizontal accuracy is better than 0.2 m and the vertical accuracy is better than 0.3 m in the midlands. The convergence time of the BDS PPP is much longer than the GPS PPP and it needs more than 60 min to achieve the accuracy better than 10 cm in both horizontal and vertical directions for dual-frequency PPP. Similar with dual-frequency PPP, the positioning accuracy of the BDS single-frequency PPP in the fringe region is worse than other regions. The positioning in the midlands can achieve 0.5 m in horizontal component and 1.0 m in the vertical component. In addition, when GPS and BDS are combined, the positioning performance of both single-frequency and dual-frequency PPP can be greatly improved.

Description: Remote Sensing, Vol. 9, Pages 833: A Study of Landfast Ice with Sentinel-1 Repeat-Pass Interferometry over the Baltic Sea Remote Sensing doi: 10.3390/rs9080833 Authors: Marjan Marbouti Jaan Praks Oleg Antropov Eero Rinne Matti Leppäranta Mapping of fast ice displacement and investigating sea ice rheological behavior is a major open topic in coastal ice engineering and sea ice modeling. This study presents first results on Sentinel-1 repeat-pass space borne synthetic aperture radar interferometry (InSAR) in the Gulf of Bothnia over the fast ice areas. An InSAR pair acquired in February 2015 with a temporal baseline of 12 days has been studied here in detail. According to our results, the surface of landfast ice in the study area was stable enough to preserve coherence over the 12-day baseline, while previous InSAR studies over the fast ice used much shorter temporal baselines. The advantage of longer temporal baseline is in separating the fast ice from drift ice and detecting long term trends in deformation maps. The interferogram showed displacement of fast ice on the order of 40 cm in the study area. Parts of the displacements were attributed to forces caused by sea level tilt, currents, and thermal expansion, but the main factor of the displacement seemed to be due to compression of the drift ice driven by southwest winds with high speed. Further interferometric phase and the coherence measurements over the fast ice are needed in the future for understanding sea ice mechanism and establishing sustainability of the presented InSAR approach for monitoring dynamics of the landfast ice with Sentinel-1 data.

Description: Bumblebees in Europe have been in steady decline since the 1900s. This decline is expected to continue with climate change as the main driver. However, at the local scale, land use and land cover (LULC) change strongly affects the occurrence of bumblebees. At present, LULC change is rarely included in models of future distributions of species. This study's objective is to compare the roles of dynamic LULC change and climate change on the projected distribution patterns of 48 European bumblebee species for three change scenarios until 2100 at the scales of Europe, and Belgium, Netherlands and Luxembourg (BENELUX). We compared three types of models: (1) only climate covariates, (2) climate and static LULC covariates and (3) climate and dynamic LULC covariates. The climate and LULC change scenarios used in the models include, extreme growth applied strategy (GRAS), business as might be usual (BAMBU) and sustainable European development goals (SEDG). We analysed model performance, range gain/loss and the shift in range limits for all bumblebees. Overall, model performance improved with the introduction of LULC covariates. Dynamic models projected less range loss and gain than climate-only projections, and greater range loss and gain than static models. Overall, there is considerable variation in species responses and effects were most pronounced at the BENELUX scale. The majority of species were predicted to lose considerable range, particularly under the extreme growth scenario (GRAS; overall mean: 64% ± 34). Model simulations project a number of local extinctions and considerable range loss at the BENELUX scale (overall mean: 56% ± 39). Therefore, we recommend species-specific modelling to understand how LULC and climate interact in future modelling. The efficacy of dynamic LULC change should improve with higher thematic and spatial resolution. Nevertheless, current broad scale representations of change in major land use classes impact modelled future distribution patterns. This article is protected by copyright. All rights reserved.

Description: Precipitation and temperature time series suffer from many problems, such as short time period, inadequate spatial coverage, missing data and biases from various causes, which is particularly critical in remote areas such as Northern Canada. The development of alternative datasets for using as proxies for inadequate/missing weather data represents a key research area. In this paper, the performance of six alternative datasets is evaluated for hydrological modeling over 12 watersheds located across Canada and the contiguous United States. The datasets can be classified into three distinct categories: (1) interpolated gridded data, (2) reanalysis data, and (3) climate model outputs. Hydrological simulations were carried out using a lumped conceptual hydrological model calibrated using standard weather data, and compared against results using a calibration specific to each alternative dataset. Prior to the hydrological simulations, the alternative datasets were all evaluated with respect to their ability to reproduce gridded daily precipitation and temperature characteristics over North America. The results show that both the reanalysis data and climate model data adequately represent the spatial pattern of daily precipitation and temperature over North America. The North American Regional Reanalysis (NARR) dataset consistently shows the best performance. With respect to hydrological modeling, the observed discharges are accurately represented by both the gridded and NARR datasets, and more so for the NARR data. The National Centers for Environmental Prediction (NCEP) dataset consistently performs worst as it is unable to even capture the seasonal pattern of observed streamflow for three out of the 12 watersheds. These results indicate that the NARR dataset could be used as a proxy for gauged precipitation and temperature for hydrological modeling over watersheds where observational datasets are deficient. The results also illustrate the ability of climate model data to be used for performing hydrological modeling when driven by reanalysis data at their boundaries, and especially so for high-resolution models.

Description: No abstract is available for this article.

Description: Flat terrain and soils with variable permeability make it difficult to assess the relative importance of surface and subsurface runoff in the Canadian Prairies, especially in riparian areas that are critical for water transmission and solute transport. The main objective of this study is therefore to determine whether patterns of hydrologic state variables, namely near-surface soil moisture (SM), soil electrical conductivity (SEC) and soil temperature (ST), can help infer riparian-to-stream soil water movement in Prairie landscapes. Focus is on the near-level Catfish Creek Watershed (south-eastern Manitoba, Canada) where three riparian sites were monitored: a naturally vegetated grassland site, a headwater forested site, and a highly impacted grassed site adjacent to an engineered drainage dyke and a man-made drainage channel. Data from nine to twelve SM, SEC and ST surveys completed at each site in 2015 using a 75-point grid are matched with riparian water table data, surface water level data from adjacent drainage channels and indicators of antecedent moisture conditions. Pattern characteristics, in the form of descriptive statistics and variogram parameters, are estimated for each state variable and then correlated to indicators of antecedent moisture conditions (AMCs), stream and subsurface water level data to infer soil water movement. Results show that potential evapotranspiration, depth to water table and antecedent precipitation have a significant yet variable impact on SM, SEC and ST patterns. A switching behaviour, between dry and wet conditions, is present in riparian areas characterized by grassland vegetation and well-drained soils. The occurrence of shallow subsurface flow is inferred during the wettest conditions. While riparian SM conditions are useful for predicting streamflow response in adjacent channels, such is not the case for riparian SEC and ST. Further investigations are however necessary to confirm the usefulness of SM spatial patterns for predicting streamflow response in other landscapes across the Canadian Prairies.

Description: Remote Sensing, Vol. 9, Pages 845: Assimilation of Sentinel-1 Derived Sea Surface Winds for Typhoon Forecasting Remote Sensing doi: 10.3390/rs9080845 Authors: Yi Yu Xiaofeng Yang Weimin Zhang Boheng Duan Xiaoqun Cao Hongze Leng High-resolution synthetic aperture radar (SAR) wind observations provide fine structural information for tropical cycles and could be assimilated into numerical weather prediction (NWP) models. However, in the conventional method assimilating the u and v components for SAR wind observations (SAR_uv), the wind direction is not a state vector and its observational error is not considered during the assimilation calculation. In this paper, an improved method for wind observation directly assimilates the SAR wind observations in the form of speed and direction (SAR_sd). This method was implemented to assimilate the sea surface wind retrieved from Sentinel-1 synthetic aperture radar (SAR) in the basic three-dimensional variational system for the Weather Research and Forecasting Model (WRF 3DVAR). Furthermore, a new quality control scheme for wind observations is also presented. Typhoon Lionrock in August 2016 is chosen as a case study to investigate and compare both assimilation methods. The experimental results show that the SAR wind observations can increase the number of the effective observations in the area of a typhoon and have a positive impact on the assimilation analysis. The numerical forecast results for this case show better results for the SAR_sd method than for the SAR_uv method. The SAR_sd method looks very promising for winds assimilation under typhoon conditions, but more cases need to be considered to draw final conclusions.

Description: Remote Sensing, Vol. 9, Pages 840: Topic Modelling for Object-Based Unsupervised Classification of VHR Panchromatic Satellite Images Based on Multiscale Image Segmentation Remote Sensing doi: 10.3390/rs9080840 Authors: Li Shen Linmei Wu Yanshuai Dai Wenfan Qiao Ying Wang Image segmentation is a key prerequisite for object-based classification. However, it is often difficult, or even impossible, to determine a unique optimal segmentation scale due to the fact that various geo-objects, and even an identical geo-object, present at multiple scales in very high resolution (VHR) satellite images. To address this problem, this paper presents a novel unsupervised object-based classification for VHR panchromatic satellite images using multiple segmentations via the latent Dirichlet allocation (LDA) model. Firstly, multiple segmentation maps of the original satellite image are produced by means of a common multiscale segmentation technique. Then, the LDA model is utilized to learn the grayscale histogram distribution for each geo-object and the mixture distribution of geo-objects within each segment. Thirdly, the histogram distribution of each segment is compared with that of each geo-object using the Kullback-Leibler (KL) divergence measure, which is weighted with a constraint specified by the mixture distribution of geo-objects. Each segment is allocated a geo-object category label with the minimum KL divergence. Finally, the final classification map is achieved by integrating the multiple classification results at different scales. Extensive experimental evaluations are designed to compare the performance of our method with those of some state-of-the-art methods for three different types of images. The experimental results over three different types of VHR panchromatic satellite images demonstrate the proposed method is able to achieve scale-adaptive classification results, and improve the ability to differentiate the geo-objects with spectral overlap, such as water and grass, and water and shadow, in terms of both spatial consistency and semantic consistency.

Description: Air temperature is correlated with precipitation oxygen isotope (δ 18 O prcp ) variability for much of the eastern and central United States, but the nature of this δ 18 O prcp -temperature relationship is largely based on data coarsely aggregated at a monthly resolution. We constructed a database of 6177 weeks of isotope and precipitation-day air temperature data from 25 sites to determine how more precise data change our understanding of this classic relationship. Because the δ 18 O prcp -temperature relationship is not perfectly linear, trends in the regression residuals suggest the influence of additional environmental factors such as moisture recycling and extratropical cyclone interactions. Additionally, the temporal relationships between δ 18 O prcp and temperature observed in the weekly data at individual sites can explain broader spatial patterns observed across the study region. For 20 of 25 sites, the δ 18 O prcp -temperature relationship slope is higher for colder precipitation than for warmer precipitation. Accordingly, northern and western sites with relatively more cold precipitation events have steeper overall relationships with higher slope values than southeastern sites that have more warm precipitation events. Although the magnitude of δ 18 O prcp variability increases to the north and west, the fraction of δ 18 O prcp variability explained by temperature increases due to wider annual temperature ranges, producing stronger relationships in these regions. When our δ 18 O prcp -temperature data is grouped by month, we observe significant variations in the relationship from month to month. This argues against a principal causative role for temperature and suggests the existence of an alternative environmental control on δ 18 O prcp values that simply co-varies seasonally with temperature.

Description: Bioretention cells, which are generally effective in controlling surface runoff and recharging groundwater, have been widely adopted as low impact development practices. However, shallow groundwater has limited their implementation in some locations due to the potential problems of a reduction in surface runoff control, groundwater pollution, and continuous groundwater drainage through the underdrain. Many guidelines have established minimum requirements for the groundwater depth below bioretention cells, but they may not be optimized for certain environmental conditions and bioretention cell designs. This study made use of a variably saturated flow model to examine the hydrologic performance of a single bioretention cell in shallow groundwater with event-based simulations, considering a wide range of initial groundwater depths, media and in situ soil types, surface runoff loads, and underdrain sizes. Performance indicators (e.g., runoff reduction, time for infiltrated water to reach the bioretention cell bottom and the groundwater table, height and dissipation time of groundwater mound) were evaluated to examine the processes of runoff generation, the formation and dissipation of groundwater mounds, and the bioretention cell's performance in a shallow groundwater environment. The most influential factors were the initial groundwater depth, the hydraulic conductivity of the media soil, and the rainfall runoff load. With a deeper initial groundwater table, infiltrated water took longer to reach the bioretention cell bottom and groundwater table. Groundwater mounds, however, took longer to dissipate even though they were smaller. The groundwater quality can be better protected if relatively less-permeable soil types (e.g., sandy loam) are used as the media, although it may compromise the performance in runoff quantity control. However, only very high surface runoff loads would cause concerns regarding a reduction in runoff quantity control and possible groundwater contamination due to the shallow groundwater. A distance of 1.5-3 m between the bioretention cell bottom and the groundwater table is generally sufficient. The results of this study could help to guide the planning and design of bioretention cells in areas of shallow groundwater.

Description: Northern peatlands are a vital component of the global carbon cycle, containing large stores of soil organic carbon and acting as a long-term carbon sink. Moss productivity is an important factor in determining whether these wetlands will retain this function under future climatic conditions. Research on unsaturated water flow in peatlands, which controls moss productivity during periods of evaporative stress, has focused on relatively deep bog systems. However, shallower peatlands and marginal connective wetlands can be essential components of many landscape mosaics. In order to better understand factors influencing moss productivity, water balance simulations using Hydrus 1-D were run for different soil profile depths, compositions and antecedent moisture conditions. Our results demonstrate a bimodal distribution of peatland realizations; either primarily conserving water by limiting evapotranspiration or, maximizing moss productivity. For sustained periods of evaporative stress, both deep water storage and a shallow initial water table delay the onset of high vegetative stress, thus maximizing moss productivity. A total depth of sand and peat of 0.8 m is identified as the threshold above which increasing peat depth has no effect on changing vegetative stress response. In contrast, wetlands with shallow peat deposits (less than 0.5 m thick) are least able to buffer prolonged periods of evaporation due to limited labile water storage, and will thus quickly experience vegetative stress and so limit evaporation and conserve water. With a predicted increase in the frequency and size of rain events in continental North America the moss productivity of shallow wetland systems may increase, but also greater moisture availability will increase the likelihood they remain as wetlands in a changing climate.

Description: Shallow aquifers typically have greater hydrologic connectivity and response to recharge and changes in surface water management practices than deeper aquifers, and are therefore often managed to reduce the risk of flooding. Quantification of the water table elevation response under different management scenarios provides valuable information in shallow aquifer systems to assess indirect influences of such modifications. The Episodic Master Recession method was applied to the 15-min water table elevation and NEXRAD rainfall data for six wells to identify water table response and individual rainfall events. The objectives of this study were to evaluate the effects of rainfall, water table elevation, canal stage, site specific characteristics, and canal structure modification/ water management practice on the fluctuations in water table elevations using multiple/stepwise multiple linear regression techniques. With the modification of canal structure and operation adjustment, significant difference existed in water table response in the southern wells due to its relative downstream position regarding the general groundwater flow direction and the structural modification locations. On average, water table response height and flood risk were lower after than before the structure modification to canals. The effect of rainfall event size on the height of water table response was greater than the effect of antecedent water table elevation and canal stage on the height of water table response. Other factors including leakance of the canal bed sediment, specific yield, and rainfall on i-1 day had significant effects on the height of water table response as well. Antecedent water table elevation and canal stage had greater and more linear effects on the height of water table response after the management changes to canals. Variation in water table response height/rainfall event size ratio was attributed to difference in S y , antecedent soil water content, hydraulic gradient, rainfall size and runoff ratio. After the structure modification, water table response height/rainfall event size ratio demonstrated more linear and proportional relationship with antecedent water table elevation and canal stage.

Description: The proper management of coastal aquifers commonly requires an understanding of regional mass flow and complete seawater-freshwater circulation. In this study, time series observations of seawater intrusion and refreshing were conducted using a column experiment based on natural flow conditions in coastal groundwater and a sampled medium from a coastal sandy aquifer without chemical treatment. Ranges of hydrodynamical and hydrochemical variables were tested and analyzed. The results showed that the zeta potential of suspended colloids in aqueous solution in an aquifer polluted with 0.5 g/kg of heavy metals exhibited an isoelectric point for pH values ranging from 5.70 to 6.07 when freshwater or seawater completely occupied the aquifer pores, which is representative of natural hydrochemical conditions. In this scenario, a high background concentration of heavy metals induced colloidal immobilization. Otherwise, seawater-freshwater circulation enabled colloid mobilization due to ionic strength and pH fluctuations. The migration of multiple heavy metals occurred at a characteristic time of approximately 1 pore volume after each intrusion stage began and when the peak rate of colloid release was reached. At these times, the colloid behavior determined the quantity and pathway of heavy metal transport. Based on the influences of seawater and freshwater interactions, the quantity of mobilized particles generally decreased and was uniformly distributed in each fraction due to particle loss and decreased porous connectivity. We speculate that the decrease in the total surface area of the migratory colloids may cause colloid-associated heavy metal transport to decrease. The experimental results provide a useful basis for testing coastal groundwater flow and mass transport models because these phenomena require full characterization to precisely evaluate the associated fluxes from the field scale to the microscopic dimension.

Description: The Natural Resources Conservation Service Curve Number model (NRCS CN) is one of the most recognizable procedures in the field of rainfall-runoff estimation. It has been widely applied for different purposes in hydrological models. In spite of its widespread use, some uncertainties have not even clarified, and must be examined for its proper application. Particularly, choosing the most representative rainfall-runoff events, and the coefficient λ which relates the parameters of the model (curve number CN and initial abstraction I a ). In this research, an advanced analysis is developed to evaluate the influence of λ for a set of representative watersheds of the ARS-USDA (Agricultural Research Service of the United Stated Department of Agriculture). They are characterized by different soil properties, land uses and climatic conditions. Finally, two novel methodologies for the selection of the most representative rainfall-runoff events and for the adaptation of coefficient λ are included, based on the pattern of rainfall distribution.

Description: s Reducing parameter uncertainty of process-based terrestrial ecosystem models (TEMs) is one of the primary targets for accurately estimating carbon budgets and predicting ecosystem responses to climate change. However, parameters in TEMs are rarely constrained by observations from Chinese forest ecosystems, which are important carbon sink over the Northern Hemispheric land. In this study, eddy-covariance data from 6 forest sites in China are used to optimize parameters of the ORCHIDEE (ORganizing Carbon and Hydrology In Dymanics EcosystEms) TEM. The model-data assimilation through parameter optimization largely reduces the prior model errors and improves the simulated seasonal cycle and summer diurnal cycle of net ecosystem exchange (NEE), latent heat (LE) fluxes as well as gross primary production (GPP) and ecosystem respiration. Climate change experiments based on the optimized model are deployed to indicate that forest net primary production (NPP) is suppressed in response to warming in the southern China but stimulated in the northeastern China. Altered precipitation has an asymmetric impact on forest NPP at sites in water-limited regions, with the optimization-induced reduction in response of NPP to precipitation decline being as large as 61% at a deciduous broadleaf forest site. We find that seasonal optimization alters forest carbon cycle responses to environmental change, with the parameter optimization consistently reducing the simulated positive response of heterotrophic respiration to warming. Evaluations from independent observations suggest that improving model structure still matters most for long term carbon stock and its changes, in particular nutrient- and age-related changes of photosynthetic rates, carbon allocation and tree mortality.

Description: Remote Sensing, Vol. 9, Pages 851: Reflectance Intensity Assisted Automatic and Accurate Extrinsic Calibration of 3D LiDAR and Panoramic Camera Using a Printed Chessboard Remote Sensing doi: 10.3390/rs9080851 Authors: Weimin Wang Ken Sakurada Nobuo Kawaguchi This paper presents a novel method for fully automatic and convenient extrinsic calibration of a 3D LiDAR and a panoramic camera with a normally printed chessboard. The proposed method is based on the 3D corner estimation of the chessboard from the sparse point cloud generated by one frame scan of the LiDAR. To estimate the corners, we formulate a full-scale model of the chessboard and fit it to the segmented 3D points of the chessboard. The model is fitted by optimizing the cost function under constraints of correlation between the reflectance intensity of laser and the color of the chessboard’s patterns. Powell’s method is introduced for resolving the discontinuity problem in optimization. The corners of the fitted model are considered as the 3D corners of the chessboard. Once the corners of the chessboard in the 3D point cloud are estimated, the extrinsic calibration of the two sensors is converted to a 3D-2D matching problem. The corresponding 3D-2D points are used to calculate the absolute pose of the two sensors with Unified Perspective-n-Point (UPnP). Further, the calculated parameters are regarded as initial values and are refined using the Levenberg-Marquardt method. The performance of the proposed corner detection method from the 3D point cloud is evaluated using simulations. The results of experiments, conducted on a Velodyne HDL-32e LiDAR and a Ladybug3 camera under the proposed re-projection error metric, qualitatively and quantitatively demonstrate the accuracy and stability of the final extrinsic calibration parameters.

Description: Remote Sensing, Vol. 9, Pages 850: What is the Direction of Land Change? A New Approach to Land-Change Analysis Remote Sensing doi: 10.3390/rs9080850 Authors: Mingde You Anthony Filippi İnci Güneralp Burak Güneralp Accurate characterization of the direction of land change is a neglected aspect of land dynamics. Knowledge on direction of historical land change can be useful information when understanding relative influence of different land-change drivers is of interest. In this study, we present a novel perspective on land-change analysis by focusing on directionality of change. To this end, we employed Maximum Cross-Correlation (MCC) approach to estimate the directional change in land cover in a dynamic river floodplain environment using Landsat 5 Thematic Mapper (TM) images. This approach has previously been used for detecting and measuring fluid and ice motions but not to study directional changes in land cover. We applied the MCC approach on land-cover class membership layers derived from fuzzy remote-sensing image classification. We tested the sensitivity of the resulting displacement vectors to three user-defined parameters—template size, search window size, and a threshold parameter to determine valid (non-noisy) displacement vectors—that directly affect the generation of change, or displacement, vectors; this has not previously been thoroughly investigated in any application domain. The results demonstrate that it is possible to quantitatively measure the rate of directional change in land cover in this floodplain environment using this particular approach. Sensitivity analyses indicate that template size and MCC threshold parameter are more influential on the displacement vectors than search window size. The results vary by land-cover class, suggesting that spatial configuration of land-cover classes should be taken into consideration in the implementation of the method.

Description: Under changing climatic conditions permafrost peatlands can play an important role in the global carbon budget through permafrost carbon feedbacks and shifts in carbon assimilation. To better predict future dynamics in these ecosystems an increased understanding of their Holocene carbon and permafrost history is needed. In Tavvavuoma, northern Sweden, we have performed detailed analyses of vegetation succession and geochemical properties at six permafrost peatland sites. Peatland initiation took place around 10 000 to 9600 cal. a BP, soon after retreat of the Fennoscandian Ice Sheet, and the peatlands have remained permafrost-free fens throughout most of the Holocene. At the four sites that showed a continuous accumulation record during the late Holocene radiocarbon dating of the shift from wet fen to dry bog vegetation, characteristic of the present permafrost peatland surface, suggests that permafrost developed at around 600–100 cal. a BP. At the other two sites peat accumulation was halted during the late Holocene, possibly due to abrasion, making it more difficult to imply the timing of permafrost aggradation. However also at these sites there are no indications of permafrost inception prior to the Little Ice Age. The mean long-term Holocene carbon accumulation rate at all six sites was 12.3±2.4 gC m −2 a −1 (±SD), and the mean soil organic carbon storage was 114±27 kg m −2 .

Description: Earth is experiencing multiple global changes that will, together, determine the fate of many species. Yet, how biological communities respond to concurrent stressors at local-to-regional scales remains largely unknown. In particular, understanding how local habitat conversion interacts with regional climate change to shape patterns in β-diversity—differences among sites in their species compositions—is critical to forecast communities in the Anthropocene. Here, we study patterns in bird β-diversity across land-use and precipitation gradients in Costa Rica. We mapped forest cover, modeled regional precipitation, and collected data on bird community composition, vegetation structure, and tree diversity across 120 sites on 20 farms to answer three questions. First, do bird communities respond more strongly to changes in land use or climate in northwest Costa Rica? Second, does habitat conversion eliminate β-diversity across climate gradients? Third, does regional climate control how communities respond to habitat conversion and, if so, how? After correcting for imperfect detection, we found that local land-use determined community shifts along the climate gradient. In forests, bird communities were distinct between sites that differed in vegetation structure or precipitation. In agriculture, however, vegetation structure was more uniform, contributing to 7%–11% less bird turnover than in forests. In addition, bird responses to agriculture and climate were linked: agricultural communities across the precipitation gradient shared more species with dry than wet forest communities. These findings suggest that habitat conversion and anticipated climate drying will act together to exacerbate biotic homogenization. Species’ abundance responses to land-use and precipitation gradients were highly correlated (a). That is, species that thrived in forests tended to be more abundant in wet regions (e.g., the pictured Barred Antshrike, Thamnophilus doliatus ), while species that benefited from agriculture tended to reach maximal abundances in drier regions (e.g., Red-winged Blackbird, Agelaius phoeniceus ). Thus, the average agricultural bird community was more dissimilar to communities found in wet forests than communities found in dry forests (b).

Description: Immiscible fluid-fluid displacement in permeable media is important in many subsurface processes, including enhanced oil recovery and geological CO 2 sequestration. Controlled by capillary and viscous forces, displacement patterns of one fluid displacing another more viscous one exhibit capillary and viscous fingering, and crossover between the two. Although extensive studies investigated viscous and capillary fingering in porous media, a few studies focused on the crossover in rough fractures, and how viscous and capillary forces affect the crossover remains unclear. Using a transparent fracture-visualization system, we studied how the two forces impact the crossover in a horizontal rough fracture. Drainage experiments of water displacing oil were conducted at seven flow rates (capillary number log 10 Ca ranging from −7.07 to −3.07) and four viscosity ratios ( M =1/1000,1/500,1/100 and 1/50). We consistently observed lower invading fluid saturations in the crossover zone. We also proposed a phase diagram for the displacement patterns in a rough fracture that is consistent with similar studies in porous media. Based on real-time imaging and statistical analysis of the invasion morphology, we showed that the competition between capillary and viscous forces is responsible for the saturation reduction in the crossover zone. In this zone, finger propagation toward the outlet (characteristic of viscous fingering) as well as void-filling in the transverse/backward directions (characteristic of capillary fingering), are both suppressed. Therefore, the invading fluid tends to occupy larger apertures with higher characteristic front velocity, promoting void-filling toward the outlet with thinner finger growth and resulting in a larger volume of defending fluid left behind.

Description: On the cover: The cover image, by Yulin Lin et al., is based on the Research Article Assessment of fracture risk in proximal tibia with tumorous bone defects by a finite element method , DOI: 10.1002/jemt.22899 .

Description: ABSTRACT Few Australian wetlands have persisted since the Last Glacial Maximum, with fewer still in existence through the entire last glacial cycle. The absence of wetlands, which itself indicates periods of moisture deficit, means there are few continuous climate and environmental change records covering this critical period. The lack of wetland persistence also raises the question of how plant and animal species that require permanent wetlands survived the last glacial cycle. Sixteen wetlands have been cored and dated on North Stradbroke Island (Minjerribah) – a large east Australian sand island – with basal dates reported from 10 sites for the first time. These wetlands range in age to over 200 000 years old, with six dating to the Last Glacial Maximum or earlier. There is no evidence of a stratigraphic discontinuity in the radiocarbon-based age–depth profiles, suggesting continuous deposition of highly organic sediment through the period covered by these ages ( c . 40 ka). The persistence of these wetlands suggests that for much of the last 40 000 years, and for perhaps much longer, the regional moisture balance has been positive. Over the last glacial cycle, persistently wet conditions may have provided a refuge from regional drying, and thus contributed to the high genetic diversity of rainforest plants in the region. Vegetation and climate records from these sites will allow hypothesis testing about the drivers of both local and regional climate variability. Copyright © 2017 John Wiley & Sons, Ltd.

Description: Honey is a valuable food produced by bees from sugary substances that they gather in nature. The transformation the nectar into honey, by bees, is long and complex. Except for honey, where heavy metals are absent or are found only in traces, the bees and their products have always been considered excellent biomarkers of such contaminants. We have assumed that the absence of heavy metals in honey is due to the presence of a detoxification system in the digestive system of bees, which involves metallothioneins, proteins that have a role in the homeostatic control of essential and non-essential metals. We have placed the beehives in three different zones: industrial, urban and rural. Investigations were carried out with ICP-MS method for the detection of heavy metals in the guts of honey bees and honey. The metallothioneins have been identified by Immunohistochemical and Western-blotting analisys. The investigations have shown the presence of heavy metals only in bees guts but not in honey, while the presence of metallothionein has been highlighted only in epithelium of the honey sac, demonstrating the existence of an efficient system of detoxification of heavy metals. In honey bees the detoxification system is primarily based on induction of Metallothioneins Metallothionein-immunopositive cells only in the honey sac.