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.     

Three new species of Batrachospermum Roth (Batrachosper-maceae,Rhodophyta) in China

1 INTRODUCTION The first published records of genus Batra-chospermum Roth in China are from Jao (1941), when eleven taxa were reported. Since then, Shi, Hu and Kumano (1993), Shi (1994a, 1994b), Hua and Shi (1996), Xie et al. (1999, 2003), Xie and Shi (2003) have described in more detail the species from China. To date, 23 taxa of Batrachospermum Roth in China have been reported, of which, B. yunnanense, B. nothocladoideum and B. transito-rium are new species of this genus, and are…     

Polyclonal antibodies to light-harvesting CHL-protein of PSII (LHC II) in marine green algaeBryopsis corticulans

Polyc·lonal antibodies raised against LHC II isolated from SDS-solubilizedBryopsis corticulans thylakiod membranes by SDS-PAGE, were characterised by double immunodiffusion, Rocket immunoelectrophoresis and antigen-antibody crossed immunoelectro-phoresis assays showed the antibodies had strong cross-reaction with allB. corticulans LHC II components (even with those which were incubated in boiling water) and showed immunological cross-reactivity with LHC II polypeptides of spinach and the marine green algaCodium fragile. The results suggested that LHC II of different species had similar antigenic determinants and also conservation of amino acid sequences of the polypeptides during evolution, and that the antibodies could cross react with apoproteins of D₂ proteins (which contain P₆₈₀) fromB. corticulans, spinach andC. fragile, but not with apoproteins of P₇₀₀ Chl-proteins. Our results indicated some similarities in primary structure between LHC II of different species, and between LHC II and D₂ proteins of marine green algae and spinach. Our finding that D₂ and P₇₀₀ Chl-proteins are not immunologically related suggested that P₇₀₀ Chl-proteins and D₂ proteins pass through independent evolutionary pathways. Contribution No 1717 from the Institute of Oceanology, Academia Sinica, This subject was supported by the National Science Foundation of China.     

Study on the relationship between speciation of heavy metals and their ecotoxicity

Heavy metal is a main pollutant in the marine ecosystem, so study on the effect of heavy metal on phytoplankton is important. Algae (Chaetoceros sp.,Dunaliella sp.,Dicrateria zhanjiangenis Hu. var. sp.) were laboratory cultured to observe the effect of heavy metals on their growth. The effect of different metal ion concentration, the detoxication effect of complexation agents and the growth of algae in different media and different nutrition levels were studied to evaluate the effect of metal speciation. It is proved that trace amount of heavy metals can stimulate the growth of algae cells but that high concentration is lethal. The sequence of toxicity is Cd²⁺>Zn²⁺>Pb²⁺. In ordinary nutrition conditions, the detoxication sequence of complexation agents toChaetoceros sp. is EDTA >sodium salicylate>sodium oxalate >sodium citrate>sulfanilic acid>O-phenanthroline. This is in good conformity with the stability constant sequence of these agents with copper and good evidence that toxicity of metal ion is related to its activity and not to its total concentration. Project supported by the Fund of Excellent Young Teachers of State Educational Commission of China.     

Harmful microalgae blooms (HAB); problematic and conditions that induce them

HAB occurrence is becoming more frequent and problematic in marine recreational waters. However, the exploitation of the coastal area for recreational use is promoting the necessary conditions for the HAB increase. In terms of the harmful effects, we can consider two types of causative organism: the toxic producers and the high-biomass producers. Toxic events can be produced by a very low concentration of the causative organism. This characteristic implies a difficulty for the monitoring programs in relation to human health. It is important to point out in the context of human health and HAB events, that in some coastal regions (e.g. the Mediterranean basin) HABs are an emerging problem. In these regions, the local population and visitors may face a health risk that is difficult to measure. The monitoring of toxic species has mainly been associated -with shellfish farming. However, the risk of intoxication could become even greater in areas not subject to legislation of aquaculture activities.     

一种单层城市冠层模式的建立及数值试验研究

