Научная статья на тему 'Impact of the Dhofar cloud forest on subsurface hydrological fluxes: modelling on the scale of a near-root zone soil continuum'

Impact of the Dhofar cloud forest on subsurface hydrological fluxes: modelling on the scale of a near-root zone soil continuum Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
96
39
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЛИСТОПАДНЫЙ ТУМАННЫЙ ЛЕС / ЗАСУШЛИВЫЙ/ПОЛУЗАСУШЛИВЫЙ РЕГИОН / ARID/SEMI-ARID REGION / СКЛОНОВАЯ РАСТИТЕЛЬНОСТЬ / ОРОГРАФИЧЕСКИЙ ДОЖДЬ / ЗАКОН ДАРСИ / МОДЕЛЬ ГРИНА-АМПТА / ЗАБОР ВЛАГИ КОРНЯМИ РАСТЕНИЙ / DECIDUOUS CLOUD FOREST / HILLSLOPE VEGETATION / OROGRAPHIC RAIN / ROOT WATER UPTAKE / INFILTRATION / DARCY LAW / GREEN-AMPT MODEL

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Касимов Анвар Рашидович, Хильдебрандт Анке, Обносов Юрий Викторович

Рассматривается склоновый лес полуаридной зоны в Омане, особенностью которого является перехват кроной деревьев жидкости (тумана и осадков слабой интенсивности) из муссонного ветрa в течение трех месяцев в году. Эта вода капает с листвы на землю по площади проекции кроны и стекает по стволу, в результате чего влага аккумулируется в прикорневой зоне почвы и используется деревом в течение длительного постмуссонного периода. За счет капиллярности и гравитации почвенная вода дрейфует к зеркалу грунтовых вод. Предложена модифицированная математическая модель Грина-Ампта для вертикального влагопереноса в зоне питания дерева, которая учитывает циклостационарность процесса «капельного самоорошения», неравномерный по глубине перехват влаги распределенными корнями дерева и возможность остановки движения воды за счет неодинаковых капиллярных давлений в зоне промачивающего и дренирующего фронтов.A slope forest in a semi-arid region of

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Oman is studied. The peculiarity of this ecosystem is an interception of the liquid (fog and drizzle) by the tree canopy from a monsoon wind during a three-month monsoon period. The intercepted water drips on the soil from the foliage and stem over the area of the projection of the canopy. The moisture accumulates in the soil profile and is used by the tree long after the monsoon. Capillarity and gravity drives moisture to the water table. A modified Green-Ampt model of vertical infiltration is suggested, with a cyclostationary «self-drip-irrigation», vertically non-uniform water uptake by the tree roots and the possibility of stopping of the fronts of a descending water slug owing to non-equal capillary pressures on the imbibition and drainage fronts of the slug.

Текст научной работы на тему «Impact of the Dhofar cloud forest on subsurface hydrological fluxes: modelling on the scale of a near-root zone soil continuum»

Y^K 630*161(581.5):517.927.4

A. Kasimov, A. Hildebrandt, Yu. Obnosov

IMPACT OF THE DHOFAR CLOUD FOREST ON SUBSURFACE HYDROLOGICAL FLUXES: MODELLING ON THE SCALE OF A NEAR-ROOT ZONE SOIL CONTINUUM

A slope forest in a semi-arid region of Oman is studied. The peculiarity of this ecosystem is an interception of the liquid (fog and drizzle) by the tree canopy from a monsoon wind during a three-month monsoon period. The intercepted water drips on the soil from the foliage and stem over the area of the projection of the canopy. The moisture accumulates in the soil profile and is used by the tree long after the monsoon. Capillarity and gravity drives moisture to the water table. A modified Green-Ampt model of vertical infiltration is suggested, with a cyclos-tationary «self-drip-irrigation», vertically non-uniform water uptake by the tree roots and the possibility of stopping of the fronts of a descending water slug owing to non-equal capillary pressures on the imbibition and drainage fronts of the slug.

Key words: deciduous cloud forest, arid/semi-arid region, hillslope vegetation, orographic rain, root water uptake, infiltration, Darcy law, Green-Ampt model.

