The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Assessing uncertainty in best management practice effectiveness under future climate scenarios

Uncertainty of best management practice (BMP) performance in future climates is an important consideration for water resources managers. The objective of this study was to quantify the level of uncertainty in performance of seven agricultural BMPs due to climate change in reducing sediment, total nitrogen, and total phosphorus loads. The Soil and Water Assessment Tool coupled with mid‐21^st century climate data from the Community Climate System Model were used to develop climate change scenarios for the Tuttle Creek Lake Watershed of Kansas and Nebraska. Uncertainty level of each BMP was determined using Latin Hypercube Sampling, a constrained Monte Carlo sampling technique. Samples were taken from distributions of several variables (monthly precipitation, temperature, CO₂, and BMP implementation parameters). Cumulative distribution functions were constructed for each BMP, pollutant, and climate scenario combination. Results demonstrated that BMP performance uncertainty is amplified in the extreme climate scenario. Among BMPs, native grass replacement generally had higher uncertainty level but also had the greatest reductions. This study highlights the importance of incorporating uncertainty analysis into mitigation strategies aiming to reduce negative impacts of climate change on water resources. Copyright © 2013 John Wiley & Sons, Ltd.     

Design and commissioning of a laminar soil container for use on small shaking tables

This paper describes the design, fabrication and commissioning of a single axis laminar shear box for use in seismic soil–structure interaction studies. A laminar shear box is a flexible container that can be placed on a shaking table to simulate vertical shear-wave propagation during earthquakes through a soil layer of finite thickness. The laminar shear box described in this paper was designed to overcome the base shear limitations of a small shaking table at The University of Western Ontario. The design details of the box are provided in addition to results of dynamic tests performed to commission the box. A synthetic clay comprising sodium bentonite mixed with diluted glycerin was used as the model soil and 1-G similitude theory was employed to maintain model to prototype similarity. The model soil was compacted into the container in lifts to achieve soil stiffness that increased with depth. A series of shaking table tests and numerical analyses that were performed to study the performance of the laminar box and non-linear seismic behavior of the model clay are described. The results of this study show that the laminar box does not impose significant boundary effects and is able to maintain 1-D soil column behavior. In addition, the dynamic behavior of the model clay during scaled model tests was found to be consistent with the behavior measured during cyclic laboratory tests.     

Complexity in astroplasmas

Ideas of complexity theory – including self-organized criticality – are being applied to the study of turbulence and astroplasmas in the solar system. An RAS MIST meeting in London brought together researchers from different disciplines to compare established paradigms with more recent ideas. Sandra Chapman , Mervyn Freeman and Sean Oughton report.     

Distorted Froude‐scaled flume analysis of large woody debris

This paper presents the results of a movable‐boundary, distorted, Froude‐scaled hydraulic model based on Abiaca Creek, a sand‐bedded channel in northern Mississippi. The model was used to examine the geomorphic and hydraulic impact of simplified large woody debris (LWD) elements. The theory of physical scale models is discussed and the method used to construct the LWD test channel is developed. The channel model had bed and banks moulded from 0·8 mm sand, and flow conditions were just below the threshold of motion so that any sediment transport and channel adjustment were the result of the debris element. Dimensions and positions of LWD elements were determined using a debris jam classification model. Elements were attached to a dynamometer to measure element drag forces, and channel adjustment was determined through detailed topographic surveys. The fluid drag force on the elements decreased asymptotically over time as the channel boundary eroded around the elements due to locally increased boundary shear stress. Total time for geomorphic adjustment computed for the prototype channel at the Q₂ discharge (discharge occurring once every two years on average) was as short as 45 hours. The size, depth and position of scour holes, bank erosion and bars created by flow acceleration past the elements were found to be related to element length and position within the channel cross‐section. Morphologies created by each debris element in the model channel were comparable with similar jams observed in the prototype channel. Published in 2001 John Wiley & Sons, Ltd.     

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.     

Graphical methods of river profile analysis to unravel drainage area change,uplift and erodibility contrasts in the Central Range of Taiwan

To reveal river channel steepness patterns and variance in settings with significant variation in rock uplift rate, rock erodibility and moving water divides, we present a series of graphical methods to interpret channel profiles. To complement Cartesian χ plots, longitudinal profiles and mapping methods, we introduce a new method based on a radial coordinate system. We map each basin onto polar coordinates in which the radial dimension is χ and the azimuthal coordinate, ?, is calculated with an increment (Δ?) scaled to the distance to neighboring channel heads. The elevation is contoured on this mapping. Average channel steepness is estimated by fitting a conical surface to the elevation. The graph simplifies the comparison of χ and elevation between channels that share a divide, and helps identify spatial changes in drainage area and patterns of erodibility. We apply this approach to derive the uplift pattern in the eastern and southern Central Range of Taiwan, where the high tectonic convergence and uplift rates combined with sub‐tropical climate and frequent typhoons results in high exhumation rate, and well‐developed, detachment‐limited river networks. Additionally, the tectonic activity leads to drainage basin reorganization. We identify examples of divide migration, discrete river capture as well as anomalous steepness that we attribute to local variability in rock erodibility. Estimated basin‐average steepness values show the highest and a near constant value from Hsinwulu basin to Liwu basin in the center of the Island. To the north and south of this region, the values gradually decrease. These estimates show good correlation with the topography of the Central Range and erosion rates derived from in situ ¹⁰Be concentrations in river‐borne quartz. We conclude that the basin steepness reflects systematic differences in rock uplift rate and erosion rate. Copyright © 2016 John Wiley & Sons, Ltd.     

