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.     

Distributed Temperature Sensing to Measure Infiltration Rates Across a Groundwater Recharge Basin

Managed aquifer recharge is used to augment groundwater resources and provide resiliency to water supplies threatened by prolonged droughts. It is important that recharge facilities operate at their maximum efficiency to increase the volume of water stored for future use. In this study, we evaluate the use of distributed temperature sensing (DTS) technology as a tool to measure high-resolution infiltration rates at a large-scale recharge facility. Fiber optic cable was laid out inside a spreading basin in a spiral pattern, at two different depths. The cables measured the propagation of diurnal surface water temperature oscillations into the basin depth. The rate of heat propagation is proportional to the velocity of the water, making it possible to estimate the infiltration rate from the temperature measurements. Our results showed that the infiltration rate calculated from DTS, averaged over the entire basin, was within 5% of the infiltration rate calculated using a conventional metering method. The high-resolution data obtained from DTS, both spatially and temporally, revealed heterogeneous infiltration rates throughout the basin; furthermore, tracking the evolution of infiltration rates over time revealed regions with consistently high infiltration rates, regions with consistently low infiltration rates, and regions that evolved from high to low rates, which suggested clogging within that region. Water utilities can take advantage of the high-resolution information obtained from DTS to better manage recharge basins and make decisions about cleaning schedule, frequency, and extent, leading to improved basin management strategies, reduced O&M costs, and increased groundwater recharge.     

Sediment transport trends and cross-sectional stability of a lagoonal tidal inlet on the Central Coast of Vietnam

Lagoonal tidal inlets are a typical morphology of the Central Coast of Vietnam. Recently, navigation channels in these inlets have become increasingly threatened by siltation. This study analyses the relations between sediment distribution and transport trends (using the technique of Sediment Trend Analysis-STA■) in the lagoonal system of the De Gi inlet and then proposes appropriate countermeasures against sand deposition in the navigation channel. The STA identified three types of transport trends in the De Gi inlet, namely dynamic equilibrium, net accretion, and net erosion. Processes associated with the tidal prism have resulted in trends of sediment transport and deposition across the flood and ebb tidal shoals, which maintain a present cross-sectional area of about 1000m^2. However, longshore sediment transport from north to south resulting from northeast waves cause additional sand deposition in the channel. In addition, the effects of refraction associated with a nearby headland and jetty also increase sedimentation. These processes provide the main reasons for sediment deposition in the De Gi inlet. Short term and regular dredging helps to maintain the navigation channel. A system comprised of three jetties (north, south, and weir) is necessary to ensure the longterm cross-sectional stability of the navigation channel.     

Determining flow patterns and emplacement dynamics from tsunami deposits with no visible sedimentary structure

In the absence of eyewitness reports or clear sedimentary structures, it can be difficult to interpret tsunami deposits or reconstruct tsunami inundation patterns. The emplacement dynamics of two historical tsunami deposits were investigated at seven transects in Okains Bay, New Zealand, using a combined geospatial, geomagnetic and sedimentological approach. The tsunami deposits are present as layers of sand and silt intercalated between soils and become finer and thinner with distance inland. The deposits are attributed to the 1960 and possibly the 1868 tsunamis, based on radiometric dating and correlation with historical records. Measurements of Magnetic Fabric (MF: Anisotropy of Magnetic Susceptibility) and particle size were used to reconstruct the evolution of flow dynamics laterally and vertically. A combination of statistical methods, including spatial autocorrelation testing, Spearman's rank order correlation, Principal Component Analysis (PCA) and K‐means cluster analysis, was applied to examine relationships between MF parameters and sediment texture, and infer depositional hydrodynamics. Flow patterns deduced from MF show the estuary channel acted as a conduit for inundation, with flow commonly aligned sub‐perpendicular to the estuary bed. MF and sediment data suggest deposition occurred from settling during laminar flow. Evidence of both uprush and backwash deposition, as well as wave reflection from infrastructure, was found. Statistical analysis of data showed significant relationships between grain size parameters and MF parameters associated with flow speed and magnetic fabric type. PCA and cluster analysis differentiated samples into two primary hydrodynamic groups: (1) samples deposited from laminar flow; and (2) samples deposited close to the limit of inundation, which includes samples deposited further inland, those affected by flow convergence, and those in the upper part of tsunami deposits. This approach has potential as a tool for reconstructing hydrodynamic conditions for palaeotsunamis and by combining spatial and statistical analyses, large‐scale investigations can be more easily performed. Copyright © 2016 John Wiley & Sons, Ltd.