Model identification for hydrological forecasting under uncertainty

Methods for the identification of models for hydrological forecasting have to consider the specific nature of these models and the uncertainties present in the modeling process. Current approaches fail to fully incorporate these two aspects. In this paper we review the nature of hydrological models and the consequences of this nature for the task of model identification. We then continue to discuss the history (“The need for more POWER‘’), the current state (“Learning from other fields”) and the future (“Towards a general framework”) of model identification. The discussion closes with a list of desirable features for an identification framework under uncertainty and open research questions in need of answers before such a framework can be implemented.     

Assimilation and Fractional Crystallization Controlled by Transport Process of Crustal Melt: Implications from an Alkali Basalt-Dacite Suite from Rishiri Volcano, Japan

Mechanisms of fractional crystallization with simultaneous crustalassimilation (AFC) are examined for the Kutsugata and Tanetomilavas, an alkali basalt–dacite suite erupted sequentiallyfrom Rishiri Volcano, northern Japan. The major element variationswithin the suite can be explained by boundary layer fractionation;that is, mixing of a magma in the main part of the magma bodywith a fractionated interstitial melt transported from the mushyboundary layer at the floor. Systematic variations in SiO₂ correlatewith variations in the Pb, Sr and Nd isotopic compositions ofthe lavas. The geochemical variations of the lavas are explainedby a constant and relatively low ratio of assimilated mass tocrystallized mass (‘r value’). In the magma chamberin which the Kutsugata and Tanetomi magmas evolved, a strongthermal gradient was present and it is suggested that the marginalpart of the reservoir was completely solidified. The assimilantwas transported by crack flow from the partially fused floorcrust to the partially crystallized floor mush zone throughfractures in the solidified margin, formed mainly by thermalstresses resulting from cooling of the solidified margin andheating of the crust. The crustal melt was then mixed with thefractionated interstitial melt in the mushy zone, and the mixedmelt was further transported by compositional convection tothe main magma, causing its geochemical evolution to be characteristicof AFC. The volume flux of the assimilant from the crust tothe magma chamber is suggested to have decreased progressivelywith time (proportional to t^–1/2), and was about 3 x 10^–2m/year at t = 10 years and 1 x 10^–2 m/year at t = 100years. It has been commonly considered that the heat balancebetween magmas and the surrounding crust controls the couplingof assimilation and fractional crystallization processes (i.e.absolute value of r). However, it is inferred from this studythat the ratio of assimilated mass to crystallized mass canbe controlled by the transport process of the assimilant fromthe crust to magma chambers. KEY WORDS: assimilation and fractional crystallization; mass balance model; magma chamber; melt transport; Pb isotope     

An alternative Boussinesq equation considering the effect of hysteresis on coastal groundwater waves

This paper presents an alternative Boussinesq equation considering hysteresis effect via a third‐order derivative term. By introducing an improved moisture–pressure retention function, this equation describes, with reasonable precision, groundwater propagation in coastal aquifers subject to Dirichlet boundary condition of different oscillation frequencies. Test results confirmed that it is necessary to consider horizontal and vertical flows in unsaturated zone, because of their variable influences on hysteresis. Hysteresis in unsaturated zone can affect the water table wave number of groundwater wave motion, such as wave damping rate and phase lag. Oscillations with different periods exert different hysteresis effect on wave propagation. Truncation/shrinkage of unsaturated zones also affects the strength of hysteresis. These impacts can be reflected in the alternative Boussinesq equation by adjusting the parameter representing the variation rate of moisture associated with pressure change, as opposed to traditional computationally expensive hysteresis algorithms. The present Boussinesq equation is simple to use and can provide feasible basis for future coupling of groundwater and surface water models. Copyright © 2016 John Wiley & Sons, Ltd.     

Glacial flood pulse effects on benthic fauna in equatorial high‐Andean streams

Equatorial glacier‐fed streams present unique hydraulic patterns when compared to glacier‐fed observed in temperate regions as the main variability in discharge occurs on a daily basis. To assess how benthic fauna respond to these specific hydraulic conditions, we investigated the relationships between flow regime, hydraulic conditions (boundary Reynolds number, Re*), and macroinvertebrate communities (taxon richness and abundance) in a tropical glacier‐fed stream located in the high Ecuadorian Andes (> 4000 m). Both physical and biotic variables were measured under four discharge conditions (base‐flow and glacial flood pulses of various intensities), at 30 random points, in two sites whose hydraulic conditions were representative to those found in other streams of the study catchment. While daily glacial flood pulses significantly increased hydraulic stress in the benthic habitats (appearance of Re* > 2000), low stress areas still persisted even during extreme flood events (Re* < 500). In contrast to previous research in temperate glacier‐fed streams, taxon richness and abundance were not significantly affected by changes in hydraulic conditions induced by daily glacial flood pulses. However, we found that a few rare taxa, in particular rare ones, preferentially occurred in highly stressed hydraulic habitats. Monte‐Carlo simulations of benthic communities under glacial flood reduction scenarios predicted that taxon richness would be significantly reduced by the loss of high hydraulic stress habitats following glacier shrinking. This pioneer study on the relationship between hydraulic conditions and benthic diversity in an equatorial glacial stream evidenced unknown effects of climate change on singular yet endangered aquatic systems. Copyright © 2013 John Wiley & Sons, Ltd.     

