On Inflow to a Tunnel in a Fractured Double-Porosity Aquifer

Inflow to a tunnel is a great public concern and is closely related to groundwater hydrology, geotechnical engineering, and mining engineering, among other disciplines. Rapid computation of inflow to a tunnel provides a timely means for quickly assessing the inflow discharge, thus is critical for safe operation of tunnels. Dewatering of tunnels is another engineering practice that should be planned. In this study, an analytical solution of the inflow to a tunnel in a fractured unconfined aquifer is obtained. The solution takes into account either the spherical or slab-shaped matrix block and the unsteady state interporosity flow. The instantaneous drainage water table and anisotropic hydraulic conductivities of the fractures network are also considered. Both uniform flux and uniform head boundary condition are considered to simulate the constant head boundary condition in the tunnel. The effects of the hydraulic parameters of the fractured aquifer on the inflow variation of the tunnel are explored. The application of the presented solution to obtain the optimum location and discharge of the well to minimize the inflow to a tunnel is illustrated.     

Lithospheric structure of NW Iran from P and S receiver functions

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b?>?5.5, epicentral distance between 30° and 95° for P receiver functions and M_b?>?5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.     

Seismic performance and damage incurred by monolithic concrete self-centering rocking walls under the effect of axial stress ratio

Self-centering rocking walls offer the possibility of minimizing repair costs and downtimes, and also nullify the residual drift after seismic events, thanks to their self-centering properties. In this paper, the effect of axial stress ratio on the behavior of monolithic self-centering rocking walls is investigated by utilizing a developed finite element model. To verify the validity of the finite element model, results and observed damage in the model are compared with those of a full-scale wall test. The axial stress ratio is varied from 0.024 to 0.30 while keeping the other structural parameters constant. For qualitative damage evaluation, the observed damage in the model compared with expected damage states of desired performance levels. In order to evaluate the incurred damage quantitatively, the amount of crushing and damage in the wall is calculated by utilizing several ratios (crushing ratio and damage ratio). Furthermore, seismic response factors (i.e., μ, R and C_d) are calculated for different axial stress ratio values. The obtained results showed that, in order to satisfy the requirements of desired performance levels, the maximum axial stress ratio should be approximately within the range of 0.10–0.15. In addition, the maximum overall damage ratio and crushing ratio are suggested to be less than 5%. For axial stress ratio higher than 0.15, the flag-shaped pattern of hysteresis curves completely disappeared and the variation of displacement ductility is less sensitive to axial stress ratio. Considering the maximum axial stress ratio limited to 0.150, values of 4 and 3.5 are conservatively proposed as a period-independent response modification factor and displacement modification factor of the investigated structural wall, respectively.     

Impacts of Methane on Carbon Dioxide Storage in Brine Formations

In the context of geological carbon sequestration (GCS), carbon dioxide (CO₂) is often injected into deep formations saturated with a brine that may contain dissolved light hydrocarbons, such as methane (CH₄). In this multicomponent multiphase displacement process, CO₂ competes with CH₄ in terms of dissolution, and CH₄ tends to exsolve from the aqueous into a gaseous phase. Because CH₄ has a lower viscosity than injected CO₂, CH₄ is swept up into a ‘bank’ of CH₄‐rich gas ahead of the CO₂ displacement front. On the one hand, this may provide a useful tracer signal of an approaching CO₂ front. On the other hand, the emergence of gaseous CH₄ is undesirable because it poses a leakage risk of a far more potent greenhouse gas than CO₂ if the cap rock is compromised. Open fractures or faults and wells could result in CH₄ contamination of overlying groundwater aquifers as well as surface emissions. We investigate this process through detailed numerical simulations for a large‐scale GCS pilot project (near Cranfield, Mississippi) for which a rich set of field data is available. An accurate cubic‐plus‐association equation‐of‐state is used to describe the non‐linear phase behavior of multiphase brine‐CH₄‐CO₂ mixtures, and breakthrough curves in two observation wells are used to constrain transport processes. Both field data and simulations indeed show the development of an extensive plume of CH₄‐rich (up to 90 mol%) gas as a consequence of CO₂ injection, with important implications for the risk assessment of future GCS projects.     

Hydrological response characteristics of Mediterranean catchments at different time scales: a meta-analysis

ABSTRACT

