CS2: A PIECEWISE-LINEAR MODEL FOR LARGE STRAIN CONSOLIDATION

This paper presents a piecewise-linear finite-difference model for one-dimensional large strain consolidation called CS2. CS2 is developed using a fixed Eulerian co-ordinate system and constitutive relationships which are defined by discrete data points. The model is dimensionless such that solutions are independent of the initial height of the compressible layer and the absolute magnitude of the hydraulic conductivity of the soil. The capability of CS2 is illustrated using four example problems involving small strain, large strain, self-weight, and non-linear constitutive relationships. In each case, the performance of the model is comparable to other available analytical and numerical solutions. Using CS2, correction factors are developed for the conventional Terzaghi theory which account for the effect of vertical strain on computed values by elapsed time and maximum excess pore pressure during consolidation. © 1997 John Wiley & Sons, Ltd.     

Flow separation on Zongo Glacier,Cordillera Real,Bolivia

Meltwaters collected from the proglacial stream escaping from Zongo Glacier (2·1 km²), Bolivia (16°S), have been monitored in order to analyse the internal drainage system of an Andean glacier. Electrical conductivity has been measured sporadically between February 1995 and March 1996, during 16 one-day field surveys, under various meteorological conditions in summer and winter. The mixing-model technique based on the electrical conductivity is used for a quantitative separation of discharge which is derived from continuous water level registration. Tracer experiments (mainly uranine dye and NaCl salt) have been carried out from March to June 1997 to obtain information about the internal drainage system. In the tropical Andes, accumulation only occurs in austral summer, whereas ablation occurs throughout the year and is higher during the accumulation season, between November and March. The assumptions involved in the use of mixing models for analysis of glacial drainage structure are applicable for tropical glaciers because glacial conduits do not suffer complete closure, and are permanently supplied by meltwaters, even in wintertime. Two components of discharge are separated: an englacial flow originating from surface meltwater which is routed without chemical enrichment, and offering low electrical conductivity; and a subglacial one routed in contact with bedrock or sediments showing high ionic concentrations. Electrical conductivity of meltwater varies diurnally, inversely to discharge fluctuations. According to this behaviour, total discharge is mainly formed by the englacial component. The drainage structures for englacial and subglacial flow have to be widely interconnected, as indicated by diurnal variations of the subglacial discharge. Comparison of hydrograph separation based on conductivity and on ¹⁸O isotope confirms that the subglacial flow is influenced by surface melting. A hydrograph separation of the subglacial flow is proposed, between a diurnal variable component, composed of water coming from the englacial network, and a base flow, which may vary seasonally. The dye tracing experiments confirm the drainage complexity of Zongo Glacier and demonstrate the interest of identifying three main drainage components. © 1998 John Wiley & Sons, Ltd.     

Assessment of observational networks with the Representer Matrix Spectra method—application to a 3D coastal model of the Bay of Biscay

The development of coastal ocean modeling in the recent years has allowed an improved representation of the associated complex physics. Such models have become more realistic, to the point that they can now be used to design observation networks in coastal areas, with the idea that a “good” network is a network that controls model state error. To test this ability without performing data assimilation, we set up a technique called Representer Matrix Spectra (RMS) technique that combines the model state and observation error covariance matrices into a single scaled representer matrix. Examination of the spectrum and the eigenvectors of that matrix informs us on which model state error modes a network can detect and constrain amidst the observation error background. We applied our technique to a 3D coastal model in the Bay of Biscay, with a focus on mesoscale activity, and tested the performance of various altimetry networks and an in situ array deployment strategy. It appears that a single nadir altimeter is not efficient enough at capturing coastal mesoscale physics, while a wide swath altimeter would do a much better job. Testing various local in situ array configurations confirms that adding a current meter to a vertical temperature measurement array improves the detection of secondary variability modes, while shifting the array higher on the shelf break would obviously enhance the model constraint along the coast. The RMS technique is easily set up and used as a “black box,” but the utility of its results is maximized by previous knowledge of model state error physics. The technique provides both quantitative (eigenvalues) and qualitative (eigenvectors) tools to study and compare various network options. The qualitative approach is essential to discard possibly inconsistent modes.     

Sensitivity-guided reduction of parametric dimensionality for multi-objective calibration of watershed models

Problem complexity for watershed model calibration is heavily dependent on the number of parameters that can be identified during model calibration. This study investigates the use of global sensitivity analysis as a screening tool to reduce the parametric dimensionality of multi-objective hydrological model calibration problems while maximizing the information extracted from hydrological response data. This study shows that by expanding calibration problem formulations beyond traditional, statistical error metrics to also include metrics that capture indices or signatures of hydrological function, it is possible to reduce the complexity of calibration while maintaining high quality model predictions. The sensitivity-guided calibration is demonstrated using the Sacramento Soil Moisture Accounting (SAC-SMA) conceptual rainfall–runoff model of moderate complexity (i.e., up to 14 freely varying parameters). Using both statistical and hydrological metrics, optimization results demonstrate that parameters controlling at least 20% of the model output variance (through individual effects and interactions) should be included in the calibration process. This threshold generally yields 30–40% reductions in the number of SAC-SMA parameters requiring calibration – setting the others to a priori values – while maintaining high quality predictions. Two parameters are recommended to be calibrated in all cases (percent impervious area and lower zone tension water storage), three parameters are needed in drier watersheds (additional impervious area, riparian zone vegetation, and percent of percolation going to tension storage), and the lower zone parameters are crucial unless the watershed is very dry. Overall, this study demonstrates that a coupled, multi-objective sensitivity and calibration analysis better captures differences between watersheds during model calibration and serves to maximize the value of available watershed response time series. These contributions are particularly important given the ongoing development of more complex integrated models, which will require new tools to address the growing discrepancy between the information content of hydrological data and the number of model parameters that have to be estimated.     

