Parallelization of the distinct lattice spring model

The distinct lattice spring model (DLSM) is a newly developed numerical tool for modeling rock dynamics problems, i.e. dynamic failure and wave propagation. In this paper, parallelization of DLSM is presented. With the development of parallel computing technologies in both hardware and software, parallelization of a code is becoming easier than before. There are many available choices now. In this paper, Open Multi‐Processing (OpenMP) with multicore personal computer (PC) and message passing interface (MPI) with cluster are selected as the environments to parallelize DLSM. Performances of these parallel DLSM codes are tested on different computers. It is found that the parallel DLSM code with OpenMP can reach a maximum speed‐up of 4.68× on a quad‐core PC. The parallel DLSM code with MPI can achieve a speed‐up of 40.886× when 256 CPUs are used on a cluster. At the end of this paper, a high‐resolution model with four million particles, which is too big to handle by the serial code, is simulated by using the parallel DLSM code on a cluster. It is concluded that the parallelization of DLSM is successful. Copyright © 2011 John Wiley & Sons, Ltd.     

The response of vegetation at the end of the last glacial period (MIS 3 and MIS 2) in littoral areas of NW Iberia

This paper describes and discusses the palaeobotanical data obtained from organic levels of two exposed deposits on the Atlantic shore of the northwestern Iberian Peninsula. Radiocarbon dating assigns these levels to a period of marine regression at the end of MIS 3 and the beginning of the Last Glacial Maximum. The pollen record shows an initial predominance of tree taxa (mainly deciduous, including the presence of Fagus far from its current limit), followed by an episode of partial forest retraction related to the end of MIS 3 and the beginning of the Last Glacial Maximum. However, the presence of numerous tree taxa in the record, even during the cold intervals associated with Heinrich events, points to the existence of sheltered refugia for these species during this period, in keeping with the conclusions of recent reviews of the vegetation dynamics of this region for the Lateglacial and the Holocene.     

Palaeoceanographic changes in the North Atlantic during the Mid‐Pleistocene Transition (MIS 31–19) as inferred from planktonic foraminiferal and calcium carbonate records

Marine sediments from the Integrated Ocean Drilling Project (IODP) Site U1314 (56.36°N, 27.88°W), in the subpolar North Atlantic, were studied for their planktonic foraminifera, calcium carbonate content, and Neogloboqudrina pachyderma sinistral (sin.) δ¹³C records in order to reconstruct surface and intermediate conditions in this region during the Mid‐Pleistocene Transition (MPT). Variations in the palaeoceanography and regional dynamics of the Arctic Front were estimated by comparing CaCO₃ content, planktonic foraminiferal species abundances, carbon isotopes and ice‐rafted debris (IRD) data from Site U1314 with published data from other North Atlantic sites. Site U1314 exhibited high abundances of the polar planktonic foraminifera N. pachyderma sin. and low CaCO₃ content until Marine Isotope Stage (MIS) 26, indicating a relatively southeastward position of the Arctic Front (AF) and penetration of colder and low‐salinity surface arctic water‐masses. Changing conditions after MIS 25, with oscillations in the position of the AF, caused an increase in the northward export of the warmer North Atlantic Current (NAC), indicated by greater abundances of non‐polar planktonic foraminifera and higher CaCO₃. The N. pachyderma sin. δ¹³C data indicate good ventilation of the upper part of the intermediate water layer in the eastern North Atlantic during both glacial and interglacial stages, except during Terminations 24/23, 22/21 and 20/1. In addition, for N. pachyderma (sin.) we distinguished two morphotypes: non‐encrusted and heavily encrusted test. Results indicate that increases in the encrusted morphotype and lower planktonic foraminiferal diversity are related to the intensification of glacial conditions (lower sea‐surface temperatures, sea‐ice formation) during MIS 22 and 20.     

One‐dimensional consolidation with a threshold gradient: a Stefan problem with rate‐dependent latent heat

During the process of one‐dimensional consolidation with a threshold gradient, the seepage front moves downward gradually, and the problem is indicated as a Stefan problem. The novel feature in this Stefan problem is a latent heat that varies inversely with the rate of the moving boundary. An exact solution for the external load that increases in proportion to the square root of time is constructed using the similarity transformation technique. Computational examples concerning the effect of different parameters on the motion of the seepage front are presented. The exact solution provides a worthwhile benchmark for verifying the accuracy of numerical and approximate methods. Copyright © 2013 John Wiley & Sons, Ltd.     

Internal erosion of chemically reinforced granular materials: a mathematical modeling approach

Internal erosion (IE) affects the stability of natural and reinforced materials by causing instability within their granular structure. The dislodgement and transport of eroded particles affect both the particulate concentration of eroding fluid and the pore network of eroded material. In this study, we examined these modifications using a transport model with a finite element code. First, IE tests on chemically reinforced sand columns were performed to obtain information about eroded material loss of mass, particulate concentration of effluent, porosity and permeability modifications, and existing IE stages. Second, based on experimental results, a mathematical one‐dimensional model has been formulated to monitor the evolution and spatial distribution of erodible solids, fluidized particles, porosity, permeability, and seepage stresses. The model consists of a set of coupled nonlinear differential equations solved in sequence. It provides valuable information about the extent and the dynamics of structural changes, which can be used to estimate an IE time for the hydraulic work to reach failure. Copyright © 2011 John Wiley & Sons, Ltd.     

