A Numerical Study for Turbulent Flow and Thermal Influence over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier–Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052–2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349–377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369–398)’s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243–258)’s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.     

Alternating skimming flow over a stepped spillway

The study of stepped spillways in laboratory scales has been essentially focused on two separated sub-regimes within skimming flow. In this paper we investigate the appearance of an unclassified alternating skimming flow regime in a 0.5 m wide stepped spillway which does not fit on these earlier definitions, and which does not occur in a 0.3 m wide spillway. Our aim is to explain the genesis of this unclassified flow which is visualised in the physical stepped spillway, by using 3D numerical modelling. Flow depths and velocities are measured using an ultrasonic sensor and Bubble Image Velocimetry in the wider flume (0.5 m). The numerical model is validated with the experimental data from the 0.5 m wide spillway. After validation, the channel width of the same numerical model is reduced to 0.3 m wide spillway in order to characterise (compare) the case without (with) alternating skimming flow. Both cases are solved using Reynolds-Averaged Navier–Stokes equations together with the Volume-of-Fluid technique and SST k-\(\omega\) turbulence model. The experimental results reveal that the alternating skimming flow regime is characterised by an evident seesaw pattern of flow properties over consecutive steps. In turn, the numerical modelling clarified that this seesaw pattern is due to the presence of a complex system of cross waves along the spillway. These cross waves are also responsible for a mass and momentum exchange in the transversal direction and for the formation of the alternating skimming flow in the spillway.     

Part I : <Emphasis Type="Italic">a priori</Emphasis> study of erosion and deposition with large eddy simulation of turbulent flow over multiple 2D sandy Gaussian hills

The micro-scale prediction of sand trapping or take-off over hilly terrains is a crucial issue in semi-arid regions for soil depletion. In this context, large eddy simulations around one or several hills are performed in order to provide statistical parameters to characterize the flow at micro-scales and provide data for mesoscale modelling. We focus on the determination of recirculation zones since they play an important role in solid particle erosion or entrapment. A new wall modeling adapted from Huang et al. (J Turbul 17:1–24, 2016) for rough boundary layers is found to improve the prediction of the recirculation zone length downstream of an isolated hill and is used for all the numerical cases presented here. A geometrical parameterization of the recirculation zones is proposed. When the recirculation region is assumed to have an ellipsoidal shape, the total surface of the recirculation can be obtained from this new parameterization and easily extrapolated to more general dune configurations. Numerical results are compared with experiments performed in our laboratory (Simoëns et al. in Procedia IUTAM 17:110–118, 2015) and good agreement is achieved. We explore general aerodynamic cases deduced from the urban canopy scheme of Oke (Energy Build 11:103–113, 1988). In this scheme the momentum and mass exchange between the upper layer and the space between hills is sorted according to the streamwise hill spacing within three basic cases of skimming, wake or isolated flow. The study of the recirculation zones, the mean velocity and Reynolds stress profiles around an isolated or two consecutive hills with different distances shows that the double hill configuration with 3H separation behaves as much as a whole to the upcoming flow. The vortex formed between the crests does not strongly affect the overall evolution of the outer flow. By an a priori prediction of the preferential zones of erosion and accumulation of fictive particles, it is shown that isolated dunes present more deposition and less erosion than two-hill configurations. The results presented in this study will be discussed in the presence of Lagrangian transport of sand particles above 2D Gaussian hills in future work.     

林窗大小和位置对丽江云杉自然更新幼苗存活和生长的影响

