Methods for the visualization of digital elevation models for landform mapping

Digital elevation models (DEMs) are increasingly used for landform mapping, particularly with the growing availability of national and global datasets. In this paper we describe a variety of techniques that can visualize a DEM. We then compare five techniques to ascertain which performs the most complete and unbiased visualization. We assess the visualization techniques by comparing landforms mapped from them against a detailed morphological map (derived from mapping of multi‐azimuth relief‐shaded DEMs cross‐checked with stereo aerial photographs). Results show that no single visualization method provides complete and unbiased mapping. The relief‐shaded visualizations are particularly prone to azimuth biasing, although they can highlight subtle landforms. We recommend curvature visualization for initial mapping as this provides a non‐illuminated (and therefore unbiased) image. Initial mapping can then be supplemented with data from relief‐shaded visualizations. Copyright © 2005 John Wiley & Sons, Ltd.     

Dynamic modelling of passive margin salt tectonics: effects of water loading, sediment properties and sedimentation patterns

We investigate the evolution of passive continental margin sedimentary basins that contain salt through two‐dimensional (2D) analytical failure analysis and plane‐strain finite‐element modelling. We expand an earlier analytical failure analysis of a sedimentary basin/salt system at a passive continental margin to include the effects of submarine water loading and pore fluid pressure. Seaward thinning sediments above a weak salt layer produce a pressure gradient that induces Poiseuille flow in the viscous salt. We determine the circumstances under which failure at the head and toe of the frictional–plastic sediment wedge occurs, resulting in translation of the wedge, landward extension and seaward contraction, accompanied by Couette flow in the underlying salt. The effects of water: (i) increase solid and fluid pressures in the sediments; (ii) reduce the head to toe differential pressure in the salt and (iii) act as a buttress to oppose failure and translation of the sediment wedge. The magnitude of the translation velocity upon failure is reduced by the effects of water. The subsequent deformation is investigated using a 2D finite‐element model that includes the effects of the submarine setting and hydrostatic pore pressures. The model quantitatively simulates a 2D approximation of the evolution of natural sedimentary basins on continental margins that are formed above salt. Sediment progradation above a viscous salt layer results in formation of landward extensional basins and listric normal growth faults as well as seaward contraction. At a later stage, an allochthonous salt nappe overthrusts the autochthonous limit of the salt. The nature and distribution of major structures depends on the sediment properties and the sedimentation pattern. Strain weakening of sediment favours landward listric growth faults with formation of asymmetric extensional depocentres. Episodes of low sediment influx, with partial infill of depocentres, produce local pressure gradients in the salt that result in diapirism. Diapirs grow passively during sediment aggradation.     

Multi‐year water balance assessment of a newly constructed wetland,Fort McMurray,Alberta

Oil sands mining in Alberta transforms the boreal landscape of forests and wetlands into open pits, tailings ponds and overburden piles. Whereas reclamation efforts have primarily focused on upland forests, rebuilding wetland systems has recently become a motivation for research. Wetland creation and sustainability in this region is complicated by the sub‐humid climate and salinity of underlying mining material. In 2012, Syncrude Canada Ltd. completed the construction of the Sandhill Fen Watershed (SFW), a 52‐ha upland‐wetland system to evaluate wetland reclamation strategies on soft tailings. SFW includes an active pumping system, upland hummocks, a fen wetland and underdrains. To evaluate the influence of management practices on the hydrology of the system, this study reports the water balance from January 2013 to December 2014, the first 2 years after commissioning. A semi‐distributed approach was taken to examine the fluxes and stores of water in uplands and lowlands. Natural and artificial inputs and outputs were measured using a series of precipitation gauges and pumps, and evapotranspiration was quantified using three eddy covariance towers. A series of near surface wells recorded water table position. Both 2013 and 2014 were normal rainfall years, with 2013 having more and 2014 less snow than normal. In 2013, inflow/outflow from pumping was the predominant hydrological fluxes, resulting in considerable variability in water table position and storage changes throughout the summer. In 2014, the artificial addition of water was negligible, yet the water table remained near the surface in lowland locations, suggesting that wetland conditions could be maintained under current conditions. Evapotranspiration rates between uplands and lowlands were similar between years and sites, ranging from 2.2 ± 1.8 to 2.5 ± 1.2 mm/day and were largely controlled by climate. These rates were less than nearby older upland systems, suggesting that water balance partitioning will change as vegetation develops. Comparison between years and with natural systems provides insight on how management practices influence hydrologic dynamics and the overall water balance of the SFW. Copyright © 2016 John Wiley & Sons, Ltd.     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

Numerical simulation of agricultural sediment and pesticide runoff: RZWQM and PRZM comparison

Agricultural sediment and pesticide runoff is a widespread ecological and human health concern. Numerical simulation models, such as Root Zone Water Quality Model (RZWQM) and Pesticide Root Zone Model (PRZM), have been increasingly used to quantify off‐site agricultural pollutant movement. However, RZWQM has been criticized for its inability to simulate sedimentation processes. The recent incorporation of the sedimentation module of Groundwater Loading Effects of Agricultural Management Systems has enabled RZWQM to simulate sediment and sediment‐associated pesticides. This study compares the sediment and pesticide transport simulation performance of the newly released RZWQM and PRZM using runoff data from 2 alfalfa fields in Davis, California. A composite metric (based on coefficient of determination, Nash–Sutcliffe efficiency, index of agreement, and percent bias) was developed and employed to ensure robust, comprehensive assessment of model performance. Results showed that surface water runoff was predicted reasonably well (absolute percent bias <31%) by RZWQM and PRZM after adjusting important hydrologic parameters. Even after calibration, underestimation bias (?89% ≤ PBIAS  ≤ ?36%) for sediment yield was observed in both models. This might be attributed to PRZM's incorrect distribution of input water and uncertainty in RZWQM's runoff erosivity coefficient. Moreover, the underestimation of sediment might be less if the origin of measured sediment was considered. Chlorpyrifos losses were simulated with reasonable accuracy especially for Field A (absolute PBIAS  ≤ 22%), whereas diuron losses were underestimated to a great extent (?98% ≤ PBIAS  ≤ ?65%) in both models. This could be attributed to the underprediction of herbicide concentration in the top soil due to the limitations of the instantaneous equilibrium sorption model as well as the high runoff potential of herbicide formulated as water‐dispersible granules. RZWQM and PRZM partitioned pesticides into the water and sediment phases similarly. According to model predictions, the majority of pesticide loads were carried via the water phase. On the basis of this study, both RZWQM and PRZM performed well in predicting runoff that carried highly adsorptive pesticides on an event basis, although the more physically based RZWQM is recommended when field‐measured soil hydraulic properties are available.     

Diagnosing snow accumulation errors in a rain‐snow transitional environment with snow board observations

Diagnosing the source of errors in snow models requires intensive observations, a flexible model framework to test competing hypotheses, and a methodology to systematically test the dominant snow processes. We present a novel process‐based approach to diagnose model errors through an example that focuses on snow accumulation processes (precipitation partitioning, new snow density, and snow compaction). Twelve years of meteorological and snow board measurements were used to identify the main source of model error on each snow accumulation day. Results show that modeled values of new snow density were outside observational uncertainties in 52% of days available for evaluation, while precipitation partitioning and compaction were in error 45% and 16% of the time, respectively. Precipitation partitioning errors mattered more for total winter accumulation during the anomalously warm winter of 2014–2015, when a higher fraction of precipitation fell within the temperature range where partition methods had the largest error. These results demonstrate how isolating individual model processes can identify the primary source(s) of model error, which helps prioritize future research.