中国北方大陆下的地幔热柱与岩石圈运动

本文首次提出中国北方大陆下存在一个地幔热柱的论证,并提出亚热柱(sub—plume)的新概念。热柱的中心与边缘部分的隙间熔浆分别为苦橄质玄武岩与碱性玄武岩。在渐新世到中新世约18.4 Ma内,北方大陆以3.26cm/a的速率向东南飘移了约600km,使日本海、渤海—华北平原等脱离热柱。导致晚第三纪日本海扩张的停止,渤海—华北平原等早第三纪火山喷发的突然中止。火山喷发期间,在热柱头部若干个亚热柱的形成,好似若干个“铆钉”穿入岩石圈,有效地阻止了岩石圈的飘移(这时的飘移速率只有0.05cm/a),我们把火山喷发称为固定岩石圈的“铆钉效应”。     

An efficient finite–discrete element method for quasi‐static nonlinear soil–structure interaction problems

An efficient finite–discrete element method applicable for the analysis of quasi‐static nonlinear soil–structure interaction problems involving large deformations in three‐dimensional space was presented in this paper. The present method differs from previous approaches in that the use of very fine mesh and small time steps was not needed to stabilize the calculation. The domain involving the large displacement was modeled using discrete elements, whereas the rest of the domain was modeled using finite elements. Forces acting on the discrete and finite elements were related by introducing interface elements at the boundary of the two domains. To improve the stability of the developed method, we used explicit time integration with different damping schemes applied to each domain to relax the system and to reach stability condition. With appropriate damping schemes, a relatively coarse finite element mesh can be used, resulting in significant savings in the computation time. The proposed algorithm was validated using three different benchmark problems, and the numerical results were compared with existing analytical and numerical solutions. The algorithm performance in solving practical soil–structure interaction problems was also investigated by simulating a large‐scale soft ground tunneling problem involving soil loss near an existing lining. Copyright © 2011 John Wiley & Sons, Ltd.     

Three dimensional elasto‐plastic interface for embedded beam elements with interaction surface for the analysis of lateral loading of piles

This paper presents a numerical formulation for a three dimensional elasto‐plastic interface, which can be coupled with an embedded beam element in order to model its non‐linear interaction with the surrounding solid medium. The formulation is herein implemented for lateral loading of piles but is able to represent soil‐pile interaction phenomena in a general manner for different types of loading conditions or ground movements. The interface is formulated in order to capture localized material plasticity in the soil surrounding the pile within the range of small to moderate lateral displacements. The interface is formulated following two different approaches: (i) in terms of beam degrees of freedoms; and (ii) considering the displacement field of the solid domain. Each of these alternatives has its own advantages and shortcomings, which are discussed in this paper. The paper presents a comparison of the results obtained by means of the present formulation and by other well‐established analysis methods and test results published in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of channel transmission losses in a dryland river reach in north‐eastern Brazil using streamflow series,groundwater level series and multi‐temporal satellite data

Scarcity of hydrological data, especially streamflow discharge and groundwater level series, restricts the understanding of channel transmission losses (TL) in drylands. Furthermore, the lack of information on spatial river dynamics encompasses high uncertainty on TL analysis in large rivers. The objective of this study was to combine the information from streamflow and groundwater level series with multi‐temporal satellite data to derive a hydrological concept of TL for a reach of the Middle Jaguaribe River (MJR) in semi‐arid north‐eastern Brazil. Based on this analysis, we proposed strategies for its modelling and simulation. TL take place in an alluvium, where river and groundwater can be considered to be hydraulically connected. Most losses certainly infiltrated only through streambed and levees and not through the flood plains, as could be shown by satellite image analysis. TL events whose input river flows were smaller than a threshold did not reach the outlet of the MJR. TL events whose input flows were higher than this threshold reached the outlet losing on average 30% of their input. During the dry seasons (DS) and at the beginning of rainy seasons (DS/BRS), no river flow is expected for pre‐events, and events have vertical infiltration into the alluvium. At the middle and the end of the rainy seasons (MRS/ERS), river flow sustained by base flow occurs before/after events, and lateral infiltration into the alluvium plays a major role. Thus, the MJR shifts from being a losing river at DS/BRS to become a losing/gaining (mostly losing) river at MRS/ERS. A model of this system has to include the coupling of river and groundwater flow processes linked by a leakage approach. Copyright © 2012 John Wiley & Sons, Ltd.     

