Effect of Mesoscale Land Use Change on Characteristics of Convective Boundary Layer:Semi-Idealized Large Eddy Simulations over Northwest China

Bangjun CAO, Shuwen ZHANG, Deqin LI, Yanlin LI, Linfan ZHOU, and Jiemin WANG

1 School of Atmospheric Sciences, Chengdu University of Information Technology, Chengdu 610225

2 Key Laboratory for Semi-Arid Climate Change of the Ministry of Education, College of Atmospheric Sciences, Lanzhou University, Lanzhou 730000

3 Shenyang Central Meteorological Observatory, Shenyang 110016

4 Cold and Arid Regions Environmental and Engineering Research Institute,Chinese Academy of Sciences, Lanzhou 730000

ABSTRACT Although large-scale topography and land use have been properly considered in weather and climate models, the effect of mesoscale and microscale heterogeneous land use on convective boundary layer (CBL) has not been fully understood yet. In this study, the influence of semi-idealized strip-like patches of oases and deserts, which resemble irrigated land use in Northwest China, on the CBL characteristics, is investigated based on the Weather Research and Forecasting (WRF)-large eddy simulation (LES) driven by observed land surface data. The influences of soil water content in oases on aloft CBL flow structure, stability, turbulent kinetic energy (TKE), and vertical fluxes are carefully examined through a group of sensitivity experiments. The results show that secondary circulation (SC)/turbulent organized structures (TOS) is the strongest/weakest when soil water content in oases is close to saturation (e.g.,when the oases are irrigated). With the decrease of soil water content in oases (i.e., after irrigation), SC (TOS) becomes weak (strong) in the lower and middle CBL, the flux induced by SC and TOS becomes small (large), which has a dramatic impact on point measurement of eddy covariance (EC) fluxes. The flux induced by SC and TOS has little influence on EC sensible heat flux, but great influence on EC latent heat flux. Under this circumstance, the area averaged heat flux cannot be represented by point measurement of flux by the EC method, especially just after irrigation in oases. Comparison of imbalance ratio (i.e., contribution of SC and TOS to the total flux) reveals that increased soil moisture in oases leads to a larger imbalance ratio as well as enhanced surface heterogeneity. Moreover,we found that the soil layer configuration at different depths has a negligible impact on the CBL flux properties.

Key words: oasis and desert, land surface heterogeneity, large eddy simulation, soil water content, secondary circulation, turbulent organized structure

1. Introduction

Land surface heterogeneities may result from topography and land use, ranging from continental scale down to mesoscale and microscale. Although large-scale topography and land use have been properly considered in weather and climate models, the effect of mesoscale and microscale heterogeneous land use on convective boundary layer (CBL) has not been fully understood yet,especially for microscale heterogeneous land use (Heinemann and Kerschgens, 2005).

Because an atmospheric general circulation model(GCM) or a mesoscale weather model is not able to correctly resolve microscale physical processes resulted from microscale topographical and land use heterogeneities, turbulence-resolving large eddy simulation (LES) has been used to investigate the interaction between heterogeneous land surface and the atmosphere aloft (Shen and Leclerc, 1995; Raasch and Harbusch, 2001; Courault et al., 2007; Huang and Margulis, 2010; Shen et al., 2016).For example, Avissar and Schmidt (1998) showed that secondary circulation (SC) was strong when the amplitude and wavelength of an ideal surface heat wave were large; Maronga and Raasch (2013) investigated the impact of real surface heterogeneity on CBL and found the developments of SC with local or roll-like patterns.

To better investigate the interaction between land surface heterogeneity and the atmosphere aloft, LES has been coupled to land surface model (LSM) (Courault et al., 2007; Huang and Margulis, 2010; Shao et al., 2013).For example, Patton et al. (2005; hereafter PSM05) reported a large contribution of SC to vertical fluxes by using an idealized one-dimensional soil moisture heterogeneity, and Shao et al. (2013) investigated the impact of real surface heterogeneities on SC and found that the heterogeneity induced fluxes were closely related to land use patterns.