通过 1997～ 2 0 0 0年对 32个林窗的连续观测调查 ,研究了滇西北玉龙雪山自然保护区云杉坪亚高山丽江云杉林林窗大小和林窗位置对自然更新幼苗存活和生长的影响 .4个生长季节的观测结果表明 :林窗与林下非林窗内的幼苗大小和幼苗存活数量差异明显 .林窗由小到大 ,单位面积内的自然更新苗木数量逐渐减少 ,小林窗中更新苗为大林窗的 5倍左右 ,而林下的更新苗很少 ,0 .2ind .10m-2 .中等林窗和小林窗内的幼苗数量在从南到中心到北的位置上几乎没有明显的差异 ;大林窗中存在由南到北的位置差异 ,更新幼苗数量逐渐减少 .从更新幼苗的生长来看 ,中等林窗内的幼苗 ,高度最大、生长最快 ,定居阶段的平均年高生长为 (7.1± 0 .5 )cma-1,小林窗次之 ,大林窗和林下幼苗个体最小 ,生长最慢 .进一步从林窗位置来看 ,中、小林窗幼苗大小几无位置差异 ,大林窗则由南到北 ,幼苗由大变小 .从幼苗存活数量和大小来看 ,中等林窗大小是丽江云杉幼苗更新的适宜面积 ,对该类型退化亚高山针叶林恢复提供了一定的参考 .图 4表 2参 2 3     

Population maintenance of the short-lived shrub Sambucus in a deciduous forest

This study quantitatively clarifies the life history of a shrub, Sambucus racemosa ssp. sieboldiana, in an old-growth forest, the Ogawa Forest Reserve, Japan, by a demographic approach using a projection matrix model that incorporates interactions between demographic parameters and canopy height dynamics. S. racemosa is a common deciduous shrub in central Japan and is known to grow predominantly at forest edges or roadsides. This indicates that it is a highly light-demanding species, and occurrence in gaps in old-growth stands suggests its "fugitive," gap-dependent life history in old-growth forests. We found that one distinctive feature of this species was that its seedlings can survive well in shaded conditions by alternating stems every year like perennial herb species. Matrix model analyses demonstrated that S. racemosa can continuously regenerate under the present disturbance regime of this forest and is highly adaptable to the structural dynamics of the old-growth forest. The maturity of S. racemosa shrubs depends on their size, and nearly all (>90%) of the mature (reproducing) individuals were found in gaps or near gaps. But wide seed dispersal by birds and the ability to form both seed banks and seedling banks, the latter of which has been regarded as a common characteristic of shade-tolerant climax species, probably increase the species' chances to encounter canopy gaps. Dynamic-canopied matrix models showed that the greatest elasticity is with shaded seedling survival. The frequent stem alternation of shaded seedlings often makes the growth rate negative, but the survival rate of seedlings in low light awaiting new gap creation is remarkably high (0.93 yr(-1)). The lower survival rate of the larger individuals and smaller minimum size to start reproduction than other canopy or subcanopy shade-tolerant species indicate that S. racemosa has the potential to reproduce before the closure of the encountered gaps and to complete its life history rapidly.     

Flows across high aspect ratio street canyons: Reynolds number independence revisited

The Reynolds number for flow in a street canyon, Re?=?U_refH/ν (where U_ref is a reference velocity, H the street canyon height, and ν the kinematic viscosity), cannot be matched between reduced-scale experiments and full-scale field measurements. This mismatch is often circumvented by satisfying the Re independence criterion, which states that above a critical Re (Re_c), the flow field remains invariant with Re. Re_c?=?11,000 is often adopted in reduced-scale experiments. In deep street canyons with height-to-width aspect ratio ≥?1.5, reduced-scale experiments have shown two recirculation vortices induced by the mean flows, but full-scale field measurements have observed only one vortex. We investigated this discrepancy by conducting water channel experiments with Re between 10⁴ and 10⁵ at three aspect ratios. The canyons with aspect ratio 1.0 have Re_c?=?11,000, the canyons with aspect ratio 1.5 have Re_c between 31,000 and 58,000, while the canyons with aspect ratio 2.0 have Re_c between 57,000 and 87,000. Therefore, the widely adopted Re_c?=?11,000 is not applicable for canyons with aspect ratio greater than 1.5. Our results also confirm that there is only one vortex in deep canyons at high Re. This single-vortex flow regime could change our fundamental understanding of deep canyons, which are often assumed to exhibit multiple-vortex flow regimes. Applications such as numerical model validation based on the multiple-vortex regime should be revisited. Our experimental data with Re up to 10⁵ could be used to validate numerical models at high Re.     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