Ocean current estimation using a Multi-Model Ensemble Kalman Filter during the Grand Lagrangian Deployment experiment (GLAD)

In the summer and fall of 2012, during the GLAD experiment in the Gulf of Mexico, the Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE) used several ocean models to assist the deployment of more than 300 surface drifters. The Navy Coastal Ocean Model (NCOM) at 1 km and 3 km resolutions, the US Navy operational NCOM at 3 km resolution (AMSEAS), and two versions of the Hybrid Coordinates Ocean Model (HYCOM) set at 4 km were running daily and delivering 72-h range forecasts. They all assimilated remote sensing and local profile data but they were not assimilating the drifter’s observations. This work presents a non-intrusive methodology named Multi-Model Ensemble Kalman Filter that allows assimilating the local drifter data into such a set of models, to produce improved ocean currents forecasts. The filter is to be used when several modeling systems or ensembles are available and/or observations are not entirely handled by the operational data assimilation process. It allows using generic in situ measurements over short time windows to improve the predictability of local ocean dynamics and associated high-resolution parameters of interest for which a forward model exists (e.g. oil spill plumes). Results can be used for operational applications or to derive enhanced background fields for other data assimilation systems, thus providing an expedite method to non-intrusively assimilate local observations of variables with complex operators. Results for the GLAD experiment show the method can improve water velocity predictions along the observed drifter trajectories, hence enhancing the skills of the models to predict individual trajectories.     

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

Revised runoff curve number for runoff prediction in the Loess Plateau of China

The soil conservation service (now Natural Resources Conservation Service) Curve Number (SCS-CN), one of the most commonly used methods for surface runoff prediction. The runoff calculated by this method was very sensitive to CN values. In this study, CN values were calculated by both arithmetic mean (CN_C) and least square fit method (CN_F) using observed rainfall-runoff data from 43 sites in the Loess Plateau region, which are considerably different from the CN₂ values obtained from the USDA-SCS handbook table (CN_T). The results showed that using CN_C instead of CN_T for each watershed produce little improvement, while replacing CN_T with CN_F improves the performance of the original SCS-CN method, but still performs poorly in most study sites. This is mainly due to the SCS-CN method using a constant CN value and discounting of the temporal variation in rainfall-runoff process. Therefore, three factors—soil moisture, rainfall depth and intensity—affecting the surface runoff variability are considered to reflect the variation of CN in each watershed, and a new CN value was developed. The reliability of the proposed method was tested with data from 38 watersheds, and then applied to the remaining five typical watersheds using the optimized parameters. The results indicated that the proposed method, which boosted the model efficiencies to 81.83% and 74.23% during calibration and validation cases, respectively, performed better than the original SCS-CN and the Shi and Wang (2020b) method, a modified SCS-CN method based on tabulated CN value. Thus, the proposed method incorporating the influence of the temporal variability of soil moisture, rainfall depth, and intensity factors suggests an accurate runoff prediction for general applications under different hydrological and climatic conditions on the Loess Plateau region.     

Halin图的色数问题

Halin图G=T∪C,其中T为每一非悬挂点(内点)度数至少为3的平面树,C为连接T的所有悬挂点的圈.文章分别讨论了Halin图的星色数、面色数及分数色数.     

天津城区春节期间大气气溶胶污染特征 和数浓度分布

利用天津城区2009-2014年春节期间大气气溶胶观测资料和相关气象资料,重点分析2013和2014年春节期间气溶胶污染特征,探求燃放烟花爆竹以及气象条件对春节期间大气气溶胶的影响。结果表明,受燃放烟花爆竹影响,春节期间PM_(2.5)质量浓度最高值均发生在除夕夜间;持续雾霾天气条件下燃放烟花爆竹,造成2013年除夕夜间PM_(2.5)质量浓度峰值达到1240μg·m～(-3),是近年来最严重的一次;2014年春节期间烟花爆竹燃放量有所减少,加之空气扩散条件较为有利,PM_(2.5)质量浓度显著低于2013年;不同天气条件下,气溶胶数浓度谱分布特征存在明显差异,燃放烟花爆竹期间气溶胶数浓度水平与严重雾-霾天气相当。