Towards medium-term (order of months) morphodynamic modelling of the Teign estuary, UK

The main objective of this paper is to address the principal mechanisms involved in the medium-term (order of months to years) morphodynamic evolution of estuaries through the application of a process-based numerical modelling. The Teign estuary (Teignmouth, UK) is the selected site. The system is forced by the macrotidal semi-diurnal tide in the English Channel and is perturbed to a minor extent by high river discharge events (freshets). Although waves have a definite influence on the adjacent coastal area, Wells (Teignmouth Quay Development Environmental Statement: Changes to Physical Processes. Report R.984c:140. ABP Marine Environmental Research Ltd., Southampton, 2002b) suggested that swell waves do not enter the estuary. Hence, wave effects are neglected in this study, as only tides and the river discharge are taken into account. The sediment grain size is highly variable, but mainly sandy. Within the frame of the COAST3D project (), four bathymetric surveys of the adjacent coastal area were carried out at a nearly weekly intervals. The outer estuary and the adjacent coastal area were also surveyed every 6 months as part of the COASTVIEW project (). Based on these data and on continuously measured parameters, such as water level, waves, wind and river discharge, numerical modelling of the morphodynamic processes can be tested. To replicate the morphological changes in the medium-term within a feasible simulation time, forcing conditions are reduced through the use of an input reduction method (called ensemble technique). In this study, simulations are based on the coupling between Telemac-2D and its non-cohesive sediment transport module, Sisyphe (version 5.3 for both modules). Three different sediment transport formulae were tested: (1) Engelund and Hansen (A monograph on sediment transport in alluvial streams, 3rd edn. Technological University of Denmark, Copenhagen, 1967) including the modifications proposed by Chollet and Cunge (J Hydraul Eng 17(1):1–13, 1979); (2) Bijker (Mechanics of sediment transport by the combination of waves and current. In: Design and reliability of coastal structures. 23rd international conference on Coastal Engineering, pp 147–173, 1968) and (3) Soulsby (Dynamics of Marine Sands. A manual for practical applications. HR Wallingford, Wallingford, p 142, 1997) modified version of van Rijn [J Hydraul Eng 110(10):1431–1456, 1984a, J Hydraul Eng 110(11):1613–1641, 1984b] formulation. Both a qualitative (i.e. visual comparison) and a quantitative tool [Brier Skill Score (BSS); described in Sutherland et al. in Coast Eng 51:917–939, 2004b] are applied to assess the similarity of simulations when compared to model predictions and observations. Tests confirmed the reliability and time efficiency of the ensemble technique, since it reproduced very well the results of a reference run, a computation based on the observed boundary conditions. For the spring-neap cycle modelled, the BSS was of 0.91 (a perfect modelling would have a BSS of 1), with a reduction in the simulation time on the order of 80%. For the 6-month-period simulation, results were also excellent: BSS=0.92 and a computer time reduction of 85%. In principle, this method has the advantage of being applied to any process-based numerical model.     

Axial poloidal electromagnetic core-mantle coupling torque: A re-examination for different conductivity and satellite supported geomagnetic field models

We investigate the temporal behaviour of the axial component of the electromagnetic core-mantle coupling torque that is associated with the poloidal part of the geomagnetic field observable at the Earth surface. For its computation, we use different models of the geomagnetic field, expanded into spherical harmonics (Wardinski and Holme, 2006; Sabaka et al., 2004), and the mantle conductivity. The geomagnetic field, which we have to know at the core-mantle boundary for the associated computations, will be inferred from the field at the Earth surface by the non-harmonic field continuation through a conducting mantle shell. The aims of this investigation are (i) to check how sensitive is the computation of the torque with respect to the different geomagnetic field models, (ii) to check its dependence on the spherical harmonic degree n, and (iii) to determine the difference between the mechanical torque derived from the observed length-of-day variations (atmospheric influence subtracted) and the poloidal electromagnetic torque in dependence on the assumed conductivity. To use the non-harmonic field continuation for the torque calculation and to obtain an insight into the influence of the different geomagnetic field models on the EM torques are the major aspects of this paper. grm@gfz-potsdam.de     