Introduction. In arid and semi-arid climates, which we, following [1], call water-limited environments (WLE), the soil moisture, salinity and temperature stresses determine the performance of plants. In these environments emerge interesting patterns of soil water, which are associated with the spatial distribution of vegetation [2]. In comparison with crops, pastures, grasslands or shrublands, forests are usually considered more water demanding [3]. In WLE higher water consumption of trees or shrubs is often explained by the ability of deep rooted vegetation (such as perennial trees or shrubs) to access water from a deep wet soil and even directly from the water table [4] (other plants with lower water use and with shallower roots can not do this). As a result of a combined water uptake by the roots from the soil, deep va-dose zone and aquifer, the experimental data [1] in Kalahari showed that on the scale of 0-60 m from the top soil to the water table the moisture content decreases in the first few meters, then remains almost constant (for approximately 5-25 m from the surface) and then increases again (25-60 m) close to the capillary fringe, with a full saturation on the phreatic surface. This «trough»-shaped moisture profile is attributed to the vertical water-hunting activity of deep tree roots.

a

b

Fig. 1. Satellite image of the study area (a) and vertical section of the mountain ridge, forest, monsoon and temperature profile (b)

© Kasimov A., Hildebrandt A., Obnosov Yu., 2010.

In this paper we consider infiltration (and water uptake) of a semi arid forests in SouthWest Oman, in the governorate of Dhofar (in Fig. 1, a a wite oval illustrates the study area on the the Arabian peninsula). The Dhofar cloud forests are extraordinary eco-hydrological systems, since their existence depends on their ability to self irrigate their root zone through interception of cloud water by the tree crowns [5, 6]. A vertical cross-section of the forest belt is shown in Fig.lb. The forests thrive in a seasonal climate with a 3 months humid (and foggy) season, which lasts from mid June to mid September, followed by a 9 months dry season. The characteristic monsoon fog is caused by orographic clouds, which are formed when moist air from the sea is pushed up against the coastal mountain range from South to North. The annual precipitation in this area is 250 mm, most of which falls during the monsoon in the form of light drizzle (on the average 2mm /day). The other major source of rainfall is cyclones, which - due to high intensity and prevailing runoff - seem, however, not to penetrate the deep soils as thoroughly as the precipitation during the monsoon [7] concluded that water, received below the tree crowns, consists of about one third of rainfall and two thirds of intercepted cloud water. This leads to a pattern where more soil water infiltrates below the tree crowns as compared to a bare soil.

A typical desiduous tree is shown in Fig. 2. Trees of this type may form a cluster with overlapping canopies.

Fig. 2. A typical tree of the cloud forest in the study area

It is noteworthy that the tree-induced drip self-irrigation in Dhofar forests guided the local engineers to install a net-mesh devises mimicking the canopy as fog-drizzle collectors with up to 58 liters/day/m2 of the mesh water condensate productivity (see, e.g., [8]).

Cloud interception is reported from many cloud forests (i.e. [9, 10]), but particularly in semiarid cloud forests, the additional water may substantially contribute to maintaining lush vegetation in a comparatively dry environment [11-14, 6]. To the best of our knowledge, in earlier works on fog forests (both humid and semi-arid) the dynamics of dripped/infiltrated water in the soil and vadose zone has not been studied in detail and only occasional measurements/modeling of the soil moisture (or capillary pressure), evaporation rate from the soil surface and soil temperature were reported. This is understandable because collecting field data from the whole soil- vadose zone - capillary fringe profile in the vicinity of a tree is instru-mentally cumbersome. Indeed, belowground hydrological sensing, e.g. with TDR, C- or Theta-probes would require a good vertical density of instruments in a borehole. Even piezome-

ters (deep in WLE) should be close to the tree trunk and well adjusted to intermittent recharge pulses from the vadose zone subjacent to the tree, with a good separation of the superjacent tree-generated recharge pulses from regional groundwater flow fluctuations or cyclicity. So, the objective of this paper is to develop a simple mathematical model, applicable to field conditions (i.e. having as few parameters as possible), that would predict the following fog-forest specific infiltration scenarios: depth-time- varying input from the soil surface, depth-time-varying water uptake by the tree roots, and ascending direction of the travel of the moisture content maximum as in [15]. The model should predict how the water content varies with time and depth, in particular, the «trough» -type moisture profile of [1].