长苞铁杉林林隙自然干扰规律

从林隙的大小结构、形成方式及形成木的数量特征等几个方面 ,对长苞铁杉纯林及混交林林隙干扰状况进行了对比研究 .结果表明 :纯林冠空隙 (CG)平均面积为 16 .4 6m2 ,扩展空隙 (EG)平均面积为 96 .98m2 ,混交林CG、EG平均面积分别为 6 9.75m2 和 2 32 .93m2 ;在长苞铁杉纯林中 ,林隙主要以枯立木形成为主 ,占总形成木的一半以上 ,且以单株形成木和多株形成木综合形成 ,形成木径级主要集中在 2 0～ 4 0cm之间 ;混交林林隙以干基折断为主 ,且多由双形成木形成 ,由双形成木形成的林隙占林隙总数的 4 4 % ,形成木以 30～ 4 0cm径级的占最大比例 ,占形成木总数的 4 3.2 4 % ;纯林及混交林形成木高度都主要集中在 2 0～ 30m之间 ;形成木物种组成的比例大小基本上与其主林层树种所占的比例相对应 ;纯林林隙大多处于林隙发育的前、中期而混交林林隙则多集中在林隙发育的中后期 .图 3表 5参 19     

Interaction between flow, transport and vegetation spatial structure

This paper summarizes recent advances in vegetation hydrodynamics and uses the new concepts to explore not only how vegetation impacts flow and transport, but also how flow feedbacks can influence vegetation spatial structure. Sparse and dense submerged canopies are defined based on the relative contribution of turbulent stress and canopy drag to the momentum balance. In sparse canopies turbulent stress remains elevated within the canopy and suspended sediment concentration is comparable to that in unvegetated regions. In dense canopies turbulent stress is reduced by canopy drag and suspended sediment concentration is also reduced. Further, for dense canopies, the length-scale of turbulence penetration into the canopy, δ _e , is shown to predict both the roughness height and the displacement height of the overflow profile. In a second case study, the relation between flow speed and spatial structure of a seagrass meadow gives insight into the stability of different spatial structures, defined by the area fraction covered by vegetation. In the last case study, a momentum balance suggests that in natural channels the total resistance is set predominantly by the area fraction occupied by vegetation, called the blockage factor, with little direct dependence on the specific canopy morphology.     

Revisiting the flow regimes for urban street canyons using the numerical Green’s function

Vortex interactions within a two-dimensional street canyon are analysed using the numerical Green’s function. On account of the inhomogeneity of the domain, vortex interactions are asymmetric: the influence of a street-level vorticity source on the roof-level shear layer differs from that of the latter on the street level. Consequently the magnitudes of the induced vertical velocities are maximised at different aspect ratios. It is argued that the transition from isolated roughness to wake interference is related to the onset of strong long-range interactions while the transition from wake interference to skimming flow is related to the weakening of these interactions. The Green’s function analysis is verified using three-dimensional large-eddy simulations.     

Scalar mixing in an urban canyon

The scalar dynamics within a unit-aspect-ratio street canyon are studied using large-eddy simulation. The key processes of ventilation and mixing are analysed with the canyon-averaged concentration, mean tracer age and variance. The results are sensitive to the source location and can be classified according to the streamline geometry. The canyon-averaged concentrations for the corner vortices, vortex sea and central vortex do not converge to the same value at large times, though the mean decay rates do. The variance measured with respect to the canyon average shows two distinct decay regimes: the early regime reflects large-scale straining and enhanced diffusion across streamlines, while the late regime is associated with escape from the canyon, i.e., ventilation. Analytical predictions for the variance-decay or mixing time scales are verified for the early regime. It is argued that the presence of an open boundary at the roof level suppresses rapid mixing of the scalar field and is responsible for differences with respect to scalar dynamics within closed domains.

     