This work examines 140 hydrological studies conducted in the Mediterranean region. It identifies key characteristics of the hydrological responses of Mediterranean catchments at various time scales and compares different methods and modelling approaches used for individual-catchment studies. The study area is divided into the northwestern (NWM), eastern (EM) and southern (SM) Mediterranean. The analysis indicates regional discrepancies in which the NWM shows the most extreme rainfall regime. A tendency for reduced water resources driven by both anthropogenic and climatic pressures and a more extreme rainfall regime are also noticeable. Catchments show very heterogeneous responses over time and space, resulting in limitations in hydrological modelling and large uncertainties in predictions. However, few models have been developed to address these issues. Additional studies are necessary to improve the knowledge of Mediterranean hydrological features and to account for regional specificities.

Editor D. Koutsoyiannis Associate editor A. Efstratiadis     

Introducing a new conceptual instantaneous unit hydrograph model based on a hydraulic approach

ABSTRACT

A novel approach is introduced for simulation of instantaneous unit hydrographs (IUHs). The model consists of a series of linear reservoirs that are connected to each other, and is referred to as the inter-connected linear reservoir model (ICLRM). By assuming the flow between two reservoirs is a linear function of the difference between the water levels in the reservoirs, a system of first-order linear differential equations is obtained as the ICLRM governing equation. By solving the equations, the discharge from the last reservoir is considered as an IUH. A small-scale laboratory device was constructed for the simulation of IUHs using the model. By studying four hydrographs extracted from the literature, and simulating them using both the ICLRM and the Nash model, it is concluded that the ICLRM can predict these hydrographs more accurately than the Nash model. Due to the simplicity of the construction and operation of the ICLRM and, more importantly, its visual aspect, the ICLRM may be considered as an effective educational tool for studying IUHs.     

Earthquake relocation in the Central Alborz region of Iran using a non-linear probabilistic method

In this study, we calculate accurate absolute locations for nearly 3,000 shallow earthquakes (≤20 km depth) that occurred from 1996 to 2010 in the Central Alborz region of northern Iran using a non-linear probabilistic relocation algorithm on a local scale. We aim to produce a consistent dataset with a realistic assessment of location errors using probabilistic hypocenter probability density functions. Our results indicate significant improvement in hypocenter locations and far less scattering than in the routine earthquake catalog. According to our results, 816 earthquakes have horizontal uncertainties in the 0.5–3.0 km range, and 981 earthquakes are relocated with focal-depth errors less than 3.0 km, even with a suboptimal network geometry. Earthquake relocated are tightly clustered in the eastern Tehran region and are mainly associated with active faults in the study area (the Mosha and Garmsar faults). Strong historical earthquakes have occurred along the Mosha and Garmsar faults, and the relocated earthquakes along these faults show clear north-dipping structures and align along east–west lineations, consistent with the predominant trend of faults within the study region. After event relocation, all seismicity lies in the upper 20 km of the crust, and no deep seismicity (>20 km depth) has been observed. In many circumstances, the seismicity at depth does not correlate with surface faulting, suggesting that the faulting at depth does not directly offset overlying sediments.     

Frequency analysis of climate extreme events in Zanjan, Iran

In this study, generalized extreme value distribution (GEV) and generalized Pareto distribution (GPD) were fitted to the maximum and minimum temperature, maximum wind speed, and maximum precipitation series of Zanjan. Maximum (minimum) daily and absolute annual observations of Zanjan station from 1961 to 2011 were used. The parameters of the distributions were estimated using the maximum likelihood estimation method. Quantiles corresponding to 2, 5, 10, 25, 50, and 100 years return periods were calculated. It was found that both candidate distributions fitted to extreme events series, were statistically reasonable. Most of the observations from 1961 to 2011 were found to fall within 1–10 years return period. Low extremal index (θ) values were found for excess maximum and minimum temperatures over a high threshold, indicating the occurrence of consecutively high peaks. For the purpose of filtering the dependent observations to obtain a set of approximately independent threshold excesses, a declustering method was performed, which separated the excesses into clusters, then the de-clustered peaks were fitted to the GPD. In both models, values of the shape parameters of extreme precipitation and extreme wind speed were close to zero. The shape parameter was less negative in the GPD than the GEV. This leads to significantly lower return period estimates for high extremes with the GPD model.