Guidelines on the use of structure-from-motion photogrammetry in geomorphic research

As a topographic modelling technique, structure-from-motion (SfM) photogrammetry combines the utility of digital photogrammetry with a flexibility and ease of use derived from multi-view computer vision methods. In conjunction with the rapidly increasing availability of imagery, particularly from unmanned aerial vehicles, SfM photogrammetry represents a powerful tool for geomorphological research. However, to fully realize this potential, its application must be carefully underpinned by photogrammetric considerations, surveys should be reported in sufficient detail to be repeatable (if practical) and results appropriately assessed to understand fully the potential errors involved. To deliver these goals, robust survey and reporting must be supported through (i) using appropriate survey design, (ii) applying suitable statistics to identify systematic error (bias) and to estimate precision within results, and (iii) propagating uncertainty estimates into the final data products. © 2019 John Wiley & Sons, Ltd.     

Mitigating systematic error in topographic models derived from UAV and ground‐based image networks

High resolution digital elevation models (DEMs) are increasingly produced from photographs acquired with consumer cameras, both from the ground and from unmanned aerial vehicles (UAVs). However, although such DEMs may achieve centimetric detail, they can also display systematic broad‐scale error that restricts their wider use. Such errors which, in typical UAV data are expressed as a vertical ‘doming’ of the surface, result from a combination of near‐parallel imaging directions and inaccurate correction of radial lens distortion. Using simulations of multi‐image networks with near‐parallel viewing directions, we show that enabling camera self‐calibration as part of the bundle adjustment process inherently leads to erroneous radial distortion estimates and associated DEM error. This effect is relevant whether a traditional photogrammetric or newer structure‐from‐motion (SfM) approach is used, but errors are expected to be more pronounced in SfM‐based DEMs, for which use of control and check point measurements are typically more limited. Systematic DEM error can be significantly reduced by the additional capture and inclusion of oblique images in the image network; we provide practical flight plan solutions for fixed wing or rotor‐based UAVs that, in the absence of control points, can reduce DEM error by up to two orders of magnitude. The magnitude of doming error shows a linear relationship with radial distortion and we show how characterization of this relationship allows an improved distortion estimate and, hence, existing datasets to be optimally reprocessed. Although focussed on UAV surveying, our results are also relevant to ground‐based image capture. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Towards a consistent geochemical model for prediction of uranium(VI) removal from groundwater by ferrihydrite

Uranium(VI), which is often elevated in granitoidic groundwaters, is known to adsorb strongly to Fe (hydr)oxides under certain conditions. This process can be used in water treatment to remove U(VI). To develop a consistent geochemical model for U(VI) adsorption to ferrihydrite, batch experiments were performed and previous data sets reviewed to optimize a set of surface complexation constants using the 3-plane CD-MUSIC model. To consider the effect of dissolved organic matter (DOM) on U(VI) speciation, new parameters for the Stockholm Humic Model (SHM) were optimized using previously published data. The model, which was constrained from available X-ray absorption fine structure (EXAFS) spectroscopy evidence, fitted the data well when the surface sites were divided into low- and high-affinity binding sites. Application of the model concept to other published data sets revealed differences in the reactivity of different ferrihydrites towards U(VI). Use of the optimized SHM parameters for U(VI)-DOM complexation showed that this process is important for U(VI) speciation at low pH. However in neutral to alkaline waters with substantial carbonate present, Ca–U–CO₃ complexes predominate. The calibrated geochemical model was used to simulate U(VI) adsorption to ferrihydrite for a hypothetical groundwater in the presence of several competitive ions. The results showed that U(VI) adsorption was strong between pH 5 and 8. Also near the calcite saturation limit, where U(VI) adsorption was weakest according to the model, the adsorption percentage was predicted to be >80%. Hence U(VI) adsorption to ferrihydrite-containing sorbents may be used as a method to bring down U(VI) concentrations to acceptable levels in groundwater.     

What is responsible for increasing flood risks? The case of Gangwon Province,Korea

We examined the anthropogenic and natural causes of flood risks in six representative cities in the Gangwon Province of Korea. Flood damage per capita is mostly explained by cumulative upper 5% summer precipitation amount and the year. The increasing flood damage is also associated with deforestation in upstream areas and intensive land use in lowlands. Human encroachment on floodplains made these urban communities more vulnerable to floods. Without changes in the current flood management systems of these cities, their vulnerability to flood risks will remain and may even increase under changing climate conditions.     