Intercomparisons of land-surface parameterizations coupled to a limited area forecast model

The goal of the Project for Intercomparison of Land-surface Parameterization Schemes (PILPS) is to improve the understanding of the interactions between the atmosphere and the continental surface in climate and weather forecast models. In PILPS Phase 4(b), selected schemes are coupled to the Limited Area Prediction System (LAPS) developed by the Australian Bureau of Meteorology. To facilitate the comparison of PILPS schemes' behavior within LAPS, a single mode of coupling is selected: explicit coupling. This type of coupling is more flexible and avoids most of the problems raised when interchanging the surface schemes. Exploratory tests are conducted. Initially, experiments are run in which the land-surface schemes use the same parameters as in their original host models. Then, in other runs, the most important surface parameters are set constant in an attempt to reduce the scatter amongst the schemes' results. In order to understand the impact of initialisation of soil moisture on the schemes' results some extreme cases (wet and dry) are performed. The partitioning between surface fluxes is studied as well as the soil moisture budget. Both regional and local results are analysed. Sensitivity between LSS is found in the precipitation field with rainfall over the Australian continent altering by about 20%, but no significant change is found in the net radiation. The scatter in the surface energy fluxes amongst the schemes is large (up to 300 W m⁻² locally, during the daytime peak) but is seldom affected by the choice of surface parameters. The dynamical range of flux partitioning between extremely dry and wet initialisation varies strongly amongst the schemes. Some major shortcoming with the BUCKET approach are seen in the re-evaporation of convective precipitation over dry land, in the very large evaporation from wet surfaces and the diurnal cycle of surface temperature.     

Internal testing of a numerical model of hillslope–channel coupling using laboratory flume experiments

The numerical model COUP 2D simulates the hydrological coupling between hillslopes and the river channel during a rainfall event. In order to test the numerical model, a 1:100 scaled laboratory flume which was modified to incorporate lateral hillslope elements, was used to run a series of experiments in which hillslope angle, channel angle, hillslope discharge and channel discharge were the varying parameters. Overall, there were 18 different experimental configurations with three replicates carried out for each condition, leading to a total of 54 experiments. These conditions were then used to parameterize and run COUP 2D. Internal model outputs of flow depth and flow velocity at four cross‐sections in the channel were compared to the measurements made in the physical model for the same parameter conditions. Statistical comparisons of the measured and modelled data were carried out for each experiment and across all experiments, using two goodness‐of‐fit measures—root mean square error and Nash–Sutcliffe coefficient of efficiency—in order to assess the performance of the model over an entire simulation as well as over all the simulations. The main effects on the goodness‐of‐fit measures for flow depth of each experimental variable, as well as the interactions between variables, were evaluated using statistical modelling. The results show that the model captures flow‐depth variations in response to changing channel and hillslope parameters. Statistical modelling suggests that the main effects on model error are cross‐section position, channel angle and channel discharge. Significant interactions also occur between all the channel variables and between the channel variables and hillslope discharge. The results of the testing procedure have significant implications for the consideration of different model components and for the interaction between data‐ and model evaluation. Copyright © 2007 John Wiley & Sons, Ltd.     

Coupled HM analysis using zero‐thickness interface elements with double nodes—Part II: Verification and application

In a companion Part I of this paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.735 ), a coupled hydro‐mechanical (HM) formulation for geomaterials with discontinuities based on the finite element method (FEM) with double‐node, zero‐thickness interface elements was developed and presented. This Part II paper includes the numerical solution of basic practical problems using both the staggered and the fully coupled approaches. A first group of simulations, based on the classical consolidation problem with an added vertical discontinuity, is used to compare both the approaches in terms of accuracy and convergence. The monolithic or fully coupled scheme is also used in an application example studying the influence of a horizontal joint in the performance of a reservoir subject to fluid extraction. Results include a comparison with other numerical solutions from the literature and a sensitivity analysis of the mechanical parameters of the discontinuity. Some simulations are also run using both a full non‐symmetric and a simplified symmetric Jacobian matrix. On top of verifying the model developed and its capability to reflect the conductivity changes of the interface with aperture changes, the results presented also lead to interesting observations of the numerical performance of the methods implemented. Copyright © 2008 John Wiley & Sons, Ltd.     