The influence of local hydrogeologic forcings on near‐stream event water recharge and retention (Upper San Pedro River,Arizona)

The rise in stream stage during high flow events (floods) can induce losing stream conditions, even along stream reaches that are gaining during baseflow conditions. The aquifer response to flood events can affect the geochemical composition of both near‐stream groundwater and post‐event streamflow, but the amount and persistence of recharged floodwater may differ as a function of local hydrogeologic forcings. As a result, this study focuses on how vertical flood recharge varies under different hydrogeologic forcings and the significance that recharge processes can have on groundwater and streamflow composition after floods. River and shallow groundwater samples were collected along three reaches of the Upper San Pedro River (Arizona, USA) before, during and after the 2009 and 2010 summer monsoon seasons. Tracer data from these samples indicate that subsurface floodwater propagation and residence times are strongly controlled by the direction and magnitude of the dominant stream–aquifer gradient. A reach that is typically strongly gaining shows minimal floodwater retention shortly after large events, whereas the moderately gaining and losing reaches can retain recharged floodwater from smaller events for longer periods. The moderately gaining reach likely returned flood recharge to the river as flow declined. These results indicate that reach‐scale differences in hydrogeologic forcing can control (i) the amount of local flood recharge during events and (ii) the duration of its subsurface retention and possible return to the stream during low‐flow periods. Our observations also suggest that the presence of floodwater in year‐round baseflow is not due to long‐term storage beneath the streambed along predominantly gaining reaches, so three alternative mechanisms are suggested: (i) repeated flooding that drives lateral redistribution of previously recharged floodwater, (ii) vertical recharge on the floodplain during overbank flow events and (iii) temporal variability in the stream–aquifer gradient due to seasonally varying water demands of riparian vegetation. Copyright © 2011 John Wiley & Sons, Ltd.     

Sedimentation of overbank floods in the confined complex channel–floodplain system of the Lower Yellow River,China

The channel boundary conditions along the Lower Yellow River (LYR) have been altered significantly since the 1950s with the continual reinforcement and construction of both main and secondary dykes and river training works. To evaluate how the confined complex channel–floodplain system of the LYR responds to floods, this study presents a detailed investigation of the relationship between the tempo‐spatial distribution of sedimentation/erosion and overbank floods occurred in the LYR. For large overbank floods, we found that when the sediment transport coefficient (ratio of sediment concentration of flow to flow discharge) is less than 0.034, the bankfull channel is subject to significant erosion, whereas the main and secondary floodplains both accumulate sediment. The amount of sediment deposited on the main and secondary floodplains is closely related to the ratio of peak discharge to bankfull discharge, volume of water flowing over the floodplains, and sediment concentration of overbank flow, whereas the degree of erosion in the bankfull channel is related to the amount of sediment deposited on the main and secondary floodplains, water volume, and sediment load in flood season. The significant increase in erosion in the bankfull channel is due to the construction of the main and secondary dykes and river training works, which are largely in a wide and narrow alternated pattern along the LYR such that the water flowing over wider floodplains returns to the channel downstream after it drops sediment. For small overbank floods, the bankfull channel is subject to erosion when the sediment transport coefficient is less than 0.028, whereas the amount of sediment deposited on the secondary floodplain is associated closely with the sediment concentration of flow. Over the entire length of the LYR, the situation of erosion in the bankfull channel and sediment deposition on the main and secondary floodplains occurred mainly in the upper reach of the LYR, in which a channel wandering in planform has been well developed.     

The life cycle of the intense IOP‐14 storm during CASPII. Part I: Analysis and simulations

Abstract

In this study, a 5‐day life‐cycle of the IOP‐14 storm during CASP II is examined using conventional observations and numerical simulations with a mesoscale version of the Canadian Regional Finite‐Element (RFE) model. Observational analysis reveals that the IOP‐14 storm forms from a lee trough, occurring along a strong baroclinic zone with an intense frontogenetic deformation, that interacts with an upper‐level travelling short‐wave trough across the Canadian Rockies. Then the storm experiences a slow, but nearly steady, growth while traversing the North American continent. It deepens explosively as it moves into the Atlantic Ocean. It appears that i) the enhanced large‐scale baroclinicity due to land‐sea temperature contrasts, ii) the tremendous latent heat release due to the transport of high‐θ_e air from the marine boundary layer, Hi) the decrease of surface drag and iv) the favourable westward tilt of the low with an amplifying trough all contribute to the explosive deepening of the storm.

