Abstract

Simulations with the WRF model have been carried out with high resolution soil data to analyze its effect on planetary boundary layer (PBL) development. The default soil texture distribution of 5′ horizontal resolution has been replaced with a 30′′ one on the basis of the Digital Kreybig Soil Information System and Hungarian Agrogeological Database in Hungary. Soil parameter values determined from HUNSODA and MARTHA Hungarian soil databases were also compared. Comparison of PBL height simulations and measurements obtained by radiometer and windprofiler shows that the impact of soil parameter differences on PBL height evolution is not negligible. The latent heat flux and PBL height daytime courses show significant ( 𝑃 < 0 . 0 1 ) differences over more than 50% of the model domain covering the Carpathian basin.

1. Introduction

The land-surface is the key climate system component in modeling of the lower atmosphere not only in the global circulation but also in the mesoscale modeling systems [1]. Since surface properties are highly variable in both space and time, the quantification of land-surface-atmosphere interaction processes is still a major challenge. The surface-atmosphere coupling is ensured by the surface energy budget equation. So, surface inhomogeneities can induce fluctuations in the spatially averaged surface energy budget terms, these in return affect the lower atmosphere [2], even up to the top of the planetary boundary layer (PBL) [3]. One major factor in the energy budget is the available soil moisture, affecting both bare soil evaporation and transpiration. However, soil moisture depends on many factors, such as soil texture, spatial variability of soil properties, and soil hydrophysical parameters. Soil parameters inherently differ between independent and local soil datasets, where the local dataset proves to be more accurate on a region’s characterization (e.g., [4, 5]). It is also known that parameter variability in soil textures can be greater than the differences in soil textures [6]. Though for soil hydraulic modeling this information is important, it arises the question whether the resulting difference in latent heat flux has substantial effect on PBL height. There had been many investigations on the role of soil in weather events (e.g., [7, 8]), but more investigation is needed (e.g., [9, 10]) as the mesoscale weather prediction systems develop.

In this study the popular WRF model is used to analyze the effects of soil database differences on PBL height changes. The simulations are organized as follows: a 25600 km2 model domain with 1 km resolution was nested into a 450000 km2 large model domain with 3 km resolution. This territory is situated approximately in center of the Carpathian-basin. In order to be able to analyze the effect of soil parameters on PBL changes, a soil texture distribution of very high resolution (30′′) has also been added. The hydrophysical parameters of soil textures are defined with the use of two Hungarian soil databases. Both contain soil samples from Hungary, but there is one order of magnitude difference in their sizes. This remarkable difference results in parameter value differences. Simulated PBL height changes are compared to radiometric and windprofiler determined PBL heights. Moreover, significance test on the simulated daytime course of latent heat flux and PBL height is also done to asses the effect of soil parameters.

2. Methods

2.1. Estimation of PBL Height

The height of PBL can be determined from various methods [11, 12]. The uncertainty range of PBL heights obtained by different instruments is 100–700 m. From many methodological options available for windprofiler, we tested separately the maximum signal-to-noise ratio (SNR) [13], the maximum windshear height, and the wind speed and wind direction change method. In most cases the SNR proved to be acceptable, clearly showing diurnal evolution of PBL; therefore, we choose to compare simulations to this method. Evaluating radiometric measurement data, we used the following PBL height determination criteria: the level of disappearing inversion during nighttime and the level of 0.4 K/100 m gradient in potential equivalent temperature during daytime [14].

2.2. Significance Test

Diurnal change of PBL height and latent heat flux is an autoregressive stochastic periodic process in a statistical sense and depends mainly on incoming radiation. In order to separate the effect of soil parameters in the diurnal courses, the natural diurnal course has been alienated from the simulated quantities with Fourier-series analysis. A more detailed description can be found in [15]. The significances were tested to 𝑃 < 0 . 0 1 probability.

3. Soil Texture

The analysis required as precise spatial distribution of soil texture as possible. The soil texture types are defined as the FAO 12 types. In the surrounding countries the available FAO-STATSGO distribution was used, which has a 5′ resolution ( 10 km). In Hungary a more detailed Digital Kreybig Soil Information System (DKSIS) [16] was applied. The database consists of soil information on a 1 : 25000 scale. Overall in Hungary, a field campaign was enacted from 1935 to 1951 to create a comprehensive and detailed soil information database. Digital processing of these maps begun in 1998, and it ended only recently. The creation of the database included a compilation of the maps [17] and soil information, and updating them with recent remotely sensed images and spatial database. The DKSIS was also compiled together with the Hungarian Agrogeological Database (HAD) which consists information of soil particle size distribution in soil samples [18]. From the available particle size distribution, the FAO classification was determined and interpolated [19] to the 30′′ soil texture grid used by the WRF model. As it can be seen in Figure 1, the distribution of soil textures is much more detailed in Hungary. Originally there were only 6 types of soil textures over the simulation area, but with the new soil information all 12 types can be found.