Physical and Mathematical Model. We consider a tree of a fog-forest and neglect the interference of the drip zones of neighbouring trees (Fig. 3). Alternatively, Fig. 3, a can be viewed as a tree cluster with an «effective canopy» (and «effective» infiltration). The real dripping intensity, id, radially decreases from the axis of the infiltration zone (stem flow «funnel») to the periphery of the canopy projection on the ground surface. Outside the throughfall zone (of radius R and area A) the intensity of background precipitation, ib, does not vary with the radial coordinate r counted from the stem axis. Both id and ib depend on time. The dependence id (r) is seldom measured and we assume in the model that i=Id(t) is an average value for both the stemflow and throughfall integrated over r.

Projection of the canopy on the ground surface, area A, infiltration ld Off-can opy M soil surface, background infiltration ih

Water table

Stem axis

Soil surface

r, radial distance from stem axis

adius of the ïajiopy projection Draina§s ar4a

Vadose zone

fr

b)

Off-tree

ont D

water

V,(t)

V Stem axis 0_

slugs^

vd(t)

Tree water slug

Imbibition front

Watertable

Fig.3. 3-D schematization of tree and subjacent soil (a), axial cross-section of the subsurface with water fronts (b), cyclostationary component of water uptake function by roots (c)

We assume the following: For one monsoon period Id continues from t=0 till t=T=const. Outside this time interval Id=0. The intensity of infiltration during 0<t<T is high enough so the soil surface is completely saturated. At t=0 a saturation (imbibition) front yi(0) appears and propagates downward as yi(t) with a vertical axis of a Cartesian coordinate system oriented upward. A purely vertical infiltration takes place and we adapt the Green and Ampt (GA) model [16] both in the under-canopy zone and outside this zone. Without any loss of generality, we suppose that during the monsoon period, transpiration is relatively small and the water uptake from the saturated zone is negligible (corresponding to [7]). Then for assumed 1-D infiltration the position of the front over the interval 0<t<T in both zones is described by the well known analytical solution of Polubarinova-Kochina [17]. The focus of our modelling exercise is in what happens after t=T, i.e. after the end of the monsoon, when infiltration ceases and transpiration starts.

When t>T another upper desaturation (drainage) front, yd(t), yd(0)=0 is formed. Obviously, yd(t)<yi(t) and this inequality should hold for any time. The four fronts are shown in Fig. 3, b

(vertical axial cross-section of Fig. 3, a) and we neglect the smearing effect on the cylindrical surface separating the zones of contrasting infiltration rate (dashed lines in the Fig. 3, b). For the sake of brevity we will study two under-tree fronts only. We track yi until the instance Ta when this advancing front arrives at the water table (depth D). After that time a groundwater mound is formed [17] and the GA model does not work.

We assume that the porous medium is homogeneous, isotropic and incompressible, water is incompressible and flow obeys the Darcy law. We introduce the pressure head p=h-y where h(y,t) is the hydraulic head within the water slug located in between the two fronts, yi(t) <y< yd(t). The basic assumptions of the GA model are: a) in the slug water is under tension but completely saturates the pore space, b)the hydraulic conductivity in the slug is k=const that coincides with the saturated conductivity, c) water moves in the slug only, i.e., ahead of the advancing front and behind the retreating front k=0.

The GA model postulates that on the front p=-pc =const. This front pressure head is one of two physical parameters of the model (another is k). Both pc and k are tabulated for common porous materials (see, e.g. [17]). In our case of a two-front slug in Fig. 1 we assume

p=-pi at y=yi, p=-pd at y=yd, (1)

where pd and pi are two given positive constants. Commonly, pd > pi that on the meniscus level of an individual soil pore corresponds to a higher contact angle of the advancing front compared to the retreating one. In our model, however, the inequality can mathematically have an opposite sign.

From the Darcy law the specific discharge q(y, t) is:

q = ~k~. (2)

dy

We assume that water is lost from the moving phase of the slug to an immobile (with respect to the soil) storage, which is prevalently due to the uptake by the roots but can also involve sorption into a stagnant film around once wetted soil particles (or water interception in dead-end pores). The intensity of this «sink» tapping the Darcian water is E(y, t). From conservation of mass the total vertical discharge through any vertical cross-section of the cylinder in Fig. 1, a is Q=A q and we have

= E (t, y). (3)

dy