Modeling the vegetation-atmosphere carbon dioxide and water vapor interactions along a controlled CO2 gradient

Ecosystem functioning is intimately linked to its physical environment by complex two-way interactions. These two-way interactions arise because vegetation both responds to the external environment and actively regulates its micro-environment. By altering stomatal aperture, and therefore the transpiration rate, plants modify soil moisture and atmospheric humidity and these same physical variables, in return, modify stomatal conductance. Relationships between biotic and abiotic components are particularly strong in closed, managed environments such as greenhouses and growth chambers, which are used extensively to investigate ecosystem responses to climatic drivers. Model-assisted designs that account for the physiological dynamics governing two-way interactions between biotic and abiotic components are absent from many ecological studies. Here, a general model of the vegetation-atmosphere system in closed environments is proposed. The model accounts for the linked carbon-water physiology, the turbulent transport processes, and the energy and radiative transfer within the vegetation. Leaf gas exchange is modeled using a carbon gain optimization approach that is coupled to leaf energy balance. The turbulent transport within the canopy is modeled in two-dimensions using first-order closure principles. The model is applied to the Lysimeter CO₂ Gradient (LYCOG) facility, wherein a continuous gradient of atmospheric CO₂ is maintained on grassland assemblages using an elongated chamber where the micro-climate is regulated by variation in air flow rates. The model is employed to investigate how species composition, climatic conditions, and the imposed air flow rate affect the CO₂ concentration gradient within the LYCOG and the canopy micro-climate. The sensitivity of the model to key physiological and climatic parameters allows it to be used not only to manage current experiments, but also to formulate novel ecological hypotheses (e.g., by modeling climatic regimes not currently employed in LYCOG) and suggest alternative experimental designs and operational strategies for such facilities.     

CO2 fluxes near a forest edge: a numerical study

In contrast with recent advances on the dynamics of the flow at a forest edge, few studies have considered its role on scalar transport and, in particular, on CO2 transfer. The present study addresses the influence of the abrupt roughness change on forest atmosphere CO2 exchange and contrasts the concentration and flux fields against those of a uniform forested surface. We use an atmospheric boundary layer two-equation closure model that accounts for the flow dynamics and vertical divergence of CO2 sources/sinks within a plant canopy. This paper characterizes the spatial variation of CO2 fluxes as a function of both sources/sinks distribution and the vertical structure of the canopy. Results suggest that the ground source plays a major role in the formation of wave-like vertical CO2 flux behavior downwind of a forest edge, despite the fact that the contribution of foliage sources/sinks changes monotonously. Such a variation is caused by scalar advection in the trunk space and reveals itself as a decrease or increase in vertical fluxes over the forest relative to carbon dioxide exchange of the underlying forest. The effect was more pronounced in model forests where the leaf area is concentrated in the upper part of the canopy. These results can be useful both for interpretation of existing measurements of net ecosystem exchange of CO2 (NEE) from flux towers in limited fetch conditions and in planning future CO2 transport experiments.     

川西高山森林生态系统林下生物量及其随林窗的变化特征