In the arid region of northwestern China, a specific type of land surface exists with a mixture of Gobi, sand desert, and oasis (Liu et al., 2011; Wang et al., 2012),where irrigation has been widely adopted (Li et al.,2013). Adding soil water into oases further enlarges the difference of land surface physical characteristics between oases and surrounding deserts. Up to now, few studies have been carried out to investigate the influence of this kind of surface heterogeneity on flow structures and vertical fluxes in the CBL with the coupled LES.

In addition, LES has been adopted to investigate point flux underestimation by the eddy covariance (EC) method(Steinfeld et al., 2007; Foken, 2008; Huang et al., 2008;Wang et al., 2009; Foken et al., 2010). Kanda et al.(2004) and Inagaki et al. (2006) demonstrated that SC with the scale between 4 and 16 km and local advection by turbulent organized structures (TOS) with the scale of 1–2 km or less were both causes for the flux underestimation by EC.

In this study, the CBL structure over semi-idealized strip-like patches of deserts and oases is first investigated; the effect of soil water content in oases on CBL flow structure is analyzed; the contribution of SC as well as TOS to energy balance is quantified; and the effect of soil layer configuration on CBL characteristics is finally examined. This paper is organized as follows. Section 2 describes the LES model and associated experimental design. Section 3 describes the analysis method. The simulation results and discussion are presented in Section 4. A summary is provided in Section 5.

2. Introduction of models and experimental design

2.1 Models

The Weather Research and Forecasting (WRF) model in its large eddy simulation (WRF-LES) coupled to the Community Noah land surface model with multi parameterization options (Noah-MP) was used (Moeng et al.,2007; Niu et al., 2011). The WRF-LES has state-of-theart numerical schemes (e.g., higher-order time and advection schemes), and has been evaluated and validated over different heterogeneous lands (Solomon et al., 2011;Talbot et al., 2012; Yamaguchi and Feingold, 2012). Talbot et al. (2012) showed that simulation with the fine resolution WRF-LES produced better surface fluxes and near-surface meteorological parameters.

In our simulations, subgrid scale (SGS) turbulence parameterization is 1.5-order TKE closure scheme(Klemp and Wilhelmson, 1978; Deardorff, 1980). Lateral boundaries are periodic condition. Surface layer scheme is based on Monin–Obukhov (MO) similarity theory. A 2-m depth of soil is divided into 4 layers with thicknesses of 0.1, 0.3, 0.6, and 1.0 m, respectively. The chosen parameterization schemes are ball-berry canopy stomatal resistance scheme, MO surface layer drag coefficient scheme, modified two-stream radiation transfer scheme, and CLASS (Canadian Land Surface Scheme)ground snow surface albedo scheme.

2.2 Experimental design

The horizontal domain in this study covers an area of 20 km × 5 km with Δx= Δy= 50 m, and vertical height is 2.6 km. The variable Δzchanges from 10 m near the surface to 40 m forz≥ 100 m. Time step is 0.4 s. A Rayleigh damping layer is set to damp the reflection of gravity waves (Klemp et al., 2008), and the thickness of the damping layer is 500 m. Initial potential temperature and moisture profiles for LES are taken from radiosonde measurements at 1200 local time (LT) 1 August 2012 at Zhangye (39°5′15″N, 100°16′39″E, 1556.06 m a.s.l.) (Li et al., 2013). Initial potential temperature is 310 K below 925 m, increases by 8 K over a strong inversion depth of 150 m; initial CBL heightziis 1000 m, and above 1000 m, stratification is 3 K km–1; initial water vapor mixing ratio is constant, with a value of 7 g kg–1below 900 m,dropping sharply to 2 g kg–1across initial inversion layer between 925 and 1075 m; initial horizontal mean wind is zero for all cases.

Incoming solar radiation is 800 W m–2and longwave solar radiation is 350 W m–2by averaging measurements at Daman (38°51′19″N, 100°22′20″E, 1556.06 m a.s.l.)in the oasis patch between 1000 and 1600 LT during the study period, respectively. Atmospheric radiation scheme is not included so any longwave radiation leaving the surface and scattered shortwave radiation are assumed to radiate into space. Albedo is diagnosed from the surface solar radiation budget.

