Inundation flow velocity of tsunami on land

An estimation of tsunami inundation flow velocity is one of the most challenging issues among tsunami research. Based on field data of inundation depth and inundation flow velocity u estimated using Bernoulli's theorem and inundation depth, fundamental characteristics of the relationship between inundation flow velocity and inundation depth are examined. Fundamental characteristics of the velocity coefficient where g is gravitational acceleration, h_f and h_r are inundation depths at the front and the back of structures such as a rectangular building with vertical walls, respectively) implicitly included in the relationship are examined through hydraulic experiments. As a result, C_v = 0.6 is recommended as its simple and practical value. It is confirmed through these examinations that the Froude number, defined by where , ranges 0.7–2.0, and when C_v = 0.6 is adopted this Froude number ranges 0.42–1.2. By using the relationship and C_v = 0.6, two simple and practical relationships are presented for two cases where inundation flow velocity exerts the largest or the smallest fluid force on structures. These relationships can be used to roughly grasp the practical side of tsunami damage, and estimate fluid force acting on individual structures, moving velocity and collision force of floating objects and sediment transport such as boulder and sand. Fundamental characteristics of the waterline (tsunami trace) distribution around/on the typical object model (square pillar, corn and column) are also examined through steady flow experiments, and it is confirmed that the maximum and the minimum values of h_f/h₀ in the full type model of the square pillar are almost the same as those of h_f/h_r obtained by field surveys where h₀ is uniform flow depth. It is also confirmed that h_r ? h₀ when the Froude number, defined by where u₀ is uniform flow velocity, is much less than 1.0. Using a newly defined velocity coefficient, tsunami inundation flow velocity on land can be estimated practically and would be useful for checking proposed sediment transport models that are now being developed by tsunami geologists.     

Mobility Enhancement of Mn Oxide during Permanganate Oxidation of TCE

Manganese (Mn) oxide precipitation during in situ permanganate oxidation of organic compounds can cause pore clogging, reduce permeability, and increase resistance to mass transfer. Stability of Mn oxide is required to enhance oxidation effectiveness. Batch tests were conducted at eight polyphosphate (PP) to permanganate () mass ratios (0 to 8) at three MnO₄⁻¹ concentrations (500, 1000, or 2000 mg/L) for identifying mass ratios to maximize stability of Mn oxide produced in the presence of trichloroethylene (TCE). In general, stability of Mn oxide was the maximum at mass ratio of approximately 4. Three column tests were designed to investigate the impact of PP on overall removal of 4.6 or 7.0 g TCE emplaced as nonaqueous phase liquid within the column porous media. Water flush, chemical flush using alone (1000 mg/L), and chemical flush using (1000 mg/L) and PP (4000 mg/L) were conducted. Mass removal of TCE and changes in media permeability were estimated over a period of 78 to 312 h (12 to 49 pore volumes [PVs]). Column tests demonstrated enhanced removal (~90%) of TCE during chemical flush with and PP in 12 PVs as compared with approximately 64% during -only flush and approximately 26% during water flush. Pressure drop changes across the column captured change in media permeability and revealed that water flush and PP and flush caused significantly lower flow resistance as compared with -only flush. These results indicate that PP was capable of mobilizing Mn oxide away from the reaction zones, thereby reducing pore clogging and enabling better and long-term contact between TCE and the aqueous phase.     

Mineral Content of Bottled and Desalinated Household Drinking Water in Kuwait

Three hundred and twenty‐two samples of desalinated household water were collected from 99 sampling locations that covered 95% of Kuwaiti's residential areas. Seventy‐one brands of bottled water were collected from Kuwaiti markets. The water quality parameters that were studied included pH, electrical conductivity (EC), total dissolved solids (TDS), F^?, Cl^?, Br^?, , , , , , , , and the major macronutrients Na⁺, K⁺, Ca²⁺, and Mg²⁺. The analysis yielded a large range of results for most of these parameters, with differences in some cases exceeding 10‐fold. With a few exceptions, the results were found to comply with US‐EPA and WHO standards. Only the water in two brands of bottled water was acidic (pH < 6.5). The TDS was found to be higher than the US‐EPA regulated value in 4 and 3% of the household samples and bottled water brands, respectively. The fluoride levels were generally higher in bottled water than in household water. However, the household water that was produced by the Doha desalination plant and some of the European brands of bottled water were the best samples studied in terms of their quantity of Ca²⁺, Mg²⁺, and Na⁺ compared with the DRI values for those substances. EC and TDS were positively correlated with , , Na⁺, K⁺, Ca²⁺, and Mg²⁺ for household water but only with Ca²⁺ and Mg²⁺ for bottled water.     