作为森林生态系统的重要组成部分,林下植被及其残体的分布受到林冠层的影响,但迄今有关林窗对林下植被和残体生物量的影响尚无研究报道.于2013 年8 月2 日至20 日,以海拔3 600 m 的川西岷江冷杉原始林林下植被为研究对象,根据区域内的坡向和林分组成等因素设置3 个100 m×100 m 的典型样地,调查其生物量及其随林窗的变化特征.在每个样地内选择3 个大林窗,在林窗、林缘和林下分别设置3 个20 m×20 m 的样方,调查粗木质残体长度或高度、大小头直径、枯立木记录胸径、腐烂等级等;在林窗、林缘和林下分别设置3 个5 m×5 m 的样方,采用“收获法”收集样方内直径在2.5-10 cm之间的细木质残体和灌木生物量;在林窗、林缘和林下分别设置3 个1 m×1 m 的样方来调查凋落物储量和草本生物量;在1m×1 m 的样方内随机选择1 个20 cm×20 cm 的小样方来调查地被植物生物量.结果表明,（1）川西高山森林生态系统总生物量为72.75 t·hm^-2,其中林下生物量为67.92 t·hm^-2,占生态系统生物量的95.17%.活体植被以灌木为主,其生物量为9.81t·hm^-2;残体部分以粗木质残体为主,其储量为53.00 t·hm^-2;（2）林窗对灌木、草本、地被植物的影响各不相同,且不同物种的灌木生物量表现出不同的分布规律;草本生物量表现出明显的“边缘效应”,在林缘显著高于林下;林窗和林缘的地被植物生物量相对较低;（3）粗木质残体储量从林下到林窗呈现减小的趋势,但总体储量仍然较大,林窗和林缘的细木质残体储量高于林下.这些结果为认识高山森林生态系统林下生物量及其格局,以及林窗在森林生态系统的重要作用提供了基础理论依据.     

The three-dimensionality of obstructed shear flows

Obstructed shear flows (i.e. those over permeable media) are common in the environment. An archetypal example, flow over a submerged vegetation canopy, is investigated here. Like any flow through complex geometry, canopy flows are characterised by strong spatial gradients. The focus of this experimental study is the three-dimensionality of aquatic canopy flow, in particular that of the coherent interfacial vortices that govern mixing into and out of the canopy. It is shown here that the vortices have a finite lateral scale that is comparable to their vertical scale; both are of the order of the drag length scale of the canopy, (C _D a)⁻¹, where a is the frontal area density and C _D is a bulk drag coefficient. The finite lateral extent of the vortices generates strong lateral hydrodynamic gradients, both instantaneously and in the long-term. The instantaneous gradients, which can contribute greatly to the dispersion of dissolved and particulate species, are far more pronounced. Finally, the potential for canopies to generate differential roughness secondary circulation is examined. In the consideration of vertical scalar transport, this circulation can be of the same order as turbulent diffusion.     

Recovery of forest structure and spectral properties after selective logging in lowland Bolivia.

Effective monitoring of selective logging from remotely sensed data requires an understanding of the spatial and temporal thresholds that constrain the utility of those data, as well as the structural and ecological characteristics of forest disturbances that are responsible for those constraints. Here we assess those thresholds and characteristics within the context of selective logging in the Bolivian Amazon. Our study combined field measurements of the spatial and temporal dynamics of felling gaps and skid trails ranging from <1 to 19 months following reduced-impact logging in a forest in lowland Bolivia with remote-sensing measurements from simultaneous monthly ASTER satellite overpasses. A probabilistic spectral mixture model (AutoMCU) was used to derive per-pixel fractional cover estimates of photosynthetic vegetation (PV), non-photosynthetic vegetation (NPV), and soil. Results were compared with the normalized difference in vegetation index (NDVI). The forest studied had considerably lower basal area and harvest volumes than logged sites in the Brazilian Amazon where similar remote-sensing analyses have been performed. Nonetheless, individual felling-gap area was positively correlated with canopy openness, percentage liana coverage, rates of vegetation regrowth, and height of remnant NPV. Both liana growth and NPV occurred primarily in the crown zone of the felling gap, whereas exposed soil was limited to the trunk zone of the gap. In felling gaps >400 m2, NDVI, and the PV and NPV fractions, were distinguishable from unlogged forest values for up to six months after logging; felling gaps <400 m2 were distinguishable for up to three months after harvest, but we were entirely unable to distinguish skid trails from our analysis of the spectral data.     

Gyre formation in open and deep lacustrine embayments: the example of Lake Geneva,Switzerland

