Characterizing Heterogeneity in Infiltration Rates During Managed Aquifer Recharge

Infiltration rate is the key parameter that describes how water moves from the surface into a groundwater aquifer during managed aquifer recharge (MAR). Characterization of infiltration rate heterogeneity in space and time is valuable information for MAR system operation. In this study, we utilized fiber optic distributed temperature sensing (FO‐DTS) observations and the phase shift of the diurnal temperature signal between two vertically co‐located fiber optic cables to characterize infiltration rate spatially and temporally in a MAR basin. The FO‐DTS measurements revealed spatial heterogeneity of infiltration rate: approximately 78% of the recharge water infiltrated through 50% of the pond bottom on average. We also introduced a metric for quantifying how the infiltration rate in a recharge pond changes over time, which enables FO‐DTS to be used as a method for monitoring MAR and informing maintenance decisions. By monitoring this metric, we found high‐spatial variability in how rapidly infiltration rate changed during the test period. We attributed this variability to biological pore clogging and found a relationship between high initial infiltration rate and the most rapid pore clogging. We found a strong relationship (R² = 0.8) between observed maximum infiltration rates and electrical resistivity measurements from electrical resistivity tomography data taken in the same basin when dry. This result shows that the combined acquisition of DTS and ERT data can improve the design and operation of a MAR pond significantly by providing the critical information needed about spatial variability in parameters controlling infiltration rates.     

Rapid nutrient load reduction during infiltration of managed aquifer recharge in an agricultural groundwater basin: Pajaro Valley,California

Artificial recharge of groundwater is an increasingly important method for augmenting groundwater supply and can have a positive or negative influence on the quality of water resources. We instrumented a managed aquifer recharge (MAR) pond in central coastal California to assess how patterns of infiltration and recharge affect the load of nitrate delivered to the underlying aquifer. The concentration of nitrate in infiltrating water consistently decreased during passage through the first metre of subsurface soils. Enrichment of ¹⁸O and ¹⁵ N in the residual nitrate in infiltrating water proceeded in a ratio of 1:2, indicating that denitrification plays a significant role in the quantitative reduction of nutrients exported during infiltration through shallow soils. The extent and rate of nitrate removal was spatially and temporally variable across the bottom of the recharge pond, with 30% to 60% of the nitrate load being removed over the first 6 weeks of managed aquifer recharge operation. During the period of highest N loading to the system, when the average infiltration rate was > 1 m/day, the recharge pond achieved a load reduction efficiency of 7 kg NO₃^?‐N/day/ha, which compares favourably to nitrate load reductions achieved by treatment wetlands. Groundwater mounding and water composition below the recharge pond suggest that recharge and subsequent lateral transport occur heterogeneously in the underlying aquifer. Nitrate concentrations in the aquifer following infiltration were lowered primarily by dilution, with little evidence for additional denitrification occurring in the aquifer in comparison to high rates documented during shallow infiltration. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupled infiltration and filtration behaviours of concrete porous pavement for stormwater management

Cementitious porous pavement (CPP) is a structural low‐impact development material for rainfall–runoff management. Both infiltration and filtration are critical functions for CPP stormwater quality and quantity control. In this study, three groups of CPP specimens exposed to rainfall–runoff for 4 years and experienced with different maintenance intervals (6, 12 and 48 months, respectively) were used to examine CPP infiltration and filtration performance. Particle mass strained on CPP surface, saturated infiltration rate I_f, temporal infiltration rate I(t), suspended sediment concentration (SSC) and turbidity (τ) were measured to evaluate the process of filtration/infiltration. I(t), SSC and τ were examined less than 50 mg/l of the suspended particle loading. It was found that the CPP surface cleaning methods used in the past 4 years, namely, high pressure wash followed by vacuuming with one atmosphere (100 kPa), were effective, and a 12‐month maintenance interval was verified suitable to maintain the pore structure an acceptable infiltration rate for stormwater management. It was also found that CPP infiltration and filtration process affect each other, and the two properties are coupled in urban stormwater quality and quantity control. On the basis of the experimental measurements, the temporal infiltration rate of the cleaned CPP under a certain particle loading could be simulated by a first‐order nonlinear rational model, and effluent turbidity–SSC relationship was found following a power law. Copyright © 2012 John Wiley & Sons, Ltd.     

