An investigation of near surface wind speed trends in arid and semiarid regions of Iran

Wind plays an important role on the ecosystems and hydrological cycles besides other meteorological parameters such as temperature, precipitation, sunshine, and relative humidity. It strongly affects evapotranspiration, especially in arid and semiarid regions where there are serious problems in regard to water resource management. Evaluating the wind speed trend can provide good information for future agricultural planning. This study was conducted in order to investigate the wind speed trends over 24 synoptic meteorological stations located in arid and semiarid regions of Iran from 1975 to 2005. Near-surface wind speed was trended by nonparametric Mann–Kendall test spatially and temporally in three time scales (annual, seasonal, and monthly). Then, Sen’s slope estimator was used to determine the amount of the changes; furthermore, 10-year moving average low-pass filter was applied to show general trends. Finally, the smoothed time series derived from the mentioned filter were classified in three clusters for each time series and then mapped to show their spatial distribution pattern. Results showed insignificant and significant, increasing and decreasing trends during the surveyed time. Wind speeds in less than 50 % of stations changed statistically in all time scales, and in most cases, the frequency of the upward trends was more than that of downward ones. The spatial distribution of significant wind speed showed that the increase mostly occurred in eastern part. Clustering gave us the turning point around 1990. Clearly, when clusters were mapped, they indicated the same pattern as the Z value maps derived from Mann–Kendall test which meant that the outputs of the mentioned method confirmed the other one. As the wind speed trends in different stations likely to follow the previous evapotranspiration (ET₀) trend results in Iran, it confirms that wind speed was an effective parameter on ET₀, even though other parameters should be considered too.     

The Ganges and Brahmaputra rivers in Bangladesh: basin denudation and sedimentation

Every year the Ganges and Brahmaputra rivers in Bangladesh transport 316 and 721 million tonnes of sediment, respectively. These high loads of suspended sediment reflect the very high rate of denudation in their drainage basins. The average mechanical denudation rate for the Ganges and Brahmaputra basins together is 365 mm 10³ yr⁻¹. However, the rate is higher in the Brahmaputra Basin than that in the Ganges Basin. Several factors, including mean trunk channel gradient, relief ratio, runoff, basin lithology and recurring earthquakes are responsible for these high denudation rates. Of the total suspended sediment load (i.e. 1037 million tonnes) transported by these rivers, only 525 million tonnes (c. 51% of the total load) are delivered to the coastal area of Bangladesh and the remaining 512 million tonnes are deposited within the lower basin, offsetting the subsidence. Of the deposited load, about 289 million tonnes (about 28% of the total load) are deposited on the floodplains of these rivers. The remaining 223 million tonnes (about 21% of the total load) are deposited within the river channels, resulting in aggradation of the channel bed at an average rate of about 3·9 cm yr⁻¹. Although the Brahmaputra transports a higher sediment load than the Ganges, the channel bed aggradation rate is much higher for the Ganges. This study also documents a wide range of interannual, seasonal and daily variation in suspended sediment transport and water discharge. Interannual variation in sediment deposition within the basin is also suggested. Copyright © 1999 John Wiley & Sons, Ltd.     

The Effect of Hydraulic Gradient and Pattern of Conduit Systems on Tracing Tests: Bench-Scale Modeling

Tracer breakthrough curves provide valuable information about the traced media, especially in inherently heterogeneous karst aquifers. In order to study the effect of variations in hydraulic gradient and conduit systems on breakthrough curves, a bench scale karst model was constructed. The bench scale karst model contains both matrix and a conduit. Eight tracing tests were conducted under a wide range of hydraulic gradients from 1 to greater than 5 for branchwork and network-conduit systems. Sampling points at varying distances from the injection point were utilized. Results demonstrate that mean tracer velocities, tracer mass recovery and linear rising slope of the breakthrough curves were directly controlled by hydraulic gradient. As hydraulic gradient increased, both one half the time for peak concentration and one fifth the time for peak concentration decreased. The results demonstrate the variations in one half the time for peak concentration and one fifth the time for peak concentration of the descending limb for different sampling points under differing hydraulic gradients are mainly controlled by the interactions of advection with dispersion. The results are discussed from three perspectives: different conduit systems, different hydraulic-gradient conditions, and different sampling points. The research confirmed the undeniable role of hydrogeological setting (i.e., hydraulic gradient and conduit system) on the shape of the breakthrough curve. The extracted parameters (mobile-fluid velocity, tracer-mass recovery, linear rising limb, one half the time for peak concentration, and one fifth the time for peak concentration) allow for differentiating hydrogeological settings and enhance interpretations the tracing tests in karst aquifers.     

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.     

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.