Numerical Simulations of the Impacts of Mountain on Oasis Effects in Arid Central Asia

The oases in the mountain-basin systems of Central Asia are extremely fragile. Investigating oasis effects and oasis-desert interactions is important for understanding the ecological stability of oases. However, previous studies have been performed only in oasis-desert environments and have not considered the impacts of mountains. In this study, oasis effects were explored in the context of mountain effects in the northern Tianshan Mountains (NTM) using the Weather Research and Forecasting (WRF) model. Four numerical simulations are performed. The def simulation uses the default terrestrial datasets provided by the WRF model. The mod simulation uses actual terrestrial datasets from satellite products. The non-oasis simulation is a scenario simulation in which oasis areas are replaced by desert conditions, while all other conditions are the same as the mod simulation. Finally, the non-mountain simulation is a scenario simulation in which the elevation values of all grids are set to a constant value of 300 m, while all other conditions are the same as in the mod simulation. The mod simulation agrees well with near-surface measurements of temperature, relative humidity and latent heat flux. The Tianshan Mountains exert a cooling and wetting effects in the NTM region. The oasis breeze circulation (OBC) between oases and the deserts is counteracted by the stronger background circulation. Thus, the self-supporting mechanism of oases originating from the OBC plays a limited role in maintaining the ecological stability of oases in this mountain-basin system. However, the mountain wind causes the “cold-wet” island effects of the oases to extend into the oasis-desert transition zone at night, which is beneficial for plants in the transition region.