The LES is run with land use conditions supplied from the LSM outputs. In the LSM simulations, the strip-like land types of oasis and desert are specified as cropland/woodland mosaic and barren as in the United States Geological Survey, respectively, with no variation in they-direction (Fig.1). Roughness length is obtained from a lookup table; the minimum and maximum background roughness lengths are 5 and 14 cm for oasis, respectively, and 1 and 1 cm for desert. The coupled LES–LSM simulations are carried out on 4 representative days of 1, 8, 15, and 21 August, representing the 3rd,10th, 17th, and 23rd days (hereafter referred as D03,D10, D17, and D23) after irrigation, and each simulation was integrated for 5 h. Initial soil water contents in oasis and desert on the above 4 days are taken from measurements at 1200 LT at Daman station and Shenshawo desert site (38°47′21″N, 100°29′35″E, 1594.00 m a.s.l.),respectively (Table 1 and Fig.2). Initial soil temperatures are based on ground-based measurements at Daman and Shenshawo at 1200 LT 1 August in all cases(Table 2).

Fig.1. Schematic diagram showing the coupling of the LES model to the LSM.

Eleven coupled simulation experiments are conducted(see Table 3). They are named as follows. D03-10S, for example, refers to a simulation on the 3rd day after irrigation with the width (i.e., λ) of one patch of oasis and desert equal to 10 km. In the first two cases (D03-OS and D03-DS, Table 3), only one type of land surface, i.e.,either oasis (OS) or desert (DS), is used with no land heterogeneity. D03-10SA is the same as D03-10S except for using a different soil layer configuration.

3. Description of methods

Under a few assumptions such as no synoptic ascending/descending motion, flat ground surface, no background horizontal wind, and periodic boundary condition, (total) vertical flux can be decomposed into the EC flux, SC flux, and TOS flux (Inagaki et al., 2006). Fur-thermore, to quantify the contribution of SC and TOS to the total flux, horizontal domain average of the imbalance ratio [I] (which can be decomposed into [Ih] for sensible heat flux and [Iq] for latent heat flux) and phase average of the imbalance ratio [I]pare defined as

Table 1. Initial soil water content (cm3 cm–3) in the four soil layers of the LSM, prescribed as the observed values at Daman (representing oasis)and Shenshawo (representing desert)

Fig.2. Temporal evolutions of soil water content at four different depths (0.04, 0.10, 0.40, and 0.80 m) from 1 to 25 August (irrigation was carried out on 29 July) at observation sites of (a) Daman and (b) Shenshawo.

Table 2. Initial soil temperature (K) in the four layers of the LSM, prescribed as the observed values at Daman (representing oasis) and Shenshawo (representing desert)

Table 3. Basic properties for various cases

wherewis vertical velocity (positive upward) andφis a scalar variable like potential temperatureθand specific humidityq; prime with subscript t or p means temporal fluctuation or temporal-plus-phase fluctuation; overbar indicates temporal averaging; [ ] and [ ]pare horizontal domain average and phase average, respectively. The variableis total flux,is EC flux,is resolved SC flux (hereafter referred to as SC), andis the resolved TOS flux (hereafter referred to as TOS). In the analysis, the EC is evaluated as the sum of computed grid-scale (GS) and subgrid-scale (SGS) components.See Appendix for details.

To acquire the quasi-steady state of the developing CBL, data analysis is carried out after 3.5-h integration.For all statistics, the averaging time is 1 h, which is a typical value used in field observations.

For normalization and comparison purposes, Table 3 lists basic properties of different cases, such as the width of one patch of oasis and desert (λ), CBL height (zi), total horizontal averaged surface sensible fluxand latent heat fluxand convective velocity scalewheregis the gravitational acceleration,θ0is reference temperature,is total surface latent heat flux,θv*is virtual potential temperature scale, andq*is humidity scale.

4. Results

4.1 CBL flow structure and atmosphere stability

To illustrate the influence of strip-like oasis and desert patches on CBL flow structure, Fig.3 shows the distributions of vertical velocity and water vapor mixing ratio averaged over 3.5–4.5-h integration for D03-10S. The ascending motion centers are near the middle of deserts atx= 7.5 and 17.5 km, confined to 30% of the wavelength,with a weak descending motion elsewhere (Fig.3a).Over oases, low water vapor zone in the upper CBL extends to the middle CBL due to the entrained dry air. In the lower CBL, water vapor mixing ratio over the oasis side is much larger than that over the desert side(Fig.3b).