新型城镇化的人地耦合异速增长分析方法

在新型城镇化背景下的城市人口规模与用地面积变化具有复杂性,研究两种要素城镇化水平的相对增长关系,有助于探索人地耦合发展的协同态势及变化动因。本文选取广州市2000—2015年的Landsat时序影像与2000、2010和2015年的人口普查与年鉴数据,在街道（乡镇）尺度上计算人口土地城镇化耦合指数识别4种耦合类型,并在此基础上提出年均增长率异速关系拟合方法。研究结果表明：① 2000—2015年广州市街道（乡镇）人口和土地城镇化呈现明显的圈层结构（核心、内、外圈层）特征,4种人地耦合类型的异速增长分别呈现正异速、负异速、负幂律3种形式,证明人地耦合分类的异速增长建模方法有助于分析人地城镇化的动态演化过程;② 基于人地耦合指数的异速标度分析表明,核心圈层—内圈层—外圈层结构上分别形成了高度集约—集约—相对粗放式扩张的发展格局,结合城市更新改造的空间分析,说明这一特色新型城镇化模式在广州具有一定成效,已初步形成了核心圈层与内圈层协同发展的态势,但需要注意防控外圈层的粗放型增长。研究结果可为新型城镇化理论内涵的认知探索和广州市案例的实证研究提供支持。     

An elastoplastic strainhardening model for soil allowing for hydraulic bonding–debonding effects

The paper presents a strainhardening constitutive model for unsaturated soil behaviour based on energy conjugated stress variables in the framework of superposed continua. The proposed constitutive law deals with hydro‐mechanical coupling phenomena. The main purpose is to develop within a consistent framework a model that can deal with possible mechanical instabilities occurring in partially saturated materials. The loss of capillary effects during wetting processes can, in fact, play a central role in unstable processes. Therefore, it will be shown that the bonding effects due to surface tensions can be described in a mathematical framework similar to that employed for bonded geomaterials to model weathering or diagenesis effects, either mechanically or chemically induced. The results of several simulations of common laboratory tests on partially saturated soil specimens are shown. The calculated behaviour appears to be in good qualitative agreement with that observed in the laboratory. In particular it is shown that volumetric collapse phenomena due to hydraulic debonding effects can be successfully described by the model. Finally, it will be highlighted the ability of the model to naturally capture the transition to a fully saturated condition and to deal with possible mechanical instabilities in the unsaturated regime. Copyright © 2008 John Wiley & Sons, Ltd.     

Microbial biomass response to different quantities and sources of organic matter in Brazilian coastal sediments

This study investigates the benthic microbial responses to organic matter (OM) variations in quantity and sources in two shallow water bays (Fortaleza and Ubatuba Bays) on the SE coast of Brazil on six occasions during the year. The pelagic and benthic compartments of the bays were evaluated by: (i) nutrients and chlorophyll a (Chl a) in the water column; (ii) quantity and sources of OM in the sediment (Chl a, total organic carbon and total nitrogen and lipid biomarker composition); and (iii) microbial biomass in sediments as an indicator of active benthic response. Although there were changes in water‐column nutrients during the year, Chl a was fairly constant, suggesting a regular supply of microalgae‐derived OM to the sea bottom. Based on the composition of lipid biomarkers in sediments, OM sources were classified as mostly marine and with high contributions of labile (microalgae‐derived) OM. Labile OM composition varied from diatoms in the summer to phytoflagellates in the winter and tended to accumulate in areas protected by physical disturbances in one of the bays. Microbial biomass followed this trend and was 160% higher in protected than in exposed areas. This study suggests that the coupling between labile OM and benthic microbial biomass occurs primarily in protected areas, irrespective of the time of the year. Since meio‐ and macrofaunal assemblages depend upon secondary microbial production within the sediments, this coupling may have an important role for the benthic food‐web.