Modelling the time‐dependent behaviour of granular material with hypoplasticity

This paper presents a constitutive model for time‐dependent behaviour of granular material. The model consists of 2 parts representing the inviscid and viscous behaviour of granular materials. The inviscid part is a rate‐independent hypoplastic constitutive model. The viscous part is represented by a rheological model, which contains a high‐order term denoting the strain acceleration. The proposed model is validated by simulating some element tests on granular soils. Our model is able to model not only the non‐isotach behaviour but also the 3 creep stages, namely, primary, secondary, and tertiary creep, in a unified way.     

The Cabo Frio upwelling overrides geographical patterns in the population dynamics of the shrimp Artemesia longinaris Spence Bate, 1888 (Decapoda: Penaeidae)

Upwelling areas are among the most productive ecosystems on the planet, influencing the biology of marine organisms. This study investigated the population dynamics of the shrimp Artemesia longinaris in two regions in southeastern Brazil, one inside (Macaé—Rio de Janeiro State) and one outside (Ubatuba—State of São Paulo) the Cabo Frio upwelling area. The aim was to verify the influence of the upwelling phenomenon on the abundance, growth, longevity, size of sexual maturity, and reproductive period of the species. In total, 188,902 individuals were captured at Macaé and 3,461 at Ubatuba. Individuals captured at Macaé showed larger maximum size, higher longevity, and slower growth rate, besides reaching sexual maturity at larger sizes than at Ubatuba. Continuous reproduction was observed in both regions, with juvenile recruitment peaks in spring and summer. Local conditions observed at Macaé were influenced by the Cabo Frio upwelling zone, characterized by productive and cooler waters that are around 20°C during most of the year. The upwelling phenomenon is probably the main factor influencing the population parameters studied here, changing the geographic patterns previously observed for the variation of these parameters in A. longinaris.     

中国区域水资源系统韧性与效率的发展协调关系评价

在界定水资源系统韧性概念的基础上，综合应用赋权法、SBM-DEA模型及发展协调度模型对2000—2016年全国31个省区的区域水资源系统的效率、韧性以及两者之间的发展协调关系进行评价。结果如下：① 中国区域水资源系统效率整体上处于非有效区，在研究期间呈现出在波动中上升的趋势。② 中国区域水资源系统韧性的平均值为0.39，总体水平较低，研究期间整体上呈波动上升趋势。③ 中国区域水资源系统效率与韧性的发展度总体呈现平稳上升—较快上升—急剧下降—上升的趋势，呈倒“U”型发展；协调度在研究期间呈现波动中上升趋势，大部分省区的发展度较好，而协调度较弱。④ 2000—2016年，水资源效率与韧性的发展协调度一直维持着“东?中?西”阶梯式递减格局，失调省区的数量降低，初级协调、中级协调、良好协调省区的数量逐渐提升，整体向协调趋势发展；空间格局上呈现由2000年倒“E”型对称式分布格局向如今北部围绕天津、中部围绕上海、南部围绕广东的三级格局演变，格局分布与中国三大经济区基本吻合，可知中国水资源系统效率与韧性的发展协调水平与经济发展水平之间有明显的关系。     

Near‐surface hydraulic conductivity of northern hemisphere glaciers

The hydrology of near‐surface glacier ice remains a neglected aspect of glacier hydrology despite its role in modulating meltwater delivery to downstream environments. To elucidate the hydrological characteristics of this near‐surface glacial weathering crust, we describe the design and operation of a capacitance‐based piezometer that enables rapid, economical deployment across multiple sites and provides an accurate, high‐resolution record of near‐surface water‐level fluctuations. Piezometers were employed at 10 northern hemisphere glaciers, and through the application of standard bail–recharge techniques, we derive hydraulic conductivity (K) values from 0.003 to 3.519 m day^?1, with a mean of 0.185 ± 0.019 m day^?1. These results are comparable to those obtained in other discrete studies of glacier near‐surface ice, and for firn, and indicate that the weathering crust represents a hydrologically inefficient aquifer. Hydraulic conductivity correlated positively with water table height but negatively with altitude and cumulative short‐wave radiation since the last synoptic period of either negative air temperatures or turbulent energy flux dominance. The large range of K observed suggests complex interactions between meteorological influences and differences arising from variability in ice structure and crystallography. Our data demonstrate a greater complexity of near‐surface ice hydrology than hitherto appreciated and support the notion that the weathering crust can regulate the supraglacial discharge response to melt production. The conductivities reported here, coupled with typical supraglacial channel spacing, suggest that meltwater can be retained within the weathering crust for at least several days. Not only does this have implications for the accuracy of predictive meltwater run‐off models, but we also argue for biogeochemical processes and transfers that are strongly conditioned by water residence time and the efficacy of the cascade of sediments, impurities, microbes, and nutrients to downstream ecosystems. Because continued atmospheric warming will incur rising snowline elevations and glacier thinning, the supraglacial hydrological system may assume greater importance in many mountainous regions, and consequently, detailing weathering crust hydraulics represents a research priority because the flow path it represents remains poorly constrained.     