Transitions to chaotic thermal convection in a rapidly rotating spherical fluid shell

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T _a =10⁶ and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R _c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R _c , R = 1.7 R _c , R = 2R _c , and R = 5 R _c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R _c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R _c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R _c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.     

Spatial and temporal infiltration dynamics during managed aquifer recharge

Natural groundwater recharge is inherently difficult to quantify and predict, largely because it comprises a series of processes that are spatially distributed and temporally variable. Infiltration ponds used for managed aquifer recharge (MAR) provide an opportunity to quantify recharge processes across multiple scales under semi-controlled conditions. We instrumented a 3-ha MAR infiltration pond to measure and compare infiltration patterns determined using whole-pond and point-specific methods. Whole-pond infiltration was determined by closing a transient water budget (accounting for inputs, outputs, and changes in storage), whereas point-specific infiltration rates were determined using heat as a tracer and time series analysis at eight locations in the base of the pond. Whole-pond infiltration, normalized for wetted area, rose rapidly to more than 1.0 m/d at the start of MAR operations (increasing as pond stage rose), was sustained at high rates for the next 40 d, and then decreased to less than 0.1 m/d by the end of the recharge season. Point-specific infiltration rates indicated high spatial and temporal variability, with the mean of measured values generally being lower than rates indicated by whole-pond calculations. Colocated measurements of head gradients within saturated soils below the pond were combined with infiltration rates to calculate soil hydraulic conductivity. Observations indicate a brief period of increasing saturated hydraulic conductivity, followed by a decrease of one to two orders of magnitude during the next 50 to 75 d. Locations indicating the most rapid infiltration shifted laterally during MAR operation, and we suggest that infiltration may function as a "variable source area" processes, conceptually similar to catchment runoff.     

Relationships between correlation lengths and integral scales for covariance models with more than two parameters

In geostatistical applications, the terms correlation length and range are often used interchangeably and refer to a characteristic covariance length ξ that normalizes the lag distance in the variogram or the covariance model. We present equations that strictly define the correlation length (r _c ) and integral range (ℓ _c ). We derive analytical expressions for r _c and ℓ _c of the Whittle–Matérn, fluctuation gradient curvature and rational quadratic covariances. For these covariances, we show that the correlation length and integral range for a given model are not fully determined by ξ. We define non-trivial covariance functions, and we formulate an ergodicity index based on ℓ _c . We propose using the ergodicity index to compare coarse-grained measures corresponding to non-trivial covariance functions with different parameters. Finally, we discuss potential applications of the proposed covariance models in stochastic subsurface hydrology.     

Field investigation of a dry detention pond with underground detention storage

Abstract

A dry pond is an urban drainage component designed to temporarily store stormwater runoff and to encourage infiltration of surface water to the subsurface layer. This paper investigates field measurement of a dry pond at Taiping Health Clinic, Perak, Malaysia that has been functioning well for five years. The pond has a surface area of 195 m², maximum depth of 32 cm, and a storage capacity of 31.88 m³. The study focused on the infiltration functionality of the constructed dry pond and the results show that it has an average infiltration rate of 125 mm/h and dries up in 330 min after being filled to a depth of 31 mm. A public-domain hydrological model was then employed to simulate hydrographs of ponding and draining, the results of which matched observations with 86–98% accuracy. These results can lead to better understanding of the system and allow duplication of such a drainage design elsewhere.

Editor D. Koutsoyiannis

Citation Lai, S.H. and Mah, D.Y.S., 2012. Field investigation of a dry detention pond with underground detention storage. Hydrological Sciences Journal, 57 (6), 1249–1255.     

Impact of Urban Surroundings on Some Limnological Characteristics of a Tropical Aquatic Ecosystem

