Anticipating the effects of groundwater withdrawal on seawater intrusion and soil settlement in urban coastal areas

Intensive pumping in urban coastal areas is a common threat to water resource quality due to seawater intrusion. In those areas where subsurface water resources are not usually used for human consumption or irrigation, intensive pumping is associated with other activities like the lowering of the water table necessary to support underground structures and building foundations. This activity also increases the likelihood of soil settlement that affects building stability and the corrosion of concrete structures due to groundwater salinity. Under these circumstances, the awareness of a certain municipality (Calonge, NE Spain) of the potential effects of groundwater withdrawal upon foundations has led to an integrated approach to anticipate seawater intrusion related to urban development. Geological mapping and correlation of borehole logs, electrical resistivity tomography, and hydrochemical data provide comprehensive knowledge of the geology and hydrogeology of the area and act as screening tools necessary to discern the influence of hydrological processes in coastal areas. Developing Strack's analytical solution, new comprehensive, dimensionless expressions are herein derived to determine the critical pumping rate necessary to prevent seawater intrusion, as well as to reproduce the evolution of the wedge toe and the water table stagnation point under different withdrawal rates. Furthermore, the Dupuit–Forchheimer well discharge formula allows the estimation of the effects of the water table lowering due to such critical pumping in the surrounding building foundations. Field data from the Calonge coastal plain illustrate this approach and provide assessment criteria for future urban development and planning. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Constraining sequence stratigraphy in north Australian basins: SHRIMP U–Pb zircon geochronology between Mt Isa and McArthur River

Fifty‐five new SHRIMP U–Pb zircon ages from samples of northern Australian ‘basement’ and its overlying Proterozoic successions are used to refine and, in places, significantly change previous lithostratigraphic correlations. In conjunction with sequence‐stratigraphic studies, the 1800–1580 Ma rock record between Mt Isa and the Roper River is now classified into three superbasin phases—the Leichhardt, Calvert and Isa. These three major depositional episodes are separated by ～20 million years gaps. The Isa Superbasin can be further subdivided into seven supersequences each 10–15 million years in duration. Gaps in the geological record between these supersequences are variable; they approach several million years in basin‐margin positions, but are much smaller in the depocentres. Arguments based on field setting, petrography, zircon morphology, and U–Pb systematics are used to interpret these U–Pb zircon ages and in most cases to demonstrate that the ages obtained are depositional. In some instances, zircon crystals are reworked and give maximum depositional ages. These give useful provenance information as they fingerprint the source(s) of basin fill. Six new ‘Barramundi’ basement ages (around 1850 Ma) were obtained from crystalline units in the Murphy Inlier (Nicholson Granite and Cliffdale Volcanics), the Urapunga Tectonic Ridge (‘Mt Reid Volcanics’ and ‘Urapunga Granite’), and the central McArthur Basin (Scrutton Volcanics). New ages were also obtained from units assigned to the Calvert Superbasin (ca 1740–1690 Ma). SHRIMP results show that the Wollogorang Formation is not one continuous unit, but two different sequences, one deposited around 1730 Ma and a younger unit deposited around 1722 Ma. Further documentation is given of a regional 1725 Ma felsic event adjacent to the Murphy Inlier (Peters Creek Volcanics and Packsaddle Microgranite) and in the Carrara Range. A younger ca 1710 Ma felsic event is indicated in the southwestern McArthur Basin (Tanumbirini Rhyolite and overlying Nyanantu Formation). Four of the seven supersequences in the Isa Superbasin (ca 1670–1580 Ma) are reasonably well‐constrained by the new SHRIMP results: the Gun Supersequence (ca 1670–1655 Ma) by Paradise Creek Formation, Moondarra Siltstone, Breakaway Shale and Urquhart Shale ages grouped between 1668 and 1652 Ma; the Loretta Supersequence (ca 1655–1645 Ma) by results from the Lady Loretta Formation, Walford Dolomite, the upper part of the Mallapunyah Formation and the Tatoola Sandstone between ca 1653 and 1647 Ma; the River Supersequence (ca 1645–1630 Ma) by ages from the Teena Dolomite, Mt Les and Riversleigh Siltstones, and Barney Creek, Lynott, St Vidgeon and Nagi Formations clustering around 1640 Ma; and the Term Supersequence (ca 1630–1615 Ma) by ages from the Stretton Sandstone, lower Doomadgee Formation and lower part of the Lawn Hill Formation, mostly around 1630–1620 Ma. The next two younger supersequences are less well‐constrained geochronologically, but comprise the Lawn Supersequence (ca 1615–1600 Ma) with ages from the lower Balbirini Dolomite, and lower Doomadgee, Amos and middle Lawn Hill Formations, clustered around 1615–1610 Ma; and the Wide Supersequence (ca 1600–1585 Ma) with only two ages around 1590 Ma, one from the upper Balbirini Dolomite and the other from the upper Lawn Hill Formation. The Doom Supersequence (<1585 Ma) at the top of the Isa Superbasin is essentially unconstrained. The integration of high‐precision SHRIMP dating from continuously analysed stratigraphic sections, within a sequence stratigraphic context, provides an enhanced chronostratigraphic framework leading to more reliable interpretations of basin architecture and evolution.     

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.