Analysis of nonstationarity in low flow in the Loess Plateau of China

An approach for nonstationary low‐flow frequency analysis is developed and demonstrated on a dataset from the rivers on the Loess Plateau of China. Nonstationary low‐flow frequency analysis has drawn significant attention in recent years by establishing relationships between low‐flow series and explanatory variables series, but few studies have tested whether the time‐varying moments of low flow can be fully described by the time‐varying moments of the explanatory variables. In this research, the low‐flow distributions are analytically derived from the 2 basic explanatory variables—the recession duration and the recession coefficient—with the assumption that the recession duration and recession coefficient variables follow exponential and gamma distributions, respectively; the derived low‐flow distributions are applied to test whether the time‐varying moments of explanatory variables can explain the nonstationarities found in the low‐flow variable. The effects of ecosystem construction measures, that is, check dam, terrace, forest, and grassland, on the recession duration and recession coefficient are further discussed. Daily flow series from 11 hydrological stations from the Loess Plateau are used and processed with a moving average technique. Low‐flow data are extracted following the pit under threshold approach. Six of the 11 low‐flow series show significant nonstationarities at the 5% significance level, and the trend curves of the moments of low flow are in close agreement with the curves estimated from the derived distribution with time‐dependent moments of the recession duration and time‐constant moments of the recession coefficient. It is indicated that the nonstationarity in the low‐flow distribution results from the nonstationarity in the recession duration in all 6 cases, and the increase in the recession duration is resulted from large‐scale ecosystem constructions rather than climate change. The large‐scale ecosystem constructions are found to have more influence on the decrease in streamflow than on the increase in watershed storage, thus resulting in the reduction of low flow. A high return period for the initial fixed design value decreases dramatically with an increasing recession duration.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

High‐T and low‐P metamorphism in the Xilingol Complex of central Inner Mongolia,China: An indicator of extension in a previous orogeny

The Xilingol Complex comprises biotite gneisses and amphibolite interlayers with extensive migmatization. Four representative samples were documented and found to record either two or three metamorphic stages. Phase modelling using thermocalc suggests that the observed assemblages represent the final stages that underwent cooling from temperature peaks, and are consistent with a fluid‐absent solidus in P–T pseudosections. Their P–T conditions are further constrained to be 5–6 kbar/680–725°C and 4–5 kbar/650–680°C for two garnet‐bearing gneiss samples, 4–5 kbar/660–730°C for a cordierite‐bearing gneiss sample, and 4–5 kbar/680–710°C for an amphibolite sample based on mineral composition isopleths, involving measured Mg content in biotite, anorthite in plagioclase, grossular and pyrope in garnet and Ti content in amphibole. The peak temperature conditions recovered are 760–790°C or >760°C at 5–6 kbar based on the composition isopleths of plagioclase, biotite, garnet and especially the comparison of melt contents between the calculated and observed. A pre‐peak heating process with slight decompression can be suggested for some samples on the basis of the core–rim increase in the plagioclase anorthite, and the stability of ilmenite. Zircon U–Pb dating using the LA‐ICP‐MS method provides systemic constraints on the metamorphic ages of the Xilingol Complex to be 348–305 Ma, interpreted to represent the post‐peak cooling stages. Moreover, metagabbroic dykes that intruded into the Xilingol Complex yield 317 ± 3 Ma from magmatic zircon, and are considered to have played a significant role for heat advection triggering the high‐T and low‐P metamorphism. Thus, the clockwise P–T paths involving pre‐peak heating, peak and post‐peak cooling recovered for the Xilingol Complex are consistent with an extensional setting in the Carboniferous that developed on a previous orogen in response to addition of mantle‐derived materials probably together with upwelling of the asthenospheric mantle.     

Neogene upper‐crustal cooling of the Olympus range (northern Aegean): major role of Hellenic back‐arc extension over propagation of the North Anatolia Fault Zone

The North Anatolian Fault Zone (NAFZ) is one of the most hazardous active faults on Earth, yet its Pliocene space‐time propagation across the north Aegean domain remains poorly constrained. We use low‐temperature multi‐thermochronology and inverse thermal modelling to quantify the cooling history of the upper crust across the Olympus range. This range is located in the footwall of a system of normal faults traditionally interpreted as resulting from superposed Middle–Late Miocene N–S stretching, related to the back‐arc extension of the Hellenic subduction zone, and a Pliocene‐Quaternary transtensional field, attributed to the south‐westward propagation of the NAFZ. We find that accelerated exhumational cooling occurred between 12 and 6 Ma at rates of 15–35 °C Ma^?1 and decreased to <3 °C Ma^?1 by 8–6 Ma. The absence of significant Plio‐Pleistocene cooling across Olympus suggests that crustal exhumation there is driven by late Miocene back‐arc extension, while the impact of the NAFZ remains limited.