To have a clear picture of vertical flow at different scales, temporal mean vertical velocityat the height of 100 m (Figs. 4a, b) is decomposed into SC part(Figs. 4c, d) and TOS part(Fig.4e, f). For a better resolution, only one oasis–desert heterogeneity cycle is depicted. Figures 4a and 4b show that ascending motion is centered over the middle of desert for D03- 10S, but is not clear for D23-10S; on the contrast, the“rib-like” pattern offor D23-10S is more obvious than that for D03-10S. Spatial structure of(Figs. 4c, d)has the same one-dimensional variation as surface heterogeneity and(Figs. 4a, b). The maximum ofoccurs at 7.5 km, the minimum is at 2.5 km, and their difference is more than 0.6 m s–1for D03-10S with an obvious SC (Figs. 4c), but less than 0.3 m s–1for D23-10S with no clear SC (Figs. 4d). As for TOSa very complicated “rib-like” pattern appears for both D03-10S and D23-10S (Figs. 4e, f), especially for D23-10S(Figs. 4f).

Fig.3. Snap shots of (a) vertical velocity and (b) water vapor mixing ratio, averaged along the y-axis of the study domain, after 3.5-h integration for D03-10S.

Fig.4. Temporal averages of vertical velocity () for (a) D03-10S and (b) D23-10S,and TOS for (e) D03-10S and (f) D23-10S.The SCaveraged along the y-axis of the domain for (c) D03-10S and (d) D23-10S. Values in all panels are at the height of 100 m, and the contours range from 0.1 to 0.6 m s–1 at an interval of 0.1 m s–1 in (a, b, e, f).

To quantitatively compare the variation of SC and TOS with the decrease of soil water content over oases from the 3rd to 23rd days after irrigation, the phase averaging standard deviation ofis calculated as

To illustrate the impact of soil water content in oases on the aloft atmosphere stability, thex-directional distribution of gradient Richardson number (Ri) at 100-m height, defined by

is shown in Fig.5c. In all cases, background Ri is equal to 0. But flow becomes dynamically stable with Ri >0.25 over oases and very unstable with Ri < 0 over deserts as in the case of D03-10S. When oases continue to become dry as in the case of D23-10S, flow becomes very dynamically unstable with Ri < 0 over both oases and deserts. In a word, the atmosphere stability over oases and deserts is different just after irrigation, but with the decrease of soil water content in oases, the stability over oases and deserts becomes similar.

Total turbulent kinetic energy (TKE) (both domain and time averaged over the CBL normalized byw*2) responds sensitively to land surface heterogeneity; for example, TKE is 0.51 just after irrigation in oases, but it rises to 0.68 on the 23rd day, due to the increased sensible heat flux.

To investigate dependence of the intensity of SC on soil water content in oases, correlation analysis between surface heat flux and vertical velocity is applied(Maronga and Raasch, 2013). A high correlation betweenandexists, especially at low levels, and becomes smaller with the decrease of soil water content in oases (Table 4). Besides, correlation coefficient decreases with the increase of height, indicating that the influence of oasis–desert patterns decreases with height.

4.2 Vertical fluxes in the CBL

4.2.1Vertical fluxes

Figure 6a shows the normalized profiles of total sensible heat fluxEC sensible heat fluxSC sensible heat fluxand TOS sensible heat fluxfor D03-10S, and the corresponding latent heat fluxes are shown in Fig.6b. Both total sensible heat fluxand total latent heat fluxdisplay classical profiles as shown by Deardorff (1972); total sensible heat flux decreases linearly with height and even becomes negative in entrainment zone (Fig.6a), and total latent heat flux increases with height up to the entrainment level and then quickly reduces to zero (Fig.6b).

In the CBL, EC sensible heat fluxdominates total sensible heat flux but decreases with height(Fig.6a). SC sensible heat fluxfirst increases with height, reaches the maximum of 0.22 atz/zi= 0.25,and then decreases. TOS sensible heat fluxhas a very slight variation with height. Relatively speaking, SC sensible heat flux is larger than TOS flux at the lower CBL.