Introduction
Oases are common in the deserts of arid areas [1][2][3][4], especially in the hinterland of arid Central Asia (CA).The formation of an oasis in this region is closely related to the geomorphic characteristics of high mountain-basin systems, because the melting of snow and glaciers and the precipitation in tall mountain ranges provide necessary water resources for oasis survival and development in these mountain-basin systems [5].Accordingly, watersheds represent a basic unit in this mountain-oasis-desert system (MODS).The Northern Tianshan Mountains (NTM) is one of the typical geomorphic regions of CA and consists of a large number of complex MODSs.Although oases account for only a small proportion of the land surface (e.g., 4-5% in Xinjiang, a typical region in the hinterlands of Central Asia), more than 90% of the local population and >95% of the socioeconomic wealth are concentrated in oases [6].Thus, oases play a vital role in social and economic development.Since the 1950s, the region has experienced distinct and intense land reclamation, characterized mainly by oasis expansion.The total oasis area has expanded to 4.23 times its original value (from 121.0 × 10 4 ha in 1949 to 512.5 ×10 4 ha in 2010).As a result, a series of ecological problems have appeared due to the limited water and soil conditions [7][8][9], such as desertification and soil salinization [10,11].These problems have hindered the sustainable development of the oases.
Oasis effects are thought to play an important role in maintaining the existence of oases over time.In the recent literature, several studies based on both field observations and numerical simulations have focused on the exchanges of water and energy between oases and the surrounding deserts [12][13][14][15].The differences in the geographical and ecological characteristics between oasis areas and desert areas result in significant differences in their energy budgets, create exchanges of momentum and water vapor and affect the amounts of heat and moisture released to the atmosphere within the convective boundary layer.During daytime, solar radiation heats the desert surface and rapidly increases the near-surface temperatures over desert areas.In contrast, the temperature over oases is lower due to the intense level of evaporation occurring in oases [16].At the same time, the humidity over oases is higher than that over the desert.The air density and pressure in the near-surface layer over the desert are lower than over the oasis, thereby causing cold moist air to flow from the oasis to the desert, resulting in a moisture inversion over the desert [17].At higher elevation, the hot dry air over the desert flows toward the oasis, resulting in the oasis becoming a cold, wet island capped by warm, dry air.This phenomenon has been observed to create a temperature inversion layer and negative sensible heat fluxes over oases [12,[18][19][20] and is known as the "oasis effects" [16].The water and heat differences between the oasis and the desert cause local atmospheric circulation patterns, called oasis breeze circulation (OBC).Recent numerical simulations of OBC have been performed using regional climate models coupled to sophisticated land surface schemes, such as the fifth-generation Penn State / National Center for Atmospheric Research (NCAR) Mesoscale Model (MM5) or the Weather Research and Forecasting (WRF) model [6,17,[21][22][23][24][25][26][27][28].These studies have improved our understanding of the mechanisms and processes associated with the oasis-desert interactions.The conception of the "self-supporting mechanism of oasis" has been proposed based on OBC.Specifically, the downdrafts over an oasis decrease evaporation within the oasis, whereas the updrafts over the desert areas adjacent to the oasis act as a wall that hinders the transfer of cold, wet air to the dry surrounding desert [16,28].Thus, the OBC helps to maintain the ecological stability of oases by reducing the exchanges of heat and moisture between oases and the surrounding desert areas.
However, an oasis cannot survive without a supply of water from the nearby large mountains in a MODS.Therefore, studies should explore the oasis effects, i.e., the oasis-desert interactions, under the impacts of these surrounding mountains rather than focusing only on the oasis and desert areas [6,17,[21][22][23][24][25][26][27][28].Failing to consider the impacts of the mountains may result in an incorrect or incomplete understanding of the oasis-desert interactions within a MODS.For example, what are the spatial patterns of water and heat throughout MODS?Do intense mountain-valley winds occur within current MODSs [29,30]?Finally, how does the mountain-valley wind, if present, impact the oasis effects and OBC?The answers to these questions are not clear and deserve additional and comprehensive investigation.Since, once OBC ceases to operate within a MODS, the self-supporting mechanism of oasis originated from OBC will disappear and the water vapor over the oasis may flow into the surrounding deserts.Consequently, the oasis will become drier and drier, perhaps resulting in a vicious cycle of water deficiency and desertification [6].Accordingly, answering the above questions is essential, both theoretically and practically, to the sustainable development and ecological stability of oases and will provide useful information for further investigating the impact of oasis expansion on regional climate changes [31] under global warming.
The WRF model coupled to the sophisticated land surface scheme Noah has been used to simulate regional climates over the past 10 years [32,33].The WRF model is sensitive to land surface properties [34,35], especially primary parameters such as land cover (LC) and secondary parameters such as the green vegetation fraction (GVF), the leaf area index (LAI) and albedo, which can be determined from the LC.However, the default values of these parameters are usually from the outdated 1992-1993 Advanced Very High Resolution Radiometer (AVHRR) dataset.Realistic representation of these land surface properties is necessary to improve WRF-based simulations over complex areas, such as MODSs [36].In this study, actual LC, albedo, LAI and GVF data from the MODerate Resolution Imaging Spectroradiometer (MODIS) products are used.
Our ultimate goals in this paper are to realistically depict the temperature, humidity and circulation patterns of a complex MODS and to describe the impacts of mountain on oasis effects within MODS in Central Asia using the WRF model.The paper is structured as follows.Section 2 presents a description of the study area, the datasets used in this work, the WRF model and the design of the experiments performed here.Section 3 presents the results of the simulations.Section 4 offers a discussion and identifies avenues for future investigation.The final section summarizes the findings of this study.

Study Area
In this study, CA includes Xinjiang Province in China and five nations of Kazakhstan, Kyrgyzstan, Tajikistan, Turkmenistan and Uzbekistan [37,38].CA is controlled by both the westerly circulation in the middle-high latitudes and the polar air mass [39] and experiences an arid continental climate with scarce and concentrated rainfall (less than 250 mm in the basin regions and 900 mm in the upwind mountains).CA contains many of the highest-altitude mountains and basins in the world.Due to the unique topography of this area, the runoff from the melting of snow, glaciers and precipitation that occurs in the mountainous areas [39] represents the main sources of water for inland basins in the mountain-basin systems.The limited abundance of runoff promotes the development of relatively high-productivity oasis systems in the mountain-basin systems by irrigation (runoff flows into the oasis then disappears in the desert areas due to strong evapotranspiration and groundwater recharge).Because of the elevation differences and the importance of oases, the mountain-basin landscapes can be artificially divided into three regions: mountainous areas, oasis areas and desert areas [5].
The Tianshan Mountains, which are located in the hinterland of CA, act as the "water tower of Central Asia" [40].The NTM is used as the study area in this paper because it includes many the special characteristics of MODSs (Figure 1a) and is a microcosm of the terrain and climate of CA.Furthermore, the periphery of the NTM region has experienced intense human activity, such as oasis expansion, urbanization and resulting desertification over the past 60 years [41,42].The oases in the NTM region developed in the groundwater overflow zone between the foothills and the surrounding Gurban Tonggut desert (Figure 1b).The simulation region, spanning a total area of 260,400 km 2 , covers the whole NTM region, thus capturing the mountain effects.

Model Description and Configuration
The WRF model is an advanced mesoscale numerical weather prediction system designed for both climate research and numerical weather prediction [43].It was developed jointly by the NCAR and the National Centers for Environmental Prediction (NCEP).The Noah land surface model was also developed through multi-institutional collaborative efforts among NCEP, NCAR, the U.S. Air Force Weather Agency, NASA and the university community [44][45][46].The Noah model, which has been implemented in the WRF and MM5 models [47], provides reasonable diurnal variations of surface heat fluxes as surface boundary conditions for coupled models, as well as correct seasonal evolution patterns of soil moisture in the context of long-term data assimilation systems [32,48].In this study, WRF version 3.6 coupled with Noah was used in our simulations.
For this study, the model was configured for fine-scale simulation with two-way nested domains (D01 and D02, Figure 1a).The outer domain (D01) and the inner domain (D02) have grid spacing of approximately 18 km and 6 km, respectively, in the horizontal direction.The inner domain (D02) is our main region of interest.It covers a total area of approximately 330.78 km × 244.23 km and contains 50 × 50 grids cells.The ERA-Interim reanalysis data provides the initial and lateral boundary conditions for simulations at 6-h intervals.The simulation results from domain D01 are stored hourly and provide hourly lateral boundary conditions for D02.Part of the grids points in each side of D02 are discarded and the remaining points are used as the analysis (blue dotted line in Figure 1a).The initial soil moisture and the soil boundary temperatures are derived from the forcing data; however, the soil moisture levels of the oases and desert areas are adjusted using observations (for details, see

Model Description and Configuration
The WRF model is an advanced mesoscale numerical weather prediction system designed for both climate research and numerical weather prediction [43].It was developed jointly by the NCAR and the National Centers for Environmental Prediction (NCEP).The Noah land surface model was also developed through multi-institutional collaborative efforts among NCEP, NCAR, the U.S. Air Force Weather Agency, NASA and the university community [44][45][46].The Noah model, which has been implemented in the WRF and MM5 models [47], provides reasonable diurnal variations of surface heat fluxes as surface boundary conditions for coupled models, as well as correct seasonal evolution patterns of soil moisture in the context of long-term data assimilation systems [32,48].In this study, WRF version 3.6 coupled with Noah was used in our simulations.
For this study, the model was configured for fine-scale simulation with two-way nested domains (D01 and D02, Figure 1a).The outer domain (D01) and the inner domain (D02) have grid spacing of approximately 18 km and 6 km, respectively, in the horizontal direction.The inner domain (D02) is our main region of interest.It covers a total area of approximately 330.78 km × 244.23 km and contains 50 × 50 grids cells.The ERA-Interim reanalysis data provides the initial and lateral boundary conditions for simulations at 6-h intervals.The simulation results from domain D01 are stored hourly and provide hourly lateral boundary conditions for D02.Part of the grids points in each side of D02 are discarded and the remaining points are used as the analysis (blue dotted line in Figure 1a).The initial soil moisture and the soil boundary temperatures are derived from the forcing data; however, the soil moisture levels of the oases and desert areas are adjusted using observations (for details, see Section 2.4).Both domains are run with 35 unevenly spaced full eta levels with 15 levels below a height of 5 km above the ground surface.The model's upper surface was fixed at 50 hPa.A series of analyses that examine the model's sensitivity to the different parameterizations of the physical atmospheric processes operating over the study region was previously performed by Qiu et al. [49].In this study, we use the optimal WRF configurations.Namely, planetary boundary layer processes are resolved using the Yonsei University scheme [50], the microphysics are simulated using the WRF Single Moment-3 scheme [51], the cumulus clouds are simulated using the Kain-Fritsch Scheme [52] and the Community Atmospheric Model scheme is used to calculate the longwave and shortwave radiation [53].

Forcing Data and In Situ Measurements
The latest global atmospheric reanalysis product ERA-Interim [54] provides the initial and lateral boundary conditions for the WRF simulation in this paper.These data were produced by the European Centre for Medium-Range Weather Forecasts.We employ the fields describing geo-potential, relative humidity, temperature, the U and V wind components at 30 pressure levels and surface forcing parameters, including the 10-m U wind, the 10-m V wind, the 2-m dew point temperature, the 2-m temperature, the mean sea level pressure, the sea surface temperature, the sea-ice cover, the skin temperature, the snow density, the snow depth, the soil temperature and the soil water content in 4 different layers.The spatial resolution of the dataset is 0.75 • × 0.75 • and the data are reported at 6-h intervals.
In situ observations from four meteorological stations within domain D02 are used to validate the simulation results and these four stations are distributed in oasis areas with diverse LC types, including irrigated crops, urban areas and combinations of both (Table 1).Hourly in situ readings of the 2-m air temperature, the 2-m relative humidity and the 10-m wind speed and direction from four meteorological stations are used to validate the simulation results.In addition, observations of latent heat and downward shortwave and longwave radiation from an eddy covariance system at station S2 are also used.Due to a damaged radiation sensor, observations of sensible heat flux are not available.

Actual Land Surface Parameters
The default LC data used in the WRF model were originally provided by the U.S. Geological Survey (USGS) and were derived from AVHRR data from 1992 to 1993 [28].In these data, the land surface is assigned to 24 different LC categories or 19 categories derived from MODIS data for January 2001 to December 2001 [55].The GVF, LAI and albedo datasets used in the default WRF configuration are based on multi-year averages of monthly maps derived from AVHRR satellite data from 1985 to 1991 [33].The outdated datasets used in the WRF model (i.e., the default LC, albedo, LAI and GVF datasets) limit its ability to accurately represent land conditions over heterogeneous and complex ground surfaces due to the underestimations of the spatial heterogeneity in LC and the associated land surface properties [34,35].
In this study, we converted an actual LC dataset with a hierarchical classification system (Table S1) into the USGS classification system used in the WRF model [56], based on the corresponding relations shown in Table S2.The actual LC data were generated by the Xinjiang Institute of Ecology and Geography (XIEG), Chinese Academy of Sciences, based on Landsat images and a 1:1,000,000 scale topographic map [57].The default albedo, LAI and GVF provided by the WRF model have been replaced with the actual MODIS products MCD43A4, MYD15A2 and a GVF dataset computed from MOD13A2 [44,47,58], respectively.These products were downloaded from https://modis.gsfc.nasa.gov/data and correspond to 3 July 2012 and the strip numbers h23v04 and h24v04.We reprocessed these products using the same coordinate systems and resolutions as in the numerical simulations (the details are provided in Sections 1-4 of the Supplementary Materials).

Experimental Design
Four simulations were conducted as follows.(i) The def simulation uses the default LC, GVF, LAI, albedo datasets provided by the WRF model; (ii) In the mod simulation, the default LC, albedo, LAI and GVF datasets are updated using the MODIS products mentioned in Section 2.3.2;(iii) The non-oasis simulation is a sensitivity simulation in which the oasis areas are replaced with desert conditions and the other conditions are the same as in the mod simulation; (iv) Finally, the reference non-mountain simulation involves setting the elevation value of the grid cells to a constant value of 300 m, while the other conditions are the same as in the mod simulation.The differences between the def simulation and the mod simulation are used to examine the performance of the WRF model using actual biophysical parameters.The differences between the mod simulation and the non-oasis simulation and the differences between the mod simulation and the non-mountain simulation were used to quantitatively investigate the effects of oases and mountains, respectively, on temperature, humidity and circulation features within the complete MODS.
All of the experimental simulations are initialized during the period of the most vigorous vegetation growth and extend from 00:00 UTC on 1 July 2012, to 18:00 UTC on 31 July 2012.Generally, in the absence of accurate, gridded initial soil moisture conditions, a spin-up period is needed to allow the soil moisture within Noah to approach equilibrium with the hydrological cycle [32].The optimal spin-up period for any particular application is uncertain and may require simulations of considerable length, as years may be required to reach equilibrium [59].In this study, the soil moisture values of the oasis and desert areas are initialized via interpolation from observed soil moistures data from similar oasis and desert regions referenced in a previous paper [28] (Table 2).In addition, also following similar previous simulations of mesoscale water, surface energy and circulation processes [6,21,26], the first 15 days of the simulations are skipped as spin-up to enable the soil moisture of the other land use types in the MODS to reach relative equilibrium.The simulation results covering the periods from 19:00 UTC on 15 July 2012, to 18:00 UTC on 20 July 2012 (i.e., from 00:00 Beijing time on 16 July to 23:00 Beijing time on 20 July), are used for analysis.This 5-day period is characterized by anticyclonic and clear-sky conditions; thus, the effects of cloud distribution are excluded from the results (not shown).All of the spatial patterns of temperature, humidity and circulation during the daytime and nighttime represent the average of the simulation results over these 5 days.

Comparison of Actual LC, Albedo, LAI and GVF Data and the Corresponding Default Datasets
We examined the differences between the actual LC dataset and the default values during the analysis period (Figure 2a,b).Large urban areas and croplands in the NTM can be clearly observed in the actual LC dataset (Figure 2b)-however, these features are absent in the default dataset (Figure 2a), which instead displays large areas of scrubland and grassland.The default LC dataset is based on AVHRR satellite data from 1992 to 1993 [28]; however, the urban areas and croplands are much larger within the current study area, due to the overexploitation of water and soil resources over the past 20 years.The actual albedo values are generally higher than the default values, which were extracted from multi-year climatological monthly albedo levels [60], with an approximate value of 0.05.Several reasons may underlie these differences [61].The larger albedo values observed over the oasis areas in the actual dataset may be attributed to severe salinization [62,63] and considerable use of plastic-mulch [64], which have developed in the study area (Figure 2c,d).The differences in the LAI and the GVF between the actual and the default datasets are shown in Figure 2e-h, respectively.Substantial differences are evident within the oasis areas.The LAI values in default dataset range from 0.2 to 1.3, whereas the actual LAI values are generally higher and range from 0.8 to 5. The GVF values range from 15% to 55% in the default dataset (Figure 2g) and from 25% to 92% in the actual image (Figure 2h).The differences within the desert areas are much less noticeable, although slight differences occur near the northern border of the oasis.

Simulation Evaluation
Figure 3 shows a comparative analysis of the observed and simulated 2-m air temperature (T2), 2-m relative humidity (RH) and latent heat flux (LE) over the oasis surface.The mod simulation is closer to the observations than the def simulation.The LE computed from the mod simulation over the oasis surface provides a much better representation of the peak values and correctly reproduces the overall shape of the observed value.
The mean LE value at the daily maximum from the mod simulation is approximately 400 W/m 2 , which is two times greater than that of the def simulation (Figure 3c).The mod LE values are closer to the observations and the mean bias (modeled values minus observations) is up to 40.6 W/m 2 smaller than those of the def simulation over the whole period (not shown).Figure 4 compares the mean bias error (MBE), root mean square error (RMSE) and the Pearson correlation coefficient (r) computed from the mod and def simulations for T2 and RH.For the four sites, the mod simulation is clearly closer to the observed values and is associated with significant reductions in the MBE and the RMSE (the p-value is less than 0.01) and a significant increase in the Pearson correlation coefficient (the p-value is less than 0.01).The performance of the WRF model in the mod simulation shows considerable improvements due to the updates of the LC, albedo, LAI and GVF datasets.In fact, the bias in T2 is corrected to approximately 0.3-2.4• C and that of RH is corrected to 1.9-12.9% in the mod simulation.
Note that the reductions in the biases associated with T2, RH and LE are larger at stations S2 and S4 than at the other two locations shown in Figure 3.This pattern occurs because the actual LAI and GVF values at S2 and S4 improved to a much greater degree than those at the other two sites.Specifically, the LAI values at S2 and S4 changed from 1.08 and 0.52 to 4.8 and 3.0, respectively and the GVF values at S2 and S4 changed from 36.6% and 13.5% to 84.7% and 82.15%, respectively (Figure 2g,h).The LC types at S2 and S4 are homogenous cropland but those at S1 and S3 include cropland combined with urban areas (S1) or desert areas (S3), respectively.As a result, more consistent trends and similar magnitudes of T2, RH and LE are observed at the four stations in the mod simulation than in the def simulation, as shown in Figure 2.This validation indicates that the simulated temperature, humidity and latent heat flux trends obtained from the mod simulation are consistent with the observations and realistically represent the water and heat conditions of the oasis surface.A comparison of observed and simulated 10-m horizontal wind at the four meteorological stations, as represented by the mod and def simulations, is shown in Figure 5. Stations S1 and S2 are located in the upper part of oasis near its southern border and stations S3 and S4 are located in the lower part of the oasis near its northern border.Most of the simulated wind speeds (WS) values are slightly higher than the observed values; 60 percent of the simulated WS values lie within a range of 2-5 m/s, whereas the observed values range from 2-4 m/s.The reason for this bias in Ws is that the uncertainty of randomized turbulence processes hinders the accurate simulation of wind [65].However, the trend in wind directions (WD) agrees well with the observations.The observed dominant WD is WNW or NW during daytime and WSW or SW at night for all four stations.The simulated WD shows direction and diurnal changes that are consistent with the observations.This validation demonstrates that the simulated WD agrees well with the observation and represents the circulation characteristics of the mountain-valley wind in the study area.A comparison of observed and simulated 10-m horizontal wind at the four meteorological stations, as represented by the mod and def simulations, is shown in Figure 5. Stations S1 and S2 are located in the upper part of oasis near its southern border and stations S3 and S4 are located in the lower part of the oasis near its northern border.Most of the simulated wind speeds (W S ) values are slightly higher than the observed values; 60 percent of the simulated W S values lie within a range of 2-5 m/s, whereas the observed values range from 2 to 4 m/s.The reason for this bias in Ws is that the uncertainty of randomized turbulence processes hinders the accurate simulation of wind [65].However, the trend in wind directions (W D ) agrees well with the observations.The observed dominant W D is WNW or NW during daytime and WSW or SW at night for all four stations.The simulated W D shows direction and diurnal changes that are consistent with the observations.This validation demonstrates that the simulated W D agrees well with the observation and represents the circulation characteristics of the mountain-valley wind in the study area.A comparison of observed and simulated 10-m horizontal wind at the four meteorological stations, as represented by the mod and def simulations, is shown in Figure 5. Stations S1 and S2 are located in the upper part of oasis near its southern border and stations S3 and S4 are located in the lower part of the oasis near its northern border.Most of the simulated wind speeds (WS) values are slightly higher than the observed values; 60 percent of the simulated WS values lie within a range of 2-5 m/s, whereas the observed values range from 2-4 m/s.The reason for this bias in Ws is that the uncertainty of randomized turbulence processes hinders the accurate simulation of wind [65].However, the trend in wind directions (WD) agrees well with the observations.The observed dominant WD is WNW or NW during daytime and WSW or SW at night for all four stations.The simulated WD shows direction and diurnal changes that are consistent with the observations.This validation demonstrates that the simulated WD agrees well with the observation and represents the circulation characteristics of the mountain-valley wind in the study area.

Spatial Patterns of Air Temperature
The 2-m air temperatures from both the mod and non-oasis simulations generally exhibit continuous stripe-like increases from the mountainous areas to the oasis areas and to the desert areas (Figure 6a-d).This pattern indicates that the spatial temperature pattern seen in the MODS reflects the lapse rate of temperature due to the elevation gradient difference within the MODS.An obvious temperature gradient line is observed along the northern boundary of the oasis in the mod simulation (Figure 6a,b) but not in the non-oasis simulation (Figure 6c,d), indicating that the water and heat differences between the oasis and desert areas obviously disturb the stripe-like increases in temperature in the MODS.

Spatial Patterns of Air Temperature
The 2-m air temperatures from both the mod and non-oasis simulations generally exhibit continuous stripe-like increases from the mountainous areas to the oasis areas and to the desert areas (Figure 6a-d).This pattern indicates that the spatial temperature pattern seen in the MODS reflects the lapse rate of temperature due to the elevation gradient difference within the MODS.An obvious temperature gradient line is observed along the northern boundary of the oasis in the mod simulation (Figure 6a,b) but not in the non-oasis simulation (Figure 6c,d), indicating that the water and heat differences between the oasis and desert areas obviously disturb the stripe-like increases in temperature in the MODS.(e,f) and the differences between the mod simulation and non-mountain simulation; and (g,h) the differences between the mod simulation and the non-oasis simulation during daytime and nighttime, respectively.The black dotted line represents the border of the oasis and the white dotted line (a) represents the location of the vertical section (longitude 85.7° E) represented in Figure 7.
The mean near-surface air temperatures over the oasis area are 30.44 °C and 23.36 °C during the daytime and nighttime, respectively (Table 3).In contrast, the near-surface air temperatures over the desert are approximately 32.5 °C and 28.0 °C, respectively (Figure 6a,b).The differences in the nearsurface air temperatures between the oasis and desert areas indeed indicate that the oasis areas represent a "cold island" compared to the surrounding desert during daytime and nighttime.The intensity of this "cold island" resulting from the oasis cooling effect is as large as −0.61 °C in daytime and as large as −3.37 °C at night, as computed using the mean oasis temperature difference between the mod simulation and non-oasis simulation (these differences are statistically significant at p < 0.01, Table 3).During the daytime, because the soil moisture and thermal capacity values are higher in the oasis areas than the desert areas, more energy is required to increase the temperatures over the oasis areas and evapotranspiration takes place in the oasis, causing cooling [28] (as can be seen in Figure S5a,c).The possible reasons for the more intense cooling effect of the oasis at night include evaporation from the surface soil (Figure S4f) and the lapse rate of temperature; the slightly higher elevation of the oasis areas (Figure 6e,f) enhances the intensity of this "cold island" effect.Note that the near-surface temperature patterns from the non-mountain simulation (Figure S1a,b) and the difference patterns between the mod simulation and the non-mountain simulation (Figure 6e,f) are also  (e,f) and the differences between the mod simulation and non-mountain simulation; and (g,h) the differences between the mod simulation and the non-oasis simulation during daytime and nighttime, respectively.The black dotted line represents the border of the oasis and the white dotted line (a) represents the location of the vertical section (longitude 85.7 • E) represented in Figure 7.
The mean near-surface air temperatures over the oasis area are 30.44 • C and 23.36 • C during the daytime and nighttime, respectively (Table 3).In contrast, the near-surface air temperatures over the desert are approximately 32.5 • C and 28.0 • C, respectively (Figure 6a,b).The differences in the near-surface air temperatures between the oasis and desert areas indeed indicate that the oasis areas represent a "cold island" compared to the surrounding desert during daytime and nighttime.The intensity of this "cold island" resulting from the oasis cooling effect is as large as −0.61 • C in daytime and as large as −3.37 • C at night, as computed using the mean oasis temperature difference between the mod simulation and non-oasis simulation (these differences are statistically significant at p < 0.01, Table 3).During the daytime, because the soil moisture and thermal capacity values are higher in the oasis areas than the desert areas, more energy is required to increase the temperatures over the oasis areas and evapotranspiration takes place in the oasis, causing cooling [28] (as can be seen in Figure S5a,c).The possible reasons for the more intense cooling effect of the oasis at night include evaporation from the surface soil (Figure S4f) and the lapse rate of temperature; the slightly higher elevation of the oasis areas (Figure 6e,f) enhances the intensity of this "cold island" effect.Note that the near-surface temperature patterns from the non-mountain simulation (Figure S1a,b) and the difference patterns between the mod simulation and the non-mountain simulation (Figure 6e,f) are also consistent with circulation patterns, illustrating that the water and heat patterns in the study area are controlled by larger-scale circulation and the climate background in the absence of high-relief terrain.
consistent with circulation patterns, illustrating that the water and heat patterns in the study area are controlled by larger-scale circulation and the climate background in the absence of high-relief terrain.The continuous stripe-like temperature patterns and the oasis cooling effect are further illustrated by the vertical section of temperature patterns from the mod and non-oasis simulations and the temperature differences between the mod and the non-oasis simulations and between the mod and non-mountain simulations (Figure 7).In the daytime, the air temperatures in the oasis area in the mod simulation are clearly lower than those over the northern desert area and the southern mountainous area at the same altitude of approximately 1500 m and this temperature pattern forms a concaveshaped vertical temperature pattern in the MODS (Figure 7a).At night, a clearly defined temperature inversion layer (TIL) is observed near the surface over the oasis area.Its thickness varies from 200 m in the southern part of the oasis with an elevation of 450 m to 600 m in the northern part with an elevation of 300 m (Figure 7b).However, the concave-shaped vertical temperature pattern observed during the daytime and the TIL that developed during the nighttime are not observed in the non-oasis simulation (Figure 7c,d), indicating the concave low-temperature center and the TIL over the oasis areas result from the oasis cooling effect.Based on the differences between the mod and non-oasis  The continuous stripe-like temperature patterns and the oasis cooling effect are further illustrated by the vertical section of temperature patterns from the mod and non-oasis simulations and the temperature differences between the mod and the non-oasis simulations and between the mod and non-mountain simulations (Figure 7).In the daytime, the air temperatures in the oasis area in the mod simulation are clearly lower than those over the northern desert area and the southern mountainous area at the same altitude of approximately 1500 m and this temperature pattern forms a concave-shaped vertical temperature pattern in the MODS (Figure 7a).At night, a clearly defined temperature inversion layer (TIL) is observed near the surface over the oasis area.Its thickness varies from 200 m in the southern part of the oasis with an elevation of 450 m to 600 m in the northern part with an elevation of 300 m (Figure 7b).However, the concave-shaped vertical temperature pattern observed during the daytime and the TIL that developed during the nighttime are not observed in the non-oasis simulation (Figure 7c,d), indicating the concave low-temperature center and the TIL over the oasis areas result from the oasis cooling effect.Based on the differences between the mod and non-oasis simulations, the "cold island" effect of the oasis extends to a height of approximately 1.5 km above the oasis surface, as shown by the n-shaped cooling center (Figure 7d,g).Moreover, the cooling effect of the oasis extends horizontally into the desert-oasis transition areas to a distance of approximately 25 km (Figures 6h and 7h) at night due to the mountain-wind.In addition, the existence of high-relief terrain (mountains) decreases the whole temperature pattern over the basin during the daytime and the nighttime, based on the differences between the non-oasis and non-mountain simulations (Figures 6 and 7e,f).

Spatial Patterns of Humidity
The mean near-surface specific humidity (Q2) values over the oasis in the mod simulation are 9.99 g kg −1 and 10.05 g kg −1 during the daytime and nighttime, respectively (Table 3).In contrast, values of approximately 8.4 g kg −1 (Figure 8a) and 9.9 g kg −1 (Figure 8b) occur at almost the same elevation in the surrounding desert areas, respectively.The oasis indeed represents a "wet island" compared with the surrounding desert at almost the same elevation but the mechanism underlying the "wet island" effect differs between daytime and nighttime.
Atmosphere 2017, 8, 212 13 of 20 simulations, the "cold island" effect of the oasis extends to a height of approximately 1.5 km above the oasis surface, as shown by the n-shaped cooling center (Figure 7d,g).Moreover, the cooling effect of the oasis extends horizontally into the desert-oasis transition areas to a distance of approximately 25 km (Figures 6h and 7h) at night due to the mountain-wind.In addition, the existence of high-relief terrain (mountains) decreases the whole temperature pattern over the basin during the daytime and the nighttime, based on the differences between the non-oasis and non-mountain simulations (Figures 6 and 7e,f).

Spatial Patterns of Humidity
The mean near-surface specific humidity (Q2) values over the oasis in the mod simulation are 9.99 g kg −1 and 10.05 g kg −1 during the daytime and nighttime, respectively (Table 3).In contrast, values of approximately 8.4 g kg −1 (Figure 8a) and 9.9 g kg −1 (Figure 8b) occur at almost the same elevation in the surrounding desert areas, respectively.The oasis indeed represents a "wet island" compared with the surrounding desert at almost the same elevation but the mechanism underlying the "wet island" effect differs between daytime and nighttime.A clearly defined humidity gradient line along the northern boundary of the oasis during the daytime is observed in the mod simulation (Figure 8a) but the gradient line is not observed at night in the mod simulation (Figure 8b) or in the non-oasis simulation (Figure 8c,d).The patterns of differences between the mod and the non-oasis simulations show that there is indeed a remarkable difference in specific humidity over the oasis areas during the daytime, whereas a slight difference occurs along the northern oasis boundary and in the desert area at night (Figure 8g,h).Based on Student's t-test, the mean humidity difference over the oasis between the mod simulation and the nonoasis simulation is extremely statistically significant (the p-value is far less than 0.01) during the daytime but statistically insignificant (the p-value is 0.075) at night (Table 3).This result indicates that the "wet island" results from the wetting effect of the oasis during the daytime but the cause during the nighttime is uncertain.The mean intensity of this "wet island" resulting from the oasis wetting effect is 1.56 g kg −1 during the daytime (Table 3).In addition, nighttime moisture differences between the mod and the non-oasis simulations are mainly located along the northern boundary of the oasis area, rather than over the oasis itself (Figure 7h).This result appears to be consistent with the circulation patterns, suggesting that the wetting effect of the oasis might be present at night but that A clearly defined humidity gradient line along the northern boundary of the oasis during the daytime is observed in the mod simulation (Figure 8a) but the gradient line is not observed at night in the mod simulation (Figure 8b) or in the non-oasis simulation (Figure 8c,d).The patterns of differences between the mod and the non-oasis simulations show that there is indeed a remarkable difference in specific humidity over the oasis areas during the daytime, whereas a slight difference occurs along the northern oasis boundary and in the desert area at night (Figure 8g,h).Based on Student's t-test, the mean humidity difference over the oasis between the mod simulation and the non-oasis simulation is extremely statistically significant (the p-value is far less than 0.01) during the daytime but statistically insignificant (the p-value is 0.075) at night (Table 3).This result indicates that the "wet island" results from the wetting effect of the oasis during the daytime but the cause during the nighttime is uncertain.The mean intensity of this "wet island" resulting from the oasis wetting effect is 1.56 g kg −1 during the daytime (Table 3).In addition, nighttime moisture differences between the mod and the non-oasis simulations are mainly located along the northern boundary of the oasis area, rather than over the oasis itself (Figure 7h).This result appears to be consistent with the circulation patterns, suggesting that the wetting effect of the oasis might be present at night but that the spatial extent of the wetting effect changes due to the mountain wind.Moreover, the surface specific humidity patterns from the non-mountain simulation (Figure S2a,b) and the differences between the mod and the non-mountain simulations (Figure 8e,f) are also consistent with the circulation patterns.Thus, the water and heat patterns in the study area would be controlled by larger-scale circulation patterns and the climate background in the absence of high-relief terrain.
Figure 9 shows the vertical section of humidity patterns from the mod and non-oasis simulations, as well as the differences between the mod and non-oasis simulations and between the mod and non-mountain simulations.The specific humidity in the mountainous area at an equivalent altitude of 1 km during the daytime and nighttime is higher than that over the basin area at the same elevation of 1 km in the non-oasis simulation (Figure 9c,d).This higher specific humidity in the mountainous areas may be due to the wetting effect of daytime evaporation from forest and grassland and the radiative cooling of mountainous areas decreases the saturated vapor pressure, leading to a higher specific humidity [66] during nighttime.
The differences between the mod and non-oasis simulations (Figure 9g,h) and between the mod and non-mountain simulations (Figure 9e,f) provide further confirmation of the wetting effects of oases and mountains based on the decrease in the cooling gradient from the oasis surface over a distance of approximately 1 km (Figure 9g,h) and from the mountainous areas to their surroundings (Figure 9e,f).Moreover, the wetting effect of the oasis areas at night extends horizontally into the desert areas to a distance of approximately 25 km (Figures 8h and 9h).This horizontal extension results primarily from the ground-hugging nighttime mountain wind, since the extension is consistent with the direction of the mountain wind (Figure 10b,d).
Atmosphere 2017, 8, 212 14 of 20 specific humidity patterns from the non-mountain simulation (Figure S2a,b) and the differences between the mod and the non-mountain simulations (Figure 8e,f) are also consistent with the circulation patterns.Thus, the water and heat patterns in the study area would be controlled by larger-scale circulation patterns and the climate background in the absence of high-relief terrain.
Figure 9 shows the vertical section of humidity patterns from the mod and non-oasis simulations, as well as the differences between the mod and non-oasis simulations and between the mod and nonmountain simulations.The specific humidity in the mountainous area at an equivalent altitude of 1 km during the daytime and nighttime is higher than that over the basin area at the same elevation of 1 km in the non-oasis simulation (Figure 9c,d).This higher specific humidity in the mountainous areas may be due to the wetting effect of daytime evaporation from forest and grassland and the radiative cooling of mountainous areas decreases the saturated vapor pressure, leading to a higher specific humidity [66] during nighttime.
The differences between the mod and non-oasis simulations (Figure 9g,h) and between the mod and non-mountain simulations (Figure 9e,f) provide further confirmation of the wetting effects of oases and mountains based on the decrease in the cooling gradient from the oasis surface over a distance of approximately 1 km (Figure 9g,h) and from the mountainous areas to their surroundings (Figure 9e,f).Moreover, the wetting effect of the oasis areas at night extends horizontally into the desert areas to a distance of approximately 25 km (Figures 8h and 9h).This horizontal extension results primarily from the ground-hugging nighttime mountain wind, since the extension is consistent with the direction of the mountain wind (Figure 10b,d).

Spatial Patterns of the 10-m Horizontal Circulation
Figure 10 shows the diurnal mean WS and WD patterns from the mod, non-oasis and non-mountain simulations, as well as the differences between the mod and non-oasis simulations.In the context of the prevailing westerly wind, a wind direction of WNW with a wind speed of approximately 3 m/s occurs during the daytime (Figure 10a,c), whereas a wind direction of SW with a wind speed of approximately 4 m/s occurs at night (Figure 9b,d), as observed in both the mod and non-oasis simulations covering the MODS.An intense southwesterly wind with a wind speed exceeding 7 m/s is observed during both daytime and nighttime in the non-mountain simulation (Figure 10e,f) and the intense southwesterly wind represents a combination of the prevailing westerly wind and the summer subtropical jet stream resulting from the Southwest Asian low-pressure system over northeastern India [39,67].However, the intense southwesterly wind is blocked in the mod and nonoasis simulations.Therefore, the diurnal change in the wind direction from WNW to SW in the mod and non-oasis simulations is produced by the high-relief terrain, further demonstrating the existence of mountain-valley winds in the MODS.
The difference between the mod simulation and the non-oasis simulation is characterized by an obvious divergent wind blowing from the center of the oasis area to the surrounding desert with a speed of less than approximately 1.5 m/s at a height of 10 m height (Figure 10g,h).This divergent wind is oriented in the opposite direction of the combination of the mountain-valley wind and the prevailing westerly wind.However, no obvious downdrafts or updrafts occur over the oasis or deserts based on the vertical component difference between the mod and non-oasis simulations (Figure S3).This divergence indicates that OBC may indeed exist between the oasis and desert areas

Spatial Patterns of the 10-m Horizontal Circulation
Figure 10 shows the diurnal mean W S and W D patterns from the mod, non-oasis and non-mountain simulations, as well as the differences between the mod and non-oasis simulations.In the context of the prevailing westerly wind, a wind direction of WNW with a wind speed of approximately 3 m/s occurs during the daytime (Figure 10a,c), whereas a wind direction of SW with a wind speed of approximately 4 m/s occurs at night (Figure 9b,d), as observed in both the mod and non-oasis simulations covering the MODS.An intense southwesterly wind with a wind speed exceeding 7 m/s is observed during both daytime and nighttime in the non-mountain simulation (Figure 10e,f) and the intense southwesterly wind represents a combination of the prevailing westerly wind and the summer subtropical jet stream resulting from the Southwest Asian low-pressure system over northeastern India [39,67].However, the intense southwesterly wind is blocked in the mod and non-oasis simulations.Therefore, the diurnal change in the wind direction from WNW to SW in the mod and non-oasis simulations is produced by the high-relief terrain, further demonstrating the existence of mountain-valley winds in the MODS.
The difference between the mod simulation and the non-oasis simulation is characterized by an obvious divergent wind blowing from the center of the oasis area to the surrounding desert with a speed of less than approximately 1.5 m/s at a height of 10 m height (Figure 10g,h).This divergent wind is oriented in the opposite direction of the combination of the mountain-valley wind and the prevailing westerly wind.However, no obvious downdrafts or updrafts occur over the oasis or deserts based on the vertical component difference between the mod and non-oasis simulations (Figure S3).This divergence indicates that OBC may indeed exist between the oasis and desert areas and the possible OBC is counteracted by the stronger background circulation produced by the combination of the mountain-valley wind and the prevailing westerly wind.

Comparison with Previous Studies and an Assessment of Model Uncertainties
In this study, the effects of oases on temperature, humidity and regional circulation patterns within an arid MODS in the NTM were investigated using the WRF model with actual LC, albedo, LAI and GVF datasets.The mod simulation provides a significantly better representation of diurnal temperature and humidity due to use of updated datasets and better reflects the exchanges of water and heat in the MODS.
The differences between the mod and non-oasis simulations show the typical "cold island" effect and "wet island" effect of an oasis during the day.Thus, a TIL develops over the oasis areas at night from the mod simulation.The oasis cooling-wetting effect reaches a vertical height of approximately 1.5 km above the oasis surface and this result is similar to findings reported in previous investigations [6,17,21,24,[26][27][28].However, our study also contributes certain new findings.The Tianshan Mountains act as the background of the oasis-desert interaction and exert a cooling and wetting effect on the whole MODS in the NTM region, as evidenced by the differences between the mod and the non-mountain simulations.Moreover, the horizontal range of the oasis wetting effect appears to extend into the desert areas, rather than being restricted over the oasis areas, due to the presence of the mountain wind at night.This horizontal extension may produce favorable conditions for desert plants in the oasis-desert transition zone.
In previous studies, researchers proposed the concept that oases have a self-supporting mechanism originating from OBC.They argued that the OBC between an oasis and the surrounding desert helps maintain the ecological stability of the oasis by reducing heat and moisture exchange [21,24,26,28].The results from the mod simulation and the differences between the mod and non-oasis simulations show that the water and heat differences between the oasis and surrounding desert do produce OBC (Figure 10e,f).However, OBC is not observed in the actual weather or climate conditions (the mod simulation in this study) (Figure 10a,b) because it is counteracted by the stronger background circulation that results from a combination of mountain-valley wind and the prevailing westerly wind in the MODS.Therefore, the effects of OBC on the ecological stability of oasis systems within an MODS may have been overestimated by previous studies.Consequently, studies related to oasis effects should consider the regional context of the MODS, if the oasis is located in a large mountain-basin system.Such investigations will help provide a better understanding of the ecological and climatic effects of oasis systems but also may result in the identification of new features of such systems.In addition, this study represents necessary preliminary work for further investigating the contributions of oasis expansion to historical and future climate changes at local or regional scales.Whether and how the typical wetting and cooling effects of oases impact local or regional climate change will be further investigated in the next stage of this work, which will include long-period simulations.
Some limitations in the modeling must be examined in future studies.For example, we found that the simulated temperatures were generally somewhat overestimated and that the relative humidity was slightly underestimated.These errors may be attributed to the evaporation from soils resulting from occasional irrigation.This process was not considered in the simulations [6,28] because the WRF model describes irrigation and irrigated areas in terms of different biophysical parameters between irrigated cropland and rain-fed cropland, instead of using a physical or mechanistic description of irrigation.An irrigation scheme will be added and a sensitivity analysis will be carried out in the future work to determine whether this bias can be corrected.

Conclusions
In this study, the effects of oasis areas were investigated within a complete MODS located in the NTM using the WRF model.Studies of oasis effects should consider the entire MODS, which will help to fully understand the ecological and climatic effects of oasis systems.The typical "cold-wet" island effects of the oasis, nighttime TIL over oasis areas and OBC are fully described within the context of the MODS in this study through the use of the WRF model with actual LC, albedo, LAI and GVF datasets.The Tianshan Mountains exert a cooling and wetting effects on the whole MODS in the NTM region.The mountain wind causes the "cold-wet" island effects of the oasis to extend into the surrounding desert area at night, which may produce conditions favorable for desert plants in the oasis-desert transition zone.However, the OBC was counteracted by stronger background resulting from a combination of the prevailing westerly winds and the mountain-valley wind in the MODS.The proposed self-supporting mechanism of the oasis associated with the OBC cannot play a significant role in maintaining the stability of the oasis in this MODS.

Figure 1 .
Figure 1.Topography (a) and land cover (b) of the study area in 2012 and the Weather Research and Forecasting (WRF) simulation domains (blue dotted line in (a) shows the horizontal extent of the analyzed area in the subsequent figures).

Figure 1 .
Figure 1.Topography (a) and land cover (b) of the study area in 2012 and the Weather Research and Forecasting (WRF) simulation domains (blue dotted line in (a) shows the horizontal extent of the analyzed area in the subsequent figures).

Figure 4 .
Figure 4.The mean bias error (MBE), root mean squared error (RMSE) and Pearson correlation coefficient (r) between the simulated and observed (a) air temperature and (b) relative humidity at 2 m at the four stations (S1, S2, S3 and S4).

Figure 5 .
Figure 5.Comparison between the observed and simulated (a) wind speed (WS) and (b) wind direction (WD) from the mod and def simulations at the four stations (S1, S2, S3 and S4) during 16 July to 20 July.

Figure 4 .
Figure 4.The mean bias error (MBE), root mean squared error (RMSE) and Pearson correlation coefficient (r) between the simulated and observed (a) air temperature and (b) relative humidity at 2 m at the four stations (S1, S2, S3 and S4).

Figure 5 .
Figure 5.Comparison between the observed and simulated (a) wind speed (WS) and (b) wind direction (WD) from the mod and def simulations at the four stations (S1, S2, S3 and S4) during 16 July to 20 July.

Figure 5 .
Figure 5.Comparison between the observed and simulated (a) wind speed (W S ) and (b) wind direction (W D ) from the mod and def simulations at the four stations (S1, S2, S3 and S4) during 16 July to 20 July.

Figure 6 .
Figure 6.Mean 2-m air temperatures (T2) from the (a,b) mod simulation and (c,d) the non-oasis simulation;(e,f) and the differences between the mod simulation and non-mountain simulation; and (g,h) the differences between the mod simulation and the non-oasis simulation during daytime and nighttime, respectively.The black dotted line represents the border of the oasis and the white dotted line (a) represents the location of the vertical section (longitude 85.7° E) represented in Figure7.

Figure 6 .
Figure 6.Mean 2-m air temperatures (T2) from the (a,b) mod simulation and (c,d) the non-oasis simulation;(e,f) and the differences between the mod simulation and non-mountain simulation; and (g,h) the differences between the mod simulation and the non-oasis simulation during daytime and nighttime, respectively.The black dotted line represents the border of the oasis and the white dotted line (a) represents the location of the vertical section (longitude 85.7 • E) represented in Figure7.

Figure 7 .
Figure 7.The same as Figure 6 but showing for a vertical cross section of air temperature along the white line shown in Figure 6a.

Figure 7 .
Figure 7.The same as Figure 6 but showing for a vertical cross section of air temperature along the white line shown in Figure 6a.

Figure 8 .
Figure 8.The same as Figure 6 but for surface specific humidity Q2 (g kg −1 ).

Figure 8 .
Figure 8.The same as Figure 6 but for surface specific humidity Q2 (g kg −1 ).

Figure 9 .
Figure 9. Same as Figure 6 but for the vertical section of specific humidity Q (g kg −1 ) along the white line (longitude 85.7° E) shown in Figure 6a.Figure 9. Same as Figure 6 but for the vertical section of specific humidity Q (g kg −1 ) along the white line (longitude 85.7 • E) shown in Figure 6a.

Figure 9 .
Figure 9. Same as Figure 6 but for the vertical section of specific humidity Q (g kg −1 ) along the white line (longitude 85.7° E) shown in Figure 6a.Figure 9. Same as Figure 6 but for the vertical section of specific humidity Q (g kg −1 ) along the white line (longitude 85.7 • E) shown in Figure 6a.

Figure 10 .
Figure 10.The mean diurnal WS and WD patterns (mean wind speed and direction at a height of 10 m during the daytime and nighttime) in the mod, non-oasis and non-mountain simulations and the differences between the mod and non-oasis simulations.(a) Daytime Ws and WD from the mod simulation; (b) nighttime Ws and WD from the mod simulation; (c) daytime Ws and WD from the nonoasis simulation; (d) nighttime Ws and WD from the non-oasis simulation; (e) daytime Ws and WD from the non-mountain simulation; (f) nighttime Ws and WD from the non-mountain simulation; (g) daytime Ws and WD difference between the mod and non-oasis simulations; (h) nighttime Ws and WD difference between the mod and non-oasis simulations.The red line represents the border of the oasis.

Figure 10 .
Figure 10.The mean diurnal W S and W D patterns (mean wind speed and direction at a height of 10 m during the daytime and nighttime) in the mod, non-oasis and non-mountain simulations and the differences between the mod and non-oasis simulations.(a) Daytime Ws and W D from the mod simulation; (b) nighttime Ws and W D from the mod simulation; (c) daytime Ws and W D from the non-oasis simulation; (d) nighttime Ws and W D from the non-oasis simulation; (e) daytime Ws and W D from the non-mountain simulation; (f) nighttime Ws and W D from the non-mountain simulation; (g) daytime Ws and W D difference between the mod and non-oasis simulations; (h) nighttime Ws and W D difference between the mod and non-oasis simulations.The red line represents the border of the oasis.

Table 1 . Station information and parameters measured at each station. Station Longitude ( • E) Latitude ( • N) Altitude (m) Land Cover Type Time Interval Parameters
D T2 represents the 2-m air temperature; RH represents the 2-m relative humidity; LE represents the latent heat flux; W S and W D represent the wind speed and direction at 10 m, respectively; SW and LW represent downward shortwave and longwave radiation.

Table 2 .
[28] moisture values for the oasis and desert areas in the four Noah soil layers[28].

Table 3 .
Statistical significance of the differences in temperature and humidity between the mod and non-oasis simulations, as assessed using Student's t-test.

Periods Statistics Temperature (°C) Specific Humidity (g kg −1 ) mod non-oasis mod non-oasis
* Indicates that the difference between the mod and the non-oasis simulations is extremely significant.

Table 3 .
Statistical significance of the differences in temperature and humidity between the mod and non-oasis simulations, as assessed using Student's t-test.

Periods Statistics Temperature ( • C) Specific Humidity (g kg −1 ) mod non-oasis mod non-oasis
* Indicates that the difference between the mod and the non-oasis simulations is extremely significant.