Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.     

Accelerating the convergence of coupled geomechanical‐reservoir simulations

The pressure variations during the production of petroleum reservoir induce stress changes in and around the reservoir. Such changes of the stress state can induce marked deformation of geological structures for stress sensitive reservoirs as chalk or unconsolidated sand reservoirs. The compaction of those reservoirs during depletion affects the pressure field and so the reservoir productivity. Therefore, the evaluation of the geomechanical effects requires to solve in a coupling way the geomechanical problem and the reservoir multiphase fluid flow problem. In this paper, we formulate the coupled geomechanical‐reservoir problem as a non‐linear fixed point problem and improve the resolution of the coupling problem by comparing in terms of robustness and convergence different algorithms. We study two accelerated algorithms which are much more robust and faster than the conventional staggered algorithm and we conclude that they should be used for the iterative resolution of coupled reservoir‐geomechanical problem. Copyright © 2006 John Wiley & Sons, Ltd.     

Water source characterization through spatiotemporal patterns of major,minor and trace element stream concentrations in a complex,mesoscale German catchment

The link between spatiotemporal patterns of stream water chemistry and catchment characteristics for the mesoscale Dill catchment (692 km²) in Germany is explored to assess the catchment scale controls on water quality and to characterize water sources. In order to record the spatiotemporal pattern, ‘snapshot sampling’ was applied during low, mean and high flow, including 73 nested sites throughout the catchment. Water samples were analysed for the elements Li, B, Na, Mg, Al, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Rb, Sr, Mo, Ba, Pb and U using inductively‐coupled‐plasma mass spectrometry, and for electric conductivity and pH. Principle component analysis and hierarchical cluster analysis were used to find typical element associations and to group water samples according to their hydrochemical fingerprints. This revealed regional hydrochemical patterns of water quality which were subsequently related to catchment attributes to draw conclusions about the controls on stream chemistry. It was found that various lithologic signals and anthropogenic point source inputs controlled the base flow hydrochemistry. During increased flows, stream waters were diluted causing additional hydrochemical variability in response to heterogeneous precipitation inputs and differences in aquifer storage capacities. The hydrochemical patterns further displayed in‐stream mixing of waters. This implied, that stream waters could be apportioned to the identified water sources throughout the catchment. The basin‐wide hydrochemical variability has the potential to outrange the tracer signatures typically inferred in studies at the hillslope scale and is able to strongly influence the complexity of the catchment output. Both have to be considered for further catchment scale tracer and modelling work. Despite the likelihood of non‐conservative behaviour, the minor and trace elements enhanced the rather qualitative discrimination of the various groundwater types, as the major cations were strongly masked by point source inputs. Copyright © 2007 John Wiley & Sons, Ltd.     

Storm‐event rainfall–runoff modelling approach for ungauged sites in Taiwan

This work develops a top‐down modelling approach for storm‐event rainfall–runoff model calibration at unmeasured sites in Taiwan. Twenty‐six storm events occurring in seven sub‐catchments in the Kao‐Ping River provided the analytical data set. Regional formulas for three important features of a streamflow hydrograph, i.e. time to peak, peak flow, and total runoff volume, were developed via the characteristics of storm event and catchment using multivariate regression analysis. Validation of the regional formulas demonstrates that they reasonably predict the three features of a streamflow hydrograph at ungauged sites. All of the sub‐catchments in the study area were then adopted as ungauged areas, and the three streamflow hydrograph features were calculated by the regional formulas and substituted into the fuzzy multi‐objective function for rainfall–runoff model calibration. Calibration results show that the proposed approach can effectively simulate the streamflow hydrographs at the ungauged sites. The simulated hydrographs more closely resemble observed hydrographs than hydrographs synthesized using the Soil Conservation Service (SCS) dimensionless unit hydrograph method, a conventional method for hydrograph estimation at ungauged sites in Taiwan. Copyright © 2008 John Wiley & Sons, Ltd.     

Coupled surface water–groundwater model and its application in the arid Shiyang River basin,China

The planning and management of water resources in the Shiyang River basin, China require a tool for assessing the impact of groundwater and stream use on water supply reliabilities and improving many environment‐related problems such as soil desertification induced by recent water‐related human activities. A coupled model, integrating rule‐based lumped surface water model and distributed three‐dimensional groundwater flow model, has been established to investigate surface water and groundwater management scenarios that may be designed to restore the deteriorated ecological environment of the downstream portion of the Shiyang River basin. More than 66% of the water level among 24 observation wells have simulation error less than 1·0 m. The overall trend of the temporal changes of simulated and observed surface runoff at the Caiqi gauging station remains almost the same. The calibration was considered satisfactory. Initial frameworks for water allocation, including agricultural water‐saving projects, water diversion within the basin and inter‐basin water transfer, reducing agricultural irrigation area and surface water use instead of groundwater exploitation at the downstream were figured out that would provide a rational use of water resources throughout the whole basin. Sixteen scenarios were modelled to find out the most appropriate management strategies. The results showed that in the two selected management options, the groundwater budget at the Minqin basin was about 1·4 × 10⁸ m³/a and the ecological environment would be improved significantly, but the deficit existed at the Wuwei basin and the number was about 0·8 × 10⁸ m³/a. Water demand for domestic, industry and urban green area would be met in the next 30 years, but the water shortage for meeting the demand of agricultural water use in the Shiyang River basin was about 2·2 × 10⁸ m³/a. It is suggested that more inter‐basin water transfer should be required to obtain sustainable water resource use in the Shiyang River basin. Copyright © 2009 John Wiley & Sons, Ltd.     

