Modeling nonlinear wave–wave interactions with the elliptic mild slope equation

An approach is developed to simulate wave–wave interactions using nonlinear elliptic mild-slope equation in domains where wave reflection, refraction, diffraction and breaking effects must also be considered. This involves the construction of an efficient solution procedure including effective boundary treatment, modification of the nonlinear equation to resolve convergence issues, and validation of the overall approach. For solving the second-order boundary-value problem, the Alternating Direction Implicit (ADI) scheme is employed, and the use of approximate boundary conditions is supplemented, for improved accuracy, with internal wave generation method and dissipative sponge layers. The performance of the nonlinear model is investigated for a range of practical wave conditions involving reflection, diffraction and shoaling in the presence of nonlinear wave–wave interactions. In addition, the transformation of a wave spectrum due to nonlinear shoaling and breaking, and nonlinear resonance inside a rectangular harbor are simulated. Numerical calculations are compared with the results from other relevant nonlinear models and experimental data available in literature. Results show that the approach developed here performs reasonably well, and has thus improved the applicability of this class of wave transformation models.     

Long-term hydrological,biogeochemical, and climatological data from Walker Branch Watershed,East Tennessee,USA

In 1967, the original Walker Branch Watershed (WBW) project was established to study elemental cycling and mass balances in a relatively unimpacted watershed. Over the next 50+ years, findings from additional experimental studies and long-term observations on WBW advanced understanding of catchment hydrology, biogeochemistry, and ecology and established WBW as a seminal site for catchment science. The 97.5-ha WBW is located in East Tennessee, USA, on the U.S. Department of Energy's Oak Ridge Reservation. Vegetation on the watershed is characteristic of an eastern deciduous, second-growth forest. The watershed is divided into two subcatchments: the West Fork (38.4 ha) and the East Fork (59.1 ha). Headwater streams draining these subcatchments are fed by multiple springs, and thus flow is perennial. Stream water is high in base cations due to weathering of dolomite bedrock and nutrient concentrations are low. Long-term observations of climate, hydrology, and biogeochemistry include daily (1969–2014) and 15-min (1994–2014) stream discharge and annual runoff (1969–2014); hourly, daily, and annual rainfall (1969–2012); daily climate and soil temperature (1993–2010); and weekly stream water chemistry (1989–2013). These long-term datasets are publicly available on the WBW website (https://walkerbranch.ornl.gov/long-term-data/ ). While collection of these data has ceased, related long-term measurements continue through the National Ecological Observatory Network (NEON), where WBW is the core terrestrial and aquatic site in the Appalachian and Cumberland Plateau region (NEON's Domain 7) of the United States. These long-term datasets have been and will continue to be important in evaluating the influence of climatic and environmental drivers on catchment processes.     

One-dimensional consolidation of soil under multistage load based on continuous drainage boundary

In engineering practice, a rapid loading rate can result in ground failure when the strength of soft soils is relatively low, and a multistage loading scheme is always utilized to deal with this situation. Firstly, under a multistage load and the continuous drainage boundary, an analytical solution of excess pore-water pressure and consolidation degree is obtained by virtue of the superposition formula of excess pore-water pressure, and a more general continuous drainage boundary under arbitrary time-dependent load is developed. Then, a comparison with existing analytical solutions is conducted to verify the present solution. A preliminary attempt on applying the continuous drainage boundary into the finite element model is made, and the feasibility of the numerical model for the one-dimensional consolidation under the continuous drainage boundary is verified by comparing the results calculated by FEM with that from present analytical solution. Finally, the consolidation behavior of soil is investigated in detail for different int erface parameters or loading scheme. The results show that, in land reclamation projects, a horizontal drain should be placed close to the boundary with a smaller interface parameter to improve the consolidation efficiency. The degree of consolidation is also related to the applied time-dependent load and interface parameters.     

Modelling landslides in unsaturated slopes subjected to rainfall infiltration using material point method