Numerical simulations were carried out to investigate gyres within open lacustrine embayments subjected to parallel-to-shore currents. In such embayments, gyre formation occurs due to flow separation at the embayment’s upstream edge. High momentum fluid from the mixing layer between the embayment and offshore flows into the embayment and produces recirculating flow. Systematic numerical experiments using different synthetic embayment configurations were used to examine the impact of embayment geometry. Geometries included embayments with different aspect ratios, depths and embayment corner angles. The magnitudes of the recirculation and turbulent kinetic energy (TKE) in the embayment vary significantly for angles in the range 40°–55°. Embayments with corner angles less than 50° have much stronger recirculation and TKE, other parameters remaining the same. The numerical findings are consistent with gyre formation observed in two embayments located in Lake Geneva, Switzerland, and thus help explain flow patterns recorded in lacustrine shoreline regions.     

Oxygen microoptodes: a new tool for oxygen measurements in aquatic animal ecology

We describe two applications of a recently introduced system for very precise, continuous measurement of water oxygen saturation. Oxygen microoptodes (based on the dynamic fluorescence quenching principle) with a tip diameter of ~50 µm, an eight-channel optode array, an intermittent flow system, and online data registration were used to perform two types of experiments. The metabolic activity of Antarctic invertebrates (sponges and scallops) was estimated in respiration experiments, and, secondly, oxygen saturation inside living sponge tissue was determined in different flow regimes. Even in long-term experiments (several days) no drift was detectable in between calibrations. Data obtained were in excellent correspondence with control measurements performed with a modified Winkler method. Antarctic invertebrates in our study showed low oxygen consumption rates, ranging from 0.03-0.19 cm³ O₂ h^-1 ind.^-1. Oxygen saturation inside living sponge specimens was affected by flow regime and culturing conditions of sponges. Our results suggest that oxygen optodes are a reliable tool for oxygen measurements beyond the methodological limits of traditional methods.     

Stratified flow interactions with a suspended canopy

Field observations of the interactions between a stratified flow and a canopy suspended from the free surface above a solid boundary are described and analysed. Data were recorded in and around the canopy formed by a large long-line mussel farm. The canopy causes a partial blockage of the water flow, reducing velocities in the upper water column. Deceleration of the approaching flow results in a deepening of isopycnals upstream of the canopy. Energy considerations show that the potential for an approaching stratified flow to be diverted beneath a porous canopy is indicated by a densimetric Froude number. Strong stratification or low-velocities inhibit vertical diversion beneath the canopy, instead favouring a horizontal diversion around the sides. The effect on vertical mixing is also considered with a shear layer generated beneath the canopy and turbulence generated from drag within the canopy. In the observations, stratification is shown to be of sufficient strength to limit the effectiveness of the first mixing process, while the turbulence within the canopy is likely to enhance vertical exchange. Velocity and temperature microstructure measurements are used to investigate the effect of the canopy on turbulent dissipation and show that dissipation is enhanced within the canopy.     

Soil drying in a tropical forest: Three distinct environments controlled by gap size

Soil water and temperature regimes in the tropical moist forest on Barro Colorado Island, Panama, were simulated directly from meteorological data using the model SWEAT. Separate field observations from root-exclusion, litter-removal and control treatments in one small and one large forest gap were used for calibration and validation. After irrigating all treatments to field capacity, soil matric potential and temperature were measured over 17 days at four depths ≤50 mm using the filter-paper technique and bead thermistors. Understorey environments were also simulated under the same initial conditions. The results suggest that three distinct scenarios, controlled by gap size, describe how the above- and below-ground processes controlling soil drying are coupled: (1) in the large gap, root water extraction by surrounding trees is negligible so soil drying is dominated by evaporation from the soil surface. Soil temperature is dominated by direct solar heating and cooling due to evaporation. (2) In the small gap, root water extraction dominates soil drying with soil evaporation playing a minor role. Soil temperature is still dominated by direct sunlight with some cooling due to evaporation. (3) In the understorey, root water extraction dominates soil drying. Soil temperature is dominated by heat conduction from deep soil layers with some evaporation and sensible heat transfer. The contrasting soil drying regimes imposed by variation in canopy structure enhance micro-environmental heterogeneity and the scope for differential germination and seedling establishment in coexisting tropical tree species.