Fig.5. (a) Standard deviation of temporal average of vertical velocity, (b) temporal average of velocity in the x-direction, and (c) Ri at the height of 100 m (dot-dashed line is background atmosphere Ri).

Normalized EC latent heat fluxreaches the maximum of 1.72 atz/zi= 0.85 (Fig.6b). SC latent heat fluxis negative with the minimum of –0.19 atz/zi= 0.30, and becomes positive in the upper CBL abovez/zi> 0.60. The contribution of TOS latent heat fluxto the total flux is positive, with the maximum of 0.34 atz/zi= 0.30. Because of the strong descending motion of SC and relative larger water vapor mixing ratio over oases than over deserts, SC latent heat flux is downward belowz/zi= 0.60 and thus water vapor at low levels is prevented from transporting into the atmosphere aloft over the oasis (Fig.3b) (Chu et al., 2005; Meng et al.,2012, 2015). Because of different experimental settings from our study, previous studies have not reported negative SC latent heat flux (e.g., Inagaki et al., 2006).

Table 4. Correlation coefficient between and at different heights

Table 4. Correlation coefficient between and at different heights

4.2.2Effect of soil water content on vertical fluxes

To show the effect of different soil water contents on vertical fluxes, the normalized profiles of the total vertical advection flux (i.e., the sum of SC and TOS), and the SC and TOS for D03-10S and D23-10S are plotted in Fig.7 for sensible heat and in Fig.8 for latent heat (water vapor). With the decrease of soil water content in oases,total vertical advection sensible heat flux decreases from the maximum of 0.26 in D03-10S to 0.20 in D23-10S atz/zi= 0.25 (Fig.7a), SC sensible heat flux atz/zi= 0.25 decreases significantly from the maximum of 0.22 in D03-10S to 0.06 in D23-10S (Fig.7b), and TOS sensible heat flux increases from the maximum of only 0.04 atz/zi= 0.20 in D03-10S to 0.14 in D23-10S (Fig.7c).In a word, SC becomes the strongest and TOS becomes the weakest when soil water content in oases is near saturation. Under this circumstance, SC is assumed to penetrate the entire CBL and exert strong influences on the CBL, and the flux induced by SC must not be ignored.

Different from sensible heat flux, total vertical advection flux for water vapor increases from 0.28 atz/zi=0.85 in D03-10S to the maximum of 0.54 in D23-10S(Fig.8a). SC latent heat flux in all cases is downward

Fig.6. Comparison of normalized total flux, EC, SC, and TOS for (a) sensible heat and (b) latent heat fluxes for case D03-10S.

Fig.7. Comparison of vertical profiles of (a) normalized vertical advection sensible heat flux (b) SC sensible heat flux and(c) TOS sensible heat flux for cases D03-10S and D23-10S.

Fig.8. As in Fig.7, but for latent heat flux.

and the absolute value of SC latent heat flux atz/zi= 0.35 decreases from 0.19 in D03-10S to 0.12 in D23-10S(Fig.8b). TOS latent heat flux in all cases is upward and increases slightly from the maximum of 0.34 in D03-10S to 0.48 in D23-10S atz/zi= 0.43 (Fig.8c). In a word, the decrease of soil water content in oases has a profound influence on SC and TOS with the magnitude of reduced SC flux relatively larger than that of increased TOS flux.

Flux patterns with λ = 20 km are similar to those with λ = 10 km, although the former flux has a slightly less magnitude (figure omitted), showing that the influence of different wave lengths on these fluxes is not significant;the reason may be related to the same atmosphere forcing used under different heterogeneity scales.

4.3 Imbalance ratio

First, we investigate the vertical variation of horizontal domain average of the imbalance ratios. The imbalance ratios for sensible heat flux [Ih] (Fig.9a) and latent heat flux [Iq] (Fig.10a) are about 6.1% and 1.2% at 10-m height, but rise to 24.1% and 3.2% at 100-m height, in D03-10S. Furthermore, the imbalance ratios for SC and TOS also increase with height, with the increase of magnitude for SC larger than that for TOS. For example, SC and TOS fractions of [Ih] are 3.4% and 2.7% at 10-m height, and increase to 18.3% (Fig.9b) and 5.2%(Fig.9c) at 100-m height, respectively.