δ13C signature of early graphite and subsequently formed microdiamond from the Saxonian Erzgebirge,Germany

Microdiamond and graphite occur in quartzofeldspathic rocks (saidenbachites) of the central Erzgebirge in Saxony. Zircon from such rocks contains not only relatively high but also variable amounts of early formed graphite and later crystallized diamond, however no carbonate. Six purified zircon concentrates were analysed by mass spectrometry for δ₁₃C(PDB) values which were between −24 and −33‰. The carbon inclusions resulted in C contents of the concentrates as high as 0.5 wt%. It is hypothesized that the observed δ₁₃C signature results from organic material in the sedimentary protoliths of the saidenbachites probably deposited not earlier than Mid‐Devonian. This material was transformed to graphite during early metamorphism and subsequently to diamond after an anatectic event at great depths (>150 km). During these processes, the δ₁₃C signature of the carbon did not (significantly) change.     

A new approximation for pore pressure accumulation in marine sediment due to water waves

The residual mechanism of wave‐induced pore water pressure accumulation in marine sediments is re‐examined. An analytical approximation is derived using a linear relation for pore pressure generation in cyclic loading, and mistakes in previous solutions (Int. J. Numer. Anal. Methods Geomech. 2001; 25 :885–907; J. Offshore Mech. Arctic Eng. (ASME) 1989; 111 (1):1–11) are corrected. A numerical scheme is then employed to solve the case with a non‐linear relation for pore pressure generation. Both analytical and numerical solutions are verified with experimental data (Laboratory and field investigation of wave– sediment interaction. Joseph H. Defrees Hydraulics Laboratory, School of Civil and Environmental Engineering, Cornell University, Ithaca, NY, 1983), and provide a better prediction of pore pressure accumulation than the previous solution (J. Offshore Mech. Arctic Eng. (ASME) 1989; 111 (1):1–11). The parametric study concludes that the pore pressure accumulation and use of full non‐linear relation of pore pressure become more important under the following conditions: (1) large wave amplitude, (2) longer wave period, (3) shallow water, (4) shallow soil and (5) softer soils with a low consolidation coefficient. Copyright © 2006 John Wiley & Sons, Ltd.     

Multifractal analysis as a tool for validating a rainfall model

A multifractal analysis was carried out in order to validate the simulation of hourly rainfall records of a local climate model for the Iberian Peninsula. Observed and simulated hourly rainfall data from four locations in Andalusia (southern Spain) were used to carry out the study. In order to detect the influence of the length of the data series on the results, two different sizes were used for the real data: 4 years, and 20 years. The results show that algebraic tails are required to fit the probability distribution of extreme rain event sizes, and rain and dry event durations for both kinds of rainfall data. Similar results are found for the extreme rain event sizes and dry event durations fits when the real and synthetic data are considered. Nevertheless, some differences appear in the cases of rain event durations. The detection of the presence of a first‐order multifractal phase transition associated with a critical moment in the empirical moment scaling exponent function and the results of the extreme rain event sizes fits, reveal that real rainfall is a self‐organized criticality (SOC) process. That behaviour is less evident in the simulated rainfall series. The same ‘synoptic maximum’ value was found for each place with both types of rainfall data. A time clustering analysis was carried out applying the count‐based periodogram and the Fano factor methods. Some periodicities have been detected in the periodograms, especially for the longest real rainfall data series. Copyright © 2007 John Wiley & Sons, Ltd.     

Long‐term trend of precipitation and runoff in Korean river basins

The spatial and temporal variations of precipitation and runoff for 139 basins in South Korea were investigated for 34 years (1968–2001). The Precipitation‐Runoff Modelling System (PRMS) was selected for the assessment of basin hydrologic response to varying climates and physiology. A non‐parametric Mann–Kendall's test and regression analysis are used to detect trends in annual, seasonal, and monthly precipitation and runoff, while Moran's I is adapted to determine the degree of spatial dependence in runoff trend among the basins. The results indicated that the long‐term trends in annual precipitation and runoff were increased in northern regions and decreased in south‐western regions of the study area during the study period. The non‐parametric Mann–Kendall test showed that spring streamflow was decreasing, while summer streamflow was increasing. April precipitation decreased between 15% and 74% for basins located in south‐western part of the Korean peninsula. June precipitation increased between 18% and 180% for the majority of the basins. Trends in seasonal and monthly streamflow show similar patterns compared to trends in precipitation. Decreases in spring runoff are associated with decreases in spring precipitation which, accompanied by rising temperatures, are responsible for reducing soil moisture. The regional patterns of precipitation and runoff changes show a strong to moderate positive spatial autocorrelation, suggesting that there is a high potential for severe spring drought and summer flooding in some parts of Korea if these trends continue in the future. Copyright © 2008 John Wiley & Sons, Ltd.     

Bayesian analysis of stage–fall–discharge models for gauging stations affected by variable backwater

Parsimonious stage–fall–discharge rating curve models for gauging stations subject to backwater complications are developed from simple hydraulic theory. The rating curve models are compounded in order to allow for possible shifts in the hydraulics when variable backwater becomes effective. The models provide a prior scientific understanding through the relationship between the rating curve parameters and the hydraulic properties of the channel section under study. This characteristic enables prior distributions for the rating curve parameters to be easily elicited according to site‐specific information and the magnitude of well‐known hydraulic quantities. Posterior results from three Norwegian and one American twin‐gauge stations affected by variable backwater are obtained using Markov chain Monte Carlo simulation techniques. The case studies demonstrate that the proposed Bayesian rating curve assessment is appropriate for developing rating procedures for gauging stations that are subject to variable backwater. Copyright © 2009 John Wiley & Sons, Ltd.