This paper presents a dynamic fully coupled formulation for saturated and unsaturated soils that undergo large deformations based on material point method. Governing equations are applied to porous material while considering it as a continuum in which the pores of the solid skeleton are filled with water and air. The accuracy of the developed method is tested with available experimental and numerical results. The developed method has been applied to investigate the failure and post‐failure behaviour of rapid landslides in unsaturated slopes subjected to rainfall infiltration using two different bedrock geometries that lie below the top soil. The models show different failure and post‐failure mechanisms depending on the bedrock geometry and highlight the negative effects of continuous rain infiltrations. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of climate indicators on continental‐scale potential groundwater recharge in Africa

In the last decades, human activity has been contributing to climate change that is closely associated with an increase in temperatures, increase in evaporation, intensification of extreme dry and wet rainfall events, and widespread melting of snow and ice. Understanding the intricate linkage between climate warming and the hydrological cycle is crucial for sustainable management of groundwater resources, especially in a vulnerable continent like Africa. This study investigates the relationship between climate‐change drivers and potential groundwater recharge (PGR) patterns across Africa for a long‐term record (1960–2010). Water‐balance components were simulated by using the PCR‐GLOBWB model and were reproduced in both gridded maps and latitudinal trends that vary in space with minima on the Tropics and maxima around the Equator. Statistical correlations between temperature, storm occurrences, drought, and PGR were examined in six climatic regions of Africa. Surprisingly, different effects of climate‐change controls on PGR were detected as a function of latitude in the last three decades (1980–2010). Temporal trends observed in the Northern Hemisphere of Africa reveal that the increase in temperature is significantly correlated to the decline of PGR, especially in the Northern Equatorial Africa. The climate indicators considered in this study were unable to explain the alarming negative trend of PGR observed in the Sahelian region, even though the Standardized Precipitation‐Evapotranspiration Index (SPEI) values report a 15% drought stress. On the other hand, increases in temperature have not been detected in the Southern Hemisphere of Africa, where increasing frequency of storm occurrences determine a rise of PGR, particularly in southern Africa. Time analysis highlights a strong seasonality effect, while PGR is in‐phase with rainfall patterns in the summer (Northern Hemisphere) and winter (Southern Hemisphere) and out‐of‐phase during the fall season. This study helps to elucidate the mechanism of the processes influencing groundwater resources in six climatic zones of Africa, even though modelling results need to be validated more extensively with direct measurements in future studies. Copyright © 2016 John Wiley & Sons, Ltd.     

Rainfall–runoff responses on Arctic hillslopes underlain by continuous permafrost,North Slope,Alaska, USA

The Arctic hydrologic cycle is intensifying, as evidenced by increased rates of precipitation, evapotranspiration, and riverine discharge. However, the controls on water fluxes from terrestrial to aquatic systems in upland Arctic landscapes are poorly understood. Upland landscapes account for one third of the Arctic land surface and are often drained by zero‐order geomorphic flowpath features called water tracks. Previous work in the region attributed rapid runoff response at larger stream orders to water tracks, but models suggest water tracks are hydrologically disconnected from the surrounding hillslope. To better understand the role of water tracks in upland landscapes, we investigated the surface and subsurface hydrologic responses of 6 water tracks and their hillslope watersheds to natural patterns of rainfall, soil thaw, and drainage. Between storms, both water track discharge and the water table in the hillslope watersheds exhibited diel fluctuations that, when lagged by 5 hr, were temporally correlated with peak evapotranspiration rate. Water track soils remained saturated for more of the summer season than soils in their surrounding hillslope watersheds. When rainfall occurred, the subsurface response was nearly instantaneous, but the water tracks took significantly longer than the hillslopes to respond to rainfall, and longer than the responses previously observed in nearby larger order Arctic streams. There was also evidence for antecedent soil water storage conditions controlling the magnitude of runoff response. Based on these observations, we used a broken stick model to test the hypothesis that runoff production in response to individual storms was primarily controlled by rainfall amount and antecedent water storage conditions near the water track outlet. We found that the relative importance of the two factors varied by site, and that water tracks with similar watershed geometries and at similar landscape positions had similar rainfall–runoff model relationships. Thus, the response of terrestrial water fluxes in the upland Arctic to climate change depends on the non‐linear interactions between rainfall patterns and subsurface water storage capacity on hillslopes. Predicting these interactions across the landscape remains an important challenge.     

Stormflow concentration–discharge dynamics of suspended sediment and dissolved phosphorus in an agricultural watershed

Concentration–discharge (C-Q) relationships are an effective tool for identifying watershed biogeochemical source and transport dynamics over short and long timescales. We examined stormflow C-Q, hysteresis, and flushing patterns of total suspended sediment (TSS) and soluble reactive phosphorus (SRP) in two stream reaches of a severely impaired agricultural watershed in northeastern Wisconsin, USA. The upper watershed reach—draining a relatively flat, row crop-dominated contributing area—showed predominantly anti-clockwise TSS hysteresis during storms, suggesting that particulate materials were mobilized more from distal upland sources than near- and in-channel areas. In contrast, the incised lower watershed reach produced strong TSS flushing responses on the rising limb of storm hydrographs and clockwise hysteresis, signalling rapid mobilization of near- and in-channel materials with increasing event flows. C-Q relationships for SRP showed complex patterns in both the upper and lower reaches, demonstrating largely non-linear chemodynamic C-Q behaviour during events. As with TSS, anti-clockwise SRP hysteresis in the upper reach suggested a delay in the hydrologic connectivity between SRP sources and the stream, with highly variable SRP concentrations during some events. A broad range of clockwise, anti-clockwise, and complex SRP hysteresis patterns occurred in the lower watershed, possibly influenced by in-channel legacy P stores and connection to tile drainage networks in the lower watershed area. Total suspended sediment and SRP responses were also strongly related to precipitation event characteristics including antecedent precipitation, recovery period, and precipitation intensity, highlighting the complexity of stormflow sediment and phosphorus responses in this severely impaired agricultural stream.     

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.     

Etude isotopique de la pluie et de la neige sur le Mont Liban: premiers résultats / Isotope study of snow and rain on Mount Lebanon: preliminary results

Résumé

Résumé Quelques analyses isotopiques préliminaires ont été réalisées sur les précipitations pluvio-neigeuses, sur un profil de neige et sur deux sources karstiques sur le Mont Liban. Elles confirment la variabilité saisonnière du signal atmosphérique et en particulier que l’excès en deutérium est en relation avec l’origine des masses d’air et avec les recharges de vapeur sur la Méditerranée. Elles montrent également une relative stabilité du signal isotopique du couvert neigeux, peu ou pas influencé par les phénomènes de sublimation, d’évaporation ou de fonte/regel. La participation progressive de la fonte du manteau neigeux à l’alimentation des sources karstiques est qualitativement observée.     

U–Th–Pb monazite dating and the timing of arc–continent collision in East Timor

The Timor Orogen represents one of the youngest arc–continent collisions exposed on the Earth. It has the potential to provide some of the key parameters about how this style of orogeny evolves. However, the metamorphic age of the highest-grade rocks formed in the collision remains controversial. Using U–Th–Pb dating of monazite from amphibolite-grade sillimanite and garnet-bearing schists we show the peak metamorphism occurred at 5.5–4.7 Ma. The young age of the monazite and the presence of significant amounts of common Pb required the development of a new protocol to simultaneously account for the ²³⁰Th disequilibrium and the ²⁰⁷Pb common Pb correction. The new estimate of metamorphic age is consistent with the estimates based on plate reconstructions for the initiation of arc–continent collision in East Timor. The metamorphic event is a result of this collision.