本文在引进先进的城市地表能量平衡方案 (Town Energy Balance, 简称TEB) 的基础上建立了一个单层城市冠层模式, 并对南京市典型居民区1 km2范围内的局地尺度地表能量平衡各分量进行离线模拟, 将模拟结果与同期观测值作了比对, 发现: TEB方案对城市地表能量平衡各分量的模拟效果良好, 而该方案的模拟性能受建筑物表面材料反照率取值的影响较大。在离线研究的基础上, 本文又将TEB方案成功耦合到南京大学区域边界层模式 (NJU－RBLM) 中, 作为该模式的地表能量平衡参数化方案之一, 分别选取该边界层模式中原有的地表能量平衡参数化方案SVAT (Soil－Vegetation－Atmosphere－Transfer model) 和新引入的TEB方案对冬夏两季不同个例进行模拟, 以常规近地面气温观测资料和Landsat卫星观测的地表反照率资料对模拟结果进行比较, 结果表明: TEB方案对原大气边界层模式的模拟效果有明显改善, 对近地面热力场的改善效果尤为明显, 可以很好地模拟出城市冠层中的“陷阱效应”。     

Relating Urban Surface-layer Structure to Upwind Terrain for the Salford Experiment (Salfex)

Profiles of wind and turbulence over an urban area evolve with fetch in response to surface characteristics. Sodar measurements, taken on 22 April 2002 during the Salford Experiment in the UK (Salfex), are here related to upstream terrain. A logarithmic layer up to z = 65m was observed in all half-hour averaged profiles. Above this height the profile showed a different vertical gradient, suggesting a change in surface cover upstream. The drag coefficient varied by a factor of two over only a 20° direction change. Turbulence intensity (σ _x ) for each wind component (x) decreased with height, but the ratio suggested an underestimate of σ _u compared to previous results. Mean urban and suburban cover fraction within the source area for each height decreased sharply between z = 20 and 50m, increasing slightly above. The near-convergence of cover fractions thus occured for source areas of minimum length ≈ 2,200 m. In comparison, the mean length scale of heterogeneity L _P was calculated from surface cover data to be 1,284 m, and the corresponding mean blending height h _b was 175 m. Finally, the mean streamline angle, α, was negative and the magnitude decreased with height. An exponential fit to α for z ≤ 65m gave an e-folding height scale of 159 m. A simple relationship between this height scale and L _P was assumed, giving L _P ≈ 1,080 m, which is in reasonable agreement with the estimate from surface cover type. The results suggest that more emphasis is required on modelling and measuring surface-layer flow over heterogeneous urban canopies.     

Numerical Study of the Response of an Atmospheric Surface Layer to a Spatially Nonuniform Plant Canopy

High-accuracy large-eddy simulations of neutral atmospheric surface-layer flow over a gapped plant canopy strip have been performed. Subgrid-scale (SGS) motions are parameterized by the Sagaut mixed length SGS model, with a modification to compute the SGS characteristic length self-adaptively. Shaw’s plant canopy model, taking the vertical variation of leaf area density into account, is applied to study the response of the atmospheric surface layer to the gapped dense forest strip. Differences in the region far away from the gap and in the middle of the gap are investigated, according to the instantaneous velocity magnitude, the zero-plane displacement, the potential temperature and the streamlines. The large-scale vortex structure, in the form of a roll vortex, is revealed in the region far away from the gap. The nonuniform spatial distribution of plants appears to cause the formation of the coherent structure. The roll vortex starts in the wake of the canopy, and results in strong fluctuations throughout the entire canopy region. Wind sweeps and ejections in the plant canopy are also attributed to the large vortex structure.     

Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies

Modelling the transfer of heat, water vapour, and CO₂ between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V _{cmax, 25}) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed.

“The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy

     

Large-Eddy Simulation of Flows over Random Urban-like Obstacles

Further to our previous large-eddy simulation (LES) of flow over a staggered array of uniform cubes, a simulation of flow over random urban-like obstacles is presented. To gain a deeper insight into the effects of randomness in the obstacle topology, the current results, e.g. spatially-averaged mean velocity, Reynolds stresses, turbulence kinetic energy and dispersive stresses, are compared with our previous LES data and direct numerical simulation data of flow over uniform cubes. Significantly different features in the turbulence statistics are observed within and immediately above the canopy, although there are some similarities in the spatially-averaged statistics. It is also found that the relatively high pressures on the tallest buildings generate contributions to the total surface drag that are far in excess of their proportionate frontal area within the array. Details of the turbulence characteristics (like the stress anisotropy) are compared with those in regular roughness arrays and attempts to find some generality in the turbulence statistics within the canopy region are discussed.