Effects of land surface inhomogeneity on convection-permitting WRF simulations over central Europe

This sensitivity study investigates the impact of the spatial scales of land use, soil moisture and orography patterns on land–atmosphere exchange fluxes, domain-wide averages and meteorologic conditions in convection-permitting RCM simulations. We perform five WRF RCM simulations, each with the same 3 km resolution atmospheric setup but different combinations of coarsely resolved (12 km) land use and soil type, initial soil moisture and orography for the heat-wave summer 2003 over central Europe. Our results indicate that a coarser-resolved orography significantly alters the flow over and around mountain ridges such as the Alps and impact the large-scale flow pattern. The smoothed mountain ridges result in weaker Föhn effects and in enhanced locally generated convective precipitation patterns, peaking earlier in the afternoon. In comparison, the impact of a coarser-resolved land use is smaller and mainly related to changes in the overall spatial fraction of a land use, rather than to the loss of heterogeneity of the different land use types on the scale analyzed here. Albeit, even small changes in the initial soil moisture (both spatial averages and local differences) have a higher potential to affect the overall simulation results, although this might also depend on the land surface model. Overall, effects induced by the coarsely resolved land surface properties are small compared to the differences between simulations with 3 km and 12 km grid spacings of the atmosphere.


Introduction
In numerical weather prediction and climate modeling, an increased model resolution generally improves the quality of simulations. While mean features on large spatio-temporal scales may be equally well-represented in coarser model resolution simulations, the added value of the higher resolution mostly emerges in the representation of mesoscale phenomena, regional-scale climate variability and daily extremes, especially for regions with complex orography, land-sea contrast and heterogeneous land use (Di Luca et al. 2012Kotlarski et al. 2014;Torma et al. 2015;Prein et al. 2016). Largest improvements occur when processes are explicitly resolved by the finer grid spacing. Especially the step from convection-parameterized to convection-permitting resolution (grid spacings below about 4 km) clearly improves the diurnal cycle, spatial distribution and intensities of precipitation (Kendon et al. 2012;Ban et al. 2014;Prein et al. 2015).
Higher resolutions better resolve the mesoscale processes in the atmosphere and more realistically represent the exchange of water and energy with the atmosphere over heterogeneous land surface including feedbacks to atmospheric processes that arise from an altered land-atmosphere coupling. Accordingly, the largest improvements of convection-permitting simulations compared to coarser resolution convection-parameterized simulations have been found over 1 3 regions with heterogeneous land use and complex orography (Prein et al. 2015).
The impact of a more highly resolved land surface can be separated into effects of land use heterogeneity and of orography. Land use types may differ largely in albedo, roughness length, leaf area index, etc., which impact the surface radiation and energy balances and the coupling with the atmosphere via the turbulent transport of momentum, heat and water vapor. Dependent on the scale and contrast of its heterogeneity, land use patterns for instance may induce sensible heat flux patterns, that generate circulations ranging from the turbulent scale, affecting the state of the atmospheric boundary layer, up to organized mesoscale circulations similar to the land-sea breeze effect (Giorgi and Avissar 1997;Avissar and Schmidt 1998).
Many studies quantified the impact of surface heterogeneity on shallow convection using idealized large eddy simulations (LES) (Raasch and Harbusch 2001;van Heerwaarden and Guerau de Arellano 2008;Brunsell et al. 2011;Garcia-Carreras et al. 2011;Wu et al. 2015). Roll-and cell-like structures create local areas of convergence associated with strong vertical upward motion, which tends to be connected to land surface discontinuities with strengths depending on the amplitude of the surface heterogeneity. Heterogeneityinduced circulations may lead to areas of higher relative humidity over warm patches with a higher potential for cloud formation (van Heerwaarden and Guerau de Arellano 2008). Huang and Margulis (2013) showed that strongly heterogeneous surfaces yield thicker and deeper shallow convective clouds. Rieck et al. (2014) report that the transition from shallow to deep convection is strongest over land surface heterogeneities of 12.8 km patch size.
Besides the static land surface characteristics, the dynamic soil moisture is a key variable for land-atmosphere interactions as it largely controls the flux partitioning into sensible and latent heat . The non-linear impact of soil moisture on the flux partitioning is strongest in transition zones between wet and dry regimes. In a wet, energy-limited regime evapotranspiration is primarily limited by the state of the atmospheric boundary layer (ABL). In a dry, moisture-limited regime, the atmospheric moisture deficit cannot be entirely compensated by moisture supply from the land surface via evapotranspiration, i.e., the latent heat flux is constrained by the soil moisture content. Besides the absolute soil moisture value, also its heterogeneity influences the ABL (Patton et al. 2004;Huang and Margulis 2013) and convective precipitation (Hohenegger et al. 2009;Adler et al. 2011;Guillod et al. 2013;Taylor et al. 2013;Froidevaux et al. 2014). Thereby the influence of both land use and soil moisture related heterogeneity depends on the background wind conditions. Avissar and Schmidt (1998) found that already moderate winds reduce almost all impacts of surface heterogeneity, while Froidevaux et al. (2014) suggest positive or negative local soil moisture-precipitation feedbacks dependent on the background wind. For real case simulations with non-uniform surface heterogeneity patterns and various weather conditions, the complex feedback processes complicate the identification of systematic patterns in the atmosphere that can clearly be related to land surface patterns.
Orography strongly impacts atmospheric flow (Smith 1979); e.g., mountains force the airflow to ascend, which may initiate the formation of clouds and precipitation via upslope condensation and upstream or downstream triggering of convection. Thus, orography impacts precipitation intensity and its spatial distribution (e.g., Barros and Lettenmaier 1994), but also changes the hydrological cycle at the surface via mountain slope runoff. Many studies that have investigated the sensitivity of simulated orographic Fig. 1 Central European WRF RCM model domain (3 kmME, orange frame, 0.0275°) nested into the EURO-CORDEX domain (EUR-11, black frame, 0.11º); grids are rotated; for the analyses, the lateral boundary relaxation zone of ten grid points each is omitted precipitation to model resolution, have found that orographic precipitation generally increases with increasing grid resolution until the important orographic features are sufficiently resolved (e.g., Leung et al. 2003;Colle et al. 2005;Pedemonte et al. 2005;Richard et al. 2007;Prein et al. 2015). This improvement is to some degree attributed to the better resolution of the atmospheric processes. Smith et al. (2015) showed, however, that the simulated orographically enhanced stratiform precipitation over the UK in winter is reduced by 10% and 23%, when the orography resolution is reduced from 1.5 km to 12 km and 40 km, respectively. For mainly convective precipitation over the Alps, Prein et al. (2013) identified, just small differences in the mean precipitation diurnal cycle between convection-permitting simulations with 3 km and 12 km orography. From the above studies, no clear conclusions on the quantitative effects of the different types of land surface heterogeneity on atmospheric motion in a non-idealized study could be drawn. Thus, we attempt in this case study to systematically investigate the impact of the spatial scales of the patterns of land use, soil moisture and orography on convection-permitting RCM simulations on land-atmosphere exchange fluxes and patterns in atmospheric state and flux variables as well as domain wide averages. To this goal, we analyze five 3 km grid size, convection-permitting seasonal simulations for central Europe, which were performed with different combinations of coarser (12 km) resolved land surface property types and allow to assess the impact of each land surface property on atmospheric processes, also in comparison to the driving 12 km simulation with parametrized convection. To our knowledge, such an analysis for central Europe and a full season does not yet exist.
In Sect. 2, we introduce the setup of the five WRF RCM simulations. In Sect. 3, the results are compared for several variables, first with a focus on seasonal means and similarities in the differences of spatial patterns, and second for their temporal evolution. Results are summarized and discussed in Sect. 4.

