Estimation of urban subtropical bahiagrass (Paspalum notatum) evapotranspiration using crop coefficients and the eddy covariance method

Eddy covariance (EC) and micro‐meteorological data were collected from May 2010 to January 2013 from urban, non‐irrigated bahiagrass (Paspalum notatum) in subtropical south Florida. The objectives were to determine monthly crop coefficients (K_c) for non‐irrigated bahiagrass by using EC evapotranspiration (ET) data and the Food and Agriculture Organization 56 Penman–Monteith reference evapotranspiration equation; compare crop ET (ET_c) calculated with new K_c values to ET_c obtained using K_c values available in the literature; and compare results and methodologies for statistical differences. New K_c values ranged from 0.62 to 0.92 and were different from K_c values found in the scientific literature for bahiagrass. Resulting ET_c calculated using literature K_c values were significantly different from EC ET data, whereas ET_c using the new K_c values was not. Specifically, literature K_c values were temporally biased to miscalculate the timing of convergence between potential and actual ET, assuming that our new K_c values calculated with EC methods were most accurate. As a consequence, ET_c calculated using the literature K_c values was either too large or too small. However, one set of literature K_c values from a similar climate and water table depth were closer to our new K_c values, indicating that climate should be considered when selecting urban non‐irrigated K_c from the literature to estimate ET. Results also indicated that more than 1 year of EC ET data was needed when establishing monthly K_c values because of annual variability in factors controlling ET, such as water availability. The new K_c values reported herein could be used as an estimate for urban non‐irrigated bahiagrass within similar climates. Copyright © 2013 John Wiley & Sons, Ltd.     

Effectiveness of introducing crop coefficient and leaf area index to enhance evapotranspiration simulations in hydrologic models

Evapotranspiration (ET) is an essential component of the hydrological cycle and plays a critical role in water resource management. However, ET is often overlooked in order to transform rainfall to runoff for better streamflow simulation. Hydrological models are commonly used to estimate areal actual evapotranspiration (AET) after calibration against observed discharge. However, classical approaches are often inadequate to appropriately simulate other hydrologic components. Hence, it is important to introduce natural heterogeneity to enhance hydrological processes and reduce water balance errors. In this study, the effectiveness of introducing a constant crop coefficient (K_c), flux tower‐based K_c, and leaf area index (LAI) to three hydrological models (Three‐Parametric Hydrologic Model [TPHM], Génie Rural à 4 paramètres Journalier [GR4J], and Catchment hydrologic cycle Assessment Tool [CAT]) is assessed for the simulation of daily streamflow and AET in a mountainous mixed forest watershed (8.54 km²) in South Korea. The results show that the streamflow simulations after introduction of K_c and LAI to the original models are quite similar. However, the effectiveness of the AET estimation was significantly enhanced after introduction of the flux tower‐based K_c and LAI. The information criterion computed to compare the models reveals that the flux tower‐based K_c‐simulated AET demonstrated better agreement with the observed AET. The Pearson's correlation coefficients (R²) of the TPHM (8%), GR4J (55%), and CAT (55%) models also showed improvements that were greater than the constant based K_c simulation. Similarly, the limitations of the three models with respect to capturing seasonal variation as well as high and low flows were enhanced after the introduction of the flux tower‐based K_c, which adequately reproduced hydrological processes with minimum water balance errors and bias. A regression analysis revealed the potential of estimating K_c as a linear function of LAI (R² = 0.84). The results of this study indicate that introduction of K_c and LAI to the conceptual rainfall–runoff models can be considered an effective approach to reduce water balance errors and uncertainties in hydrological models and improve the reliability of climate change studies and water resource management.     

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.     

Evaluation of actual evapotranspiration measured by large-scale weighing lysimeters in a humid alpine meadow,northeastern Qinghai-Tibetan Plateau