MODELING HEADCUT DEVELOPMENT AND MIGRATION IN UPLAND CONCENTRATED FLOWS

On hillslopes and agricultural fields, discrete areas of intense, localized soil erosion commonly take place in the form of migrating headcuts. These erosional features significantly increase soil loss and landscape degradation, yet the unsteady, transient, and migratory habits of headcuts complicate their phenomenological and erosional characterization. Here a unique experimental facility was constructed to examine actively migrating headcuts typical of upland concentrated flows. Essential components of the facility include a deep soil cavity with external drainage, rainfall simulator, capacity for overland flow, and a video recording technique for data collection. Results from these experiments show that： （1） after a short period of adjustment, headcut migration attained a steady-state condition, where the rate of migration, scour hole geometry, and sediment discharge remain constant with time; （2） boundary conditions of higher rates of overland flow, steeper bed slopes, and larger initial headcut heights produced systematically larger scour holes with higher rates of soil erosion; and （3） during migration, the turbulent flow structure within the scour hole remained unchanged, consisting of an overfall nappe at the brink transitioning into a reattached wall jet with two recirculation eddies within the plunge pool. The systematic behavior of headcut development and migration enabled the application of modified jet impingement theory to predict with good success the characteristics of the impinging jet, the depth of maximum scour, the rate of headcut migration, and the rate of sediment erosion. These laboratory data and the analytical formulation can be used in conjunction with soil erosion prediction technology to improve the management of agricultural areas impacted by headcut development and ephemeral gully erosion.     

The impact of storm-triggered landslides on sediment dynamics and catchment-wide denudation rates in the southern Central Range of Taiwan following the extreme rainfall event of Typhoon Morakot

Extreme erosion events can produce large short-term sediment fluxes. Such events complicate erosion rates estimated from cosmogenic nuclide concentrations in river sediment by providing sediment with a concentration different from the long-term basin average. We present a detrital ¹⁰Be study in southern Taiwan, with multiple samples obtained in a time sequence bracketing the 2009 Typhoon Morakot, to assess the impact of landslide sediment on ¹⁰Be concentrations (N¹⁰Be) in river sediment. Sediment samples were collected from 13 major basins, two or three times over the last decade, to observe the temporal variation of N¹⁰Be. Landslide inventories with time intervals of 5–6 years were used to quantify sediment flux changes. A negative correlation between N¹⁰Be and landslide areal density indicates dilution of N¹⁰Be by landslide sediment. Denudation rates estimated from the diluted N¹⁰Be can be up to three times higher than the lowest rate derived from the same basins. Observed increases imply that, 3 years after the passage of Typhoon Morakot, fluvial channels still contain a considerable amount of sediment provided by hillslope landslides during the event. However, higher N¹⁰Be in 2016 samples indicate that the contribution from landslide sediment at the sampled grain size has decreased with time. The correlation between changes in N¹⁰Be and landslide area and volume is not strong, likely resulting from the stochastic and complex nature of sediment transport. To simultaneously evaluate the volume of landslide-derived sediment and estimate the background denudation rate, associated with less impulsive sediment supply, we constructed a sediment-mixing model with the time series of N¹⁰Be and landslide inventories. The spatial pattern of background erosion rate in southern Taiwan is consistent with the regional tectonic framework, indicating that the landscape is evolving mainly in response to the tectonic forcing, and this signal is modified, but not obscured by impulsive sediment supply. © 2019 John Wiley & Sons, Ltd.     

Incipient weathering rind development on introduced machine‐polished granite discs in an Arctic alpine environment,northern Scandinavia

Weathering rinds, zones of alteration on the exterior surfaces of rock outcrops and coarse unconsolidated surficial debris are widely used by geomorphologists and Quaternary geologists as indicators of the relative age of landforms and landscapes. Additionally they provide unique insights into the earliest stages of rock and mineral weathering, yet the origin of these alteration zones is relatively poorly understood. This lack of understanding applies especially to the initial stages of rind formation. The study reported in this paper has two principal objectives. The first is to use lightly polished granite discs inserted in soil profiles under several different plant communities in an Arctic alpine environment for a period of four or five years to investigate the nature of incipient weathering rind development. The second is to investigate the factors responsible for spatial variability in the nature and rates of rind formation. Incipient weathering rind development on the outer edges of the granite discs is observable and measurable over a period of time as short as four years in the mild Arctic alpine environment of Swedish Lapland. The earliest stages of rind development involve the development of a porous structure consisting of a combination of pits and fractures which have been solutionally enlarged and modified. Solution appears to be preferentially concentrated on the surfaces of feldspars and, to a lesser extent, quartz. In addition, iron oxides are present along grain boundaries and in grain interiors and are interpreted to have been derived from the oxidation of ferromagnesian minerals. Spatial variability in weathering rind development appears to be particularly driven by differences in moisture but is not related to soil pH. Copyright © 2005 John Wiley & Sons, Ltd.