4. Soil Parameters

In this study, two Hungarian soil databases and the parameters derived from them are compared. The first one is the HUNSODA [20, 21]. It contains information on 476 soil samples. Each sample refers to one soil horizon. The samples were taken only from plains. Measured soil moisture is available for soil moisture retention curves in the whole suction range (pF 0, 0.4, 1, 1.5, 2, 2.3, 2.7, 3.4, 4.2, and 6.1). The other soil database is the MARTHA [22]. The number of soil samples is about 15 times larger ( 7500 soil samples) than in the former case. The spatial extent of the collected information is about 5 times larger ( 2300 sample sites) with respect to the former database. The samples originate from the whole country. Information on soil moisture retention curves is available at least for standard pF values: 0, 2.5, 4.2, and 6.2 for every sample. Soil parameters were determined using Campbell’s [23] method, where the field capacity is defined at the pF = 2.5 value. The greatest differences in soil hydraulic parameters were found for sand, loamy sand, clay loam, and silty clay (Table 1). For loamy sand the field capacity differs with almost 200 mm/m, and the saturated soil moisture with 170 mm/m. The least differences were found for silt loam, sandy clay loam, and loam.

5. Measurements

For this study, vertical profiling measurements were conducted at the observatory of Szeged (46.25572N, 20.09023E) by the Hungarian Meteorological Service. The site is outside the city surrounded by agricultural cultivations (corn, barley). During the measurement period (August 1–31, 2011), the sum of precipitation was only 1.9 mm; this is about 1/30 of monthly average. Two instruments were used to verify the simulations, a radiometer and a windprofiler. The Radiometrics MP-3000A is a ground-based microwave radiometer which can continuously measure temperature, humidity, liquid water profiles to 10 km height. The instrument has 21 calibrated channels in 22–30 GHz (K-band) and 14 in 51–59 GHz (V-band); it is also equipped with sensors for surface temperature, relative humidity, and pressure. Measurements are available from 0–500 m height at every 50 meters, from 500–2000 m height at every 100 meters, and from 2000–10000 m at every 250 meters. At the site a Vaisala LAP-3000 is also used which is a lower atmospheric wind profiler with radio acoustic sounding system (RASS). Measurements of vertical profiles of horizontal wind speed and direction are available up to an altitude of 4 km. The instrument’s operating frequency is 1290 MHz; measurements are available at every 15 minutes. The vertical resolution of the measurements is 220 m starting at around 150 m.

6. Model Setup

Simulations were carried out with the WRF-ARW v3.1 model [24]. The number of vertical levels was 44, the mother domain’s horizontal resolution was set to 3 km, while the nest’s was 1 km. The respective domain size was 2 3 5 × 2 1 1 grid points covering the Carpathian-basin and 1 6 0 × 1 6 0 grid points centered by a synoptic station around the south border of Hungary, on the bank of the Tisza River. Simulation step size was set to 15 and 5 seconds, and the 6th order diffusion integration option was chosen. From 00 UTC the model run was 18 hour long and no measurements were assimilated. Both domains were run with the same physical parameterizations: RRTM [25] for radiation, Bougeault-Lacarrere [26] for PBL, Thompson [27] for microphysics and Noah [28] for land-surface processes. From the point of view of the analysis it has to be noted that the PBL scheme predicts the turbulent kinetic energy (TKE) to estimate the PBL height. Also in the Noah scheme the latent heat flux mainly depends on the difference of soil moisture and wilting point, and the difference in field capacity and wilting point. Initial and boundary conditions were taken from ECMWF MARS (European Centre for Medium-Range Weather Forecasts Meteorological Archive and Retrieval System) database, with 0.125° horizontal resolution and updated hourly.

7. Case Studies

From the measurement period, two days had been chosen for simulations. On 16 August 2011, a cold front passed the region a day before and a high pressure system was to be found to the west. High and mid-level clouds remained until late night. The minimum temperature was 16–19°C, while the maximum was only 25–30°C. On 26 August 2011, the synoptic conditions showed a high pressure system to the northeast and an approaching cold front from west which resulted strong southerly winds on lower levels. The day was cloud free, with 2 m temperatures between 17–22°C in the morning, reaching 32–37°C during afternoon.

8. Results

The difference between field capacity and wilting point is on average 30 mm/m greater for HUNSODA derived soil parameters compared to MARTHA ones. Since the soil moisture initial condition is the same, on average it results in lower latent heat flux and higher surface temperatures, which should enhance PBL growth. The result of this effect is first looked closely on comparing measurements with the model simulations and then the spatial distribution of affected area is discussed.

8.1. Model Validation