Characterizing the hydraulic properties of unsaturated sand considering various porosities and drying‐wetting paths of the retention curve

Previous studies have shown that water retention curve (WRC) and the hydraulic conductivity vary because of changes of the void ratio or porosity of soil. However, limited documents pointed out the change of hydraulic properties of soil when compacted to different porosities while considering both of the drying and wetting processes of the WRC. This information is sometimes necessary for research like finger flow analysis or the occurrence of wetting and drying cycles as what would be seen in the field. Therefore, this study aims to examine the change of WRC characteristics with varied porosity considering both of the drying and wetting path in WRC by conducting a sand box experiment. Results show that the same type of sand compacted to various porosities have different hydraulic parameters. Hydraulic conductivities generally decrease with reduced porosities; shape parameter α of the van Genuchten equation (1980) linearly decreases with declining porosity and shape parameter n in a reversal manner for the sands of interest whether in the drying process or wetting process. The unsaturated properties of sand are further characterized by inspecting the variations of moisture content, matric suction and vertical displacement of soil body subject to periodic changes of the water level by another sand box experiment. The outcomes suggest that the saturated water content and residual water content are changing during the wetting–drying process, which can be an implication of the changed properties of WRC. The characteristics of volumetric deformation might be varied as well because of the observation of the dissimilar patterns of the changing vertical displacements among each wetting–drying process. Infiltration patterns of the sands also are identified through numerical modelling by introducing a constant infiltration flux from the surface followed by a no‐influx condition. Results indicate that less water accumulates in the sand near the surface for the sand compacted to higher porosity, but water can move deeper. Hydraulic conductivity is found as the prime factor dominating the evolvement of wetting fronts. However, shape parameters of water retention curves also affect the infiltration pattern to some extent. In addition, different sands with similar porosities can have quite different infiltrating characteristics. Copyright © 2012 John Wiley & Sons, Ltd.     

Abrupt behaviours of streamflow and sediment load variations of the Yangtze River basin,China

Monthly sediment load and streamflow series spanning 1963–2004 from four hydrological stations situation in the main stem of the Yangtze River, China, are analysed using scanning t‐test and the simple two‐phase linear regression scheme. Results indicate significant changes in the sediment load and streamflow from the upper reach to the lower reach of the Yangtze River. Relatively consistent positive coherency relations can be detected between streamflow and sediment load in the upper reach and negative coherency in the middle and lower reaches. Interestingly, negative coherency is found mainly for larger time scales. Changes in sediment load are the result mainly of human influence; specifically, the construction of water reservoirs may be the major cause of negative coherency. Accentuating the human influence from the upper to the lower reach results in inconsistent correlations between sediment load and streamflow. Decreasing sediment load being observed in recent years has the potential to alter the topographical properties of the river channel and the consequent development and recession of the Yangtze Delta. Results of this study are of practical significance for river channel management and evaluation of the influence of human activities on the hydrological regimes of large rivers. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling combined open channel flow by artificial neural networks

Combined open channel flow is encountered in many hydraulic engineering structures and processes, such as irrigation ditches and wastewater treatment facilities. Extensive experimental studies have conducted to investigate combined flow characteristics. Nevertheless, there is no simple relationship that can fully describe the velocity profiles in a turbulent flow. The artificial neural network (ANN) has great computational capability for solving various complex problems, such as function approximation. The main objective of this study is to evaluate the applicability of the ANN for simulating velocity profiles, velocity contours and estimating the discharges accordingly. The velocity profiles measured by an acoustic doppler velocimeter in the open channel of the Chihtan purification plant, Taipei, with different discharges at fixed measuring section and different depths are presented. The total number of data sets is 640 and the data sets are split into two subsets, i.e. training and validation sets. The backpropagation algorithm is used to construct the neural network. The results demonstrate that the velocity profiles can be modelled by the ANN, and the ANN constructed can nicely fit the velocity profiles and can precisely predict the discharges for the conditions investigated. Copyright © 2005 John Wiley & Sons, Ltd.     

A comparison of fibre‐optic distributed temperature sensing to traditional methods of evaluating groundwater inflow to streams

There are several methods for determining the spatial distribution and magnitude of groundwater inputs to streams. We compared the results of conventional methods [dye dilution gauging, acoustic Doppler velocimeter (ADV) differential gauging, and geochemical end‐member mixing] to distributed temperature sensing (DTS) using a fibre‐optic cable installed along 900 m of Ninemile Creek in Syracuse, New York, USA, during low‐flow conditions (discharge of 1·4 m³ s^?1). With the exception of differential gauging, all methods identified a focused, contaminated groundwater inflow and produced similar groundwater discharge estimates for that point, with a mean of 66·8 l s^?1 between all methods although the precision of these estimates varied. ADV discharge measurement accuracy was reduced by non‐ideal conditions and failed to identify, much less quantify, the modest groundwater input, which was only 5% of total stream flow. These results indicate ambient tracers, such as heat and geochemical mixing, can yield spatially and quantitatively refined estimates of relatively modest groundwater inflow even in large rivers. DTS heat tracing, in particular, provided the finest spatial characterization of groundwater inflow, and may be more universally applicable than geochemical methods, for which a distinct and consistent groundwater end member may be more difficult to identify. Copyright © 2011 John Wiley & Sons, Ltd.     

Garnet crystal plasticity in the continental crust,new example from south Madagascar

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.