Automated forecasting of volcanic ash dispersion utilizing Virtual Globes

There are over 100 active volcanoes in the North Pacific (NOPAC) region, most of which are located in sparsely populated areas. Dispersion models play an important role in forecasting the movement of volcanic ash clouds by complementing both remote sensing data and visual observations from the ground and aircraft. Puff is a three-dimensional dispersion model, primarily designed for forecasting volcanic ash dispersion, used by the Alaska Volcano Observatory and other agencies. Since early 2007, the model is in an automated mode to predict the movement of airborne volcanic ash at multiple elevated alert status volcanoes worldwide to provide immediate information when an eruption occurs. Twelve of the predictions are within the NOPAC region, nine more within the southern section of the Pacific ring of fire and the others are in Europe and the Caribbean. Model forecasts are made for initial ash plumes ranging from 4 to 20 km altitude above sea level and for a 24-h forecast period. This information is made available via the Puff model website. Model results can be displayed in Virtual Globes for three-dimensional visualization. Here, we show operational Puff predictions in two and three-dimensions in Google Earth^®, both as iso-surfaces and particles, and study past eruptions to illustrate the capabilities that the Virtual Globes can provide. In addition, we show the opportunity that Google Maps^® provides in displaying Puff operational predictions via an application programming web interface and how radiosonde data (vertical soundings) and numerical weather prediction vertical profiles can be displayed in Virtual Globes for assisting in estimating ash cloud heights.     

Identifying subsoil sediment sources with carbon and nitrogen stable isotope ratios

Increased sediment loads from accelerated catchment erosion significantly degrade waterways worldwide. In the South East Queensland region of Australia, sediment loads are degrading Moreton Bay, a Ramsar listed wetland of international significance. In this region, like most parts of coastal Australia, sediment is predominantly derived from gully and channel bank erosion processes. A novel approach is presented that uses carbon and nitrogen stable isotope ratios and elemental composition to discriminate between these often indistinguishable subsoil sediment sources. The conservativeness of these sediment properties is first tested by examining the effect of particle size separation (testing for consistency during transport) and the effect of sampling at different times (testing for temporal source consistency). The discrimination potential of these sediment properties is then assessed with the conservative properties, based on the particle size and temporal analyses, modelled to determine sediment provenance in three catchments. Nitrogen sediment properties were found to have significant particle size enrichment and high temporal variance indicative of non‐conservative behaviour. Conversely, carbon stable isotopes had very limited particle size and temporal variability highlighting their suitability for sediment tracing. Channel erosion was modelled to be a significant source of sediment (μ 51%, σ 9%) contrasting desktop modelling research that estimated gully erosion is the predominant sediment source. To limit the supply of sediment to Moreton Bay, channel bank and gully erosion must both be targeted by sediment management programs. By distinguishing between subsoil sediment sources, this approach has the potential to enhance the management of sediment loads degrading waterways worldwide. Copyright © 2014 John Wiley & Sons, Ltd.     

Petrological and experimental evidence for differentiation of water-rich magmas beneath St. Kitts,Lesser Antilles

St. Kitts lies in the northern Lesser Antilles, a subduction-related intraoceanic volcanic arc known for its magmatic diversity and unusually abundant cognate xenoliths. We combine the geochemistry of xenoliths, melt inclusions and lavas with high pressure–temperature experiments to explore magma differentiation processes beneath St. Kitts. Lavas range from basalt to rhyolite, with predominant andesites and basaltic andesites. Xenoliths, dominated by calcic plagioclase and amphibole, typically in reaction relationship with pyroxenes and olivine, can be divided into plutonic and cumulate varieties based on mineral textures and compositions. Cumulate varieties, formed primarily by the accumulation of liquidus phases, comprise ensembles that represent instantaneous solid compositions from one or more magma batches; plutonic varieties have mineralogy and textures consistent with protracted solidification of magmatic mush. Mineral chemistry in lavas and xenoliths is subtly different. For example, plagioclase with unusually high anorthite content (An_≤100) occurs in some plutonic xenoliths, whereas the most calcic plagioclase in cumulate xenoliths and lavas are An₉₇ and An₉₅, respectively. Fluid-saturated, equilibrium crystallisation experiments were performed on a St. Kitts basaltic andesite, with three different fluid compositions (XH₂O = 1.0, 0.66 and 0.33) at 2.4 kbar, 950–1025 °C, and fO₂ = NNO ? 0.6 to NNO + 1.2 log units. Experiments reproduce lava liquid lines of descent and many xenolith assemblages, but fail to match xenolith and lava phenocryst mineral compositions, notably the very An-rich plagioclase. The strong positive correlation between experimentally determined plagioclase-melt Kd_Ca–Na and dissolved H₂O in the melt, together with the occurrence of Al-rich mafic lavas, suggests that parental magmas were water-rich (> 9 wt% H₂O) basaltic andesites that crystallised over a wide pressure range (1.5–6 kbar). Comparison of experimental and natural (lava, xenolith) mafic mineral composition reveals that whereas olivine in lavas is predominantly primocrysts precipitated at low-pressure, pyroxenes and spinel are predominantly xenocrysts formed by disaggregation of plutonic mushes. Overall, St. Kitts xenoliths and lavas testify to mid-crustal differentiation of low-MgO basalt and basaltic andesite magmas within a trans-crustal, magmatic mush system. Lower crustal ultramafic cumulates that relate parental low-MgO basalts to primary, mantle -derived melts are absent on St. Kitts.     