The water quality of an urban pond in the thickly populated area of Varanasi city (5 km apart) was studied and compared with a rural pond in the Banaras Hindu University campus for transparency, conductivity and nutrient richness (Cl^?, SO, PO? P, NO? N, organic carbon, Ca²⁺, Mg²⁺, K⁺, Na⁺) at three depths (surface, 1.5 m, 3 m) at monthly intervals between February 1982 and February 1983. This was done to assess the effects of urban surroundings of a very ancient city sector on pond water quality in reference to that of a rural pond. The rural pond had a lush growth of 12 macrophytic species, whereas the urban one had only such a growth with many phytoplanktonic species. Transparency was maximum in the winter season and the rural pond water was more transparent, while the electrolytical conductance was maximum in the rainy season, being higher in the urban pond. Electrolytical conductivity was negatively correlated to transparency: urban: EC = 1081.612–6.575 T, r² = 0.897, F_1,11 = 96, P <0.005; rural: EC = 728.981–4.328 T, r₂ = 0.892, F_1,11 = 91, P <0.005. Chloride and sulphate concentrations were highest in summer months, but the former was much higher in the urban pond while the latter in the rural pond. NO₃–N was highest in the rainy season in the rural pond and in early winter in the urban one and showed a definite trend with change in depth. PO₄–P also varied with depth and time and it was higher in late summer and the early rainy season in the rural pond and in early winter in the urban pond. But both these nutrients were much higher in the urban pond. The maximum organic carbon concentration was found in the rainy season in the rural pond and in summer months in the urban pond. The variation of organic carbon with depth was distinct. Both summer and winter seasons showed almost similar values of calcium concentration in the rural pond, but in the urban pond it was maximum in summer. Organic carbon and calcium were higher in the urban pond. The magnesium concentration was highest in rainy months in both the ponds, but the periodicity of the minimum differed. The distribution of calcium with depth was not well defined. The highest concentration of potassium was found in the winter season in both the ponds. The sodium concentration in the rural pond was observed maximum in summer and minimum in the rainy season, but in the urban pond the trend was different. The variation of potassium and sodium with depth was not well defined. Magnesium and sodium were also higher in the urban pond but potassium was almost at the same concentration in both the ponds. The effect of urbanisation may be one of the factors which might be responsible to the shift of the species composition towards phytoplanktonic flora.     

Sediment sources and transport pathways in a rural catchment,Herefordshire, UK

The relative significance of primary and secondary suspended sediment sources has been identified in a small (15 km²) agricultural catchment through the application of a linear unmixing model using mineral magnetic, geochemical and radionuclide signatures of fine sediments. Roads appear to be important as a secondary source of suspended sediment and in the conveyance of topsoils to the river channel. Roads also alter the source signatures of sediment in transport, especially that derived from agricultural topsoils. Modelling suggests that c. 30% of the suspended sediment collected from suspended sediment traps in the River Leadon were derived from, or transported via, roads. Relative contributions from primary sources of c. 43% and c. 27% were estimated from a combined (subsoil and channel bank) subsurface source, and a combined (arable and grassland) topsoil source, respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

Infiltration model in sloping layered soils and guidelines for model parameter estimation

A Green-Ampt type model for sloping layered soils (GASLS) was developed to investigate infiltration processes. We introduced a factor c, which is the same for all layers and represents the ratio of effective hydraulic conductivity over saturated hydraulic conductivity. Guidelines to estimate the factor c were established based on 234 scenarios under various conditions. The model with the estimated factor c can describe infiltration processes better than that with c = 1. For fine soils, or layered formations with finer soils on the top, c is smaller than 1. The factor c for coarse soils, or layer formations with coarse soils on the top is close to 1. Comparison with laboratory experiments on a sloping surface indicated that the GASLS model with a slope factor that is adjusted by the sine of the slope angle can represent the sloping surface effects. The GASLS model can incorporate any slope factor.     

Variations of crop coefficient and its influencing factors in an arid advective cropland of northwest China

Understanding the variation and magnitude of crop coefficient (K_c) is important for accurate determination of crop evapotranspiration and water use. In this study, we calculated K_c in an irrigated maize field with ground mulching by eddy covariance evapotranspiration measurements during the whole growing periods in 2009 and 2010 in an arid region of northwest China. A semi‐empirical practical approach for estimating K_c was proposed by introducing the dynamic fraction of canopy cover and incorporating the effect of leaf senescence as a function of days after sowing. The contribution of arid advection of sensible heat resulting from irrigation to K_c and the response of K_c to canopy conductance (G_c) were investigated. The averaged values of daily K_c were lower than typical values obtained previously without mulching due to decreasing effect of mulching on K_c, with 0.82 and 0.80 for the 2 years, respectively. The maximum average K_c occurred at the heading stage, with 1.21 and 1.04 for the 2 years, respectively. The difference of K_c was attributed to the difference of leaf area index. The semi‐empirical practical approach could well estimate the variations of K_c, thus could be a robust and useful tool for the practical users and water managers. The contributions to daily K_c from the arid advection were 4.4–28.0% of the measured K_c. The G_c had stronger control on daily K_c at the early and later stages than at the middle stage. When G_c, leaf area index and relative soil extractable water were lower than the respective threshold values of 20 mm s^?1, 3.0 m² m^?2 and 0.5, the daily K_c increased significantly with the increase of the three factors, and almost remained constant when the three factors were beyond the threshold values. These results are helpful for quantifying contributions of individual factors to K_c, and subsequently improving water management practices according to K_c. Copyright © 2014 John Wiley & Sons, Ltd.     