Thus E in (3) is the intensity of interception of vertical descending seepage into an immobile liquid (root) compartment. Within the roots the moisture keeps moving (upward, laterally or downward, see [1] but this in-root motion does not affect the GA motion in the soil. A similar model with a distributed sink term is used in the Dupuit-Forchheimer approximation of prevalently horizontal groundwater flow in an unconfined aquifer [17]. In our case E also should be small enough for the basic model assumption of vertical infiltration to hold. Integration of (3) yields:

J E(t, y )dy = Qs (t),

- D

where Qs is the total liquid flow rate, diverted from seepage to roots. We will assume that all Qs are converted into the sap flow, i.e. thus ignoring tree root internal storage and evaporation from the upper front. Moreover, this assumption ignores a potentially complex topology of water motion within the root system [1], which one can hardly measure.

Combining (2) and (3) we get the governing equation

d2h

77 = -e(t, y), (4)

dy

where e=E/(A k). This equation is similar to one for the water table dynamics in an uncon-fined aquifer with a spatially nonuniform and transient recharge [17].

We assume

e(t, y) = e0f(t )exp[oy], (5)

where a is a constant, characterizing the root system, and f(t) is a function, characterizing a season-dependent uptake of the whole root continuum (one example of this function is shown in Fig. 3, c). Instead of an exponential decrease of the uptake with depth, in (5) any other function can be taken (see, e.g. [18]). If f=const and a=0 then the water uptake is steady and uniform within the travelling slug. We recall that the GA limitations require that neither capillarity nor uptake create unsaturated conditions in the slug.

Integration of (4) gives:

h(t, y) = CjO) + c2(t)y - e0f (t)exp[ay]/a2, (6)

where c1 and c2 are two functions which are expressed from the boundary conditions (1) applied on both fronts:

h(t,-y (t)) = c (t) - C2 (t)yt (t) - ej(t) exp[ -ay, (t)] / a2 = -y (t) - p,

(7)

h(t,-yd(t)) = C1(t) - C2(t)yd(t) - e0f (t) exp[-ayd(t)] /a2 = -yd(t) - Pd.

Eqns (7) make a system of two linear equations with respect to c1 and c2 from which we immediately obtain by eliminating c2 the following expression:

c (t) = 1 + Pi - Pd - e0f (t)(exp[-ay, ] - exp[-ayd ]) (8)

2 yt - yd a2(yt - yd) '

The expression for c1 and elimination of this function follows from the first equation in (7).

On the fronts the linear average velocity and Darcian velocity are related through the so-called kinematic boundary conditions [17]:

m ^ = k dh«,-y, (>)) = k [c2 - ej-(,) exp[-ay, (t)] / a] dt dy

(9)

md^ = k dh(t,~ yd (t)) = k h - ej (t) exp[-ayd (t)] / a], dt dy

where md and mi are constant effective porosities on the corresponding fronts. Normally md < mi (if the slug propagates through a dry soil) but mathematically the sign of this inequality can be opposite. Eqns. (9) constitute a system of two ordinary nonlinear differential equations [19] with respect to yi(t) and yd(t), c1 and c2 taken from eqns (7)-(8). Without any loss of generality we rescale time as tr=t-T and drop «r» as a subscript. Then the initial conditions for eqn. (9) are:

yd (0) = 0, yt (0) = y0, (10)

where y0 is a given monsoon-wetting depth, which we retrieve from the monsoon stage of infiltration.

If both fronts asymptotically (effectively) stop before reaching the depth D then the slug «hangs» in the vadose zone, waiting for the next-season stronger slug to come, overtake the stopped one and resume the descending motion as a «merged» slug. If f(t) is high enough then the drainage front eventually catches the imbibition one. Mathematically, the Cauchy problem can be solved for any |yi(t)|<|yd(t)| but physically this inequality is unrealistic. If it appears from integration of (9)-(10) this means that the GA model becomes inapplicable, i.e. a tension-saturated flow can not be realized. Then an unsaturated flow model should be used as in e.g. [15].

Modeling Results. For a free standing tree 3 m of the canopy projection radius is a large crown. In Gogub site of Dhofar the average distance between trunks is 2.5m. The trees there are arranged in clusters, one cluster containing five and more stems in a very small area (1 m radius). The average distance between the clusters is much larger (about 5 m). In computations below we selected 1.5 m as the crown radius i.e. a tree density of 0.14.