Two consecutive simulations covering a total of 102 h during the storm development are carried out with a grid size of 50 km. The RFE model reproduces very well the formation of the surface low on the lee side of the Rockies, the track and deepening rates, the explosive development and decay of the storm, and various mesoscale phenomena (e.g., a “bent‐back” warm front, a “T‐bone” thermal pattern, a cold frontal “fracture”, an upper‐level “eye” and warm‐core structures), as verified by conventional observations, satellite imagery, flight‐level and dropsonde data from a research aircraft. It is found from potential vorticity (PV) analysis that the storm reaches its peak intensity as the upper‐level dry PV anomaly, the low‐level moist PV anomaly and surface thermal warmth are vertically superposed. PV inversions reveal that these anomalies contribute about 60%, 30% and 10%, respectively, to the 900‐hPa negative height perturbation. It is shown that the warm‐core structure near the cyclone centre is produced by advection of warmer air ahead of the cold front, rather than by adiabatic warming associated with subsidence.     

On the testing of fully integrated surface–subsurface hydrological models

Studies employing integrated surface–subsurface hydrological models (ISSHMs) have utilized a variety of test cases to demonstrate model accuracy and consistency between codes. Here, we review the current state of ISSHM testing and evaluate the most popular ISSHM test cases by comparing the hydrodynamic processes simulated in each case to the processes found in well‐characterized, real‐world catchments and by comparing their general attributes to those of successful benchmark problems from other fields of hydrogeology. The review reveals that (1) ISSHM testing and intercode comparison have not adopted specific test cases consistently; (2) despite the wide range of ISSHM metrics available for model testing, only two model performance diagnostics are typically adopted: the catchment outflow hydrograph and the catchment water balance; (3) in intercode comparisons, model performance is usually judged by evaluating only one performance diagnostic: the catchment outflow hydrograph; and (4) ISSHM test cases evaluate a small number of hydrodynamic processes that are largely uniform across the model domain, representing a limited selection of the processes of interest in well‐characterized, real‐world catchments. ISSHM testing would benefit from more intercode comparisons using a consistent set of test cases, aimed at evaluating more catchment processes (e.g. flooding) and using a wider range of simulation diagnostics (e.g. pressure head distributions). To achieve this, a suite of test case variations is required to capture the relevant catchment processes. Finally, there is a need for additional ISSHM test problems that compare model predictions with hydrological observations from intensively monitored field sites and controlled laboratory experiments. Copyright © 2012 John Wiley & Sons, Ltd.     

The importance of wave-tide-surge interaction to coastal material transport and diffusion

This study is to combine a coastal high-resolution (2′×2′) two-way coupled wave-tide-surge numerical model (including 3 main physical mechanisms) with a material transport/diffusion model for understanding the law of material transport/diffusion. Results show that the law of material trans- port/diffusion driven by background current field simulated by the coupled wave-tide-surge model is dif- ferent from that simulated by pure tide-surge, and more different from traditional ones driven by tidal current. The coupled background current should be taken into account for the simulation.     

Response evaluation of axially loaded fixed‐head pile groups in clayey soils

The aim of this paper is to investigate the interaction between the piles in a group with a rigid head and correlate the response of a group of piles to that of a single pile. For this purpose, a computationally intensive study using 3‐D nonlinear numerical analysis was carried out for different pile group arrangements in clayey soils. The responses of the groups of piles were compared with that of a single pile and the variation of the settlement amplification factor R_a was then quantified. The influence of the number of piles, the spacing, and the settlement level on the group response is discussed. A previously proposed relationship for predicting the response of a pile group, based on its configuration and the response of a single pile, has been modified to extend its applicability for any pile spacing. The modified relationship provides a reasonable prediction for various group configurations in clayey soils. Copyright © 2009 John Wiley & Sons, Ltd.