The spatial heterogeneity of riverbed saturated permeability coefficient in the lower reaches of the Heihe River Basin,Northwest China

In arid region, direct infiltration from rainfall contributes little to groundwater compared with localized recharge from streams. How to quantify riverbed infiltration to groundwater systems is an important area of research in hydrology. In this study, saturated permeability coefficient of a riverbed in an arid inland river basin located in the northwest of China was obtained by Guelph Permeameter and laboratory analysis methods. The characteristics of riverbed infiltration and its spatial patterns were analysed using geostatistical method and kriging method. The results showed that the saturated permeability coefficient varied from 0.089 to 2.802 m/d, indicating moderate degree of variability. The Guelph Permeameter and laboratory test methods provided consistent estimates of saturated permeability coefficient. There was a strong spatial correlation for K_fs of the riverbed in this study area when Range (A) was less than 0.276°, suggesting that the maximum sampling distance for saturated permeability coefficient of the riverbed was 0.276° under isotropic conditions. The K_fs near the centre of the riverbed was higher than the value near riverbank. The K_fs values decreased in the direction of upstream to downstream in the Heihe River Basin. The riverbed mechanical composition, initial soil water content and bulk density have significant influence up on the riverbed infiltration. Besides, the topographical factors including the width, altitude and distance factors of the riverbed together impacted the riverbed infiltration and the slope of the riverbed and also influenced the riverbed infiltration. Copyright © 2015 John Wiley & Sons, Ltd.     

Statistical Detection and Characterization of a Deviation from the Gutenberg-Richter Distribution above Magnitude 8

— We present a quantitative statistical test for the presence of a crossover c₀ in the Gutenberg-Richter distribution of earthquake seismic moments, separating the usual power-law regime for seismic moments less than c₀ from another faster decaying regime beyond c₀. Our method is based on the transformation of the ordered sample of seismic moments into a series with uniform distribution under condition of no crossover. A simulation method allows us to estimate the statistical significance of the null hypothesis H₀ of an absence of crossover (c₀=infinity). When H₀ is rejected, we estimate the crossover c₀ using two different competing models for the second regime beyond c₀ and the simulation method. For the catalog obtained by aggregating 14 subduction zones of the Circum-Pacific Seismic Belt, our estimate of the crossover point is log(c₀)=28.14 ± 0.40 (c₀ in dyne-cm), corresponding to a crossover magnitude m_W=8.1 ± 0.3. For separate subduction zones, the corresponding estimates are substantially more uncertain, so that the null hypothesis of an identical crossover for all subduction zones cannot be rejected. Such a large value of the crossover magnitude makes it difficult to associate it directly with a seismogenic thickness as proposed by many different authors. Our measure of c₀ may substantiate the concept that the localization of strong shear deformation could propagate significantly in the lower crust and upper mantle, thus increasing the effective size beyond which one should expect a change of regime.     

Biological and Physical Clogging in Infiltration Wells: Effects of Well Diameter and Gravel Pack

Gravity-driven infiltration into the shallow subsurface via small-diameter wells (SDWs), i.e., wells with an inner diameter smaller than 7.5 cm (3 inches) and no gravel pack) has proven to be a cost-efficient and flexible tool for managed aquifer recharge (MAR), as it provides relatively high recharge rates with minimal construction effort. SDWs have a significantly smaller open filter area than larger diameter wells with gravel pack, making the infiltration of low-quality waters through these wells more at risk clogging. To investigate their susceptibility for biological and physical clogging, 24 physical models with different well setups were evaluated by infiltrating either nutrient-poor but turbid water or nutrient-rich but clear water. The experiments showed that smaller diameters and the lack of a gravel pack increase the well's susceptibility to both kinds of clogging. However, this effect was observed to be much more pronounced for physical than for biological clogging. Our conclusion is that SDWs show severe disadvantages with respect to the infiltration of highly turbid waters in comparison to large diameter wells with a gravel pack. Nevertheless, this disadvantage is much less severe when it comes to the infiltration of clear but nutrient-rich waters (e.g., treated wastewater). Depending on the economic and geological circumstances of a MAR-project, this disadvantage could be outweighed by the significantly lower construction costs of SDWs.     