Removal of Bromate from Water Using De‐Acidite FF‐IP Resin and Determination by Ultra‐Performance Liquid Chromatography‐Tandem Mass Spectrometry

Contamination of water due to bromate is a severe health hazard. The aim of the present study was to remove bromate from water using a crosslinked polystyrene based strongly basic anion exchange resin De‐Acidite FF‐IP. Batch experiments were performed to study the influence of various experimental parameters such as effect of pH, contact time, temperature, and effect of competing anions on bromate removal by De‐Acidite FF‐IP resin. At optimum parameters, the removal rate of bromate was very fast and 90% removal took place in 5 min and equilibrium was established within 10 min. The presence of competitive anions reduced the bromate adsorption in the order of Cl^? > F^? > CO > SO > NO > PO. The practical utility of this resin has been demonstrated by removing bromate in some of the commercial bottled water from Saudi Arabia. The level of bromate was determined using a very sensitive, precise and rapid method based on ultra‐performance liquid chromatography‐tendem mass spectrometry (UPLC‐MS/MS).     

From an initial transient-state to a steady-state in metamorphic reactions: An experimental approach in the system dolomite-quartz-H2O

Abstract We carried out hydrothermal experiments in the system dolomite‐quartz‐H₂O to track the temporal change in reaction rates of simultaneous reactions during the development of reaction zones. Two types of configurations for the starting materials were prepared: dolomite single crystals + quartz powder + water and quartz single crystals + dolomite powder + water, both sealed separately in gold capsules. Runs at 0.1GPa and 600°C with cold seal pressure vessels gave the following results. (i) In short duration (45–71 h) runs metastable layer sequences involving wollastonite and talc occur in the reaction zone, whereas they disappear in longer duration (168–336 h) runs. (ii) The layer sequence of the reaction zones in short duration runs differs from place to place on the dolomite crystal even in the same run. (iii) The diversity of layer sequences in the short duration runs merges into a unique layer sequence in the longer duration runs. (iv) The reaction zone develops locally on the dolomite crystal, but no reaction zone was observed on the quartz crystal in any of the runs. The lines of evidence (i)–(iii) show that the system evolves from an initial transient‐ to a steady‐state and that the kinetic effect is important in the development of reaction zones. A steady diffusion model for the unique layer sequence Qtz/Di/Fo + Cal/Dol + Cal/Dol shows that the Dol + Cal layer cannot be formed by diffusion‐controlled process and that the stability of the layer sequence Qtz/Di/Fo + Cal/Dol depends not only on L‐ratios (a = /L_CaOCaO and b = /L_MgOMgO) but also on the relative rate P = (−2ξ₁ − ξ₂)/(–ξ₁ − 2ξ₂) of competing reactions: Dol + 2Qtz = Di + 2CO₂ (ξ₁) and 2Dol + Qtz = Fo + 2Cal + 2CO₂ (ξ₂). For smaller P the stability field will shift to higher values of a and b. The steady diffusion model also shows that the apparent‐non‐reactivity on the quartz surface can be attributed to void formation in a large volume fraction in the diopside layer.     

VSP traveltime inversion for anisotropy in a buried layer