Soil properties under the trees varied and systematic hydropedological studies were not undertaken in the study area. For the modeling below we assumed that the soil is a laveen loam (43% sand, 44% silt, 13% clay according to the USDA taxonomy [20]). From experiments in [20] k=0.216 m/day, pd=0.47 m, pi=0.22 m. Without any loss of generality we assumed D=20 m, md=0.2, mi=0.3.

The annual rainfall outside the canopy is about 100 mm. Annual precipitation over the canopy projection is 300 mm/year [7]. In computations based on [7] we assumed that 85% of precipitation in both zones (canopy-covered and bare) goes for transpiration. For the selected crown radius this gives 1750 liters/year/tree under the canopy. For comparison, the most common fruiting tree in the region, date palm, consumes 30-50 thousand liters of water through regular irrigation. Most of the water uptake in the Dhofar forest occurs after the end of the monsoon until a certain instance when water becomes limiting and transpiration decays [7]. It is noteworthy that irrigation of date palms and other crops requires thorough scheduling over season with the range of irrigation 50-150 l/day (winter-summer) that is time, labour and instrument consuming. The natural forest does self-irrigation spontaneously.

The following balance equation holds:

365

{ Q(t)dy = 1.75 m3

0

wherefrom we calculated e0 assuming different values of a. Experimental data for a are not available for our trees. Along with the exponentially decreasing uptake - depth function other «integrable» functions can be included into the model, for example, linear ones as in [21].

y,

14 1.2 1 0.8 0.6 0.4 0.2

Yd

2

4

50 100 150 200 250 300 350

Fig. 4. Dynamics offronts for two different a

In Fig.4 the fronts yi and yd are shown as functions of t in the interval 0<t<365 days for a cyclostationary sink with f(t)=sin2(rat), ra=2 tc/182.5 1/day. We solve eqns. (9)-(10) as a Cauchy problem using a standard Runge-Kutta NDSolve routine of computer algebra [23]. We fixed yo=0.6 m and selected a=0.1 (eo=0.00064) and a=2 1/m, (eo=0.013). Curves 1 and 2 correspond to the imbibition fronts and curves 3-4 indicate the drainage fronts for these two a values. As we can see from the graphs, the fronts decelerate with fluctuations (the slugs

1

t

shrink). Our computations also showed that with the increase of the initial slug thickness its migration rate increases. Hydrological implications of this dynamics are clear: in Darcian flows inertia is neglected [17] and therefore the motion of two consequitive slugs can be studied by specifying their initial positions only. If a larger «new» slug moves faster it may catch up a smaller one. Tehn two slugs merge and move even faster to the water table.

Discussion and Conclusion

Soil moisture distribution and its dynamics in the soil profile are vitally important for the flora in WLE. Cloud forests in arid/semi-arid climates of the Arabian peninsula in the Dhofar region of Oman live in unique hydrological niches. The forest in this region is highly susceptible to variations of the precipitation fluxes, which are characterized by a 3-month long monsoon duration and annual cyclicity. Precipitation consists of two components: (1) Orographic low-intensity rain (drizzle) and (2) fog collection (turbulent deposition of cloud droplets). Fog collection leads to heterogeneous infiltration: Under tree canopies infiltration is estimated to be twice as high, compared to the off-canopy surroundings. The latter leads to temporary and radially varying compartmentalization of water in the root zone.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

The forest-enhanced precipitation stimulates the biota and possibly contributes to the recharge of the water table. Net precipitation under the tree canopies has been estimated to be around 300 mm over the entire 2004 monsoon season. The transient subsurface storage accumulated during the monsoon is consumed later by evapotranspiration (net ascending flux) and is composed of the gravitational, capillary and viscous (Darcian) forces (net descending flux). In this paper we predict how much out of the 300 mm net precipitation might reach the aquifer, provided the water uptake by the tree roots varies with depth and time that results in a «trough-shaped» moisture profiles with depth.

We considered the monsoon-generated soil moisture volume as a shrinking parcel, which travels through the vadose zone, and implemented the Green-Ampt type model on the scale of a near-trunk cylindrical porous column. Water motion is characterized by two fronts (drainage and imbibition) and soil water removal by the roots. The advancing front flux is greater than the retreating one. The ordinary differential equations for the front loci are solved numerically, and the possibility of water table recharge is estimated.

If the slugs reach the water table, then groundwater mounds (axisymmetric with respect to the tree roots) emerge and decay cyclostationary. The possibility for the deep tree roots to be immersed into the capillary fringe should be assessed in the future. If the forest is approximated as a regular (e.g., chess-board) pattern of trees without interference of adjacent trunk-centered moisture parcels, then the catchment-scale recharge of the aquifer can be evaluated, in particular, the effect of a partial uphill deforestation (for example due to cattle browsing or infections) on the drop of the groundwater level and hence capillary fringe. The latter may induce further stress on phreatophytic species in the downhill zone.

References

1. Lubczynski, M. W. The hydrogeological role of trees in water-limited environments / M. W. Lubczynski // Hydrogeology J. - 2008. - V. 17. - P. 247-259.

2. Brown, G. Community composition and population dynamics in response to artificial rainfall in an undisturbed desert annual community in Kuwait / G. Brown // Basic Appl. Ecol. - 2002. - V. 3. - P. 145-156.

3. Van Dijk, A. I. J. M. Planted forests and water in perspective / A. I. J. M. Van Dijk // Forest Ecology and Management. - 2007. - V.251 - P. 1-9.

4. Steward, D. R. An analytic solution for groundwater uptake by phreatophytes spanning spatial scales from plant to field to regional / D. R. Steward, T. S. Ahring // Journal of Engineering Mathematics. - 2009. -V. 64. - P.85-103.

5. Hildebrandt, A. Forest on the edge: Seasonal cloud forest in Oman creates its own ecological niche / A. Hildebrandt, E. A. B. Eltahir// Geophys. Res. Lett. - 2006. - V. 33. - L11401.

6. Hildebrandt, A. Using a horizontal precipitation model to investigate the role of turbulent cloud deposition in survival of a seasonal cloud forest in Dhofar / A. Hildebrandt, E. A. B. Eltahir // J. Geophys. Res. - 2008. - V. 113. - G04028.

7. Hildebrandt, A. Ecohydrology of a seasonal cloud forest in Dhofar: 1. Field experiment / A. Hildebrandt, M. A. Al-Aufi, M. Amerjeed and othes //Water Resources Research . - 2007. - V. 43. -W10411.

8. Abdul-Wahab, S. A. Total fog and rainwater collection in the Dhofar region of the Sultanate of Oman during the monsoon season / S. A. Abdul-Wahab, A. M. Al-Damkhi, H. Al-Hinai and othes // Water International. - 2010. - V. 35. - P.100-109.

9. Bruijnzeel, L. A. Hydrology of tropical montane cloud forests: A reassessment / L. A. Bruijnzeel // Land Use and Water Resources Research. - 2001. - V. 1. - P. 1.8-1.18.

10. Foster, P. The potential negative impacts of global climate change on tropical montane cloud forests / P. Foster // Earth-Science Reviews. - 2001. - V. 55. - P. 73-106.

11. Hursh, C. R. Field moisture balance in the Shimba Hills, Kenya / C. R. Hursh, H. C. Pereira // The East African Agricultural Journal. - 1953. - V. 18. - P. 139-148.

12. Juvik, J. O. Relationship between rainfall, cloud-interception, and canopy throughfall in a Hawaiian montane forest / J. O. Juvik, D. Nullet // In Tropical Montane Cloud Forests. - 1995, edited by L. S. Hamilton, et al., Springer Verlag. - New York, 1995.

13. Hutley, L. B.Water balance of an Australian subtropical rainforest at altitude: The ecological and physiological significance of intercepted cloud and fog / L. B. Hutley, D. Doley, D. J. Yates, A. Boonsaner // Australian Journal of Botany. - 1997. - V. 45. - P. 311.

14. del-Val, E. Rain forest islands in the chilean semiarid region: Fog-dependency, ecosystem persistence and tree regeneration / E. del-Val, J. J. Armesto, O. Barbosa and othes // Ecosystems. - 2006. - V. 9. - P. 598608.

15. Philip, J. R. Redistribution of soil water from plane, line, and point sources / J. R. Philip, J. H. Knight // Irrigation Science. - 1989. - V. 12. - P. 169-180.

16. Kacimov, A. R. The Green-Ampt 1-D infiltration from a ponded surface into a heterogeneous soil / A. R. Kacimov, S. Al-Ismaily, A. Al-Maktoumi // J. Irrigation and Drainage ASCE. - 2010. -V. 136(1). - P. 6872.

17. Polubarinova-Kochina, P. Ya. Theory of Ground-water Movement / P. Ya. Polubarinova-Kochina. -Princeton: Princepton Univ. Press, 1962.

18. Jarvis, N. J. A simple empirical model of root water uptake / N. J. Jarvis // J. Hydrology. - 1989. -V. 107. - P. 57-72.

19. Ince, E. L. Ordinary Differential Equations / E. L. Ince // Longmans, Green and Co. - London, 1926.

20. Bouwer, H. Rapid field measurements of air entry value and hydraulic conductivity of soil as significant parameters in flow system analysis / H. Bouwer // Water Resources Research. - 1966. - V.2. - P. 729-738.

21. Molz, F. J. Extraction models of soil moisture use by transpiring plants / F. J. Molz, I. Remson //Water Resources Research. - 1970. - V.6. - P. 1346-1356.

22. Prada, S. Fog precipitation and rainfall interception in the natural forests of Madeira Island (Portugal) / S. Prada, M. M. de Sequeira, C. Figueira, M. A. da Silva // Agricultural and Forest Meteorology. - 2009. -V. 149. - P. 1179-1187.

23. Wolfram, S. Mathematica. A System for Doing Mathematics by Computer / S.Wolfram. - Addison-Wesley, Redwood City, 1991.

Статья поступила в редакцию 03.08.10.

А. Р. Касимов, А. Хильдебрандт, Ю. В. Обносов

ВЛИЯНИЕ СКЛОНОВОГО ЛЕСА ОМАНА (РЕГИОН ДОФАР), ПИТАЕМОГО ОСАДКАМИ ИЗ ТУМАНА, НА ПОДЗЕМНЫЙ ГИДРОЛОГИЧЕСКИЙ РЕЖИМ: МОДЕЛИРОВАНИЕ В МАСШТАБЕ КОРНЕВОЙ ЗОНЫ

Рассматривается склоновый лес полуаридной зоны в Омане, особенностью которого является перехват кроной деревьев жидкости (тумана и осадков слабой интенсивности) из муссонного ветрa в течение трех месяцев в году. Эта вода капает с листвы на землю по площади проекции кроны и стекает по стволу, в результате чего влага аккумулируется в прикорневой зоне почвы и используется деревом в течение длительного постмуссонного периода. За счет капиллярности и гравитации почвенная вода дрейфует к зеркалу грунто-

вых вод. Предложена модифицированная математическая модель Грина-Ампта для вертикального влагопереноса в зоне питания дерева, которая учитывает циклостационар-ность процесса «капельного самоорошения», неравномерный по глубине перехват влаги распределенными корнями дерева и возможность остановки движения воды за счет неодинаковых капиллярных давлений в зоне промачивающего и дренирующего фронтов.

Ключевые слова: листопадный туманный лес, засушливый/полузасушливый регион, склоновая растительность, орографический дождь, забор влаги корнями растений, закон Дарси, модель Грина-Ампта.

КАСИМОВ Анвар Рашидович - кандидат физико-математических наук, профессор, декан факультета гидрологии, почвоведения и сельхозинжиниринга Университета Султана Кабуса (Оман). Область научных интересов - течения в пористых средах. Автор более 130 публикаций.

E-mail: anvar@squ.edu.om

ХИЛЬДЕБРАНДТ Анке - младший профессор по экологическому моделированию Института наук о Земле, Университет Фридриха Шиллера, Йена, Германия. Область научных интересов - измерения и моделирование забора влаги корнями растений, интерпретация полей почвенной влаги, обратные методы моделирования движения влаги в почве, взаимоотношения между структурой экосистем и вертикальным потоком воды из атмосферы в почву и обратно, динамика растительности в полуаридном климате. Автор семи публикаций.

E-mail: anke.hildebrandt@ufz.de

ОБНОСОВ Юрий Викторович - доктор физико-математических наук, профессор, заведующий кафедрой дифференциальных уравнений Казанского Федерального университета. Область научных интересов - краевые задачи теории голоморфных функций и их приложения в задачах механики сплошных сред. Автор более 130 публикаций.

E-mail: yobnosov@ksu.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.