The accurate estimation of evapotranspiration (ET) is essential for assessing water availability and requirements of regional-scale terrestrial ecosystems, and for understanding the hydrological cycle in alpine ecosystems. In this study, two large-scale weighing lysimeters were employed to estimate the magnitude and dynamics of actual evapotranspiration in a humid alpine Kobresia meadow from January 2018 to December 2019 on the northeastern Qinghai-Tibetan Plateau (QTP). The results showed that daily ET_a averaged 2.24 ± 0.10 mm day ⁻¹ throughout the study period, with values of 3.89 ± 0.14 and 0.81 ± 0.06 mm day⁻¹ during the growing season and non-growing season, respectively. The cumulative ET_a during the study period was 937.39 mm, exceeding precipitation (684.20 mm) received at the site during the same period by 37%, suggesting that almost all precipitation in the lysimeters was returned to the atmosphere by evapotranspiration. Furthermore, the cumulative ET_a (805.04 mm) was almost equal to the maximum potential evapotranspiration estimated by the FAO-56 reference evapotranspiration (ET₀) (801.94 mm) during the growing season, but the cumulative ET_a (132.25 mm) was 113.72% less than the minimum equilibrium ET_eq) (282.86 mm) during the non-growing season due to the limited surface moisture in frozen soil. The crop coefficient (K_c) also showed a distinct seasonal pattern, with a monthly average of 1.01 during the growing season. Structural equation model (SEM) and boosted regression tree (BRT) show that net radiation and air temperature were the most important factors affecting daily ET_a during the whole study period and growing season, but that non-growing season ET_a was dominated by soil water content and net radiation. The daily K_c was dominated by net radiation. Furthermore, both ET_a and K_c were also affected by aboveground biomass.     

A meteorological-based crop coefficient model for estimation of daily evapotranspiration

Estimation of daily evapotranspiration (ET) over cloudy regions highly desires models which rely on meteorological data only. Notwithstanding, the conventional crop coefficient (K_c) method requires detailed knowledge of geo/biophysical properties of the coupled land-vegetation system, precipitation, and soil moisture. Six Eddy Covariance (EC) towers in Iowa, California and New Hampshire of the USA (covering corn, soybeans, prairie, and deciduous forest) were selected. Investigation on 6 years (2007–2012) 15-min micrometeorological records of these sites revealed that there is an indubitable strong interaction between relative humidity (RH), reference ET (ET_o), and actual ET at different timescales. This allowed to bypass the need for the non-meteorological inputs and express K_c as a second-order polynomial function of RH and ET_o, the ambient regression evapotranspiration model (AREM). The coefficients of the empirical function are crop-specific and may require calibration over different soil types. The mean absolute percentage error (MAPE) of the regression against daily EC observations was 17% during the growing season, and 32% throughout the year with root mean square error (RMSE) of 0.74 mm day⁻¹ and coefficient of determination of 0.71. The model was fully operational (MAPE of 34% and RMSE of 0.82 mm day⁻¹) over the four Iowan sites based on inputs from local weather stations and NLDAS-2 forcing data of NASA. AREM was capable of capturing the dynamics of ET at 15-min and daily timescales irrespective of varying complexities associated with biophysical, geophysical and climatological states.     

Direct and indirect methods to simulate the actual evapotranspiration of an irrigated overhead table grape vineyard under Mediterranean conditions

Two methods, indirect and direct, for simulating the actual evapotranspiration (E) were applied to an irrigated overhead table grape vineyard during summer, situated in the Mediterranean region (south Italy), over two successive years. The first method, indirect but more practical, uses the crop coefficient (K_c) approach and requires determination of the reference evapotranspiration E₀ (FAO (Food and Agriculture Organization) method). This method underestimated on average by 17% the daily values of the actual evapotranspiration E. The analysis in this paper shows that the values of K_c for the table grapes determined by the FAO method seem to not be valid in our experimental conditions. Similar conclusions can be found in the literature for the table grape cultivated under different experimental conditions and using different training systems. The second method, is a direct method for estimating the evapotranspiration. It requires development of a model for the overhead table grape vineyard E, following the Penman–Monteith one‐step approach, and using standard meteorological variables as inputs for the determination of the canopy resistance. This method, which needs a particularly simple calibration, provided a better simulation of the hourly and daily evapotranspiration than the indirect method. In additon, the standard error of the daily values for the direct method ( ± 0 · 41 mm) was about 50% lower than that obtained for the indirect method, also when the indirect method used a locally calibrated coefficient K_c instead of a generic K_c. Both, for practical application and theoretical issues, the advantages and disadvantages linked to the use of each tested method are discussed in detail. Copyright © 2007 John Wiley & Sons, Ltd.     

Accounting for the influence of runoff on event soil erodibilities associated with the EI30 index in RUSLE2