The surface boundary condition in nonhydrostatic ocean models

With increasing resolution in numerical ocean models, nonhydrostatic pressure effects have to be accounted for. In sigma-coordinate mode split ocean models, this pressure may be regarded as a pressure correction. An elliptic equation must be solved for the nonhydrostatic pressure, and the gradients are used to correct the provisional hydrostatic velocity components in each time step. The focus in the present work is on the surface boundary condition for the elliptic equation. In the literature, both Dirichlet and Neumann boundary conditions are suggested and applied. To investigate the sensitivity of the numerical results to the choice of boundary condition, three numerical experiments are performed. The first and second experiments are studies of the propagation and steepening of nonlinear internal waves. The first study is on tank scale and the second experiment is on ocean scale. In the tank-scale experiment, the density and the flow fields are very robust to the choice of boundary condition. In the ocean-scale experiment, the waves produced with a Dirichlet boundary condition become more damped than the waves produced with a Neumann boundary condition. The third study involves a surface buoyant jet. It is shown that well-known characteristics of the plume front are reproduced with a Neumann boundary condition, but the rotating turbulent core of this front is lost with a Dirichlet condition. It is accordingly argued that the appropriate surface boundary condition in mode split nonhydrostatic ocean models is the Neumann condition.     

Sediment fingerprinting in fluvial systems： review of tracers,sediment sources and mixing models

Suspended sediments in fluvial systems originate from a myriad of diffuse and point sources, with the relative contribution from each source varying over time and space. The process of sediment fingerprinting focuses on developing methods that enable discrete sediment sources to be identified from a composite sample of suspended material. This review identifies existing methodological steps for sediment fingerprinting including fluvial and source sampling, and critically compares biogeochemical and physical tracers used in fingerprinting studies. Implications of applying different mixing models to the same source data are explored using data from 41 catchments across Europe, Africa, Australia, Asia, and North and South America. The application of seven commonly used mixing models to two case studies from the US （North Fork Broad River watershed） and France （Bldone watershed） with local and global （genetic algorithm） optimization methods identified all outputs remained in the acceptable range of error defined by the original authors. We propose future sediment fingerprinting studies use models that combine the best explanatory parameters provided by the modified Collins （using correction factors） and Hughes （relying on iterations involving all data, and not only their mean values） models with optimization using genetic algorithms to best predict the relative contribution of sediment sources to fluvial systems.     

Quantifying urban land cover change between 2001 and 2006 in the Gulf of Mexico region

We estimated urbanization rates (2001–2006) in the Gulf of Mexico region using the National Land Cover Database (NLCD) 2001 and 2006 impervious surface products. An improved method was used to update the NLCD impervious surface product in 2006 and associated land cover transition between 2001 and 2006. Our estimation reveals that impervious surface increased 416 km² with a growth rate of 5.8% between 2001 and 2006. Approximately 1110.1 km² of non-urban lands were converted into urban land, resulting in a 3.2% increase in the region. Hay/pasture, woody wetland, and evergreen forest represented the three most common land cover classes that transitioned to urban. Among these land cover transitions, more than 50% of the urbanization occurred within 50 km of the coast. Our analysis shows that the close-to-coast land cover transition trend, especially within 10 km off the coast, potentially imposes substantial long-term impacts on regional landscape and ecological conditions.