The effects of rock fragments on hydrologic and hydraulic responses along a slope

A model describing the three‐dimensional matrix flow along a slope with rock fragments or impermeable blocks was developed. The model was combined with modified Picard's iteration to ensure mass conservation in the unsaturated flow. We found that rock fragments obstruct water flow along the slope. The groundwater table must be raised to provide a sufficient pore water pressure gradient to facilitate water flow, but higher pore water pressure may induce slope failure. We also conducted a bench‐scale laboratory flume experiment to examine the effects of impermeable blocks on downstream seepage flow. In addition, a numerical experiment was conducted to examine how different arrangements of impermeable blocks affect downstream seepage flow and pore water pressure. This research demonstrated that the hydraulic phenomena were affected when impermeable blocks were present, and pore water pressure increased as the position of impermeable blocks was lowered. Copyright © 2006 John Wiley & Sons, Ltd.     

Non‐local damage modelling of concrete: a procedure for the determination of model parameters

This paper presents a procedure for the determination of parameters of non‐local damage models. This is to assure a consistent response of a non‐local damage model, as choice of the internal length and other parameters of the model are varied. Correlations between the internal length and other parameters governing the local constitutive behaviour of the model are addressed and exploited. Focus is put on the relationship between the internal length of the non‐local model and the width of the fracture process zone. Numerical examples are used to demonstrate the rigour of the proposed method. Copyright © 2006 John Wiley & Sons, Ltd.     

Mobilization of arsenic and iron from Red River floodplain sediments, Vietnam

Sediments from the Red River and from an adjacent floodplain aquifer were investigated with respect to the speciation of Fe and As in the solid phase, to trace the diagenetic changes in the river sediment upon burial into young aquifers, and the related mechanisms of arsenic release to the groundwater. Goethite with subordinate amounts of hematite were, using Mössbauer spectroscopy, identified as the iron oxide minerals present in both types of sediment. The release kinetics of Fe, As, Mn and PO₄ from the sediment were investigated in leaching experiments with HCl and 10 mM ascorbic acid, both at pH 3. From the river sediments, most of the Fe and As was mobilized by reductive dissolution with ascorbic acid while HCl released very little Fe and As. This suggests As to be associated with an Fe-oxide phase. For oxidized aquifer sediment most Fe was mobilized by ascorbic acid but here not much As was released. However, the reduced aquifer sediments contained a large pool of Fe(II) and As that is readily leached by HCl, probably derived from an unidentified authigenic Fe(II)-containing mineral which incorporates As as well. Extraction with ascorbic acid indicates that the river sediments contain both As(V) and As(III), while the reduced aquifer sediment almost exclusively releases As(III). The difference in the amount of Fe(II) leached from river and oxidized aquifer sediments by ascorbic acid and HCl, was attributed to reductive dissolution of Fe(III). The reactivity of this pool of Fe(III) was quantified by a rate law and compared to that of synthetic iron oxides. In the river mud, Fe(III) had a reactivity close to that of ferrihydrite, while the river sand and oxidized aquifer sediment exhibited a reactivity ranging from lepidocrocite or poorly crystalline goethite to hematite. Mineralogy by itself appears to be a poor predictor of the iron oxide reactivity in natural samples using the reactivity of synthetic Fe-oxides as a reference. Sediments were incubated, both unamended and with acetate added, and monitored for up to 2 months. The river mud showed the fastest release of both Fe and As, while the effect of acetate addition was minor. This suggests that the presence of reactive organic carbon is not rate limiting. In the case of the river and aquifer sediments, the release of Fe and As was always stimulated by acetate addition and here reactive organic carbon was clearly the rate limiting factor. The reduced aquifer sediment apparently can sustain slower but prolonged microbially-driven release of As. The highly reactive pools of Fe(III) and As in the river mud could be due to reoxidation of As and Fe contained in the reducing groundwater from the floodplain aquifers that are discharging into the river. Deposition of the suspended mud on the floodplain during high river stages is proposed to be a major flux of As onto the floodplain and into the underlying aquifers.     

Geochemistry and petrogenesis of Permian ultramafic-mafic complexes of the Jinping–Song Da rift (southeastern Asia)