Soil erodibilities (K) associated with the EI₃₀ index vary not only with soil properties but also with soil moisture as it varies in time and space. In Revised Universal Soil Loss Equation Version 2 (RUSLE2), temporal variations in soil erodibility in the USA are calculated using monthly precipitation and temperature as independent variables. K_UM, the soil erodibility factor associated with the Q_REI₃₀ index, varies independently of runoff and the product of K_UM and the runoff ratio for the unit plot (Q_R1) provides an alternative to the temporally varying Ks currently used in predicting storm soil loss in RUSLE2. Comparisons were made between the product of Q_R1 and K_UM and RUSLE2 Ks for representative storms at four locations representing the north to south variation in climate in the USA. Peak erosion associated with the current approach used in RUSLE2 was slightly higher at two locations and slightly lower at the other two locations. One other location, Morris, MN, provided an exception with the peak loss predicted by using the product of Q_R1 and K_UM being 1.7 times that obtained using RUSLE2 Ks. In theory, average annual K_UM values should be better related to soil properties than the average annual values of K frequently used when the average annual values of EI₃₀ are used to predict soil loss. However, work has yet to be performed to determine how K_UM varies directly with soil properties and in space and time. Copyright © 2014 John Wiley & Sons, Ltd.     

Relationships between riparian evapotranspiration and groundwater depth along a semiarid irrigated river valley

Evapotranspiration (ET) from riparian vegetation can be difficult to estimate due to relatively abundant water supply, spatial vegetation heterogeneity, and interactions with anthropogenic influences such as shallower groundwater tables, increased salinity, and nonpoint source pollution induced by irrigation. In semiarid south-eastern Colorado, reliable ET estimates are scarce for the riparian corridor that borders the Arkansas River. This work investigates relationships between the riparian ecosystem along the Arkansas River and an underlying alluvial aquifer using ET estimates from remotely sensed data and modelled water table depths. Results from a calibrated, finite-difference groundwater model are used to estimate weekly water table fluctuations in the riparian ecosystem from 1999 to 2009, and estimates of ET are calculated using the Operational Simplified Surface Energy Balance (SSEBop) model with over 200 Landsat scenes covering over 30 km² of riparian ecosystem along a 70-km stretch of the river. Comparison of calculated monthly SSEBop ET to estimated alfalfa reference ET from local micrometeorological station data indicated statistically significant high linear correspondence (R² = .87). Daily calculated SSEBop ET showed statistically significant moderate linear correspondence with data from a local weighing lysimeter (R² = .59). Simulated monthly SSEBop ET values were larger in drier years compared with wetter years, and ET variability was also larger in drier years. Peak ET most commonly occurred during the month of June for all 11 years of analysis. Relationships between ET and water table depth showed that peak monthly ET was highest when groundwater depths were less than about 3 m, and ET values were significantly lower for groundwater depths greater than 3 m. Negative sample Spearman correlation highlighted riparian corridor locations where ET increased as a result of decreased groundwater depths across years with different hydroclimatic conditions. This study shows how a combination of remotely sensed riparian ET estimates and a regional groundwater model can improve our understanding of linkages between riparian consumptive use and near-river groundwater conditions influenced by irrigation return flow and different climatic drivers.     

Field-aligned conductance values estimated from Maxwellian and kappa distributions in quiet and disturbed events using Freja electron data

We study the question of what difference it makes for the derived field-aligned conductance (K) values if one uses Maxwellian or kappa distributions for the fitting of low-orbiting satellite electron flux spectra in the auroral region. This question has arisen because sometimes a high-energy tail is seen in the spectra. In principle, the kappa fits should always be better, because the kappa distribution is a generalization of the Maxwellian. However, the physical meaning of the parameters appearing in the Maxwellian is clearer. It therefore makes sense to study under which circumstances it is appropriate to use a Maxwellian. We use Freja electron data (TESP and MATE) from two events. One of the events represents quiet magnetospheric conditions (stable arc) and the other represents disturbed conditions (surge). In these Freja events, at least, using kappa rather than Maxwellian fitting gives a better fit to the observed distribution, but the difference in K values is not large (usually less than 20%). The difference can be of either sign. However, sometimes even the kappa distribution does not provide a good fit, and one needs a more complicated distribution such as two Maxwellians. We investigate the relative contributions of the two Maxwellians to the total field-aligned conductance value in these cases. We find that the contribution of the high-energy population is insignificant (usually much less than 20%). This is because K is proportional to n/E_c E_c, where n is the source plasma density and E_c is the characteristic energy.     