RCM simulation experimental setup
This study is based on convection-permitting regional climate simulations using the non-hydrostatic Weather Research and Forecasting model (WRF, v.3.6.1) (Skamarock et al. 2008) in a one-way double-nesting setup (see Fig. 1). The central European domain with a 0.0275° (~ 3 km) grid spacing (3kmME, 480 × 456 grid points) is nested into the pan-European 0.11° (~ 12 km) EURO-CORDEX domain (EUR-11, 448x436 grid points) with 50 vertical levels up to 20 hPa. For both domains the same model settings and parametrizations are used, except for the Grell-Freitas convection scheme (Grell and Dévényi 2002) used in the EUR-11 domain and switched off in the 3kmME domain. Both use WSM-5 microphysics (Hong et al. 2004), the RRTMG radiation scheme (Iacono et al. 2008), the YSU PBL scheme (Hong et al. 2006) and the NOAH land surface model (LSM) (Ek et al. 2003).
Five simulations at 3 km grid size have been carried out each with the same parametrization setup and identical initial conditions of the atmosphere. The 3 km simulations only differ from each other in the respective resolution (either 12 km or 3 km) of the land surface characteristics (a) land use and soil type, (b) initial soil moisture, and (c) orography. See Table 1 for an overview of how the surface characteristics have been changed and Fig. 2 for the differences in the spatial distribution of the static fields.
Since the soil type classification base data have a coarser than 3 km resolution, no additional simulations with 12 km soil type and, e.g., 3 km land use have been carried out. Soil type is varied together with the land use.
The setup of the driving EUR-11 simulation is identical for all five simulations, and since there is no feedback from the inner 3kmME to the EUR-11 domain (1-way nesting), identical lateral boundary conditions are provided for the 3kmME domain for each of the five simulations. Thus,  -level soil moisture (JJA mean, first row) and its temporal evolution (15-daily means, second-last row) in the REF simulation (first column) and in the setups A-D (second-fifth column) displayed as difference to REF. Domain averages and differences are shown in the upper right corner ◂ differences in the individual runs only result from the differently resolved land surface properties. Lateral boundary forcing data for the EUR-11 domain simulations are provided by 6-hourly ERA-Interim reanalysis data at 0.75° grid resolution.
The experiment encompasses the summer months June, July, August (JJA) in the year 2003 as this season was characterized by strong land-atmosphere coupling conditions (Fischer et al. 2007). It had mainly convective precipitation events in June and July as well as longer dry and hot periods in July peaking in a severe heat wave in the first half of August.
The land use data are based on IGBP-MODIS data for all simulations but interpolated either to 3 km or 12 km resolution. As all setups use the same 3 km model grid, the coarse 12 km pattern of the land surface properties is generated by nearest neighbor resampling from the 12 km grid of the parent EUR-11 simulation. Thus, for setups A, B and C the land use and soil types of 4 × 4 3 km-wide grid boxes are identical (see Fig. 2d). Accordingly, the percentages of individual land use types slightly change between the 3 km and 12 km resolutions (Fig. 3). The LSM used in our simulations does not use a mosaic approach but the dominant land use type of one grid box; thus subgrid-scale heterogeneity is ignored. The coarsely resolved orography setup for A and D is generated by a weighted bilinear interpolation from an orography resolved at 12 km; to avoid very steep slopes between the neighboring grid cells, which would generate unrealistic wind patterns.
All simulations have been initialized with interpolated atmospheric state variables from ERA-Interim at 1 st of June 2003. To reduce the spin-up time for the soil variables (soil temperature and moisture) their initial values are taken from previous long-term WRF climate simulations started in 2001 with the same EUR-11/3kmME nesting setup; for an in depth information on these simulations, the reader if referred to Knist et al. (2018). For setups A and B the initial soil variables are taken from the 12 km resolved EUR-11 climate simulation and for setups C, D and REF from the 3kmME climate simulations. The mean soil moisture values are slightly different between the EUR-11 and 3kmME climate simulation used for the initialization, mainly because of precipitation differences; overall the 3kmME climate simulations are slightly drier than the EUR-11 climate simulations. For that reason, an additional setup C that uses 12 km land use and soil type (like setup B) is initialized with 3 km soil moisture from the former 3kmME climate simulation from Knist et al. (2018).