We present a method for calculating the anisotropy parameter of a buried layer by inverting the total traveltimes of direct arrivals travelling from a surface source to a well‐bore receiver in a vertical seismic profiling (VSP) geometry. The method assumes two‐dimensional media. The medium above the layer of interest (and separated from it by a horizontal interface) can exhibit both anisotropy and inhomogeneity. Both the depth of the interface as well as the velocity field of the overburden are assumed to be known. We assume the layer of interest to be homogeneous and elliptically anisotropic, with the anisotropy described by a single parameter χ. We solve the function describing the traveltime between source and receiver explicitly for χ. The solution is expressed in terms of known quantities, such as the source and receiver locations, and in terms of quantities expressed as functions of the single argument x_r, which is the horizontal coordinate of the refraction point on the interface. In view of Fermat's principle, the measured traveltime T possesses a stationary value or, considering direct arrivals, a minimum value, . This gives rise to a key result ‐‐ the condition that the actual anisotropy parameter . Owing to the explicit expression , this result allows a direct calculation of in the layer of interest. We perform an error analysis and show this inverse method to be stable. In particular, for horizontally layered media, a traveltime error of one millisecond results in a typical error of about 20% in the anisotropy parameter. This is almost one order of magnitude less than the error inherent in the slowness method, which uses a similar set of experimental data. We conclude by detailing possible extensions to non‐elliptical anisotropy and a non‐planar interface.     

Experimental Designs Applied to Desorption of Dichromate Ions after Separation and Preconcentration from Natural and Industrial Water by Modified Nano‐Alumina

Nano‐alumina modified by 9‐aminoacridine was used as a sorbent for separation and determination of dichromate ions from water. Statistical method, based on surface response design, has been used for the optimization of dichromate ions elution from 9‐aminoacridine nano‐alumina. The adsorbed dichromate ions were found to be eluted quantitatively with 0.8 mol L^?1 KCl in 1.6 mol L^?1 NaOH which optimized by response surface design. Under optimum conditions, the accuracy, precision (relative standard deviation, RSD%) and R‐square of the method were calculated as >98, <3, and >94%, respectively. Remarkable agreement between experimental and theoretical data was confirmed the predicted assumption. The method was applied to the simultaneous determination of dichromate in natural and industrial water samples. We also examined the retention of dichromate anions in the presence of Cl^?, , and anions at pH 3.     

Spatial and temporal variations in summer precipitation in Japanese mountain areas

We examined spatial and temporal variations in precipitation measured during summer season between 1976 and 2007 for 28 stations located in mountain areas across Japan using the amount of precipitation (Pr), the mean depth of precipitation events (η), and the inverse of the mean interval times (λ). We obtained positive correlations between the period mean Pr (Pr ) and the period mean η ( ) and between Pr and the period mean λ ( ) for the 28 stations. Pr was more strongly related to than to , indicating the spatial variations in Pr that are primarily related to the variations in . In addition, Pr was more strongly related to η than to λ for most stations on the basis of data for 1976–2007, indicating that the year‐to‐year variations in Pr are primarily related to η. We also examined temporal trends in Pr, η and λ for 1976–2007 and found no systematic trends for 23 of the 28 stations, suggesting long‐term trends that are not common in mountain areas of Japan. The relationships between Pr and and between Pr and η presented in this study enable us to generate a temporal precipitation distribution pattern based on only Pr and Pr data, respectively. Furthermore, probabilistic stochastic hydrological models require precipitation characteristics as input; thus, this study contributes to the determination of hydrological cycles and their possible future changes in Japanese mountain areas and therefore to water resource management. Copyright © 2010 John Wiley & Sons, Ltd.     

Theoretical investigation of the time variation of drainage density