A Single-Station Spectral Model of the Monthly Median foF2 and M(3000)F2

A new single-station model (SSM) for monthly median values of the ionospheric parameters foF2 and M(3000)F2 has been developed. Fourier analysis provides a tool for decomposing the time-varying ionospheric parameters. The 12–month smoothed sunspot number R ₁₂ was used as an external solar characteristic because of its availability and predictability. However, for the first time, the solar activity is described not only by R ₁₂ , but also by the linear coefficient K _R representing the tendency of the change of solar activity. A general non-linear approximation of the influence of the solar-cycle characteristics R ₁₂ and K _R and ionospheric parameters foF2 and M(3000)F2 was accepted. The new SSM is applied to several European stations and its statistical evaluation shows better results than the other two SSMs used in the paper. The approach described in the paper does not contradict the use of different synthetic ionospheric indices (as the T-index, MF2–index); the basic aim is to show only that using one additional new characteristic of the solar-cycle variations, such as K _R , improves the monthly median model.     

Impacts of different vegetation restoration options on gully head soil resistance and soil erosion in loess tablelands

Vegetation restoration is identified as an effective approach to control soil erosion and affects soil detachment and resistance to concentrated flow on the Loess Plateau. However, the effects of vegetation restoration at gully heads in loess-tableland remains unclear. This study was performed to investigate the effects of nine vegetation restoration types at gully heads on soil detachment rate (Dr) and soil resistance to concentrated flow (i.e. soil erodibility, K_r and critical shear stress, τ_c). Undisturbed soil samples were collected from nine vegetation-restored lands and one slope cropland (as the control) and were subjected to a hydraulic flume to obtain Dr values of gully heads under six inflow discharges (0.5–3.5 L s^-1). The results showed that the Dr values of nine revegetated gully heads were 77.11% to 95.81% less than that of slope cropland, and the grassland dominated by Cleistogenes caespitosa and the shrubland dominated by Hippophae rhamnoides had a relatively greater decrease in Dr than those of other seven restoration types. The Dr value of nine revegetated gully heads could be better simulated by stream power than by flow velocity and shear stress and was also significantly affected by soil disintegration rate (positively), soil bulk density, saturated hydraulic conductivity, organic matter content, and water-stable aggregate stability (negatively). Additionally, roots with diameters of 0 to 0.5 mm showed a greater effect on Dr than those with larger diameters. Compared to cropland, the nine restored types reduced K_r by 76.26% to 94.26% and improved τ_c by 1.51 to 4.68 times. The decrease in K_r and the increase in τ_c were significantly affected by organic matter content, water-stable aggregate, mean weight diameter of aggregate and root mass density. The combination of grass species (Cleistogenes caespitosa) and shrub (Hippophae rhamnoides) could be considered the best vegetation restoration types for improving soil resistance of gully heads to concentrated flow. © 2019 John Wiley & Sons, Ltd.     

Spatial variation in soil resistance to flowing water erosion along a regional transect in the Loess Plateau

The factors influencing soil erosion may vary with scale. It remains unclear whether the spatial variation in soil erosion resistance is controlled by regional variables (e.g. precipitation, temperature, and vegetation zone) or by local specific variables (e.g. soil properties, root traits, land use, and farming operations) when the study area enlarges from a hillslope or catchment to the regional scale. This study was performed to quantify the spatial variations in soil erosion resistance to flowing water under three typical land uses along a regional transect on the Loess Plateau and to identify whether regional or local specific variables are responsible for these changes. The results indicated that the measured soil detachment capacities (D_c) of cropland exhibited an irregular trend along the regional transect. The D_c of grassland increased with mean annual precipitation, except for two sites (Yijun and Erdos). The measured D_c of woodland displayed an inverted ‘U’ shape. The changes in rill erodibility (K_r) of three land uses were similar to D_c, whereas no distinguishable trend was found for critical shear stress (τ_c). No significant correlation was detected between D_c, K_r and τ_c, and the regional variables. The spatial variation in soil erosion resistance could be explained reasonably by changes in soil properties, root traits, land use, and farming operations, rather than regional variables. The adjustment coefficient of K_r for grassland and woodland could be well simulated by soil cohesion and root mass density (R² = 0.70, P < 0.01), and the adjustment coefficient of critical shear stress could be estimated with aggregate stability (R² = 0.57, P < 0.01). The results are helpful for quantifying the spatial variation in soil detachment processes by overland flow and to develop process‐based erosion model at a regional scale. Copyright © 2015 John Wiley & Sons, Ltd.     