Surface properties impact on the evolution of atmospheric states
We first analyze the mean differences of variables between the reference run REF with all land surface characteristics resolved at 3 km, and the results of setups with all (A) or part (B-D) of the land surface properties resolved at 12 km. As there are at least two setups that have one land surface property in coarse resolution (12 km) in common, similarities in the difference patterns of these setups with reference to the REF run allows to separate the impacts of the individual land surface characteristics. The REF run is used as the reference in this sensitivity study for the other setups, which have at least one land surface property at a coarse resolution; thus it is not to be seen as the truth which could only be provided by observations or-to some degree-by assimilating the observations into the simulations. We investigate also the temporal evolution of the difference patterns between the individual setups for several variables using 15-day averages. That analysis is based on the hypothesis, that some difference patterns may depend on specific weather conditions and the large-scale synoptic conditions during JJA 2003, and that similarities in the difference patterns may remain consistent throughout the simulated time period despite different synoptic conditions.

Latent and sensible heat fluxes
Summer 2003 in central Europe was anomalously warm and dry and culminated in a severe heat wave beginning of August (e.g., Fischer et al. 2007). The domain spatial mean soil moisture impacting on flux partitioning between the latent and sensible heat flux continuously decreases in the reference simulation from 0.23 m 3 /m 3 beginning of June to 0.13 m 3 /m 3 end of August (Fig. 4, leftmost column).
While the latent heat flux dominates the turbulent fluxes in June, the system undergoes a transition from energylimited to soil moisture-limited conditions accompanied with a strong decrease of latent heat flux and an increase in sensible heat flux until the end of August (Figs. 5 and 6, leftmost column). The spatial patterns of latent and sensible heat flux are very heterogeneous due to their dependence on soil moisture, land use, and soil type. Urban areas are, e.g., clearly identified by high sensible heat fluxes and very low latent heat fluxes. High latent heat fluxes and low sensible heat fluxes are also simulated for the forested areas in the Different temporal turbulent flux evolutions between the reference and the different land surface characteristics setups occur instantly for coarse land use simulations (setups A, B and C) at grid points that changed their land use type with respect to the original 3 km reference run ( Fig. 5 and 6, second row). Besides the noisy pattern with alternating positive and negative differences on the local scale, the setups A and B simulate for large areas higher mean latent and lower mean sensible heat fluxes. These differences relate to the slightly wetter initial soil moisture conditions of the 12 km resolution (Fig. 4), which lead to 5.4 W/m 2 higher (4.1 W/ m 2 lower) domain average latent (sensible) heat fluxes. Interestingly, the excess in latent (and deficit in sensible heat) fluxes in setups A and B peak in mid July and not in June when the soil moisture difference is largest. This results from the non-linear relationship between soil moisture and latent heat flux: the transition from the energy-limited state (when soil moisture does not impact evapotranspiration) to the soil moisture-limited state (when evapotranspiration decreases with decreasing soil moisture) is reached earlier in the initially drier REF simulation than in the initially moister A and B setups. Until the end of August differences decrease as the latent heat flux is low in all simulations.
For setup C with an initial soil moisture identical to REF, strong differences mainly occur just for grid points with disparate land use types caused by the coarser resolution of the land use type maps. The domain-averaged latent heat flux for setup C is consistently lower than for setups A and B, but 3-4 W/m 2 higher compared to the REF run because of the slightly different percentages of land use types (see Figs. 3 and 5).
The individual land use type contributions to the mean latent heat flux difference between REF and setup C are illustrated in Fig. 7, which shows time series of daily latent heat flux differences averaged over groups of grid points that undergo land use changes from REF to setup C. Accordingly, the initially higher domain average latent heat flux in setup C mainly results from the increase of the number of grid points with land use type crop in setup C that are urban (low evapotranspiration) in REF. A similar effect holds for land use changes from mixed forest in REF to crop in setup C, and vice versa: in June crop grid points typically evapotranspirate more than mixed forest land use grid points, which reverses with the drying of the soil conditions in July and August. The slightly different fluxes between setup C and REF even for the same land use types at the beginning of the simulation relate to differences in mean vegetation fraction between the 3 km and 12 km resolution setups.
For setup D (like REF but coarse orography) initially only small differences compared to REF occur (Figs. 5 and  6). The evolving noisy latent and sensible heat flux difference patterns are mainly related to different soil moisture evolutions as a result of the differently evolving precipitation (see below). Overall, the differences in turbulent fluxes of all setups compared to REF are primarily related to the flux partitioning at the surface and less to differences in the incoming solar radiation, and hence cloudiness, which is discussed in the next section.