Seasonal variation of infiltration rates through pond bed in a managed aquifer recharge system in St-André, Belgium

In Belgium, IWVA uses managed aquifer recharge (MAR) to recharge the aquifer with treated wastewater generated from the communities to sustain the potable water supply on the Belgian coast. This MAR facility is faced with a challenge of reduced infiltration rates during the winter season when pond water temperatures near 4°C. This study involves the identification of the predominant factor influencing the rate of infiltration through the pond bed. Several factors, including pumping rates, natural recharge, tidal influences of the North Sea and pond-water temperature, were identified as potential causes for variation of the recharge rate. Correlation statistics and linear regression analysis were used to determine the sensitivity of the infiltration rate to the aforementioned factors. Two groundwater flow models were developed in visual MODFLOW to simulate the water movement under the pond bed and to obtain the differences in flux to track the effects of variation of hydraulic conductivity during the two seasons. A 32% reduction in vertical hydraulic gradient in the top portion of the aquifer was observed in winter, causing the recharge rates to fluctuate. Results showed that water temperature caused a 30% increase in hydraulic conductivity in summer as compared with winter and has the maximum impact on infiltration rate. Cyclic variations in water viscosity, occurring because of seasonal temperature changes, influence the saturated hydraulic conductivity of the pond bed. Results from the models confirm the impact on infiltration rate by temperature-influenced hydraulic conductivity.     

Calibration and validation of a general infiltration model

A general infiltration model proposed by Singh and Yu (1990) was calibrated and validated using a split sampling approach for 191 sets of infiltration data observed in the states of Minnesota and Georgia in the USA. Of the five model parameters, f_c (the final infiltration rate), S_o (the available storage space) and exponent ‘n’ were found to be more predictable than the other two parameters: m (exponent) and a (proportionality factor). A critical examination of the general model revealed that it is related to the Soil Conservation Service (1956) curve number (SCS‐CN) method and its parameter S_o is equivalent to the potential maximum retention of the SCS‐CN method and is, in turn, found to be a function of soil sorptivity and hydraulic conductivity. The general model was found to describe infiltration rate with time varying curve number. Copyright © 1999 John Wiley & Sons, Ltd.     

The vertical heterogeneity of soil detachment by overland flow on the water-level fluctuation zone slope in the Three Gorges Reservoir,China

Periodic submersion and exposure due to the operation of the Three Gorges Reservoir (TGR) alter the soil properties and plant characteristics at different elevations within the water level fluctuation zone (WLFZ), possibly influencing the soil detachment capacity (D_c), but the vertical heterogeneity of this effect is uncertain. Soil samples were taken from 6 elevation segments (5 m per segment) along a slope profile in the WLFZ of the TGR to clarify the vertical heterogeneity of D_c. Scouring experiments were conducted at 5 slope gradients (17.6%, 26.8%, 36.4%, 46.6%, and 57.7%) and 5 flow rates (10, 15, 20, 25, and 30 L min⁻¹) to determine D_c. The results indicate that the soil properties and biomass parameters of the WLFZ exhibit strongly vertical heterogeneity. D_c fluctuates with increasing elevation, with maximum and minimum average values at elevations of 145–150 m and 165–170 m, respectively. Linear equations accurately describe the relationships between D_c and hydrodynamic parameters, for which the shear stress (τ), stream power (ω), and unit energy of water-carrying section (E) perform much better than the unit stream power (U). Furthermore, a clear improvement is achieved when using a general index of flow intensity to estimate D_c. Furthermore, D_c is significantly and negatively correlated with the mean weight diameter (MWD, p < 0.05) and organic matter content (p < 0.01) but not significantly correlated with other soil properties (p > 0.05). The rill erodibility at elevations of 145–150 m and 170–175 m is greater than that at other elevations. The critical hydraulic parameters were highest in the 165–170 m segments. Both the rill erodibility and the critical parameters fluctuate vertically along the sloping surface. This research highlights the vertical heterogeneity of D_c and is helpful for better understanding the mechanisms responsible for soil detachment in the WLFZ of the TGR.     

Numerical Modeling of Unsaturated Flow in Wastewater Soil Absorption Systems