A theoretical equation was developed to express the time variation of drainage density in a basin or geomorphic surface: D_i(t, T) is the drainage density at time T on the i-th basin or geomorphic surface, which was formed at time t; β(τ) is a factor related to the erosional force causing the development of the rivers of the basin or surface at time τ; δ_i is the maximum drainage density; and D_i is the initial drainage density on the i-th geomorphic surface or basin. The equation is based on the assumption that the drainage density increases with time until it reaches a specific upper limit δ_i(t)), the maximum drainage density, which is related to certain physical properties of the basin. The equations for various dated basins or geomorphic surfaces can be combined into one modified equation if the same relative erosional forces have acted on those basins or surfaces (β(t) = β(t) and if the basins or surfaces have the same physical properties δ_i(t) = δ_i(t), (D_i = D₀). The application of this equation to coastal terraces and glacial tills shows that the model is compatible with observed drainage densities on various dated basins or surfaces.     

Chemical dynamics of N-containing ionic species in a boreal forest snowcover during the spring melt period

A study of the changes in the ionic loads of NO, NH, SO and H⁺ in a boreal forest snowpack at Lake Laflamme, Québec was carried out using hydrological and chemical data from field lysimeters. The results showed that depletion of the N-containing species occurs periodically in the snowpack during meltwater discharge. Rain-on-snow events led to in-pack losses of NO and NH at a rate of 130 μeq m^?2 day^?1 and 101·3 μeq m^?2day^?1 respectively. On dry days, however, dry deposition and deposition of organic debris from the canopy resulted in increases of 183·3 μeq m^?2day^?1 for NO and 4·5 μeq m^?2day^?1 for NH in the pack. In contrast, SO₄^2? showed continual in-pack increases due to deposition of 5·0 μeq m^?2day^?1 for wet days and 92·6 μeq m^?2day^?1 for dry days. The depletion of NO and NH is due to microbiological uptake of these nutrients during periods when the free water content of the pack is high. Controlled melts in a laboratory snowmelt simulator containing snow and organic matter from the forest canopy at Lake Laflamme showed losses of NO and NH similar to those observed in the field. As the microbiological uptake proceeds at a rate comparable to that of ionic load increases in the pack by dry deposition, models of the chemical dynamics of snowmelt should take the former into account in any system where organic content of the snowpack is appreciable.     

Average sand particle trajectory examined by the Raindrop Detachment and Wind‐driven Transport (RD‐WDT) process

Recent studies of soil loss by the integrated action of raindrop impact and wind transport have demonstrated the significance of this mechanism. This paper presents data obtained during wind‐tunnel experiments examining the ‘Raindrop Detachment and Wind‐driven Transport’ (RD‐WDT) process to investigate average sand particle trajectory and the spatial extent at which the process operates. In the experimental design, at the same time as the horizontal wind velocities of 6·4, 10, and 12 m s^–1 passed through the tunnel, rainfall was simulated falling on very well sorted dune sand. The aspect and slope of the sand bed was varied to reproduce both windward (Ww) and leeward (Lw) slopes of 4º and 9º with respect to the prevailing wind direction. The average sand particle trajectories by the RD‐WDT process ( ) were estimated by a mass‐distribution function, which was integrated over a 7‐m uniform slope segment. The results showed that depended statistically upon the wind shear velocity (u*), and the effect of the slope gradient (θ) was insignificant on . This was different from that of the windless rain process ( ), ‘Raindrop Detachment and Splash‐driven Transport’ (RD‐ST), the spatial range of which relies strongly on θ. Additionally, was approximately 2·27 ± 2·2 times greater than the average path of a typical saltating sand particle of the rainless wind ( ), ‘Wind Erosion Saltation Transport’ (WE‐ST). Copyright © 2009 John Wiley & Sons, Ltd.     

Diel patterns in coastal‐stream nitrate concentrations linked to evapotranspiration in the riparian zone of a low‐relief,agricultural catchment

Evapotranspiration (ET) can cause diel fluctuations in the elevation of the water table and the stage in adjacent streams. The diel fluctuations of water levels change head gradients throughout the day, causing specific discharge through near‐stream sediment to fluctuate at the same time scale. In a previous study, we showed that specific discharge controls the residence time of groundwater in streambed sediment that, in turn, exerted the primary control on removal from groundwater passing through the streambed. In this study, we examine the magnitude of diel specific discharge patterns through the streambed driven by ET in the riparian zone with a transient numerical saturated–unsaturated groundwater flow model. On the basis of a first‐order kinetic model for removal, we predicted diel fluctuations in stream concentrations. Model results indicated that ET drove a diel pattern in specific discharge through the streambed and riparian zone (the removal zones). Because specific discharge is inversely proportional to groundwater travel time through the removal zones and travel time determines the extent of removal, diel changes in ET can result in a diel pattern in concentration in the stream. The model predictions generally matched observations made during summertime base‐flow conditions in a small coastal plain stream in Virginia. A more complicated pattern was observed following a seasonal drawdown period, where source components to the stream changed during the receding limb of the hydrograph and resulted in diel fluctuations being superimposed over a multi‐day trend in concentrations. Copyright © 2013 John Wiley & Sons, Ltd.     

The Impact of Geology of Recharge Areas on Groundwater Quality: A Case Study of Zhob River Basin,Pakistan

Groundwater is a major source of water supply for domestic and irrigation uses in semiarid, remote but rapidly developing Kilasaifullah district part of Zhob River Basin, located at Pakistan–Afghanistan Border. Zhob River is among few major rivers of perennial nature in Balochistan, which flows from WSW to ENE and falls in Gomal River, a tributary of Indus River. Keeping in view the important geopolitical position and rapid development of the region, this study is primarily focused on groundwater chemistry for contamination sources as well as agriculture development. Water samples from open and tube wells are analyzed and calculated for electrical conductivity (EC), total dissolved solids (TDS), turbidity, pH, K⁺, Na⁺, Ca²⁺, Mg²⁺, HCO, Cl^?, NO, SO, PO, sodium percent (Na%), sodium adsorption ratio (SAR), Kelly's index (KI), and heavy metals (Fe, Cu, Cr, Zn, Pb, and Mn). On the basis of the chemical constituents two zones within the study area are identified and possible causes of the contaminants are pointed out. Two recharge areas were responsible for the different chemical results in groundwater, e.g., zone A was recharged from NNW saline geological formations (Nisai, Khojak, Multana, Bostan formations, and Muslim Bagh ophiolites), which are concentrated with high sodium and chloride. On the other hand Zone B was sourced from SSW from carbonate rich rocks (Alozai, Loralai, Parh formations, and Muslim Bagh ophiolites). The groundwater is classified as C2–S1, C3–S1, C3–S2, C4–S2 on the basis of EC and SAR values which indicate that most of the water of both zones can be used for irrigation safely except the samples plotted in C3–S2 and C4–S2 categories which could be dangerous for soil and crops. Groundwater samples are plotted in good to permissible limits with some samples excellent to good and few samples belong to doubtful category based on sodium percent. Groundwater of zone A is unsuitable for irrigation use due to higher values of KI (more than one) but water of zone B are good for irrigation based on KI. In general, water of both zones is suitable for irrigation but care should be taken during the selection of crops which are sensitive to alkalinity or sodium hazards particularly in zone A.     

3D angle-domain common-image gathers for migration velocity analysis

Angle‐domain common‐image gathers (ADCIGs) are an essential tool for migration velocity analysis (MVA). We present a method for computing ADCIGs in 3D from the results of wavefield‐continuation migration. The proposed methodology can be applied before or after the imaging step in a migration procedure. When computed before imaging, 3D ADCIGs are functions of the offset ray parameters (p, p) ; we derive the geometric relationship that links the offset ray parameters to the aperture angle γ and the reflection azimuth φ. When computed after imaging, 3D ADCIGs are directly produced as functions of γ and φ. The mapping of the offset ray parameters (p, p) into the angles (γ, φ) depends on both the local dips and the local interval velocity; therefore, the transformation of ADCIGs computed before imaging into ADCIGs that are functions of the actual angles is difficult in complex structure. By contrast, the computation of ADCIGs after imaging is efficient and accurate even in the presence of complex structure and a heterogeneous velocity function. On the other hand, the estimation of the offset ray parameters (p, p) is less sensitive to velocity errors than the estimation of the angles (γ, φ). When ADCIGs that are functions of the offset ray parameters (p, p) are adequate for the application of interest (e.g. ray‐based tomography), the computation of ADCIGs before imaging might be preferable. Errors in the migration velocity cause the image point in the angle domain to shift along the normal to the apparent geological dip. By assuming stationary rays (i.e. small velocity errors), we derive a quantitative relationship between this normal shift and the traveltime perturbation caused by velocity errors. This relationship can be directly used in an MVA procedure to invert depth errors measured from ADCIGs into migration velocity updates. In this paper, we use it to derive an approximate 3D residual moveout (RMO) function for measuring inconsistencies between the migrated images at different γ and φ. We tested the accuracy of our kinematic analysis on a 3D synthetic data set with steeply dipping reflectors and a vertically varying propagation velocity. The tests confirm the accuracy of our analysis and illustrate the limitations of the straight‐ray approximation underlying our derivation of the 3D RMO function.     

Factors controlling groundwater chemistry in an agricultural area with complex topographic and land use patterns in mid‐western South Korea

Shallow and bedrock groundwater from granitic aquifers were investigated for the hydrogeochemistry of major and minor constituents in an agricultural area. Nitrate concentrations were observed up to 49 mg/l as NO₃‐N, with 22% of samples exceeding the drinking water standard, which could pose a significant threat because most residents rely on groundwater as their drinking water source. Principal component analysis revealed three principal components (PCs): (1) nitrate contamination, contributed by major cations, Cl^?, SO and NO , (2) reduction processes positively involving Fe, Mn and B, and negatively involving dissolved oxygen and NO and (3) natural mineralization, involving HCO and F^?. Cluster analysis, performed on the PC scores, resulted in seven sample groups, which were successfully identified by total depth, elevation and land use. The nitrate‐contaminated groups had mixed land uses, with locally concentrated residential areas. Uncontaminated groundwater groups were found in the natural environment, including high‐altitude spring water and bedrock groundwater with a higher degree of natural mineralization. Shallow groundwater groups in paddy fields in lowlands were affected by reducing environments, of which one group was characterized by high Fe, Mn and B, and negligible nitrate. Groundwater with intermediate nitrate and lower Cl^? and SO was found primarily in hilly terrains with orchards and vegetable gardens, indicating lower contaminant loadings than lowland areas. Higher concentrations of F^? and nitrate were observed in the nitrate‐contaminated water, which seemed unlikely to be explained by groundwater mixing. The strong acidity generated from nitrification may infiltrate deeper into the aquifer, induce accelerated weathering of bedrock and result in the coexistence of F^? and nitrate, which may be an evidence of intense nitrate loading, leading to soil acidification. Multivariate statistical analysis successfully delineated hydrochemical characteristics of groundwater attained by natural and anthropogenic processes in an agriculturally stressed area with complex topographic land use patterns. Copyright © 2009 John Wiley & Sons, Ltd.     

Untersuchungen zur Phosphataufnahme durch Unterwasserpflanzen

The uptake of phosphorus by Potamogeton pectinatus, P. lucens, Ranunculus fluitans, Elodea canadensis, Myriophyllum spicatum and Ceratophyllum demersum is investigated in batch tests at 0.02, 0.04, 0.4 and 4.0 mg/l PO in the medium for 4, 28 and 52 h. Uptake is strictly dependent on concentration, a saturation function being discernible already after 52 h. At rising concentration the mechanism of uptake is determined by a decreasing rate of incorporation. The rates of uptake are dependent also on the actual and maximum cell nutrient contents. Species with a low phosphorus content of <5 mg PO/g dry weight as Elodea canadensis, Potamogeton pectinatus and Myriophyllum spicatum show their maximum phosphorus incorporation at low substrate concentrations. At phosphorus contents >5 mg PO/g dry weight the phosphorus incorporation is increased only at substrate concentrations of 0.04… 0.4 mg/l PO.     

Capturing temporal variability for estimates of annual hydrochemical export from a first‐order agricultural catchment in southern Ontario,Canada

This 2‐year study (2000, 2001) reports annual nutrient (phosphorus, nitrate) export from a first‐order agricultural watershed in southern Ontario based on an intensive monitoring programme. The importance of storm and melt events in annual export estimates is demonstrated and the temporal variability in nutrient loading during events is related to processes occurring within the catchment. The feasibility of predicting event‐related nutrient export from hydrometric data is explored. The importance of sampling frequency throughout events is also shown. Export of total phosphorus (TP), soluble reactive phosphorus (SRP) and nitrate ( ) for 2000 and 2001 averaged 0·35 kg ha^?1 year^?1, 0·09 kg ha^?1 year^?1, and 35 kg ha^?1 year^?1 (as N) respectively. Approximately 75% of annual TP export, 80% of annual SRP export and 70% of annual export occurred during 28 events per year. A small number of large‐magnitude events (>34 mm) accounted for 18–42% of annual TP export, 0–61% of annual SRP export and 13–33% of annual NO export. Our results show that temporal variability in nutrient export is largely governed by discharge in this basin, and export can be predicted from discharge. SRP and TP export can also be predicted from discharge, but only for events that are not large in magnitude. The sampling interval throughout events is important in obtaining precise estimates of nutrient export, as infrequent sampling intervals may over‐ or under‐estimate nutrient export by ± 45% per event for P. This study improves our understanding of and P export patterns and our ability to predict or model them by relating temporal variability in event nutrient export to discharge and processes occurring within the basin, and also by exploring the significance of sampling interval in the context of the importance of individual events, season and temporal variability during events. Copyright © 2006 John Wiley & Sons, Ltd.     

Water Quality Management Based on Division of Dry and Wet Seasons in Pearl River Delta,China

In the Pearl River Delta (PRD), river water quality deteriorates continually due to the population increase and ongoing industrialization and urbanization. In this study, a water quality management paradigm based on the seasonal variation is proposed. For better exploring the seasonal change of water quality, wavelet analysis was used to analyze the division of dry and wet seasons in the PRD during 1952–2009. Then water quality seasonal variation in 2008 and relevant impact factors were analyzed by multivariate statistic methods as a case to make some management measures. The results show that there are some differences of dry and wet seasons division among different years. Wet season mainly appear from April to September, which occupy the largest proportion among the 58 years (about 70%) and then followed by the wet season from May to October (about 13.8% of the total years). As to the water quality of 2008, significant differences exist between dry and wet seasons for 17 water quality parameters except TP, , Fe²⁺, and Zn²⁺. Levels of parameters pH, EC, COD_Mn, BOD₅, , , and Cl^? in dry season are much higher than those in wet season. In dry season the variations of river water quality are mainly influenced by domestic sewage, industrial effluents, and salt water intrusion. While in wet season, except the aforementioned pollution sources, drainages from cultivated land and livestock farm are also the main factors influencing water pollution. Thus, water quality management measures are proposed in dry and wet seasons, respectively. The results obtained from this study would further facilitate water quality protection and water resources management in the PRD.     

Spatial average aquifer recharge through atmospheric chloride mass balance and its uncertainty in continental Spain

The atmospheric chloride mass balance (CMB) method allows spatial evaluations of the average diffuse aquifer recharge by rainfall () in large and varied territories when long‐term steady conditions can be assumed. Often, the distributed average CMB variables necessary to calculate have to be estimated from the available variable‐length data series, which may be of suboptimal quality and spatial coverage. This paper explains the use of these data and the reliability of the results in continental Spain, chosen as a large and varied territory. The CMB variables have been regionalized by ordinary kriging at the same 4976 nodes of a 10 km × 10 km grid. Nodal values vary from 14 to 810 mm year^–1, 90% ranging from 30 to 300 mm year^–1. The recharge‐to‐precipitation ratios vary from 0.03 in low‐permeability formations and semiarid areas to 0.65 in some carbonate massifs. Integrated average results for the whole of continental Spain yield a potential aquifer recharge of 64 km³ year^?1, the net recharge over permeable formations (40% of the territory) being 32 km³ year^?1. Two main sources of uncertainty affecting (given by the coefficient of variation, CV), induced by the inherent natural variability of the variables (CV_R) and from mapping (), have been segregated. The average CV_R is 0.13 and could be improved with longer data series. The average is 0.07 and may be decreased with better data coverage. The estimates were compared with other regional and local recharge estimates, being 4% and 1% higher, respectively. Copyright © 2012 John Wiley & Sons, Ltd.