Breeze effects at a large artificial lake : summer case study

Natural lakes and big artificial reservoirs can affect the weather regime of surrounding areas but, usually, consideration of all aspects of this impact and their quantification is a difficult task. The Alqueva reservoir, the largest artificial lake in western Europe, located on the south-east of Portugal, was filled in 2004. It is a large natural laboratory that allows the study of changes in surface and in landscape and how they affect the weather in the region. This paper is focused on a 3-day case study, 22–24 July 2014, during which an intensive observation campaign was carried out. In order to quantify the breeze effects induced by the Alqueva reservoir, two simulations with the mesoscale atmospheric model Meso-NH coupled to the FLake freshwater lake model has been performed. The difference between the two simulations lies in the presence or absence of the reservoir on the model surface. Comparing the two simulation datasets, with and without the reservoir, net results of the lake impact were obtained. Magnitude of the impact on air temperature, relative humidity, and other atmospheric variables are shown. The clear effect of a lake breeze (5–7 m s−1) can be observed during daytime on distances up to 6 km away from the shores and up to 300 m above the surface. The lake breeze system starts to form at 09:00 UTC and dissipates at 18:00– 19:00 UTC with the arrival of a larger-scale Atlantic breeze. The descending branch of the lake breeze circulation brings dry air from higher atmospheric layers (2–2.5 km) and redistributes it over the lake. It is also shown that despite its significant intensity the effect is limited to a couple of kilometres away from the lake borders.