The Jinping–Song Da rift structure in the Emeishan Large Igneous Province is composed of Permian high- and low-Ti volcanic and volcanoplutonic ultramafic-mafic associations of different compositions and genesis. High-Ti picrites, picrobasalts, basalts, and dolerites are enriched in LREE and depleted in HREE and show low Al₂O₃/TiO₂ ratios (~4), commensurate ε_Nd(T) values (+0.5 to +1.1), and low (Th/Nb)_PM ratios similar to those of OIB-enriched mantle source. The established geochemical characteristics evidence that the parental melts of these rocks were generated from garnet lherzolite at the depths of garnet stability (~200 to 400 km). Later, high-Mg low-Ti volcanics (komatiites, komatiitic basalts, and basalts) and associating small peridotite-gabbro massifs and komatiite-basalt dikes were produced as a result of ~20% partial melting of depleted water-poor (≤0.03 wt.% H₂O) peridotite substratum from the hottest upper part of mantle plume at relatively shallow depths (100–120 km). The LREE-depleted komatiites and komatiitic basalts are characterized by low (Ce/Yb)_CH values, ¹⁸⁷Re/¹⁸⁸Os = 0.05–1.2, ⁸⁷Sr/⁸⁶Sr = 0.704–0.706, positive ε_Nd(T) values (+3 to +8), γ_Os = –0.5 to +0.9, and strong negative anomalies of Ba, K, and Sr on the spidergrams. The scarcer LREE-enriched komatiites, komatiitic basalts, and basalts vary greatly in chemical composition and values of ε_Nd(t) (+6.4 to –10.2), ⁸⁷Sr/⁸⁶Sr (0.706–0.712), and γ_Os (+14.8 to +56), which is due to the different degrees of crustal contamination of parental magmas. The Rb-Sr isotopic age of basaltic komatiite is 257 ± 24 Ma. The Re-Os age determined by analysis of 12 komatiite samples is 270 ± 21 Ma. These data agree with the age of flood basalts of the Emeishan Large Igneous Province. The komatiite-basalt complex of the Song Da rift is still the only Phanerozoic PGE-Cu-Ni-complex of this composition. The geochemistry of accompanying Cu-Ni-PGE-ores confirms their relationship with komatiite-basaltic magmatism.     

Monitoring severe aquifer-system compaction and land subsidence in Taiwan using multiple sensors: Yunlin,the southern Choushui River Alluvial Fan

During 1992–2007, excessive pumping of groundwater caused large-scale aquifer-system compaction and land subsidence in the Choshui River Alluvial Fan, especially in the area of Yunlin county. The subsidence impedes surface-water runoff and endangers the operation of Taiwan High Speed Rail. Leveling, Global Positioning System (GPS), multi-level compaction monitoring well, and Differential Interferometric Synthetic Aperture Radar (DInSAR) are used to study the extent of subsidence in Yunlin and its mechanism. These sensors complement each other in spatial and temporal resolutions. A leveling network totaling 434 km in length was deployed to derive subsidence at every 1.5 km along the routes, and the result is accurate to few mm and shows a basin-like subsidence pattern centering at Tuku Township. Four multi-level compaction monitoring wells, co-located with GPS pillars, detect compactions at different depths, showing that the aquifer-system compaction (the cause of subsidence) occurs mostly below depths >200 m, where reduction of groundwater pumping is most needed. The vertical displacements from GPS and leveling agree to within 1 cm, and are larger than the cumulative compaction detected by the compaction-monitoring wells, suggesting that compaction also occurs below 300 m (the depth of the wells). The vertical displacements derived using DInSAR and 8 ENVISAT SAR images agree with the leveling result to 1–2 cm.     

New Methodology to Characterize Shear Behavior of Joints by Combination of Direct Shear Box Testing and Numerical Simulations