横向非均匀地球负荷问题的CLFE有限元算法的有效性

本文设计了三个数值算例分别模拟径向分层黏弹地球对二阶球谐负荷、均厚圆盘负荷和全球末次冰期真实冰负荷的响应.利用耦合位扰动Laplace方程的三维有限元法(CLFE),计算出地球表面的位移、大地水准面的变化、冰后海平面变化和这些量的现今变化速率,并将其结果与传统谱方法的结果进行对比.结果表明,两种方法的结果总体上表现出良好的一致性,因此,CLFE方法被证实是正确的、有效的、能被用于研究地幔横向非均匀性的影响.值得注意的是,新的方法完全没有横向非均匀小扰动假设的限制.     

Caldera rim collapse: A hidden volcanic hazard

Following the emblematic flank collapse of Mount St Helens in 1981, numerous models of flank sliding have been proposed. These models have allowed to largely improve the understanding of mechanisms involved in such landslides, which represent a tremendous risk for populations living around volcanoes. In this article, a new mode of landslide formation, related to buried calderas, is described. The model emphasizes the paramount importance of the hidden ring fault that, even when the caldera is buried, still remains a plane of weakness in the core of the edifice. Under certain conditions, this plane of weakness becomes activated as the upper part of a pre-existing critical slip surface and is used in the emplacement of huge landslides which travel downslope at a very high velocity. A natural example is taken from Piton de la Fournaise Volcano (La Réunion Island, Indian Ocean). It reveals that the primary cause triggering caldera rim collapse is partial unbuttressing of the flank of the volcano. In the natural example, this occurs through regressive erosion that excavates deep canyon in the direction of the buried caldera but other mechanisms may exist. On account of the large volumes of material involved in caldera rim collapse as well as their long runout distances, such a volcanic hazard should be taken into account on every volcano where buried calderas are suspected.     

Hypoxia impairs embryo development and survival in black bream (Acanthopagrus butcheri)

Coastal environments are threatened by the increasing frequency, extent and severity of hypoxic events. Hypoxia affects vast areas around the world and often causes fish kills, reduced abundance, altered distribution, low benthic biomass and declines in fisheries. In Australia, many fisheries are based on sparid fishes and in the southern states black bream (Acanthopagrus butcheri) is important to both the recreational and commercial sectors. This species completes its entire life cycle in estuaries and annual recruitment is highly variable and very likely influenced by environmental conditions during the spawning season.In a laboratory-based experiment, fertilised black bream eggs (embryos) were exposed to five different levels of dissolved oxygen (DO). The DO levels were maintained in small test wells using nitrogen gas in a novel chamber design. Embryo development was assessed over a 2-day period and hatched larvae were observed until Day 2 post-hatch. Significant differences (p < 0.05) were observed in embryonic development and survival as a function of DO level. In severely hypoxic conditions (30% saturation) survival to 1 day was reduced and no hatching occurred. In moderately hypoxic conditions (45–55%S), both precocious and delayed hatching was observed and hatch rates were reduced, whilst the number of hatched larvae with deformities increased, resulting in reduced larval lengths. No larvae survived to Day 2 post-hatch when held in hypoxic conditions (<55%S). This study demonstrates the detrimental effect that severe hypoxia can have on the early development of black bream which could result in reduced recruitment and lowered abundance. Other species that share similar early life histories may also be at risk.     

Geomorphological changes in an arid transgressive coastal dune field due to natural processes and human impacts

Geomorphological changes in recent decades in an arid transgressive coastal dune system (Maspalomas, Gran Canaria, Canary Islands) are analyzed. The methodology used is based on the generation of two geomorphological maps (1961 and 2003) by interpretation of digital orthophotos. The overlay of both maps in a geographic information system (GIS) enabled the spatial and surface changes of the landforms to be determined, and the processes that generated these changes. Twelve cultural and geomorphological processes were identified from highest to lowest importance, namely: anthropization by urban occupation (114 ha changed), stabilization (92.5 ha), barchanization (37 ha), salinization/halophytication (15 ha), anthropization (12.4 ha), deflation (11.8 ha), dune loss/beach gain (11.3 ha), dune formation (9.6 ha), progradation (8 ha), retrogradation (7.7 ha), destabilization (2.7 ha) and flooding (0.7 ha). Geomorphological changes are associated with a combination of five main factors, three of anthropogenic origin and two natural ones. The natural factors are: (1) the arid climate, which favors changes occurring at high speed; (2) the existence of a progressive sedimentary deficit. Anthropogenic factors are: (3) construction of tourist urbanizations, infrastructures and facilities; (4) installation of equipment or infrastructure on the beaches; (5) the activities carried out by users. These human factors have altered the aeolian dynamics and reduced the area occupied by vegetation in some areas, causing changes in aeolian sedimentary processes. The geomorphological processes identified can be used as indicators of environmental change, allowing us to synthesize the changes in landforms detected, and group all combinations derived from the analysis by GIS and analyze them spatially. Thus, the environmental changes in the transgressive coastal dune systems could be interpreted more effectively. Copyright © 2018 John Wiley & Sons, Ltd.