Shortwave radiation
Incoming shortwave radiation (SW_d) controls the variability of surface net radiation, which is the main driver of the turbulent fluxes. Accordingly, SW_d is inversely related to daytime cloudiness. Since there is no direct dependency of SW_d on land surface properties, differences in SW_d between the setups (Fig. 8) are smaller compared to the turbulent fluxes (Figs. 5 and 6). Due to the complex and non-linear land-atmosphere feedback processes, systematic difference patterns are difficult to identify not only due to the random-like internal variability evolving during the time integrations of the RCM. Nevertheless, some similarities are noticeable for the different surface characteristics setups.
The largest similarity in the difference patterns-also in terms of spatial correlation (Table 2)-is observed for the setups with the same coarse orography (A and D). Similar to the temporal evolution of the area-mean values, less shortwave radiation reaches the surface over the Alpine ridge because of the consistently higher convective activity and hence cloudiness (data not shown). Other regional patterns only sporadically occur, indicating most likely a dependency on weather conditions. Setups B and C, which differ from REF by their coarse land use type resolution, show only scattered differences in incoming solar radiation, which mostly cancel each other in the mean. Setup B, which additionally differs from REF by an initially slightly wetter soil, generates somewhat less incoming overall shortwave radiation, due to more cloudiness (data not shown). We see very similar differences between both setups with coarse orography but different initial soil moisture, A and D. These differences most likely result from the overall higher evapotranspiration and thus more clouds in the moister setups A (and B). This causality is corroborated by the anti-correlated temporal evolution of the difference in average shortwave radiation and latent heat flux except for setup D, which suggests that in that case, the differences are dominated by different flow pattern due to the coarse orography. The difference patterns for setups B and C show, despite the noisy structure, some similarities on the regional scale (average daily correlation 0.49), that suggest a systematic impact of the coarsely resolved land use. For setups that have only one surface property in common, the highest spatial correlation exists between setup A and D (same 12 km orography) while the lowest correlation between setups that differ by only one surface property is found for setup A and B, i.e., different orography (Table 2).