A new procedure is presented, which combines big shear box tests on rocks and corresponding numerical simulations with explicit consideration of joint roughness to get deeper insight into the shear behavior of rock joints. The procedure consists of three parts: (1) constant normal load- or CNS-shear box tests with registration of shear- and normal-components of stress and displacements and deduction of basis rock mechanical parameters; (2) high resolution 3D-scanning of joint surface to deduce joint topography; and (3) set-up, run and evaluation of 3-dimensional numerical model with explicit duplication of joint roughness as back-analysis of shear box tests. The numerical back-analysis provides deeper insight into the joint behavior at the micro-scale. Several parameters can be deduced, like micro-slope angle distribution, aperture size distribution, local normal stress distribution and detailed analysis of dilation in relation to shear direction. The potential of the new procedure is illustrated exemplary by shear box tests on slate.     

Statistical data fusion of multi-sensor AOD over the Continental United States

This article illustrates two techniques for merging daily aerosol optical depth (AOD) measurements from satellite and ground-based data sources to achieve optimal data quality and spatial coverage. The first technique is a traditional Universal Kriging (UK) approach employed to predict AOD from multi-sensor aerosol products that are aggregated on a reference grid with AERONET as ground truth. The second technique is spatial statistical data fusion (SSDF); a method designed for massive satellite data interpolation. Traditional kriging has computational complexity O(N³), making it impractical for large datasets. Our version of UK accommodates massive data inputs by performing kriging locally, while SSDF accommodates massive data inputs by modelling their covariance structure with a low-rank linear model. In this study, we use aerosol data products from two satellite instruments: the moderate resolution imaging spectrometer and the geostationary operational environmental satellite, covering the Continental United States.     

Floodplain hydrology of the Mekong Delta,Vietnam

The Mekong Delta is one of the largest and most intensively used estuaries in the world. Each year it witnesses widespread flooding which is both the basis of the livelihood for more than 17 million people but also the major hazard. Therefore, a thorough understanding of the hydrologic and hydraulic features is urgently required for various planning purposes. While the general causes and characteristics of the annual floods are understood, the inundation dynamics in the floodplains in Vietnam which are highly controlled by dikes and other control structures have not been investigated in depth. Especially, quantitative analyses are lacking, mainly due to scarce data about the inundation processes in the floodplains. Therefore, a comprehensive monitoring scheme for channel and floodplain inundation was established in a study area in the Plain of Reeds in the northeastern part of the Vietnamese Delta. This in situ data collection was complemented by a series of high‐resolution inundation maps derived from the TerraSAR‐X satellite for the flood seasons 2008 and 2009. Hence, the inundation dynamics in the channels and floodplains, and the interaction between channels and floodplains, could be quantified for the first time. The study identifies the strong human interference which is governed by flood protection levels, cropping patterns and communal water management. In addition, we examine the tidal influence on the inundation in various parts of the Delta, since it is expected that climate change‐induced sea level rise will increase the tidal contribution to floodplain inundation. Copyright © 2011 John Wiley & Sons, Ltd.     

Ensemble rainfall forecasting with numerical weather prediction and radar‐based nowcasting models

The overall objective of this study is to improve the forecasting accuracy of the precipitation in the Singapore region by means of both rainfall forecasting and nowcasting. Numerical Weather Predication (NWP) and radar‐based rainfall nowcasting are two important sources for quantitative precipitation forecast. In this paper, an attempt to combine rainfall prediction from a high‐resolution mesoscale weather model and a radar‐based rainfall model was performed. Two rainfall forecasting methods were selected and examined: (i) the weather research and forecasting model (WRF); and (ii) a translation model (TM). The WRF model, at a high spatial resolution, was run over the domain of interest using the Global Forecast System data as initializing fields. Some heavy rainfall events were selected from data record and used to test the forecast capability of WRF and TM. Results obtained from TM and WRF were then combined together to form an ensemble rainfall forecasting model, by assigning weights of 0.7 and 0.3 weights to TM and WRF, respectively. This paper presented results from WRF and TM, and the resulting ensemble rainfall forecasting; comparisons with station data were conducted as well. It was shown that results from WRF are very useful as advisory of anticipated heavy rainfall events, whereas those from TM, which used information of rain cells already appearing on the radar screen, were more accurate for rainfall nowcasting as expected. The ensemble rainfall forecasting compares reasonably well with the station observation data. Copyright © 2012 John Wiley & Sons, Ltd.     

Simultaneous identification of parameter,initial condition,and boundary condition in groundwater modelling

In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd.