Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar

Wetlands represent one of the world's most biodiverse and threatened ecosystem types and were diminished globally by about two‐thirds in the 20th century. There is continuing decline in wetland quantity and function due to infilling and other human activities. In addition, with climate change, warmer temperatures and changes in precipitation and evapotranspiration are reducing wetland surface and groundwater supplies, further altering wetland hydrology and vegetation. There is a need to automate inventory and monitoring of wetlands, and as a study system, we investigated the Shepard Slough wetlands complex, which includes numerous wetlands in urban, suburban, and agricultural zones in the prairie pothole region of southern Alberta, Canada. Here, wetlands are generally confined to depressions in the undulating terrain, challenging wetlands inventory and monitoring. This study applied threshold and frequency analysis routines for high‐resolution, single‐polarization (HH) RADARSAT‐2, synthetic aperture radar mapping. This enabled a growing season surface water extent hyroperiod‐based wetland classification, which can support water and wetland resource monitoring. This 3‐year study demonstrated synthetic aperture radar‐derived multitemporal open‐water masks provided an effective index of wetland permanence class, with overall accuracies of 89% to 95% compared with optical validation data, and RMSE between 0.2 and 0.7 m between model and field validation data. This allowed for characterizing the distribution and dynamics of 4 marsh wetlands hydroperiod classes, temporary, seasonal, semipermanent, and permanent, and mapping of the sequential vegetation bands that included emergent, obligate wetland, facultative wetland, and upland plant communities. Hydroperiod variation and surface water extent were found to be influenced by short‐term rainfall events in both wet and dry years. Seasonal hydroperiods in wetlands were particularly variable if there was a decrease in the temporary or semipermanent hydroperiod classes. In years with extreme rain events, the temporary wetlands especially increased relative to longer lasting wetlands (84% in 2015 with significant rainfall events, compared with 42% otherwise).     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Study of the flood control scheduling scheme for the Three Gorges Reservoir in a catastrophic flood

The Three Gorges Project is the world's largest water conservancy project. According to the design standards for the 1,000‐year flood, flood diversion areas in the Jingjiang reach of the Yangtze River must be utilized to ensure the safety of the Jingjiang area and the city of Wuhan. However, once these areas are used, the economic and life loss in these areas may be very great. Therefore, it is vital to reduce this loss by developing a scheme that reduces the use of the flood diversion areas through flood regulation by the Three Gorges Reservoir (TGR), under the premise of ensuring the safety of the Three Gorges Dam. For a 1,000‐year flood on the basis of a highly destructive flood in 1954, this paper evaluates scheduling schemes in which flood diversion areas are or are not used. The schemes are simulated based on 2.5‐m resolution reservoir topography and an optimized model of dynamic capacity flood regulation. The simulation results show the following. (a) In accord with the normal flood‐control regulation discharge, the maximum water level above the dam should be not more than 175 m, which ensures the safety of the dam and reservoir area. However, it is necessary to utilize the flood diversion areas within the Jingjiang area, and flood discharge can reach 2.81 billion m³. (b) In the case of relying on the TGR to impound floodwaters independently rather than using the flood diversion areas, the maximum water level above the dam reaches 177.35 m, which is less than the flood check level of 180.4 m to ensure the safety of the Three Gorges Dam. The average increase of the TGR water level in the Chongqing area is not more than 0.11 m, which indicates no significant effect on the upstream reservoir area. Comparing the various scheduling schemes, when the flood diversion areas are not used, it is believed that the TGR can execute safe flood control for a 1,000‐year flood, thereby greatly reducing flood damage.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Petrochronology of oxidized granulites from southern Peru

Sapphirine–quartz granulites from the Cocachacra region of the Arequipa Massif in southern Peru record early Neoproterozoic ultrahigh‐temperature metamorphism. Phase equilibrium modelling and zircon petrochronology are used to quantify timing and pressure–temperature (P–T) conditions of metamorphism. Modelling of three magnetite‐bearing sapphirine–quartz samples indicates peak temperatures of >950°C at ~0.7 GPa and a clockwise P–T evolution. Elevated concentrations of Al in orthopyroxene are also consistent with ultrahigh‐temperature conditions. Neoblastic zircon records ages of c. 1.0–0.9 Ga that are interpreted to record protracted ultrahigh‐temperature metamorphism. Th/U ratios of zircon of up to 100 reflect U‐depleted whole‐rock compositions. Concentrations of heavy rare earth elements in zircon do not show systematic trends with U–Pb age but do correlate with variable whole‐rock compositions. Very large positive Ce anomalies in zircon from two samples probably relate to strongly oxidizing conditions during neoblastic zircon crystallization. Low concentrations of Ti‐in‐zircon (<10 ppm) are interpreted to result from reduced titania activities due to the strongly oxidized nature of the granulites and the sequestration of titanium‐rich minerals away from the reaction volume. Whole‐rock compositions and oxidation state have a strong influence on the trace element composition of metamorphic zircon, which has implications for interpreting the geological significance of ages retrieved from zircon in oxidized metamorphic rocks.     

Determination of the regional deformation rates of Shanghai and Kashima VLBI stations based on ITRF97

The vertical deformation rates (VDRs) and horizontal deformation rates (HDRs) of Shanghai VLBI station in China and Kashima and Kashima34 VLBI stations in Japan were re-analysed using the baseline length change rates from Shanghai to 13 global VLBI stations, and from Kashima to 27 stations and from Kashima34 to 12 stations, based on the NASA VLBI global solution glb1123 (Ma, 1999). The velocity vectors of the global VLBI stations were referred to the ITRF97 reference frame, and the Eulerian vectors of different models of plate motion were used for comparative solutions. The VDR of Shanghai station is estimated to be −1.91±0.56 mm/yr, and those of Kashima and Kashima34 stations, −3.72±0.74mm/yr and −8.81±0.84mm/yr, respectively. The difference between the last two was verified by further analysis. Similar estimates were also made for the Kokee, Kauai and MK_VLBA VLBI stations in mid-Pacific.     

On the stability of particular solutions of the singly averaged Hill problem

We consider the particular solutions of the evolutionary system of equations in elements that correspond to planar and spatial circular orbits of the singly averaged Hill problem. We analyze the stability of planar and spatial circular orbits to inclination and eccentricity, respectively. We construct the instability regions of both particular solutions in the plane of parameters of the problem.     

冬半年青藏高原低涡发生发展的天气学分析

本文针对冬半年高原低涡几个不同源地的低涡生成背景,着重分析500hpa环流场,指出:地形对低涡生成是非常重要的条件,但一定范围内有利的环境流场也是低涡生成不可缺少的外部条件。并给出了不同涡源区有利于低涡发生,发展的500hpa环流形势。