A Practical Measurement Strategy to Estimate Nonlinear Chlorinated Solvent Sorption in Low foc Sediments

In this study, we tested a practical strategy useful for accurate chlorinated volatile organic compound (cVOC) sorption prediction. Corresponding to the feature of the superposition of adsorption due to thermally altered carbonaceous matter (TACM) with organic carbon‐water partitioning, a nonlinear Freundlich sorption isotherm covering a wide range of aqueous concentrations was defined by equilibrium sorption measurement at one or a few low concentration points with extrapolation to the empirical organic carbon‐water partition coefficient (K_oc,e) near compound solubility. We applied this approach to obtain perchloroethene equilibrium sorption isotherm parameters for TACM‐containing glacial sand and gravel subsoil samples from a field site in New York. Sorption and associated K_oc,c applicable to low (5–500 µg/L) and high (>100,000 µg/L) aqueous concentrations were determined in batch experiments. (The K_oc,c is the organic carbon‐normalized sorption partition coefficient corresponding to aqueous concentration C_w.) The K_oc,c measurements at low concentration (~5 µg/L) were 6 to 34 times greater than the K_oc,e. The importance of this type of data is illustrated through presentation of its substantial impact on the site remedy. In so doing, we provide an approach that is broadly applicable to cVOC field sites with similar circumstances (low carbon content glacial sand and gravel with TACM).     

Feasibility of grain-size analysis methods for determination of vertical hydraulic conductivity of streambeds

Accurate estimation of streambed vertical hydraulic conductivity (K_v) is of great importance in the analysis of water quantity exchange and solute transfer between a stream and its sediments. The paper analyzed the inaccuracy of hydraulic conductivity values of sediments derived from grain-size distribution (K_g), which were determined from eight empirical grain-size methods to represent streambed K_v. In this study, the values of K_v for a streambed were derived using falling-head standpipe permeameter tests conducted at eight study sites in the Elkhorn River, Nebraska, and the tested streambed columns were then collected for grain-size analysis by sieving. These empirical methods were used to calculate the K_g values of the streambed from grain-size distribution data of sediments. Unlike many other studies, this study verifies K_g from grain-size distribution with K_v from permeameter tests on the basis of the same samples of streambed sediments. The K_g values derived from the eight empirical methods were larger than the K_v from permeameter tests; there are five methods that give K_g values of about 3–6 times larger than these K_v. The K_g values from the Kozeny formula followed by the Hazen formula give the largest overestimation error if they are used to represent the K_v of the streambed. The USBR and Shepherd formulas generated K_g values close to K_v, but these K_g values are still larger in general than the K_v values. Moreover, the new values of coefficient C for the empirical formulas were revised so that they can be used to calculate the approximate K_v of a streambed. Among the eight methods, the ratios of the original C values to the average new C range from 1.3 to 5.9. It can be hypothesized that smaller C values must be used in the estimation of K_v for general soil samples if these empirical formulas are used to calculate K_v.     

The timing and magnitude of changes to Hortonian overland flow at the watershed scale during the post-fire recovery process

Extreme hydrologic responses following wildfires can lead to floods and debris flows with costly economic and societal impacts. Process-based hydrologic and geomorphic models used to predict the downstream impacts of wildfire must account for temporal changes in hydrologic parameters related to the generation and subsequent routing of infiltration-excess overland flow across the landscape. However, we lack quantitative relationships showing how parameters change with time-since-burning, particularly at the watershed scale. To assess variations in best-fit hydrologic parameters with time, we used the KINEROS2 hydrological model to explore temporal changes in hillslope saturated hydraulic conductivity (K_sh) and channel hydraulic roughness (n_c) following a wildfire in the upper Arroyo Seco watershed (41.5 km²), which burned during the 2009 Station fire in the San Gabriel Mountains, California, USA. This study explored runoff-producing storms between 2008 and 2014 to infer watershed hydraulic properties by calibrating the model to observations at the watershed outlet. Modelling indicates K_sh is lowest in the first year following the fire and then increases at an average rate of approximately 4.2 mm/h/year during the first 5 years of recovery. The estimated values for K_sh in the first year following the fire are similar to those obtained in previous studies on smaller watersheds (<1.5 km²) following the Station fire, suggesting hydrologic changes detected here can be applied to lower-order watersheds. Hydraulic roughness, n_c, was lowest in the first year following the fire, but increased by a factor of 2 after 1 year of recovery. Post-fire observations suggest changes in n_c are due to changes in grain roughness and vegetation in channels. These results provide quantitative constraints on the magnitude of fire-induced hydrologic changes following severe wildfires in chaparral-dominated ecosystems as well as the timing of hydrologic recovery.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