It is common practice in the United States to use wastewater soil absorption systems (WSAS) to treat domestic wastewater. WSAS are expected to provide efficient, long-term removal of wastewater contaminants prior to ground water recharge. Soil clogging at the infiltrative surface of WSAS occurs due to the accumulation of suspended solids, organic matter, and chemical precipitates during continued wastewater infiltration. This clogging zone (CZ) creates an impedance to flow, restricting the hydraulic conductivity and rate of infiltration. A certain degree of clogging may improve the treatment of wastewater by enhancing purification processes, in part because unsaturated flow is induced and residence times are significantly increased. However, if clogging becomes excessive, the wastewater pond height at the infiltrative surface can rise to a level where system failure occurs. The numerical model HYDRUS-2D is used to simulate unsaturated flow within WSAS to better understand the effect of CZs on unsaturated flow behavior and hydraulic retention times in sandy and silty soil. The simulations indicate that sand-based WSAS with mature CZs are characterized by a more widely distributed flow regime and longer hydraulic retention times. The impact of clogging on water flow within the silt is not as substantial. For sand, increasing the hydraulic resistance of the CZ by a factor of three to four requires an increase in the pond height by as much as a factor of five to achieve the same wastewater loading. Because the degree of CZ resistance directly influences the pond height within a system, understanding the influence of the CZ on flow regimes in WSAS is critical in optimizing system design to achieve the desired pollutant-treatment efficiency and to prolong system life.     

A simple parameterization of bulk canopy resistance from climatic variables for estimating hourly evapotranspiration

This paper examines a model for estimating canopy resistance r_c and reference evapotranspiration ET_o on an hourly basis. The experimental data refer to grass at two sites in Spain with semiarid and windy conditions in a typical Mediterranean climate. Measured hourly ET_o values were obtained over grass during a 4 year period between 1997 and 2000 using a weighing lysimeter (Zaragoza, northeastern Spain) and an eddy covariance system (Córdoba, southern Spain). The present model is based on the Penman–Monteith (PM) approach, but incorporates a variable canopy resistance r_c as an empirical function of the square root of a climatic resistance r* that depends on climatic variables. Values for the variable r_c were also computed according to two other approaches: with the r_c variable as a straight‐line function of r* (Katerji and Perrier, 1983, Agronomie 3 (6): 513–521) and as a mechanistic function of weather variables as proposed by Todorovic (1999, Journal of Irrigation and Drainage Engineering, ASCE 125 (5): 235–245). In the proposed model, the results show that r_c/r_a (where r_a is the aerodynamic resistance) presents a dependence on the square root of r*/r_a, as the best approach with empirically derived global parameters. When estimating hourly ET_o values, we compared the performance of the PM equation using those estimated variable r_c values with the PM equation as proposed by the Food and Agriculture Organization, with a constant r_c = 70 s m^?1. The results confirmed the relative robustness of the PM method with constant r_c, but also revealed a tendency to underestimate the measured values when ET_o is high. Under the semiarid conditions of the two experimental sites, slightly better estimates of ET_o were obtained when an estimated variable r_c was used. Although the improvement was limited, the best estimates were provided by the Todorovic and the proposed methods. The proposed approach for r_c as a function of the square root of r* may be considered as an alternative for modelling r_c, since the results suggest that the global coefficients of this locally calibrated relationship might be generalized to other climatic regions. It may also be useful to incorporate the effects of variable canopy resistances into other climatic and hydrological models. Copyright © 2005 John Wiley & Sons, Ltd.     

Growing Spirodela polyrrhiza in Swine Wastewater for the Production of Animal Feed and Fuel Ethanol: A Pilot Study

To evaluate the performance of Spirodela polyrrhiza grown in swine wastewater for protein and starch production under field conditions, a pilot‐scale duckweed culture pond was installed at Barham Farm, Zebulon, North Carolina and operated from May to November 2010. The anaerobically treated swine wastewater was fed to the duckweed pond intermittently to provide nutrients for the growth of duckweed, and the duckweed biomass was harvested regularly from the pond and prepared as a protein‐ or starch‐rich feedstock for the production of animal feed or fuel ethanol. Over the experimental period, the duckweed pond produced protein and starch at rates of 2.68 and 1.88 g m^?2 day^?1, respectively. During the same time, NH₄–N and o‐PO₄–P in the wastewater were, respectively, removed at rates of 92.9 and 2.90 mmol m^?2 day^?1.