The daytime course of PBL height estimations can be seen in Figure 2. The PBL height estimated from SNR shows a gradual building-up mixing layer during the day, while PBL height inferred from the radiometer changes quickly. In the latter case the inversion diminishes at around 4 : 00 UTC in both cases. Afterwards the potential temperature profile still shows the typical form of a stable atmosphere [29], but the heat from the incoming solar radiation enhances turbulent flow in the lower few hundred meters of the atmosphere, slowly making it unstable. This results in a decreasing estimated PBL height, as the potential temperature profile changes. After the mixing layer became unstable, the critical gradient in potential temperature was reached quickly. In the first case study the large scale stable environment slowed the change from stable to unstable stratification. At sunset the PBL height change was sudden, because the sign of potential temperature gradient in the lowest 100 m changed quickly. During night the windprofiler usually shows a higher PBL height, than the radiometric measurement. This mainly relates to the fact that the windprofiler only has 3 measurement levels in the lowest 750 m, while the radiometer has 12. The difference in the vertical resolution suggests that the windprofiler is not able to show the nighttime stable PBL layer correctly. However after reviewing a month of windprofiler measurements, it was found that during nighttime the maximum of SNR can reach up to 1000 m without any atmospheric trigger. This means however that during night the maximum of SNR indicates the top of the residual layer rather than the top of the stable PBL. PBL heights calculated by the model are the average of the 16 grid points in the nest which are nearest to the synoptic station. Simulations show good correlation with PBL height derived from both windprofiler and radiometer measurements. In the first case study, the maximum of the model simulated PBL height is closer to the one estimated with radiometer. The passing of the cold front resulted a synoptic scale downdraft which has a reducing effect on turbulent kinetic energy (TKE). Though the difference in soil parameters caused a difference in latent heat flux in average of about 25 W/m2 and as a result a 0.5°C temperature difference occurred at 2 m, the PBL height difference was small. It has to be noted that the northwestern wind was only between 1–3 m/s through the whole day. On the second day, the effect of soil parameter difference is more obvious, since clouds were not present. Furthermore, as opposed to the previous case, strong (6-7 m/s) southeastern flow characterized the lower 2 km atmosphere. The soil texture is clay loam southeast to Szeged (1) over a vast area in the model domain, for which the wilting point and field capacity soil moisture differs with −41 mm/m and −86 mm/m, respectively. Over this type of soil texture the latent heat flux for HUNSODA parameters was reduced with about 40 W/m2 compared to MARTHA simulations. As a result, at the measurement site the temperature was about 1°C higher during most of the day in case of HUNSODA. Along with the added heat from southeasterly flow, the TKE at 1000 m height was doubled compared to MARTHA simulations causing higher PBL height. The difference between model runs between 10 and 14 UTC is around 250 m (1 model level), and late afternoon is around 500 m.

8.2. Significance Test

With the significance test, first the daytime (from 5 UTC to 17 UTC) course of latent heat flux has been analyzed (Figure 3.). On 16.08.2011, significant differences were found on the 54% of the mother domain. Usually over cites and lakes, no diurnal course was observed; in these cases the method cannot be applied, which is denoted as errors. The nest remarkably showed that in case of soil textures for which the parameters did not differ too much, the diurnal course did not change as well. For instance, the north-south directed patch of sandy loam along the Tisza River in Serbia is such a territory of nonsignificant differences. On the other day, 9% more significant grid points occurred with respect to the former case. The territory of significant differences increased in Croatia and Slovenia and also on the Plain of Hungary. The sensible heat flux daytime course (not shown) was significant over 75% and 89% of the domain, respectively.

The spatial distribution of significant differences in PBL height (Figure 4) shows little correspondence to latent heat flux change. On the first day, 63% of the grid point PBL heights showed notable difference, while this was 71% on the second day. The nest showed microstructure of the PBL height changes, parallel to the average wind direction. Daytime PBL height differences varied from 10–500 m, approximately 150 m on average.

9. Conclusions

Comparing high resolution nested WRF simulated PBL heights to measurements, the effect of soil parameters shows up in PBL height 5 hours after dawn. In some cases, the differences in PBL height can reach up to 500 m. The sudden collapse of PBL around sunset corresponds to measurements obtained by both the windprofiler and radiometer. According to the significance test applied to diurnal courses of surface heat fluxes and PBL height, the spatial distribution of significant heat flux differences usually does not correspond to significant PBL changes prompted by soil parameter value differences. It has to be also noted that the latent flux daytime courses show less response to soil parameter differences than PBL height, primary due to sensible heat flux changes.

Acknowledgments

The project is supported by the European Union and co-financed by the European Social Fund (Grant Agreement no. TÁMOP 4.2.1/B-09/1/KMR-2010-0003 and TÁMOP 4.2.2/B-10/1-2010-0030) and by the Hungarian Scientific Research Fund (OTKA K-81432). The authors also thank the help of István Aszalos with the radiometric data.