Precipitation
Precipitation in the JJA 2003 simulations is characterized by high amounts in the central and eastern Alps and in the central German uplands (Fig. 9, left column). The noisy pattern with large variability on the local scale results from single convective events, whose tracks also cause lines of high precipitation amounts as seen, for example, in northern Germany-oriented WSW-ENE.
Since the location of single convective cells and its evolution is mostly random also the difference pattern of the individual setup runs with respect to the REF simulation shows large variability on the local scale. However, there are some similarities in the difference patterns on the regional scale and between the mean values over the whole domain. Consistent with the results for incoming shortwave radiation, we observe more precipitation in the setups with less shortwave radiation (more clouds). For example, setups A and D with the coarse orography are consistently wetter over the Alpine ridge where up to 30% more precipitation is simulated. Difference patterns in other regions are more heterogeneous and vary with time, while some similarities exist between the setups with coarsely resolved land surface properties. Similar to shortwave radiation, setups A and D are correlated highest (0.61) followed by setups B and C (0.48) ( Table 2). Again, positive and negative differences mostly cancel each  other and lead to overall precipitation differences < 10% with the higher initial soil moistures leading to higher precipitation for all 15-day-averaged domain averages. Setups A and B (coarse land use and higher initial soil moisture) lead to about 7% more mean precipitation compared to their counterparts D and REF, respectively, with the same orographic resolution. Setups C and REF, which only differ in the resolution of land use and soil type, are very similar with 3% higher mean precipitation in REF (higher resolved land use distribution compared to setup D).