As for the impact of soil water content in oases on horizontal domain average of the imbalance ratios, the impact at low levels is very small for sensible and latent heat fluxes (Figs. 9a, 10a), SC flux (Figs. 9b, 10b), and TOS flux (Figs. 9c, 10c). At high levels (60–100 m), horizontal domain average of the imbalance ratios for sensible heat flux [Ih] and the magnitude of SC sensible heat and latent heat flux all decrease with the decrease of soil water content in oases. On the contrary, the imbalance ratios for latent heat flux [Iq], TOS sensible heat flux, and TOS latent heat flux increase. In a word, the decrease of soil water content has a significant influence on the imbalance ratio at high levels, but the influence becomes weak with the decreased height.

Fig.9. Vertical variations of horizontal domain average of the imbalance ratios for (a) sensible heat flux, (b) SC sensible heat flux, and (c) TOS sensible heat flux.

Fig.10. As in Fig.9, but for latent heat flux.

Second, to investigate the variation of phase average of the imbalance ratios ([Ih]pand [Iq]p), Table 5 lists their values at threex/λlocations and four levels under two soil water conditions (i.e., wet soil on the 3rd day and dry soil on the 23rd day after irrigation), respectively. The magnitude of phase average of the imbalance ratios increases with height and becomes small with oasis drying.The phase average of the imbalance ratios has a very different performance at threex/λlocations: over the center of oasis (x/λ= 0.25), the ratios always have relatively large absolute values with [Ih]pnegative and [Iq]ppositive, indicating that EC flux near the center of oasis has a large deviation from horizontal averaged total flux; over the joint region of oasis and desert (x/λ= 0.5), the ratio is always positive but its magnitude is relatively small; over the center of desert, [Ih]pis always positive no matter whether the oasis is wet or dry, and [Iq]pis negative but becomes relatively small with oasis drying.

Based on the above results, the mechanism for the increase of soil water content in oases on the CBL characteristics is summarized as follows. With increase of soil water content in oases, both updraft over deserts and downdraft over oases become stronger, shear (buoyant)forcing becomes weaker (stronger), and TOS (SC) becomes weaker (stronger). Meanwhile, the imbalance ratio for sensible (latent) heat flux becomes larger(smaller), and vice versa. When soil water content in oases decreases, surface heterogeneity becomes weaker and thus the difference of turbulent structure parameters between oasis and desert becomes small.

4.4 Sensitivity experiment

In the above analysis, the 2-m depth soil with fixed 4-layer thicknesses of 0.1, 0.3, 0.6, and 1.0 m was used.However, previous studies have shown that soil layer configuration had an influence on the simulation of surface heat flux and soil moisture (Zhang et al., 2009,2010; Liu and Shao, 2013). Therefore, further experiment by using different soil layer configurations was conducted (hereafter D03-10SA), in which the default 4-layer thicknesses were respectively reduced to 0.01, 0.03,0.06, and 0.10 m (Liu and Shao, 2013), with the rest configuration the same as D03-10S.

Compared to D03-10S, the near-surfaceis 35 W m–2(approximately 20%) smaller andis larger by the same magnitude over oases in D03-10SA, while bothandover deserts are almost similar (Table 3).Compared with measurements at Daman and Shenshawo observation sites, the fluxes from D03-10S are closer to the measurements than those from D03-10SA. In addition, with the decreased thickness of soil layers, soil temperature at 0.1-m depth in D03-10S is larger than that in D03-10SA (Fig.11a), and both of them are larger than the measurements; soil water content at 0.1-m depth in D03-10S is larger than that in D03-10SA (Fig.11b), and both of them are smaller than the measurements. Furthermore, total soil water in soil column in D03-10SA is smaller than that in D03-10S. Therefore, the increased sensible heat and decreased latent heat for D03-10SA are due to the following two factors: first, increase in soil temperature for D03-10SA causes sensible heat increase and latent heat decrease, and second, the decreased soil water in D03-10SA further makes latent heat decrease.Moreover, the soil layer configuration has a negligible impact on total normalized sensible heat flux profile (Fig.12a) and a small impact on total normalized latent heat flux profiles with the maximum of 1.92 in D03-10S reduced to 1.71 in D03-10SA atz/zi= 0.85(Fig.12b). In addition, the impact of soil layer configuration on EC, SC, and TOS sensible heat flux and SC latent heat flux is almost negligible and that on EC and TOS latent heat flux is not large (figure omitted).