Introduction
Human activities, such as urbanization, deforestation, or water reservoir building, change the properties of the surface (vegetation cover, emissivity, albedo) that rule the surface energy fluxes (Cotton and Pielke, 2007).As a consequence, changes in surface energy fluxes affect local weather and climate.Lakes and reservoirs contain about 0.35 % of global freshwater storage (Hartmann, 1994) and cover only 2 % of the continental surface area (Segal et al., 1997).Thermal circulations triggered by lake-land thermal contrast have an impact on dispersion of air pollution and lake catchment transport (Lee et al., 2014).Big lakes, being a significant source of atmospheric moisture, can intensify storm formation (Samuelsson et al., 2006;Lin et al., 2012).Lakes and reservoirs, compared to land surfaces, have higher thermal inertia and heat capacity, and lower albedo and roughness length (Bonan, 1995).They can affect meteorological conditions and atmospheric processes at mesoscales and synoptic scales (Pielke, 1974;Bates et al., 1993;Pielke, 2013).
Normally, near-surface relative humidity is increased while daily air temperature is decreased above lake and shore areas.During the warm summer periods the relatively colder lake surface interacts with the atmosphere above, which leads to a reduction in clouds and precipitation.Formation of the local high-pressure areas over the lake surface in the summer season supports atmospheric circulation, which can be observed as a lake breeze (Bates et al., 1993).In autumn and winter it has the opposite effect: due to the fact that water is warmer than the air above, increases in evaporation and cloud formation can be observed (Ekhtiari et al., 2017).These regional lake effects have been seen in previous studies, e.g.Elqui Valley reservoir in Chile (Bischoff-Gauß et al., 2006) and the great African lakes (Thiery et al., 2014).
Published by Copernicus Publications on behalf of the European Geosciences Union.
The theoretical aspects of lake breeze formations are well known; however, this phenomenon remains unexplored.Difficulties in studies of lake breeze are due to the diversity and complexity of lake shapes and surrounding landscapes, and the lack of observational data at sufficiently fine spatial resolution (Segal et al., 1997).
Lake breezes are mainly determined by the landscape and weather conditions.Formation and intensity of the breeze depend on the set of parameters such as large-scale winds, sensible heat flux, geometry of the lake, and terrain types of the surrounding area (Segal et al., 1997;Drobinski and Dubos, 2009;Crosman and Horel, 2012).
The focus of this work is on the study of the lake breeze at the Alqueva reservoir and its impact on atmospheric parameters of the surrounding area.This large artificial reservoir was filled in 2004, which makes it a big natural laboratory for studying physical, chemical, and biological effects.Few studies about the influence of Alqueva on atmosphere and climate were published.The first report, in Portuguese, was published even before the construction of the dam by Miranda et al. (1995), as a part of the environmental impact study of the reservoir on the basis of numerical simulations performed with the NH3D (non-hydrostatic threedimensional) mesoscale model from Miranda and James (1992).It was concluded that the climate impact of the multipurpose Alqueva project should be mainly due to the irrigation of surrounding area.The influence of the reservoir itself was unclear as at that time it was not possible to perform high-resolution simulations.The studies were continued and improved by Salgado (2006) who made the first attempt to quantify the direct effect of the reservoir on the local climate, in particular on winter fog.Using the Meso-NH (nonhydrostatic mesoscale atmospheric model) model, the author concluded that the introduction of the reservoir should increase the winter fog slightly in the surrounding area, but decrease it over the filled area.Later on, Policarpo et al. (2017) used observations from two periods of 10 years (before and after the Alqueva reservoir) combined with Meso-NH simulations and showed a slight increase in the average number of days with fog during the winter (about 4 days per winter after 2003 in a downwind site).
Mesoscale atmospheric models, such as Meso-NH, allow obtainment of results with adequate horizontal resolution (250 m in present study) for studying the local effects of air temperature changes and the generation of smallscale circulations under different large-scale atmospheric situations.In this work simulations have been run for the Intensive Observation Period (IOP) of the ALEX project (ALqueva hydro-meteorological EXperiment, http://www.alex2014.cge.uevora.pt/,28 September 2018).Data collected during this experiment were used to validate the numerical simulations.
The article outline is as follows.Section 2 provides a brief description of the Alqueva reservoir.Section 3 introduces information about the ALEX experiment and the measure-ments used in this paper.Section 4 contains a brief description of the numerical models used in this work: Meso-NH and FLake.Sections 5 and 6 are dedicated to the case study on 22-24 July 2014: validation of simulation results using in situ measurements and the studies of the lake effects respectively, with an illustration and discussion of the magnitude of the impact and intensity of a lake breeze.Section 7 summarizes the results and conclusions.

Object of study
The Alqueva reservoir, established in 2002, is an artificial lake located in the south-east part of Portugal.It spreads along 83 km over the Guadiana river valley covering, when completely filled, an area of 250 km 2 with a capacity of 4.15 km 3 , which makes it the largest artificial lake in western Europe (Fig. 1a).The maximum and average depths of the reservoir are 92 and 16.6 m respectively.The dam is located in the southern part of the reservoir (Fig. 1c).
The Alqueva reservoir is mainly used to provide water supply, irrigation, and hydroelectric power.The region where it is located is known for the irregularity of its hydrological resources, with long periods of drought that could last for more than one consecutive year (Silva et al., 2014).The region has a Mediterranean climate with dry and hot summers (Csa according to the Köppen climate classification), with a small area within of the mid-latitude steppe (BSk) category.During summer, the maximum air temperature ranges between 31 and 35 • C on average (July and August), often reaching values close to 40 • C, or even higher.The incident solar radiation at the surface is one of the highest in Europe, with mean daily values (integrated over 24 h) of about 300 W m −2 and the daily maximum in July often reaching 1000 W m −2 .Rainfall periods are seasonal, normally from October to April.The annual average of accumulated precipitation in the climatological station of Beja, located 40 km away from the Alqueva reservoir, is 558 mm (1981-2010 normals from http://www.ipma.pt/pt/index.html,28 September 2018).
Two major factors determine synoptic circulations over the region during the summer period: the shape and location of the Azores anticyclone, and the frequent establishment of a low-pressure system over the Iberian Peninsula, induced by the land-ocean thermal contrasts.The later phenomenon is known as the Iberian thermal low (Hoinka and Castro, 2003), which organizes the sea breezes generated at the Iberian coasts on a peninsular scale.At the west coast, the sea breeze controls the transport of the maritime air masses from the Atlantic to Iberian Peninsula, on distances more than 100 km reaching the Alqueva region in the late afternoon.This effect is observed in the local increase in wind intensity and in its rotation (Salgado et al., 2015).

Measurement data
The measured data used in this work were obtained during the ALEX campaign -a multidisciplinary observational experiment at the Alqueva reservoir that lasted from June to October 2014.One of the aims of this project was to perform a wide set of measurements of chemical, physical, and biological parameters in the water, air columns, and over the water-atmosphere interface.To reach this goal the project operated the following facilities: seven sites with meteorological measurements: two platforms (Montante and Mourão); one permanent weather station located in a small island nearby the dam (Alquilha), two dedicated weather stations (Barbosa and Cid Almeida), two compact weather stations in Solar Park and Amieira; four floating platforms where water quality and biological sampling was carried out: Montante, Mourão, Captação, and Alcarrache; three weather stations of the Institute of Earth Sciences (ICT), located in Mitra, Portel, and at the University of Évora; two air quality mobile units: Amieira and Solar Park; three atmospheric electricity stations: Amieira, Solar Park, and Beja.
Also, data from 42 IPMA (Portuguese Institute of Sea and Atmosphere) meteorological stations located in the region were integrated into the ALEX database.1c).These locations allowed the characterization of the lake effects on a fine scale.Land stations stored data at 1 min resolution including horizontal wind speed, relative humidity, air temperature, and downwelling short-wave radiation.Montante floating platform was the principal experimental site inside the reservoir.The following equipment was installed there on 2 June 2014 and data have been collected until the end of the campaign: an eddy-covariance system that provides data for pressure, temperature, water vapour and carbon dioxide concentrations, 3-D wind components, momentum flux, sensible and latent heat fluxes (with 30 min time step), carbon dioxide flux, and evaporation; one albedometer and one pirradiometer in order to measure upwelling and downwelling shortwave and total radiative fluxes; nine thermistors to measure water temperature profile.
The intensive observation period of the ALEX project from 22 to 24 July and included the launch of meteorological balloons every 3 h.In total, 18 radiosondes were launched: 2 from the boat over the lake and 16 from the land.Atmospheric profiles of air temperature, relative humidity, wind, and pressure were obtained.This period was chosen for a case study as it was well documented with typical anticyclonic conditions: hot, dry, and with low near-surface wind speed.
Data collected during the ALEX field campaign have already been used to study lake-atmosphere interactions, including the heat and mass (H 2 O and CO 2 ) fluxes in the water-air interface (Potes et al., 2017); the effects of inland water bodies on the atmospheric electrical field (Lopes et al., 2016); and the evolution of the vertical electrical charge profiles and its relation with the boundary layer transport of moisture, momentum, and particulate matter (Nicoll et al., 2018).

Meso-NH atmospheric model
To study the lake breeze effects in the Alqueva reservoir, the Meso-NH model (Lac et al., 2018) was used.Meso-NH is a non-hydrostatic mesoscale atmospheric research model.It can simulate the evolution of the atmosphere on scales ranging from large (synoptic) to small (large eddy) and has a complete set of physical parametrizations.Meso-NH is coupled with the SURFEX (Surface Externalisée, Masson et al., 2013) platform of models for the representation of surface-atmosphere interactions by considering different surface types (vegetation, city, ocean, inland waters).
Meso-NH allows a multi-scale approach through a gridnesting technique (Stein et al., 2000).In this work, three nesting domains were used: a 400 × 432 km 2 domain with 4 km horizontal resolution to take into account the large-scale circulations, namely the influence of the sea breeze (Fig. 1a); an intermediate 96 × 72 km 2 domain with 1 km horizontal resolution centred at the Alqueva reservoir (Fig. 1b); and a finer 40 × 40 km 2 domain with 250 m spatial resolution to track the small-scale effects of the lake (Fig. 1c).Hereinafter we denote these three domains A-C correspondingly.The two-way nesting technique used in Meso-NH allows to conduct simulations on different horizontal resolutions at the same time.Domain A is the "father" domain for B, which means that simulation results in domain A are interpolated and used as initial and boundary conditions for domain B. The same scheme applies for domains B and C. European Centre for Medium-Range Weather Forecast (ECMWF) operational analyses, updated every 6 h, were used for Meso-NH initialization and domain A boundary forcing.
For surface and orography, ECOCLIMAP II (Faroux et al., 2013) and SRTM (Shuttle Radar Topography Mission, Jarvis et al., 2008) databases were used, respectively, both updated with the inclusion of the Alqueva reservoir by Policarpo et al. (2017).All model domains had 68 vertical levels starting with 20 m up to 22 km, including 36 levels for the lower atmospheric level (2 km).The model configuration included a turbulent scheme based on a one-dimensional 1.5 closure (Bougeault and Lacarrere, 1989), and a mixed-phase microphysical scheme for stratiform clouds and explicit precipitation (Cohard and Pinty, 2000;Cuxart et al., 2000), which distinguishes six classes of hydrometeors (water vapour, cloud water droplets, liquid water, ice, snow, and graupel), was used.Longwave and shortwave radiative transfer equations are solved for independent air columns (Fouquart and Bonnel, 1980;Morcrette, 1991).Atmosphere-surface exchanges are taken into account through physical parametrizations: the surface soil and vegetation are described by the Interface Soil Biosphere Atmosphere (ISBA) model (Noilhan and Mahfouf, 1996); the town energy balance was handled according to Masson (2000).Basic parameters for each model domain are shown in the Table 1.Deep and shallow convection parametrization schemes were activated in the coarser domain A. The 1 km and 250 m resolutions of domains B and C are fine enough for the deep and shallow convection to be represented explicitly.The following schemes were used (see Table 1): KAFR (Kain and Fritsch, 1990;Bechtold et al., 2001), EDKF (Pergaud et al., 2009), WENO (Lunet et al., 2017), and ICE3 (Pinty and Jabouille, 1998).
To track the impact of the reservoir on the weather conditions, two numerical simulation were performed: one with the surface input files updated to include the Alqueva reservoir (ECOCLIMAP database version updated by Policarpo et al., 2017) and another with the previous version of this database where the reservoir is absent.In order to distinguish these simulations hereinafter we denote them LAKE1 and LAKE0, correspondingly.Both simulations covered the case study period, 22-24 July 2014, with 1 h output.To reproduce the atmospheric conditions more realistically the simulations included the previous 24 h (21 July), so the model was integrated for 96 h.The differences between these two simulations were then computed, with the aim of evaluating the direct influence of the lake on the atmosphere.

FLake model
In order to better represent the evolution of the lake surface temperature and therefore the water-air heat fluxes, the FLake (Freshwater Lake) model (Mironov, 2008) was used.FLake is a bulk-type model capable of predicting the evolution of the lake water temperature at different depth on timescales from a few hours to many years.For an unfrozen lake it uses a two-layer approach: upper mixing layer, with a constant water temperature, and the thermocline beneath it where the temperature decreases with depth.Parametrization of the thermocline profile is based on the concept of selfsimilarity, assuming that such an approach could be applied to all natural and artificial freshwater lakes.The FLake model requires at least the following sets of variables and parameters to run: four initial parameters to describe the lake temperature structure, six atmospheric input variables for each time step, and two parameters -lake depth and the attenuation coefficient of light in the water.This coefficient is used to compute the penetration of the solar radiation in the water body.In this work, the attenuation coefficient was set to 0.85 based on in situ measurements carried out in Alqueva by Potes et al. (2017).
The FLake prognostic variables that need to be initialized are water temperature at the bottom, temperature and depth of the mixing layer, and shape factor C f -a parameter that describes the shape of the thermocline curve.In the parametrization proposed by Mironov (2008) for the normalized temperature profile it varies from 0.5 to 0.8.The initial values of the shape factor C f , a water mixing layer temperature, and depth were determined using a fitting technique applied to the observed water temperature profile at a given time in Montante platform (Fig. 2a).Short-term FLake model runs are very sensitive to initial parameters.The fitting technique is based on the assumption that the bottom temperature is fixed and given by the value of the lowermost sensor.Thereby, the other three parameters could vary within some range until the best set is found.The initial conditions for our simulations were obtained following this technique: C f = 0.8, mixed layer temperature is 23.8 • C, and depth is 3.4 m.Test simulation with these set of inputs was carried out using a stand-alone version of the FLake model (not coupled to Meso-NH).The results of the comparison of water mixing layer temperature is shown in Fig. 2b.The maximum difference does not exceed 0.8 • C, which is a very good result for such a short-term simulation.
The observed daytime temperature profiles showed strong skin effects (higher temperatures in the first 10 cm) and could not be well fitted by a FLake type temperature profile, which assumes a constant temperature in the mixed layer.Thus, the midnight profile was used initially and the simulation started at midnight, 21 July 2014.
The required atmospheric variables were taken interactively from the Meso-NH simulation since FLake was implemented in the SURFEX model by Salgado and Le Moigne (2010).
The depth of artificial lakes decreases rapidly from the centre to the shore, because the bottom of the reservoirs used to be valleys.In case of Alqueva, when completely filled, the mean depth is of about 17 m (http://www.edia.pt/pt/,28 September 2018).However, the local depth at Montante platform can reach 70 m.As a 1-D bulk model, FLake has only one depth value, which should be seen as an effective depth and is not easy to assess.Moreover, the FLake model is not capable of representing deep lakes; it works well for depths from 20 to 50 m with the sediments routine switched off.After a series of sensitivity tests of short-term (2-4 days) and long-term (2-4 months) simulation it was found that the best simulations results can be obtained with the bottom depth value of 20-30 m.Thus, since the deepest temperature probe was installed at the depth of 27 m, this value was chosen for the effective lake depth in this work.The comparison between measurements of water temperature near the surface (at 1 m depth) and FLake-simulated values of mixed layer temperature is shown in Fig. 2b.The sensor at 1 m depth was chosen because it always stays in the mixed layer and is not affected by surface "skin" effects.Modelled values are close to measurements, which indicates that the initial conditions were realistically imposed.

Validation
The simulation LAKE1 results were validated against radiosonde data (vertical profiles) and meteorological data from ALEX and IPMA stations located on land and in the floating platforms.All three domains were considered in the validation.The size of the domain A was enough to consider 12 synoptic stations located in the region, domain B was used to track the radiosondes trajectory, and domain C results were validated against stations installed on the lake shores and on the Montante floating platform.The variables under analysis were air temperature, relative humidity, wind speed, and sensible and latent heat fluxes.

Comparison with radiosonde data
The ALEX IOP took place between 22 and 24 July 2014 at the Alqueva reservoir.It included the launch of 18 meteorological balloons every 3 h.The radiosondes took measurements of air temperature, humidity, pressure, and wind speed.As the balloons did not ascend vertically and flew several kilometres away from the launching point, a trajectory profile comparison was performed.Each balloon had a GPS tracker to register its coordinates every 2 s, which was used to build a corresponding numerical trajectory on the simulation domain.Radiosondes reached the altitude of the top of the model (about 22 km in about 2.5 h).Therefore in order to build the simulated profile, three consecutive hourly outputs from the model were used.Figure 3a-c represent examples of the daytime profiles of air temperature, relative humidity, and horizontal wind speed.Examples of night profiles can be found in Fig. 3d-f.In general, the night-simulated profiles show slightly better accordance with model results.
Figure 4 shows the same profiles, but for the lower troposphere level (3000 m altitude).Simulations are in good accordance with observations.The simulation results are within the confidence interval of the measurements as given by the radiosonde accuracy -±0.5 • C, ±5 % relative humidity, and ±0.15 m s −1 wind speed with 2σ confidence level (95.5 %).
The principal features of the profiles are well represented by the model.During daytime, air temperature and relative humidity curves indicate that the model tends to represent the height of the boundary layer at 2-2.5 km altitude well (around 2 km in Fig. 4a and b).Overall, Meso-NH reproduces the air temperature above the surface layer (over 500 m) very well.Near the surface, the Meso-NH tends to anticipate the development of the unstable boundary layer in the morning (09:00 and 12:00 UTC), simulating higher temperatures in the lower levels.In the late afternoon (18:00 and 21:00 UTC) the model also tends to anticipate the decrease in the temperature in the surface layer (see the Supplement, Fig. S1).
The vertical patterns of relative humidity and wind speed are good, as observed and modelled curves look similar; nevertheless simulations tend to be more conservative and their values do not represent all changes in altitude recorded by the radiosonde.Nocturnal low-level jets at the edge of the boundary layer are represented by the model, but their magnitude is slightly weaker than that observed.The complete set of observed versus simulated profiles can be found in the Supplement (Fig. S2).The moisture vertical profile in the boundary layer is well reproduced by the model, as can be seen in the graphs of the relative humidity (Figs.4b, e, and S3).Above the boundary layer, the radiosondes show a dry layer, which is also well simulated.From 23 July the observations show the appearance of a moist layer close to the troposphere, the magnitude of which is overestimated by the model.At dawn on 24 July the radiosondes and the model indicate the existence of a very moist layer close to the surface, with the formation of low clouds that were not formed in the simulations.Statistical results for them are the following: temperature average bias is −0.13 • C, RMSE is 1.49 • C, and the correlation coefficient is 0.99; relative humidity average bias is 0.59 %, RMSE is 11.26 %, and the correlation coefficient is 0.87; and for the wind speed average bias is 0.05 m s −1 , RMSE is 2.07 m s −1 , and the correlation coefficient is 0.90.These values testify that the simulation is in a good accordance with the observations, in line with similar studies of Meso-NH validation against radiosonde data (e.g.Masciadri et al., 2013).

Comparison with IPMA stations data
The model was also validated against 12 IPMA automatic meteorological stations.For this comparison the output of the bigger domain A was used.Geographical positions of the stations can be found in Fig. 1a.Scatter plots of air temperature, relative humidity, and wind speed are shown in Fig. 5.It should be mentioned that not all stations provided the same set of variables.The scatter plots show the intercomparison of the model data (x axis) and the measured values (y axis) over all stations in the simulated period.The model tends to overestimate lower values of air temperature (14-24 • C) and slightly underestimate higher values (> 30 • C), as can be seen from Fig. 5a.
In some stations and for several times, the model simulates lower values of relative humidity within the range from 40 to 100 % (Fig. 5b). Figure 5c indicates that the wind speed is slightly underestimated by the model.
Statistical parameters (bias, mean absolute error, root mean square, and the correlation coefficient) for each station are shown in Table 2.
Simulated and observed air temperature are highly correlated (the correlation coefficient is higher that 0.91) with bias always less than 1 • .The worst results are observed in Portalegre (square points in relative humidity plot in Fig. 5b), and a possible reason for this relies on the location of the station, which is installed in a complex mountain area, possible not well represented in the model orography.Regarding wind speed, Meso-NH provides values for 10 m height while the measurements were made at 2 m.For proper comparison modelled values were interpolated to the height of the sensors using a known logarithmic approach and a roughness length (which was also obtained from the model).Table 2 shows small biases, in general lower than 1 m s −1 , and relatively high correlation coefficients for wind simulations (0.68-0.92).The lowest correlation coefficient is also obtained for Portalegre data.Overall, simulation results are in good agreement with synoptic station data, and the statis-   tical parameters are similar to other published works that use Meso-NH (e.g.Lascaux et al., 2013Lascaux et al., , 2015)).

Comparison with data from the ALEX database
In addition to the validation against the IPMA synoptic stations, comparisons were made with data obtained at ALEX meteorological stations (Barbosa, Cid Almeida, and Montante platform).Their coordinates were used to locate corresponding grid points on the C domain output.For land stations with grid points associated with water fraction, the nearest land grid point was chosen.
Figure 6 shows the time evolution of simulated and observed air temperature and wind speed at Cid Almeida, Bar-   bosa, and Montante sites.Overall, the simulation results are slightly more conservative (except wind speed over the Montante platform), but in general, the patterns are well represented.The model could not represent the maximum and minimum temperatures well, especially in land stations where the temperature range is larger.Regarding wind speed, the model underestimates the maximum values at land stations (Fig. 6e and f) and, on the contrary, overestimates the values at Montante platform (Fig. 6d), but the principal features of the curves are represented.Statistical values for this validation are as follows.For Barbosa, air temperature average bias is 0.23 • C, RMSE is 1.37 • C and the correlation coefficient is 0.98; for wind speed these values are as follows: average bias is 0.55 m s −1 , RMSE is 1.08 m s −1 , and the correlation coefficient is 0.73.For Cid Almeida, temperature average bias is 0.5 • C, RMSE is 1.57• C, and the correlation coefficient is 0.98, wind speed average bias is −0.36 m s −1 , RMSE is 1.24 m s −1 , and the correlation coefficient is 0.69.
For the Montante platform, air temperature average bias is −0.1 • C, RMSE is 1.22 • C, correlation is 0.98, wind speed average bias is −0.97 m s −1 , RMSE is 1.55 m s −1 , and the correlation coefficient is 0.61.The temporal evolution of simulated and observed latent and sensible heat fluxes at Montante platform is shown in Fig. 7a and b.Overall patterns of the curves are similar but the simulated one is more smooth.Comparison demonstrates that for latent heat the RMSE is 57.34 W m −2 with correlation coefficient of 0.47, and for sensible heat the RMSE is 13.39 W m −2 with the correlation of 0.82.
Both observed and simulated curves of Fig. 7a reveal that sensible heat flux has two different periods during the day, positive when air-water temperature difference is negative, and vice versa, showing that during daytime the water receives energy from the air and during nighttime the water warms the nearby air.The transition between the two regimes is well captured by the model.Figure 7b also shows that the magnitude of the sensible heat flux is relatively small when compared with the other terms of the energy balance.Daily maximum positive (negative) fluxes between 30 and 60 W m −2 (−15 and −30 W m −2 ) are well reproduced by the model.An apparently strange behaviour appears on the afternoon of 22 July, with the sensible heat flux being almost zero.This effect, unfortunately not documented due to the lack of data, will be discussed later and is linked to the fact that the wind is very weak during this period (see Fig. 6d).
More detailed analysis of Fig. 7a shows that the lowest heat flux values usually occur during the afternoon (12:00-18:00 UTC), under windless conditions, and high peaks in the early evening (20:00-21:00 UTC).The simulation reproduces these peaks with 1-2 h delays that are related to the delay on the simulated wind speed.The magnitude of the latent heat flux daily maximum (order of 200-250 W m −2 ) is well captured by the model.The delay in the simulation of the peaks reduces the value of the correlation coefficient and is a manifestation of the so-called double-penalty that penalize high-resolution model scores.As seen in Fig. 7b the simulated latent heat flux is almost zero between 14:00 and 16:00 UTC of 22 July.As pointed out before, there is a gap in the flux measurements during this period, but data from the day before indicate that the results are realistic.This effect of almost zero evaporation from water on a very hot day is contrary to common sense and will be discussed later.
Wind direction at ALEX stations is represented in Fig. 8. Different behaviour in wind direction between the two stations from 21 to 23 July is clearly seen from measurement data (green dots).In Barbosa station the wind changes from a north-west to a south regime during daytime while in Cid Almeida this effect is not observed.In the simulations this difference is not so clear, but is still visible during the afternoon on 22 July.Barbosa station, located on the north-west shore of the lake, indicates the presence of the lake breeze because its direction is the opposite to the dominant wind.At Cid Almeida station, on the south-east shore, the lake breeze is co-directed with the dominant wind.

Lake impact
To analyse the impact of the Alqueva reservoir on the local area the changes in the following atmospheric variables were considered: air temperature and potential temperature, relative humidity and water mixing ratio, and vertical and horizontal wind speed.In this section only B and C domain datasets were used.

Impact on air temperature and relative humidity
The lower layers of air are the first to be affected by the presence of the water surface.Differences in air temperature at 2 m during 22 July are shown in Fig. 9, the warmest day of the IOP and thus with a stronger lake breeze.The positive anomaly (up to 3-4 • C) can be traced during the period from 1 h after the sunset (21:00 UTC) to 1 h after the sunrise (07:00 UTC).By positive and negative anomalies here we mean differences between LAKE1 and LAKE0 simula- Hydrol.Earth Syst.Sci., 22, 5191-5210, 2018 www.hydrol-earth-syst-sci.net/22/5191/2018/ 2.5 m s -1 tions.Examples of positive night anomalies are illustrated in Fig. 9a and b.Night north-west wind transports warm air from the lake to the south-east part of the reservoir for up to 2 km away from the shore.The daytime period is characterized by negative temperature anomalies up to 7 • C (Fig. 9cf).Initially, the effect is essentially limited by the lake borders.When the large-scale sea breeze system arrives, temperature trace of the lake impact is advected by the wind and can be found in 10-12 km away from the south-east part of the reservoir (Fig. 9f).Vertical cross sections help to illustrate the processes at different altitudes.Two different cross sections Cs1 and Cs2 are shown to provide a better visualization of the threedimensional structure of air circulation above the lake.The first one crosses the lake near Montante platform and the sec-ond in the middle, exact locations are indicated in Fig. 9a.Cross sections Cs1 along 38.215 • N (Fig. 10) show the evolution of wind and potential temperature during 22 July in the experiment with Alqueva (simulation LAKE1).The highest impact on the air temperature can be observed in the early afternoon (12:00-14:00 UTC).The boundary layer is cooling down and its height decreases from more than 2 km above the land outside to values close to 1 km over the lake surface (Fig. 10a).The thermal anomaly induced by the presence of the reservoir seems to affect an area greater than what was identified at the surface, especially in the middle of the boundary layer.Later on, at 19:00-20:00 UTC the powerful ocean breeze system reaches the area and cools the lower (1 km) layer of air by 6-7 K.The progression of the sea breeze front is impressively well shown in Fig.  (20:00 UTC), when it reaches the border of the reservoir, and in Fig. 10e and f (21:00-22:00 UTC), when it is already beyond the east bank of the Alqueva reservoir.More cross sec-tions of 22 July 2014 with potential temperature can be found in Fig. S4.
Alqueva causes a similar anomaly on 2 m relative humidity, which is shown in Fig. 11.At night when the temperature impact is negative some small negative differences in relative humidity can be seen over the lake surface (Fig. 11a).There are also traces of daytime positive impact, essentially due to the decrease in air temperature, advected by the sea breeze in the south-east direction.In the morning, however, the difference of relative humidity cannot be detected because the thermal impact is not strong enough (Fig. 11b).Figure 11c-e show how relative humidity increases during the daytime over the water surface.The peak of the difference can reach 50 % in the af-ternoon (Fig. 11f).In general, lake impact on relative humidity is limited by the area of the reservoir and does not spread over the surrounding land.

Breeze effects
Differences in near-surface sensible heat fluxes and contrast of the air temperature over the land and water surfaces during the daytime induce the formation of the lake breeze system.The development of the lake breeze is illustrated in Fig. 9 (arrows that correspond to the wind speed lesser than 0.5 m s −1 www.hydrol-earth-syst-sci.net/22/5191/2018/ Hydrol.Earth Syst.Sci., 22, 5191-5210, 2018 Figure 13.Observed and simulated water vapour mixing ratio at the Montante platform. are not plotted).During the night, the large-scale circulation (Fig. 9a and b) is dominant in the area.After the sunrise (07:00-08:00 UTC) the air temperature over the water surface becomes lower than the air temperature over the surrounding areas, which induce a thermal circulation directed from the centre of the lake to its shores.The breeze intensifies during the afternoon reaching 6 m s −1 in some areas (Fig. 9d  and e).Daytime cross sections Cs1 in Fig. 10a-c indicate that the direct lake breeze can be found on altitudes up to 300 m above the lake, with a divergent flow over the water surface.The lake breeze intensity and pattern depends on the local orography, but usually the traces can be found 4-6 km away from the lake shores (Fig. 10c).In altitude, a return flow is visible in the eastward wind component over the west shore and a westward component in the east of the reservoir, which causes an upper-level convergence that can be seen in Fig. 10a-c.This will be discuss together with the effects of the reservoir on the moisture field, in which the structure of the lake breeze system is more visible.
In the late afternoon (18:00 UTC) the negative temperature anomaly due to the presence of the lake is getting weaker, and the breeze system starts to wane and dissipate.At 19:00-20:00 UTC the ocean breeze arrives to the area and overlaps the local circulations (Fig. 10d-f).
Cross sections Cs1 and Cs2 presented in Fig. 12 show that the lake breeze system includes a descending branch over the reservoir that carries dry air from a height of about 2-2.5 km and redistributes it over the lake surface.Additional cross sections illustrating this process can be found in the Supplement (Figs.S5 and S6).
This dry downstream is confirmed by the measurements of water vapour mixing ratio at the Montante platform.As can be seen in Fig. 13 the observed and the simulated mixing ratio of water vapour have a daily minimum with average val-ues of about 8-8.5 g kg −1 around 14:00-16:00 UTC.During the afternoon of 22 July, the day with a strong lake breeze, the minimum reached a value lower than 6 g kg −1 .Out of the period in which the air over the lake subsides, the water vapour mixing ratio returns back to 9-10.5 g kg −1 .The presence of this dry downstream was proposed as a hypothesis by Potes et al. (2017) and is proved through the performed simulations.In the same figure it is clearly seen that the model tends to overestimate the mixing ratio, except during the afternoon of 22 July.
However, Fig. 12 also shows that outside the reservoir there are zones of low-level convergence and upward motion that increase the moisture of the boundary layer and form some kind of lake breeze fronts.The complex shape of the reservoir implies also an complex 3-D structure of the breeze system.Towards the southernmost part, near the dam, the low-level divergent breeze circulation is very clear, but the convergence upper-level return current is weaker (Fig. 12ac).In contrast, near the middle of the reservoir (cross section Cs2 in Fig. 12d-f) where two water branches exist, the circulation near the surface is more complex due to the presence of a land area in between, but the subsidence motion is more prominent, inducing a decrease in mixing ratio through the boundary layer, which reaches a magnitude of about 4 g kg −1 at 16:00 UTC (Fig. 12f).
Figure 14 illustrates this process in a horizontal plane.At midnight (Fig. 14b) the reservoir does not directly affect vapour mixing ratio in the air.In the morning hours, when the sun rises, but the breeze system is not yet formed, a positive impact on the moisture over the lake can be seen due to the increase in the evaporation.This anomaly affects the air above the central and southern parts of the reservoir and is advected to nearby areas (Fig. 14c).Later in the afternoon, with the formation of the lake breeze, a negative impact can be traced over the water surface due to the descending branches of the local circulation (Fig. 14d and e).This explains the afternoon decrease in the water vapour mixing ratio observed at the Montante platform, as seen in Fig. 13.The localization of the area of this negative anomaly changes in time, but predominantly it is over the larger southern part of the reservoir.With the dissipation of the local lake breeze system and the arrival of the stronger large-scale north-western wind, the negative moisture anomaly over the reservoir disappears and a positive effect is visible in the downwind region (Fig. 14a and f), due to the increase in evaporation (note that Fig. 14a corresponds to the night of 21 to 22 July, when the effect was more noticeable).
During daytime, water temperature is lower than air temperature, which is associated with a very weak air circulation over the water surface, which leads to very low evaporation from the lake (refer to low latent heat flux values in Fig. 7a).This period of day is characterized by higher evaporation over the land than over water.By late afternoon when the dominant sea breeze system reaches the region, the northwestern wind accelerates significantly when passing over the smooth surface of the lake.As result, evaporation from the lake becomes very intense.

Conclusions
In this work the authors studied the formation and magnitude of the summer lake breeze in the Alqueva reservoir, southern Portugal, and its impacts on local weather.The study was based on Meso-NH simulations of a well documented case www.hydrol-earth-syst-sci.net/22/5191/2018/ Hydrol.Earth Syst.Sci., 22, 5191-5210, 2018 study of 22-24 July 2014.This period was used for several reasons.First, a large volume of meteorological data was collected during these days, which allowed for a validation of the simulation results.Secondly, this period was hot and dry, which is typical for most summer days in the region.
The model allowed simulations with horizontal resolution of 250 m, which is fine enough to resolve such relatively small-scale lake breeze and to spot the impact of the reservoir on the detailed local boundary layer structure.Due to the "youth" of the Alqueva reservoir it is possible to run an atmospheric model with the surface conditions prevailing before the filling of the reservoir.Two simulations, one with Alqueva and another one without it, allow the raw impact of the lake on the local weather regime to be evaluated.
Formation and dissipation of the daytime breeze system induced by the reservoir are described in the work.On hot summer mornings the difference between air temperatures above water and neighbouring land surfaces induces the radial movement of air from the lake.The breeze system starts to form in the morning and the peak of the wind speed reaches 6 m s −1 during the afternoon.Simulation results show that the lake breeze could be detected at a distance of more than 6 km away from the shores and on altitudes up to 300 m above the water surface.In late afternoon the dissipation stage of the lake breeze system is anticipated with the arrival of the larger-scale sea breeze from the Portuguese west Atlantic coast.In early evening (19:00-20:00 UTC) the local lake breeze system cannot be detected anymore.No reverse land breeze is detected during the night.
Lake breeze system brings dry air from upper atmospheric layers (2-2.5 km) to near-surface levels above the reservoir.This effect leads to the fact that the air above the surface of the lake becomes more dry in terms of water vapour mixing ratio, in spite of the fact that its relative humidity can increase up to 50 % due to the decrease in air temperature.
The simulations testify to the observed very low evaporation from water surface during the daytime (0-120 W m −2 in terms of latent heat flux), due to weak winds and the stable stratification of the internal atmospheric surface layer.At nighttime, the strong wind associated with the peninsular larger-scale circulation induced by the sea-land contrasts, originates a very high evaporation rate (200-250 W m −2 ).
The cooling effect of the reservoir can decrease the air temperature by up to 7 • C; nevertheless it is limited by the lake borders and normally cannot be seen farther than few kilometres away from the shore, mostly in the south-east direction.The cooling can be found up to 1200 m above the lake surface.
Further work suggests two options: first, tuning the lake model and its initialization in order to obtain more accurate results and reduce validation biases; second, carry out a longer experiment, which would cover a 12-month period.Such a simulation could reveal seasonal aspects of the impact of Alqueva on local weather.

Figure 1 .
Figure 1.Nested domains used in the simulations: (a) the "father" domain at 4 km horizontal resolution with 100 × 108 grid points, with location of the 12 IPMA synoptic stations used for validation; (b) intermediate 1 km horizontal resolution domain, 96 × 72 grid points; (c) finer 250 m resolution domain comprising 160 × 160 grid points, together with the location of the ALEX land stations, the Montante floating platform and the dam.

Figure 2 .
Figure 2. Water temperature observed and fitted profiles on 21 July, 00:00 UTC at the Montante platform (a) and comparison of upper level water temperature between measurements and FLake results (b).

Figure 3 .
Figure 3. Examples of observed and simulated vertical profiles for 22 July at 12:00 UTC and 23 July at 03:00 UTC of air temperature (a, d), relative humidity (b, e), and wind speed (c, f).

Figure 4 .
Figure 4. Profiles of air temperature, relative humidity, and wind speed in low atmosphere on 22 July at 12:00 UTC and 23 July at 03:00 UTC.

Figure 5 .
Figure 5. Scatter plots of the comparison between Meso-NH simulation LAKE1 and measured values at synoptic stations.Air temperature (a), relative humidity (b), and horizontal wind speed (c).

Figure 6 .
Figure 6.Observed and simulated air temperature and wind speed at 2 m for ALEX stations: Montante platform (a, d), Barbosa (b, e), and Cid Almeida (c, f) sites.

Figure 7 .
Figure 7. Observed and simulated latent (a) and sensible (b) heat fluxes at Montante floating platform.

Figure 8 .
Figure 8. Observed and simulated wind direction at Barbosa and Cid Almeida stations.

MCs2Figure 9 .
Figure 9. Anomalies in air temperature at 2 m (in filled contours) and horizontal wind in LAKE1 experiment (arrays, the scale is indicated in the upper right corner of each figure) of the reservoir on 22 July 2014 in domain C. Horizontal lines in (a) indicate the location of cross sections Cs1 (southern) and Cs2 (northern).

Figure 10 .
Figure 10.East-west direction cross sections along 38.215 • N (Cs1, crosses the lake near Montante platform, Fig. 9a) of potential temperature (filled contours) and wind vectors in the plane of the cross section (arrows), at different hours (indicated in the top of each figure) in the LAKE1 experiment at 250 m horizontal resolution.The wind vertical and horizontal scales are indicated in the upper right corner of each figure.The blue line on the surface level indicates the location of the reservoir. 10d

Figure 11 .
Figure 11. 2 m relative humidity anomalies (in filled contours) on 22 July 2014 in domain C.

Figure 12 .
Figure 12.East-west direction cross sections Cs1 along 38.215 • N (a-c) and Cs2 along 38.274 • N (d-f) with the difference (LAKE1 minus LAKE0 simulations) of water mixing ratio (filled contours), and wind vectors in the plane of the cross section (arrows) in the LAKE1 experiment at 250 m horizontal resolution at different times (indicated in the top of each figure).The blue line on the surface level indicates the location of the reservoir.

Figure 14 .
Figure 14.Observed and simulated water vapour mixing ratio anomalies in filled contours on 22 July 2014 in domain C for selected hours.

Table 1 .
Summary of the Meso-NH physical schemes used in the simulations.
M. Iakunin et al.: Breeze effects at a large artificial lake: summer case study

Table 2 .
Statistics for the hourly values of the station validation.