CAPE and 850 hPa level patterns
Summer precipitation is mostly convective; thus, the precipitation response should be related to changes in the convective available potential energy (CAPE, Fig. 10a-e).
High CAPE not necessarily results in stronger precipitation on average, but indicates the potential for strong convective events. The runs with coarse land use (setups A, B and C) show slightly higher overall CAPE, while the runs with coarse orography (setups A and D) simulate lower CAPE compared to their respective counterparts (setup B and REF). Moreover, the coarse orography runs (setups A and D) result in consistently lower CAPE south of the Alps, which we relate to different flow patterns (see below) over and around the Alpine ridge with a smoother orography. The coarsely resolved orography (setups A and D) influences the flow pattern in 850 hPa over the whole period much stronger than in the highly resolved setups B and C ( Fig. 10f-j). The strongest differences occur south of the Alps. On any occurrence of north-westerly flows towards the Alps (which is the case for most of the simulated period, Fig. 10f), lee cyclogenesis results in lower surface pressure and the 850 hPa geopotential height, which is, however, weaker in the coarse orography setups A and D (Fig. 10g,  j). A possible reason is the overall lower height of the main Alpine ridge in the smoothed 12 km orography, which leads to a weaker Föhn effect resulting in slightly lower temperatures at the downwind side and less precipitation at the windward side of the mountains (Fig. 10 l and o). The rest of the domain does not exhibit significant differences in geopotential and flow pattern between the individual setups. Runs with coarsely resolved land use but highly resolved orography (setups B and C) do not differ much in mean 850 hPa geopotential height compared to REF, thus orography dominates the observed changes in the flow pattern.
In general, all strong local temperature differences between setups with coarse and fine resolved land use vanish above the boundary layer and no significant effect of land use heterogeneity is seen at 850 hPa ( Fig. 10k-o). The overall slightly colder 850 hPa mean temperature (0.1-0.5 K) in setups A and B are likely related to the overall higher latent heat flux and more evaporative cooling due to slightly wetter initial soil in those setups (see Fig. 5 and 6). This is consistent with difference patterns of 850 hPa specific humidity that is slightly higher in setups A and B compared to REF (Fig. 10p-t). Again, the strong heterogeneity at the surface completely vanishes at the 850 hPa level.

Mean diurnal cycle of orographic precipitation
Besides the larger scale difference pattern caused by the altered flow pattern, the coarsely resolved orography also leads to overall more precipitation and clouds especially over the Alps. The mostly convective summer precipitation follows a pronounced mean diurnal cycle with maximum precipitation in the afternoon. While nighttime precipitation amounts are quite similar, the coarse orography runs simulate higher (+ 40%) and earlier (− 1 h) afternoon maxima compared to REF (Fig. 11). The mean hourly precipitation patterns suggest that the coarse orography facilitates the initiation and evolution of more organized, longer lasting, and larger convective cells. The higher convective activity is consistent with the higher mean CAPE values reached in these runs (see Fig. 10a-e).