A comparative study on a simple two-parameter monthly water balance model and the Kajiyama formula for monthly runoff estimation

ABSTRACT

A two-parameter monthly water balance model to simulate runoff can be used for a water resources planning programme and climate impact studies. However, the model estimates two parameters of transformation of time scale (c) and of the field capacity (SC) by a trial-and-error method. This study suggests a modified methodology to estimate the parameters c and SC using the meteorological and geological conditions. The modified model is compared with the Kajiyama formula to simulate the runoff in the Han River and International Hydrological Programme representative basins in South Korea. We show that the estimated c and SC can be used as the initial or optimal values for the monthly runoff simulation study in the model.

EDITOR M.C. Acreman; ASSOCIATE EDITOR S. Kanae     

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.     

Determination of Partition Coefficients for Selected PAHs between Water and Dissolved Organic Matter

The fate and transport of highly hydrophobic chemicals are affected by the partitioning between water and dissolved organic carbon. Large variation in the partition coefficient (K_DOCw) is often found, due to the selection of model organic matter or potential experimental artifacts. To investigate the roles of the type of organic matter on the partitioning of highly hydrophobic compounds, the partition coefficients of eight selected polycyclic aromatic hydrocarbons (PAHs), with 3–6 aromatic rings, were determined using a passive dosing/extraction method between water and model dissolved organic matter (humic acid, fulvic acid, sodium dodecyl sulfate micelle (SDS), and 2‐hydroxypropyl‐β‐cyclodextrin). Although the K_DOCw values for 3–4 ring PAHs in this study were close to those reported in the literature, experimental K_DOCw values between Aldrich humic acid (AHA) and water were higher than values reported in the literature for 5–6 ring PAHs. The K_DOCw values were highest for AHA, followed by SDS and Suwannee river fulvic acid (SFA). The slopes of the linear regression between log K_DOCw and log K_ow were 1.23 (± 0.13), 0.82 (± 0.09), and 0.59 (± 0.13) for AHA, SDS, and SFA, respectively. The differences in the K_DOCw values between AHA and the other organic matter (SDS, SFA, and CD) increased with increasing hydrophobicity of the PAHs, showing that the sorption of highly hydrophobic chemicals to the humic acid fraction may be important in the presence of mixed organic matter.     

Application of an S‐shaped curve model to the temporal development of tafoni of salt‐weathering origin

To describe temporal change in tafone development, an S‐shaped curve equation is proposed: Z = Z_c [1 ? (n + 1) exp (? β t ) + n exp (? (1 + 1/n) β t )] , where Z is observed tafone depth, Z_c is ultimate tafone depth, t is time, and n and β are constants. The applicability of this model is examined using tafone data selected from seven sites, which are categorized into three different salt‐weathering environments: a spray/splash‐dominant (occasionally wave‐affected) supra‐tidal zone, aerosol‐affected coastal regions, and inland desert areas. The results indicate that the equation can well describe tafone development in each of these environments. An investigation based on the values of n and β, determined through a best fit of the equation to the data, suggests that n characterizes site‐specific environmental conditions and β reflects the magnitude of factors controlling the recession mechanism of tafone surfaces. It is found that (1) the maximum rate of tafone growth dramatically decreases from supra‐tidal, through coastal, to desert environments, and (2) the growing mode of tafoni is different depending on the environmental settings. The erosional force to facilitate the development of tafoni at supra‐tidal sites is estimated to be about 400 times greater than that in the general coastal area. Copyright © 2011 John Wiley & Sons, Ltd.