5. Summary and discussion

To study the influence of strip-like patches of desertsand oases on CBL flow structures such as SC and TOS,four experiments with different soil water contents in oases using observed data were carried out, representing real-time monthly timescale variation of soil water content in oases due to irrigation in a semi-idealized mode.Phase averaging and temporal averaging were used to decompose the total vertical flux into EC turbulent flux, SC flux, and TOS flux.

Table 5. Phase average of the imbalance ratios

Fig.11. Comparison of the simulated (a) soil temperature and (b) soil water content over oases for D03-10S (dashed line) and D03-10SA (solid line) with the observations (dots) at depth of 0.1 m.

Fig.12. Normalized vertical profiles of (a) total heat flux and (b) total latent heat flux from D03-10S and D03-10SA.

SC flux is larger than TOS flux in the middle and lower CBL for both sensible and latent heat fluxes; SC latent heat flux moves downward while the other fluxes are upward; SC is the strongest and TOS is the weakest when soil water content in oases is close to saturation.

With the decrease of soil water content in oases, SC becomes weaker, TOS becomes stronger, and the atmospheric stratification over oases gradually becomes unstable. Moreover, normalized total vertical advection heat flux remains almost unchanged with the SC fraction decreasing and TOS fraction increasing, but normalized total vertical advection flux for water vapor increases significantly in the upper CBL with the SC fraction decreasing in magnitude and TOS fraction increasing.

Due to SC and TOS, horizontal domain averages of the imbalance ratios for both sensible and latent heat fluxes are large at high levels, and decreases with the decreased height. Furthermore, with soil water content in oases increasing (decreasing), the imbalance ratio becomes larger (smaller) at high levels. Because of co-existence of SC and TOS, the EC flux measured at one point may not represent the total flux, especially just after irrigation in oases.

The sensitivity experiments demonstrate that soil temperature becomes slightly higher and soil water content becomes smaller for thin soil layer, and thus, surface sensible heat flux increases while latent heat flux decreases. The soil layer configuration has only a small impact on total, EC, and TOS latent heat fluxes, but almost no impact on the rest fluxes.

Finally, different horizontal grid resolutions, horizontal winds, and averaging periods used might have influ-enced SC and TOS, so further investigation is needed by using different horizontal wind speeds, domain sizes, and the degree of heterogeneity. In addition, although SC and TOS have been estimated by LES over the semi-idealized strip-like patches of desert and oasis in this paper,further study is also needed over real heterogeneous land surfaces.

Appendix

The variable [ φ]pis phase averaged as follows (Hadfield et al., 1991; Inagaki et al., 2006):

whereis the new coordinate aligned withx, which represents the position in relation to the wave cycle,=xmodλ.NxandNyare the total number ofx- andy-grids,respectively,Npis the total number of the waves in thexdirection,Np=NxΔx/λ, Δxis the grid spacing in thex-direction. The horizontal heterogeneity in thex-direction as a step function makes it possible to use phase averaging and temporal averaging to describe the vertical heat transport process of the SC. Thus,denotes the phase fluctuation and represents the deviation from the phase mean at a particular model grid.

The temporal and phase average of vertical velocity(w) (positive upward) and scalar (φ) can be expressed by(Kanda et al., 2004; Steinfeld et al., 2007; Huang et al.,2008)

Accordingly, the temporal and phase average results are written as

Thus, to obtain the horizontal domain average of imbalance ratio [I], ([Ih] for sensible heat flux, and [Iq] for latent heat flux), Eq. (A3) is substituted into the advection term

Acknowledgments. The observation data used in this paper were kindly provided by the MUlti-Scale Observation EXpEriment on Evapotranspiration over heterogeneous land surfaces 2012 (MUSOEXE-12) under the Heihe Watershed Allied Telemetry Experimental Research (HiWATER). See their website http://westdc.westgis.ac.cn/archives/news/sciencenews/archive-144.html for details.