Discussion and conclusion
Five WRF simulations each with the same atmospheric setup at 3 km resolution but different combinations of coarsely resolved (12 km) land use and soil type, initial soil moisture, and orography have been performed for JJA 2003. Our results indicate that a coarser-resolved orography significantly alters the flow over and around extensive mountain ridges such as the Alps and hence impacts the large-scale flow pattern. This could be a so far undiscovered added value of convection-permitting simulations on the lee side of the western Alps, that could improve the simulated climate in this region, especially the Po valley, which is known to be error-prone in climate simulations (Nikulin et al. 2011;Kotlarski et al. 2014;Prein et al. 2016). The coarser resolution leads to weaker Föhn effects and might also influence lee cyclogenesis, which could be investigated in future studies with a larger domain that would also have to cover larger parts of the northern Mediterranean Sea.
The smoothed orography furthermore results in enhanced locally generated convective precipitation patterns peaking earlier in the afternoon. This is at variance with other studies, which suggest lower precipitation over smoothed orography (e.g. Smith et al. 2015). These seemingly conflicting results can likely be attributed to the character of precipitation, which was mostly convective in our summer 2003 study and predominantly stratiform in the winter cases analyzed in Smith et al. (2015). Also, Prein et al. (2013) identified only small differences in simulated convective summer precipitation over the eastern Alpine region between 3 km and 12 km orography runs. They concluded that the typically positive precipitation bias in convection-parameterized simulations mainly results from the parameterization scheme and the less resolved atmospheric dynamics rather than from the coarsely resolved orography. Our study does not confirm this conclusion. The mean diurnal cycle of precipitation over the Alpine region in the setups REF and D and in the driving EUR-11 simulation (Fig. 11) suggests that the much higher precipitation amount in EUR-11 (+80%) is at least partly attributed to the coarsely resolved orography that leads to an increase of 40% in setup D compared to the REF simulation. The steeper slopes and higher differences between valleys and hill tops in the 3 km orography may initialize stronger vertical wind speeds and lift air more easily to the levels of condensation and free convection, which would suggest overall less convective precipitation with coarser orography, contrary to our results where the smoothed orography seems to facilitate the initiation and evolution of more organized and overall larger convective cells that last longer. This reasoning is also consistent with the higher mean CAPE values observed for the simulations with smooth orography. However, differences to previous studies may also result from different model setups and the simulated season.
Besides local inevitable effects on near-surface variables caused by differing land use types, the reduction of heterogeneity in coarsely resolved land use distributions, going from 3 km to 12 km resolution, does not systematically affect the atmospheric flow, temperature or precipitation patterns. However, domain average differences in latent heat flux (about +3% between setup C and REF) and consequently differences in mean cloudiness and precipitation (about +3%) may result from the slightly changing the percentages of land use types upon aggregation. No clear differences in the spatial distribution of hourly precipitation intensities and the mean diurnal cycle are, however, observed between the different setups (not shown). Positive spatial correlations between precipitation and shortwave radiation difference patterns of setups with coarse land use also suggest a nonrandom effect on cloud formation. However, additional perturbed ensemble simulations with each setup and an extension to several summer seasons would be necessary to better assess systematic influences on these difference patterns.
The initial soil moisture distribution more strongly alters the regional-scale atmospheric patterns compared to the scale of the land use type heterogeneity. In our setup, the overall higher (+2.3% domain average) initial soil moisture taken from transient EUR-11 climate runs (setup A and B) compared to the 3kmME climate runs (setups REF, C, and D) leads to systematic differences in domain-averaged latent and sensible heat flux and consequently also in cloudiness and precipitation in the order of 6-8% domain average. The impact of soil moisture differences clearly depends on the preceding weather conditions and is largest for transitions from energy to soil moisture-limited land-atmosphere coupling regimes, that widely occurred during the simulated summer season. In line with previous studies, our results highlight the sensitivity of land-atmosphere coupling to soil moisture dynamics in particular for this anomalous dry summer 2003 (Fischer et al. 2007;Seneviratne et al. 2010;Keune et al. 2016;Knist et al. 2017) and confirms a contrasting response of different land use types (e.g., crop and forest) to heat wave conditions (Teuling et al. 2010).
In summary, when decreasing the spatial resolution of individual land surface properties, orography has the largest impact on the overall simulation results as it alters the flow over and around mountain ridges and largely influences the local precipitation pattern and intensity over regions with complex orography. The effect of a coarser-resolved land use distribution is overall smaller and mainly related to changes in the overall percentages of the different land use types, rather than to the loss of heterogeneity in the surface pattern on the scale analyzed here. But even small changes in soil moisture (both spatial averages and local differences) have a higher potential to affect the overall simulation results, although this might also depend on the particular land surface model used, the treatment of subsurface processes and the investigated season. Overall, differences caused by coarsely resolved land surface properties are much smaller compared to differences between simulations with 3 km and 12 km grid spacings of the atmosphere (e.g., Ban et al. 2014;Kendon et al. 2012;Knist et al. 2018;Prein et al. 2013).