Investigations of boundary layer structure , cloud characteristics and vertical mixing of aerosols at Barbados with large eddy simulations

Large eddy simulations (LESs) are performed for the area of the Caribbean island Barbados to investigate island effects on boundary layer modification, cloud generation and vertical mixing of aerosols. Due to the presence of a topographically structured island surface in the domain center, the model setup has to be designed with open lateral boundaries. In order to generate inflow turbulence consistent with the upstream marine boundary layer forcing, we use the cell perturbation method based on finite amplitude potential temperature perturbations. In this work, this method is for the first time tested and validated for moist boundary layer simulations with open lateral boundary conditions. Observational data obtained from the SALTRACE field campaign is used for both model initialization and a comparison with Doppler wind and Raman lidar data. Several numerical sensitivity tests are carried out to demonstrate the problems related to “gray zone modeling” when using coarser spatial grid spacings beyond the inertial subrange of three-dimensional turbulence or when the turbulent marine boundary layer flow is replaced by laminar winds. Especially cloud properties in the downwind area west of Barbados are markedly affected in these kinds of simulations. Results of an additional simulation with a strong trade-wind inversion reveal its effect on cloud layer depth and location. Saharan dust layers that reach Barbados via long-range transport over the North Atlantic are included as passive tracers in the model. Effects of layer thinning, subsidence and turbulent downward transport near the layer bottom at z≈ 1800 m become apparent. The exact position of these layers and strength of downward mixing is found to be mainly controlled atmospheric stability (especially inversion strength) and wind shear. Comparisons of LES model output with wind lidar data show similarities in the downwind vertical wind structure. Additionally, the model results accurately reproduce the development of the daytime convective boundary layer measured by the Raman lidar.

Abstract.Large eddy simulations (LESs) are performed for the area of the Caribbean island Barbados to investigate island effects on boundary layer modification, cloud generation and vertical mixing of aerosols.Due to the presence of a topographically structured island surface in the domain center, the model setup has to be designed with open lateral boundaries.In order to generate inflow turbulence consistent with the upstream marine boundary layer forcing, we use the cell perturbation method based on finite amplitude potential temperature perturbations.In this work, this method is for the first time tested and validated for moist boundary layer simulations with open lateral boundary conditions.Observational data obtained from the SALTRACE field campaign is used for both model initialization and a comparison with Doppler wind and Raman lidar data.Several numerical sensitivity tests are carried out to demonstrate the problems related to "gray zone modeling" when using coarser spatial grid spacings beyond the inertial subrange of three-dimensional turbulence or when the turbulent marine boundary layer flow is replaced by laminar winds.Especially cloud properties in the downwind area west of Barbados are markedly affected in these kinds of simulations.Results of an additional simulation with a strong trade-wind inversion reveal its effect on cloud layer depth and location.Saharan dust layers that reach Barbados via long-range transport over the North Atlantic are included as passive tracers in the model.Effects of layer thinning, subsidence and turbulent downward transport near the layer bottom at z ≈ 1800 m become apparent.The exact position of these layers and strength of downward mixing is found to be mainly controlled atmospheric stability (especially inversion strength) and wind shear.Comparisons of LES model output with wind lidar data show similarities in the downwind vertical wind structure.Additionally, the model results accurately reproduce the development of the daytime convective boundary layer measured by the Raman lidar.

Introduction
A series of ground-based and airborne remote sensing measurements took place at and around Barbados during the SALTRACE (Saharan Aerosol Long-range Transport and Aerosol-Cloud-Interaction Experiment) 2013 summer campaign.Since Barbados is the easternmost island in the Caribbean and steady easterly trade winds are present, it is not affected by other surrounding islands.For that reason, Barbados is suitable for island effect studies both from the measurement and the modeling point of view.First of all, mineral dust emitted from the Saharan region is transported for more than 4000 km over the Atlantic Ocean with almost no anthropogenic influence.Dust layers arriving at Barbados can be detected with respect to layer height and thickness as well as aerosol composition.Secondly, cloud studies are possible due to persistent trade-wind circulation at the eastern Caribbean.For example, extensive investigations on shallow cumulus cloud properties and their response to different ambient cloud condensation nuclei (CCN) num-ber concentrations took place during the CARRIBA (Cloud, Aerosol, Radiation and tuRbulence in the trade wInd regime over BArbados) project in 2010/2011 (Siebert et al., 2013).Within CARRIBA, airborne in situ measurements were conducted east of Barbados.The field site of the Max Planck Institute for Meteorology (MPI-M), Hamburg, Germany, with ground-based instruments is located at the east coast as well.The choice of these locations ensures that the island itself has very little to no influence on the measurements and thus marine boundary layer properties can be accurately investigated.During SALTRACE, the TROPOS (Leibniz Institute for Tropospheric Research, Leipzig) and LMU (Ludwig-Maximilians-Universität Munich) field sites were located at the area of the local Caribbean Institute for Meteorology and Hydrology (CIMH) near the west coast of Barbados (see Fig. 1), whereas the DLR (Deutsches Zentrum für Luftund Raumfahrt) research aircraft Falcon was stationed at the international airport of Barbados.Regarding the measurement field site, incoming air masses at these sites are already influenced by the island due to surface roughness change, different energy fluxes and topographical features.Whereas the first two properties primarily influence the atmospheric boundary layer (ABL), gravity waves caused by the latter also propagate within the free troposphere.
There are several works regarding the understanding of airflow and thermodynamic quantities around Barbados.A first detailed observational study using pilot balloon measurements was done by DeSouza (1972) and further interpreted by Garstang et al. (1975).DeSouza's calculated vertical wind velocity fields showed a daytime divergence and nighttime convergence over the island.Mahrer and Pielke (1976) did a series of two-and three-dimensional numerical studies and found that DeSouza's calculations only hold for a flat island, because he neglected significant effects of terrain slope in his divergence calculations.Heat island effects on vertical mixing of aerosols at Cape Verde islands were studied by Engelmann et al. (2011) using aircraft lidar measurements and idealized large eddy simulations (LESs) with flat island surfaces.They found indications that the differential heating and the orographic impact control downward mixing of African aerosols, which results in a complex vertical layering over the Cape Verde region.Taking the topographical structures into account, Mahrer and Pielke (1976) pointed out some main characteristics, e.g., diurnal changes in the vertical wind velocity fields downwind (i.e., west coast of Barbados) with sinking motions over the center and western part of the island and an upwind cell off the west coast.Considering numerical sensitivity studies by Savijärvi and Matthews (2004, SM04 hereafter), the general conclusion was that these forced rising and sinking motions and their consecutive effects can only be explained if island orography is included in the numerical models.In their 2-D study, SM04 added a 200 m high central mountain to a 20 km wide island and showed that sea-breeze circulations are enhanced by upslope winds during the day.These topographically forced components will dominate if the large-scale mean wind is in the order of magnitude of at least 10 m s −1 , which is the case for Barbados.Smith et al. (1997) assigned different island structures to different mountain wake types.Since the highest elevation of Barbados (Mt.Hillaby, with a height of 340 m a.s.l.) does not exceed the critical height for wave breaking, no wind wake can develop.However, a long trail of cumulus clouds extending westwards can evolve during the afternoon hours (cumulus cloud street).Kirshbaum and Fairman (2015) found that surface fluxes control the downwind circulation strength and the trade inversion controls precipitation and thus the disruption of cloud trails.Other influence factors like terrain height, wind speed and their interactions have multiple impacts on flow regimes, turbulence, cloud trail lengths etc.Another study on island effects with similar topographical heights compared to Barbados was done by Minda et al. (2010).They investigated the evolution of the convective boundary layer (CBL) above Okinawa Island, Japan.It was found that for a flat island simulation, the warmed land already induces a distinct roll cloud that is in agreement with the observations.However, the inclusion of island terrain leads to reinforced moisture uplifts, which in turn induce strong convection that can penetrate into the free atmosphere.Idealized numerical studies were conducted by Kirshbaum and Grant (2012) to investigate the impact of mesoscale ascent (with an island height of 500 m) on cumulus convection.There, a particularly important process with regard to the mean horizontal cloud size has been found.The broader the clouds are, the lower is the fractional entrainment rate in these clouds, which in the end leads to an increase in precipitation rates downstream.A key result from another combined theoretical and numerical study by Kirshbaum and Wang (2014) was that nonlinear interactions between mechanical and thermal flow over taller mountains were significant and thus lead to a strengthening of the leeside convergence band.
There are also many studies where the focus lies on the orographic influence of tall islands (e.g., Hawaii Island or Dominica with mountain heights above 1 km) on the leeward flow and precipitation patterns.Esteban and Chen (2008) state that for a strong trade-wind flow, the daily rainfall totals at the windward side of the island of Hawaii show a nocturnal maximum due to the convergence of katabatic flow, whereas for weak trades (≤ 5 m s −1 ) the rainfall amounts have their maximum in the late afternoon due to anabatic winds.In a work by Smith et al. (2009), orographic precipitation for the Caribbean island Dominica was studied.There, the conditionally unstable trade-wind layer together with terrainforced lifting leads to convective triggering over the windward slope.The reduced instability on the lee side destroys convective clouds and creates a rainless area.A complementary study with airborne observations and cloud-resolving modeling for the same island was performed by Minder et al. (2013).The comparison showed that the dynamical structures are very well reproduced but that it was difficult to reproduce the observed rainfall using the model.Overall, mesoscale flow controls convection and rainfall over Dominica.At lower wind speeds, the circulations seem to be more thermally driven by solar heating.

Atmos
The main objective of this work is to study local island effects on the modification of the boundary layer structure, microphysical properties and downwind vertical mixing of aerosols for selected days during the first SALTRACE field campaign.Regarding aerosols, especially Saharan dust, it is known from several studies that notable amounts of mineral dust reach Barbados via long-range transport over the North Atlantic, e.g., from first observations at the end of the 1960s (Prospero et al., 1970;Prospero and Carlson, 1970) or from back-trajectory calculations by Ellis and Merrill (1995).
Within this work, the following questions are addressed: -How does the model setup have to be chosen to get an as realistic as possible representation of an island-ocean system in the trade-wind regime through the example of Barbados?
-How do turbulent inflow characteristics and grid spacing affect the simulation results?
-Can the daytime convective island boundary layer explain downward mixing of low-altitude Saharan dust layers?
-Are the simulation results comparable with lidar measurements over and in the lee of the island?
This paper is structured as follows.Section 2 deals with the general model setup.There, the numerical method, model physics, computational domain, boundary conditions, initial data and forcings are described.To generate a turbulent marine boundary layer, a novel method based on potential temperature perturbations is adopted, verified and applied to our numerical model and particular setup.Results for two case studies in June 2013 and two sensitivity tests are presented and discussed in Sect.3. In Sect.4, the simulation results are compared with stationary and airborne lidar data.Section 5 provides a summary and concluding remarks.

Model setup
All LESs are performed with the latest version of the non-hydrostatic, fully compressible All Scale Atmospheric Model (ASAM).An extensive model description is presented in Jähn et al. (2015), both covering numerical discretization methods and physical parameterizations.A special feature of ASAM is the usage of so-called cut cells for the orography.There, a grid box is cut by the intersection of the orographical structure.This method can handle steep terrain gradients and prevents discretization errors compared to traditional methods like terrain-following coordinates, also conserving the original shape of the topography to a high degree.The dynamical core solves the flux-form tendency equations for mass, momentum, energy (in terms of density potential temperature) and other scalars.The most important physical parameterizations include a Smagorinsky subgrid-scale model and a two-moment cloud microphysics scheme.Further details on the model are described in Appendix A.
In the next subsections, the computational domain, boundary conditions (BC), data initialization and forcings for the cases of study are described, followed by a novel method to generate inflow turbulence.

Domain and boundary conditions
To simulate atmospheric flow for the island-ocean system, the size of the model domain has to have appropriate values dependent on the island size.The main criterion in this case is that a marine boundary layer has to develop at least several kilometers before it interferes with the island area.Also, the downwind area should approximately be twice of the island width so that resulting structures induced by the island can be properly represented.Since Barbados is a 24 km wide (west-east) and 34 km long (south-north) island, a model domain with a spatial extent of 102.4 km × 102.4 km is chosen.The island is located at the domain center.The model top is set to 5 km altitude.Because of the required domain size and for computational reasons, the horizontal grid spacing is set to x = y = 200 m.Such a resolution can be considered a "coarse" LES; however, it is sufficient to resolve some portion of inertial range scales, as will be shown later on through spectral analysis.
Due to the presence of the island area, non-cyclic lateral boundary conditions have to be used.Within the finite volumes/differences discretization strategy adopted herein, a "zero-gradient" boundary condition is applied to all scalars and velocity components at each lateral boundary (north, east, south, west).This means that the boundaryperpendicular flux for these quantities is set to zero, which leads to a simple radiation condition near the outlets with minimal wave reflection.A pressure correction for sound waves is applied to each actual normal velocity component and not to the initial wind profile, which also suppresses artificial wave reflection near the inflow boundary.This setup ensures stability for the whole simulation time and works appropriately with the turbulence generation method, as shown at the end of this section.
For the top boundary, a free-slip condition is applied, i.e., the gradient of the tangential velocity component is zero.In order to prevent gravity wave reflection, an additional relaxation term is applied on the right-hand side of the momentum equations: with a damping function depending on the distance to the top boundary d: (2) This damping layer is applied above d w = 4 km model height (20 vertical layers) with a damping parameter d f = 1 × 10 −3 .
Surface boundary conditions are represented by a momentum flux parameterization based on the Monin-Obukhov similarity theory (Monin and Obukhov, 1954): C m is the drag coefficient for momentum, which is defined as follows: and φ m representing the integrated similarity function.L stands for the Obukhov length and k is the von Kármán constant.The topographic data are obtained from the Consortium for Spatial Information (CGIAR-CSI) Shuttle Radar Topography Mission (SRTM) data set (http://srtm.csi.cgiar.org)at 200 m resolution.A simple smoothing algorithm is applied to guarantee a proper grid pre-processing.In the smoothed data set, the maximum elevation is lowered by about 15 m compared to the raw topography data, which is an acceptable level.
Table 1 summarizes the model configuration for the Barbados LESs performed in Sect.3.

Initial data
The two cases examined (22 and 27 June 2013) mainly differ in their atmospheric state and geostrophic forcing.Measured nighttime radiosonde profiles of temperature and humidity are directly used for model initialization (Fig. 2), which reduces the complexity of the simulations due to the absence of horizontal inhomogeneities and a time-varying background state.There are two reasons behind the choice of using single profiles instead of averaging multiple profiles.Firstly, a single initial profile is better for comparing the LES results with lidar data (cf.Sect.4), which are obtained for a few selected cases during SALTRACE.Secondly, trade-wind inversions are only poorly represented when the soundings are averaged over many cases.This becomes apparent when considering the sharply defined inversion at the 27 June case, which is shown later on.Air density and pressure profiles are obtained by vertical integration with respect to hydrostatic equilibrium.Some simplifications are assumed for the geostrophic forcing.The wind direction is purely east (i.e., d = 90 • and v g = 0), which is also for simplicity and to make it easier to define upwind and downwind regimes later on.The vertical wind profiles are expressed as piecewise linear functions for both cases.For the 22 June case, the initial wind at first linearly decreases above z = 1600 m altitude and then increases 3.0 km < z ≤ 5.0 km, with a roughness length z 0 = 0.01 m.A change in wind direction to southwest is observed within the layer where the wind speed decreases.However, this is not captured by the LES due to the simplifications and assumptions mentioned above.Therefore, the effect of wind directional shear might be underestimated in the model for this case.The change in wind direction (±15 • ) is rather small at other altitudes, so the LES input profile can be considered a good approximation.For the 27 June case, the initial wind linearly decreases above z = 3000 m altitude: In this profile there is no distinct change in wind direction.Figure 3 visualizes the measured (green lines) and parameterized (red lines) velocity profiles for both cases.The LES background wind profiles are parameterized to closely match the soundings.Within the boundary layer, the LES profile should be near the nighttime measurements (dark-green line) because this is mainly a representation of the marine boundary layer.For the free troposphere, the LES profile should roughly be a mean of all three soundings, since no large-scale advection term is applied on the wind components during the simulation time.
Table 2 shows a comparison of the two simulated cases with respect to mean flow properties, trade inversion strength, moisture load (all derived from radiosonde profiles), CCN concentrations (obtained by ground-based measurements at Ragged Point station) and the location of the Saharan dust layer (estimated from BERTHA lidar measurements at CIMH).The differences in the geostrophic forcing are already discussed.Regarding the atmospheric stability, there is a much stronger trade inversion for the 27 June case with a local virtual potential temperature gradient of 14 K km −1 .As mentioned in the introduction, the trade inversion controls the amount of precipitation and the lifetime of cloud streets.Furthermore, there is a 18 % stronger moisture load for the 27 June case, where a faster cloud development  CCN number concentrations, a mean value of 300 cm −3 has been chosen for both cases, which is a typical magnitude for days with a moderate dust load, where aerosol optical depths between 0.2 and 0.4 are observed.

Forcings
Surface sensible and latent heat fluxes over the island and the ocean are obtained by separate 1-D simulations with full model physics.The parameterizations there include the radiation scheme (Fu and Liou, 1993) as well as land-use and soil models.The soil class "loam" was chosen to represent the average island soil type.Hydraulic and thermal parameters of this soil type can be found in Doms et al. (2011) and Jähn et al. (2015).For land surface parameterization, "shrubland" appears to be a good compromise between coastal beach areas and forest in the island interior.The roughness length of this land type is z R,island = 0.2 m, whereas the ocean roughness length is set to z R,ocean = 0.01 m.The usage of direct (compared to interactive) fluxes reduces computational costs for the LES runs and makes it easier to potentially reproduce these simulations by other models, especially due to a large number of existing radiation and land-use models.Figure 4 shows the diurnal variation in sensible and latent heat fluxes over the island area.The maximum sensible heat flux over the island is Q SHF max ,island = 425 W m −2 and the corresponding maximum latent heat flux is Q LHF max ,island = 105 W m −2 .Surface heat fluxes over the ocean are constant during the whole simulation time with Q SHF,ocean = 6 W m −2 and Q LHF,ocean = 56 W m −2 .Sunrise is at 05:36 LT and sunset is at 18:29 LT, whereby the fluxes are shifted by 30 min to represent the delay due the fact that the soil has to be heated first before energy exchange with the lower atmosphere can take place.
Later on, reference simulations with periodic boundary conditions are performed to obtain information of marine boundary layer characteristics.For these simulations, large-scale forcings from the BOMEX LES study of trade-wind cumulus convection (Siebesma et al., 2003) are applied.They include a piecewise linear subsidence velocity profile with an absolute peak value of −560 m day −1 , radiative cooling of −2 K day −1 and large-scale advection of dry air into the lower boundary layer of −1 g kg −1 day −1 :

Turbulence generation -the cell perturbation method
The LES modeling technique has the advantage of allowing explicit resolution of turbulent production and part of the inertial range scales, and is today the most accurate and computationally feasible modeling approach in the context of high Reynolds number flows.LES results are strongly dependent on boundary conditions, therefore requiring specification of realistic inflow turbulence characteristics that propagate through the domain into the area of interest.In order to ensure that the incoming boundary layer characteristics at Barbados correspond to fully developed turbulence consistent with the imposed marine boundary layer forcing, we use the cell perturbation method recently proposed by Muñoz-Esparza et al. (2014).The cell perturbation method uses a novel stochastic approach based upon finite amplitude perturbations of the potential temperature field applied within a region near the inflow boundaries of the LES domain.This method has demonstrated superior performance when compared to a state-of-the-art synthetic turbulence generator and is computationally inexpensive (Muñoz-Esparza et al., 2015).
Previous studies where the cell perturbation method was developed and validated dealt with transitions from smooth mesoscale flow to nested LES (Muñoz-Esparza et al., 2014, 2015).In these idealized cases, boundary conditions at the LES domain boundaries were imposed from the mesoscale model instantaneous solution (Dirichlet boundary conditions), in which moisture effects were not considered.Herein, we further extend the application of the cell perturbation method to turbulence inflow generation for cloud modeling including terrain effects.As explained in earlier sections, zero-gradient open radiative lateral boundary conditions need to be used in order to minimize wave reflections at the boundaries that do develop in fully compressible codes like the ASAM LES model when the domain includes terrain features.In order to test the best configuration for the cell perturbation method in this particular context, we perform a series of calculations where only the upstream region of the Barbados island is considered (i.e., incoming marine boundary layer).The reduced subset of the domain consists of a 51.2 km× 51.2 km area in the horizontal, with the same vertical extent and large-scale forcing described in Sect.2.3 for the 22 June 2013 case study.To represent the marine boundary layer conditions that are going to be imposed through the entire simulation period, constant sensible and latent heat fluxes of Q SHF,ocean = 6 W m −2 and Q LHF,ocean = 56 W m −2 are used (see Fig. 4).
We explore the sensitivity of the generated turbulence by the cell perturbation method to the optimum perturbation Eckert number, Ec = U 2 g /c p θ pm = 0.2, where θ pm is the maximum potential temperature perturbation, and the perturbations are random and uniformly distributed in the interval − θ pm , + θ pm .Three square cells adjacent to the east boundary are used, which were found to provide the fastest transition to a fully developed turbulent state (Muñoz-Esparza et al., 2015).The cell size is set to 4 × 4 grid points to ensure that the cell wavelength falls within the inertial range of three-dimensional turbulence.The perturbation timescale, t p , was obtained from = t p U 1 (4dx) −1 = 1 (Muñoz-Esparza et al., 2015), with U 1 being the horizontal wind speed in the first vertical layer, resulting in a frequency to seed instantaneous perturbations of t p = 145 s. Figure 5 shows instantaneous contours of vertical velocity at z = z i /2 = 375 m for different perturbation Eckert numbers, Ec = 0.2, 0.33 and 0.4, and for the periodic reference run.The cell perturbation method for the three Ec numbers considerably accelerates the formation of three-dimensional turbulent structures that agree with the ones obtained in the reference simulation using periodic lateral boundary conditions.As the perturbation Eckert number increases (maximum perturbation amplitude decreases), the strength of the vertical velocities induced by the temperature perturbations is progressively reduced, and the onset of forcing-consistent turbulence seems to qualitatively occur at earlier distances from the inflow boundary.
In order to have a better understanding of the turbulence initiation and development processes, the energy spectrum evolution in the streamwise direction for the three velocity components is presented in Fig. 6.The cell perturbation method causes a rapid development of the upperwavenumber portion of the energy spectrum for the u and v components.The larger scales (lower wavenumbers) require longer distances to be established due to large buoyant plumes having to emerge from the surface and populate across the entire extent of the boundary layer.This flow development pattern is consistent with the findings from Muñoz-Esparza et al. ( 2014) for convective conditions.In contrast, the energy spectrum for the vertical velocity reveals a rapid growth of turbulent energy that reaches levels 10 times greater than the periodic quasi-equilibrium solution (dashed black line) and that progressively dissipates as the flow transitions through the domain.We attribute this behavior to the cell size, 4dx, which for the resolution employed in this study may fall in the vicinity of the limit of the inertial range.Smaller cell sizes were not considered due to the energy dissipation at high wavenumbers present in finite differences/volumes discretizations.There, an interaction with fully resolved scales and triggering of an accelerated transition to a developed turbulence state would not have taken place.In addition, the use of zero-gradient lateral boundary conditions helps to maintain the signature of the perturbations more than in the case of Dirichlet boundary conditions, hence contributing to strengthen the periodically seeded perturbations.By increasing the perturbation Eckert number from 0.2 to 0.4 (first row vs. third row in Fig. 6), the energy overestimation is damped, and results after a fetch of 40 km for Ec = 0.4 are in close agreement with the periodic simulation used as a reference and have reached quasiequilibrium converged statistics.The Ec = 0.2 case results in an energy deficit at wavenumbers close to the integral length scale, and also at the highest wavenumbers for the w component.When the cell perturbation method is not used (NOCP panels, bottom row in Fig. 6), dramatic energy deficits are found, together with an unrealistic spiky energy distribution in which the expected energy production and cascade processes are not present.
Finally, we examine the vertical distribution of relevant boundary layer quantities at a downstream distance of 40 km from the east boundary (i.e., x = 11.2 km).Vertical profiles (Fig. 7) show the best agreement with the periodic simulation for the Ec = 0.4 and 0.33 cases, in particular for the turbulent kinetic energy levels and boundary layer structure.Momentum flux profiles exhibit slightly larger values in the first 250 m, due to the differences in the horizontal wind speed distribution near the surface.However, the boundary layer structure is similar, with the differences being related to distinct quasi-equilibrium solutions for the periodic and the open boundary condition simulations.Similar conclusions are found for the sensible and latent heat fluxes.The cell perturbation method was originally developed and tested in the context of dry boundary layers (Muñoz-Esparza et al., 2014, 2015).It is worth emphasizing that we have herein demonstrated for the first time, as can be seen from the latent heat flux profile, that the cell perturbation method has the ability to develop turbulent moisture features that are in agreement with the imposed forcing.The Ec = 0.2 case fails to produce a boundary layer structure that is similar to the reference periodic case, with excessive mixing attributed to an enhanced effect of the perturbations for the reasons mentioned above.Also, the NOCP case does not provide realistic turbulent boundary layer features corresponding to a strongly underdeveloped turbulent state.Therefore, we select the Ec = 0.4 setup as the inflow to be used for the island simulations presented in the remainder of the paper since it produces the most rapid development and stabilization of forcing-consistent turbulence.For the island cases, we use a domain with horizontal extent of 102.4 km × 102.4 km, which leaves sufficient fetch for the marine boundary layer to develop prior to the start of interaction with the topography of Barbados and its local stability effects.

Results
To investigate the effects of the Barbados island area on boundary layer properties, cloud generation and vertical mixing of aerosols, we define two subdomains that are considered to be representative of the upwind and downwind area, respectively.Figure 8 shows the position of these two subdomains.They both cover a base area of 10 km × 20 km and are used for averaging of vertical profiles and time series of the relevant quantities.The upwind domain east of Barbados (representing the marine boundary layer) is approximately 15 km away from the eastern boundary to avoid contamination from the inflow boundaries where turbulence has to be generated first.Looking into the model data, it becomes apparent that the flow has to pass at least half of the island area (≈ 12 km) before a well-mixed convective layer can fully de-velop.For that reason, the downwind subdomain is located between 35 km < x < 45 km and thus covers the west coast island area and the marine offshore area in equal parts.The following analysis mainly consists of comparisons between these two regimes to investigate island effects on various parameters.

Overview of simulations performed
Besides the two mentioned case studies, two additional sensitivity studies are part of the island effect analysis.Here, the 22 June 2013 case serves as a reference case (REF).For the first sensitivity case (NOCP), the cell perturbation method is disabled so that the upwind flow is strongly underdeveloped.With this setup, the effect of having a realistic turbulent boundary layer around the island rather than idealized constant winds is investigated.In the next sensitivity case (DX400), the grid resolution is halved from 200 to 400 m horizontally and from 50 to 100 m vertically to point out the deficiencies in the use of coarser resolution without appropriate resolved turbulence and gray zone modeling (Wyngaard, 2004) for particular aspects of interest in boundary layers and cloud modeling.In this simulation, the cell perturbation method is also put off since the usage of a turbulent inflow in coarse resolution studies has not been utilized before and, moreover, appears to be questionable because the inertial subrange of the turbulence spectrum is not resolved anymore.The simulation ensemble is completed by the 27  3 summarizes the settings for all simulations that deviate from the standard configuration in Tables 1 and 2.

Boundary layer and cloud characteristics
To get a qualitative impression of the local situation simulated by the LES model, Fig. 9 shows a three-dimensional snapshot of the temperature and humidity field as well as cumulus clouds with up-and downdrafts visualized by isosurface fields at 12:00 LT for the reference case.The daytime convection is clearly visible by multiple updraft cells distributed over the whole island area, which subsequently leads to the development of non-precipitating shallow cumulus clouds.Advection of heated air from the central and southern part of the island towards the west can be seen in the surface temperature field (which is meant as temperature of the lowest model layer in this context), whereas the cooler marine flow narrows the thermal wake toward the meridional center of the domain up to 40 km downwind.This effect is connected with an island-induced change in wind speed and direction.The change in the humidity profile can be observed in the vertical cut plane at the western model boundary.A large amount of moisture is transported vertically upwards in the central region where also occasional cumulus clouds are present.A few tens of kilometers away in the y direction, dryer air from heights of 500-1000 m is mixed downward.
For further insight into flow dynamics, especially for the downwind region, Fig. 10   of about 10 km in the zonal direction and vertical velocities up to 2.5 m s −1 develop all over the island area.However, one main band at y ≈ 52 km remains persistent, even at higher altitudes.This updraft band is a result of the dynamic and thermal instability over the island, forming quasi twodimensional horizontal vortex rolls with their axes aligned in the downwind direction (e.g.Etling and Brown, 1993).Toward the evening, as the surface sensible heat flux is not positive anymore and convection fades away, the band decouples from the island and vanishes (not shown).Turbulent updraft cells within the marine boundary layer with vertical velocities between 0.5 and 1 m s −1 are also visible since a turbulent inflow is generated with the cell perturbation method described in Sect.2.5.In the INV case, these updrafts are a bit weaker, which is most likely due to the stronger mean horizontal wind speed.Wave-like structures in the upwind verti-cal velocity field are observed in the NOCP case.There, the flow remains laminar in this region, and since no perturbation is applied but surface fluxes are present, these artificial convergence lines are forming.Note that this effect is not seen in the REF and INV case.This underscores the importance of having an explicit inflow turbulence generation when working on LES scales.Just by visibly comparing the "coarse" simulation DX400 with the other cases, it becomes apparent that there is a lot of structure loss in the vertical wind field.All up-and downdraft bands -even the main updraft band downwind -are almost perfectly aligned in the x direction.This shows the importance of using a grid spacing that resolves the inertial subrange of the velocity spectrum (cf.Bryan et al., 2003).Note that with coarser grid spacings the orographical structures of the island are also less represented.Figure 11 shows the surface wind fields and liquid water path for all simulated cases at 14:00 LT.In all these cases, the island convection affects both the strength (up to 4 m s −1 stronger wind speeds compared to marine surface winds) and direction (±30 • ) of the wind in the downwind area of Barbados, thus leading to strong surface convergence and subsequently forming the updraft band as seen in Fig. 10 at y = 52 km.Despite having this elongated band, very little cloud formation is observed in this area for the REF, NOCP and INV cases, which is also the case for other times of the day (not shown).This means that no continuous cloud street is modeled on 22 and 27 June 2013, respectively.While cloud streets occur on around 60 % of undisturbed days, there are several effects that suppress cloud street generation (Kirsh-  Besides the distinct cloud bands, the DX400 case shows further very notable differences in the cloud field.First of all, clouds are broader because of the coarser grid spacing.In addition to that, a continuous cloud street is modeled, which can be considered an artifact since such a cloud band is not seen in either other simulations or satellite observations.Furthermore, the downwind horizontal velocity field is slightly stronger compared to the other cases.We attribute this behavior to the lack of resolved small scales that cannot extract energy from the large eddies and therefore grow and become more coherent.This effect is also observed to a lesser extent for the NOCP case. In the following, the diurnal development of the convective island boundary layer is investigated.Figure 12 3).
Figure 11.Horizontal cut planes (xy) of surface wind vectors and contours of liquid water path for all four simulation cases (see Table 3).where θ v0 is the virtual potential temperature at the surface.Being relatively similar in boundary layer characteristics and cloud cover, there is a clear distinction between REF and INV with regard to cloud microphysical properties.Due to higher cloud base height in the REF case, the cumulus clouds tend to grow deeper, which is seen in the mean liquid water path (LWP) values as well as in the cloud base and top heights (see Table 4).The cloud cover, however, is fairly comparable for these two cases.Peak updraft values also show a diurnal variation, starting from approximately w max = 1 m s −1 (which is equivalent to the upwind area value) up to w max = 7 m s −1 around noon.The DX400 case has overall weaker peak updrafts; however, this does not mean that there is less vertical transport of energy, moisture, momentum etc. Due to the coarser grid spacing in every spatial direction there is a higher net upward transport.
To further investigate daytime-dependent vertical mixing and layering, Fig. 13 shows hourly averaged vertical profiles boundary layer and cloud parameters for the downwind domain compared to the daily upwind average.Comparing again the REF and INV cases, daytime-dependent differences in the density potential temperature and specific humidity profiles can be noticed (not shown).Lower levels at z < 700 m are warmer and dryer compared to the marine background.The vertical turbulent transport is evidenced by the profiles of sensible and latent heat fluxes.The sensible heat flux is linearly decreasing within the mixing layer up to heights between 700 m < z < 900 m, depending on the time of day.The maximum latent heat fluxes occur between 600 m < z < 800 m.Above that layer, the cloud water content reaches its maximum, which is connected to latent heat release and thus to an increase in the sensible heat flux and a decrease in the latent heat flux.In the INV case, the trade inversion around z ≈ 2000 m inhibits further cloud development above this height, whereas in the REF simulation there is also a notable amount of cloud water above 2000 m.The presence of wind shear above 1500 m height leads to a secondary maximum of TKE around z = 2000 m, which is not the case in the shearless INV simulation.
In the NOCP case, i.e., without a turbulent inflow, persistent updraft bands form over the island area, which is consistent with the modeled cloud field from Fig. 11.The inflow characteristics have little effect on boundary layer properties like TKE, vertical velocity variance (not shown), sensible and latent heat flux (cf.Fig. 13).However, the values for the NOCP case tend to be a bit higher than in REF, especially between 500 and 700 m.One reason for this could be that the upwind marine boundary layer already vertically transported some amount of energy, which is missing in the NOCP case.More remarkable differences are noticeable with regard to cloud development.The LWC around z ≈ 1000 m during noon is nearly doubled for the NOCP case.There is also a particularly pronounced secondary maximum of LWC around z ≈ 2100 m in the same order of magnitude.With the average over the whole daytime period taken, the mean LWP is more than doubled in the NOCP case compared to the REF case, which is in agreement with higher LWC values and also higher cloud cover (≈ +3 %).
More undesired effects become apparent when using a coarser spatial resolution as in the DX400 case, which is most noticeable in the averaged vertical profiles.First of all, there is less variability in the potential temperature and specific humidity fields for altitudes z > 1000 m, which can be explained by the lack of turbulent vertical transport within the boundary layer (this effect can be seen in the vertical profile of TKE in Fig. 13 and in the profile of the vertical velocity variance).The LWC, however, has maximum val-Table 4. Diagnostics the four sensitivity simulations (see Table 3), including cloud cover, cloud base height z cb , cloud top height z ct , boundary layer height z i , water vapor path (WVP), liquid water path (LWP), and maximum updraft values w max .All quantities are spatially averaged for the downwind area and temporally averaged between 06:00 and 18:00 LT.

Case
Cloud cover ues of 0.028 g m −3 around z ≈ 2000 m, which is a factor of 2 higher compared to the REF case.This is accompanied by strong latent heat fluxes in these layers.Cloud growth is also more inhibited at finer resolutions due to explicit entrainment of dryer environmental air (e.g.Bryan et al., 2003).Having a distinct and quite symmetric diurnal variation in boundary layer and cloud properties in the other cases, the evening transition in the DX400 case is poorly represented, where still a notable number of clouds exist and a deeper boundary layer is modeled around 20:00 LT.

Vertical mixing of aerosols
After the long-range transport of Saharan dust into the Caribbean region, these dust layers arrive at Barbados with mean base heights of about 1.5-2 km a.s.l.Due to a possible interaction with the convective island boundary layer, vertical mixing of aerosols is investigated in this subsection.
As already shown in Sect.2, these aerosol layers are represented by passive tracers in the model.They are initialized with a relative concentration of 1 within the layer where the aerosol is detected and 0 otherwise.This approach has already been used for heat island effect studies in Engelmann et al. (2011).These relative concentrations can be related to mass concentrations of Saharan dust, e.g., 180 µg m −3 .This mass concentration and the Saharan dust layer heights are estimated from ground-based multi-wavelength aerosol lidar measurements and provide a rough idea of the magnitude of these quantities.Especially the particle depolarization ratio indicates that the pure dust layer begins at around 1.5 km altitude.A detailed analysis of dust layers during the campaign can be found in Groß et al. (2015).Figure 14 displays height-distance profiles of the boundary layer tracer φ BLT and the Saharan dust tracer φ SDT near the west coast of Barbados.In both cases, the turbulent character over the island section is visible as a vertical distribution of the passive tracer within the higher boundary layer and the corresponding decrease in tracer concentration.It is more pronounced for the southern part of the island, which is due to the broader land area width (20 km in the south compared to 10 km in the north).It also indicates the wind shear at the island boundaries (e.g., at y = 43 km), which causes the advection of air masses from z > 700 m into the boundary layer.This effect is more pronounced further west (not shown).In the REF case, the mean boundary layer height around noon was calculated to be z i ≈ 1400 m.The passive tracer analysis additionally shows some local overshoots at heights over 2 km a.s.l.The Saharan dust tracers do have a different vertical structure.
For the REF simulation, the tracer is thinned out, with maximum concentrations between 1.9 km < z < 2.4 km, whereas in the INV case it is between 1.7 km < z < 2.5 km.There are also no overshoots visible beyond z = 1.7 km.The stronger turbulent mixing in the REF case can be explained by the presence of wind shear around z ≈ 1.5 km height, whereas in the INV case the strong trade-wind inversion suppresses further development of turbulence in higher altitudes.The logarithmic scale indicates the tracer diffusion, which shows that about 1 % of the maximum concentration is present at altitudes between 1.3 and 1.5 km (depending on the case) and locally down to 1.0 km altitude due to the island effect.
Although there are already some indications of downward aerosol transport, a better quantification of these effects is still needed to achieve a better understanding of the processes behind it.For that reason, mean vertical profiles are calculated in the same manner as the boundary layer and cloud analysis.Figure 15 shows vertical profiles of tracerrelated quantities for the REF and the INV case.The tracer profile is being deformed by vertical transport processes, which can be originated by larger-scale or turbulent processes.To distinguish between those two, the total vertical tracer flux wφ SDT and the turbulent vertical tracer flux w φ SDT are computed.According to the model data, there is already a persistent downward movement east of Barbados between 1.5 km < z < 3.0 km, which corresponds to a subsidence velocity of w = −407 m day −1 for the REF case and w = −345 m day −1 for the INV case.There is a total relative downward flux for the downwind area in the REF case at 1.2 km height, with a maximum magnitude of wφ SDT = −0.03m s −1 .The height of this local extremum depends on the time of day, with higher altitudes in the morning and around noon and lower altitudes toward the evening.A net positive upward flux is always present above z = 2 km, becoming zero at z = 3 km.About one-third to one-half of the total downward flux is caused by turbulent mixing, having its local maxima at z = 1.8 km, which is the lower bound of the Saharan dust layer.The daytime turbulent downward mixing is about 5 times stronger in the INV case but is shifted approximately 400 m toward the surface compared to REF.Due to the tracer subsidence, the layer also reaches these altitudes of maximum turbulent downward mixing, which would not be the case if the dust layer remained at its initial height of 1.7 km.The downwind total tracer flux, however, is positive (upwards) for all altitudes and is almost 1 order of magnitude stronger compared to REF.Since the mean wind speed only slightly differs between the two cases and surface fluxes are the same, it can be concluded that atmospheric stability (especially inversion strength) and the presence of wind shear in the sub-inversion layer have a major impact on altitude and strength of layers with preferably upward or downward mixing properties.The model data suggest that the net effect, i.e., effects from both mean transport and turbulent mixing, is downward transport around the dust layer base for the downwind region in REF and a net upward transport in INV. because no nighttime radiosonde data were available on this particular day.Although the measurements and the simulation correspond to different days, the comparison of the radiosonde profiles used for the LES initialization launched on 27 June and the dropsonde measurements obtained during the measurement flight on 20 June show a good agreement in the altitude of the trade inversion, relative humidity levels and temperature profile (Fig. 16).Especially the measured horizontal wind speed profile from the dropsondes matches better to the 27 June than the 22 June simulation case.The airborne DWL used for this comparison was deployed onboard the DLR Falcon 20 research aircraft during the SALTRACE campaign.The system, based on an instrument developed by Lockheed Martin Coherent Technologies and enhanced by DLR to provide airborne measurement capabilities, can be operated in either nadir-pointing mode or scanning mode (Reitebuch, 2012;Chouza et al., 2015).The nadirpointing mode allows the retrieval of vertical wind speeds with a vertical resolution of 100 m and a horizontal resolution of approximately 200 m with a random error lower than 0.15 m s −1 and a systematic error lower than 0.05 m s −1 .Note that these resolutions are almost identical to the grid spacings used in the LES.The overflight took place between 10:36 and 10:44 LT at a flight altitude of 2900 m.
Figure 17 shows the flight track corresponding to the overflight (plane 2 in the figure) together with the LES results of the vertical wind speed for the 27 June case at 10:30 LT.For this overflight, vertical wind measurements on the lee side and over the island are available.
A comparison between the measured and the simulated vertical wind speed profiles is displayed in Fig. 18, where some main structures can be recognized in both profiles.Strong vertical winds associated with convective activity over Barbados can be observed in both the simulation and the measurements.In the case of the measurements, the presence of convective clouds limits the lidar coverage over the island (between 47 and 63 km on the x axis).For the measurements performed on the lee side of Barbados and for altitudes above 1 km a series of waves with a wavelength of approximately 10 km and an amplitude of 2 m s −1 can be recognized.A similar feature can be seen in the LES data, but with a slightly weaker amplitude.
To get a better idea of the horizontal distribution of these waves at a certain level, Fig. 19 shows the vertical wind speed at the trade-wind inversion height for different times during the day.It seems that the strength of daytime convective activity plays a minor role since the overall pattern looks very similar in all three snapshots, with the exception of minimally stronger fluctuations toward the afternoon hours.A marked wave structure in the lee of the island is visible, which is the result of trapped gravity waves due to the strong inversion.Again, a similar wavelength amplitude of about 10 km can be seen and stronger amplitudes (comparable to the DWL measurements) originating from the northern part of Barbados are also visible.
Figure 20 shows a comparison between the mean and the variance of the measured and simulated vertical wind speed on the lee side (between 0 and 35 km on the x axis) of Barbados.The mean vertical wind profiles of the measurements (black line) show a reasonable agreement with the LES results regarding the overall structure (e.g., cut planes 1, 2 and 5, which are located at the southern and northern edge of Barbados, covering a similar spatial area as in the flight route).The mean vertical velocity below 1 km height indicates a downdraft region on the lee side of Barbados.This can be explained by the daytime circulation pattern in the lee of Barbados since the updraft band located between the center and southern part of the island causes downdrafts at its lateral flanks.Planes 3 and 4 reflect this behavior with positive vertical velocities below 1 km height.The variance of the vertical wind field shows an overall increase with altitude both for measurement and model results.However, the strength of these fluctuations is increased in the measurement data, which could result from differences in the forcing, e.g., stronger surface fluxes than the prescribed ones in the LES model.by a second, weaker dust layer between 2.8 and 3.7 km which was more humid (40-50 % RH).Temperature inversions set the limits of the total dust layer or Saharan Air Layer.The lidar signal shows clouds all over the day around 1 km height and close to 1800 m.On average, 20 clouds per hour cross the lidar beam with an increasing number in the afternoon, leading to a cloud cover (in this context defined on a temporal scale) of roughly 50 %.At nighttime only a single cloud was detected within 3 h of measurement, which confirms the convective character of the daytime clouds at the investigated day over the west coast of Barbados.The bottom panel in Fig. 21 shows the LES model output of the corresponding heighttime profile of relative humidity for this day and the same location.The diurnal variation in convective vertical moisture transport is clearly pronounced.Clouds (white/red coloring) develop at the same altitudes as in the BERTHA lidar measurements.Even higher and deeper clouds up to the top of the trade inversion at z = 2 km are resolved.Higher cloud activity can be noticed in the measurements, whereas the LES model results show fewer clouds toward the late afternoon.This effect can be attributed to the relatively low ocean latent heat flux from the marine boundary layer forcing, which tends to dry out the lower troposphere over a longer period of time (i.e., toward the afternoon and evening hours).Additionally, there is a notable increase in moisture during the afternoon (see the 16:01 LT radiosonde launch compared to the other ones) within the first 800 m of the lower troposphere.

Multi-wavelength
Because of the fact that no additional large-scale forcing (in this case, advection of moist air) is applied during the simulation time, this effect expectedly cannot be captured by the LES.
Further comparisons regarding the cloud base height are conducted, revealing an even better agreement between the BERTHA lidar measurements and the model results.Figure 22 shows the temporal evolution of the cloud base height derived from the REF case LES output and the cloud base height derived from BERTHA lidar measurements during 22 June 2013.The cloud base was detected from the lidar signal with an accuracy of ±50 m.The algorithm sets the cloud base if the 52.5 m vertically smoothed 532 nm total signal increases by a factor of 2 within 50 m.To get an overview of the day, the temporal resolution was set to 30 s.This procedure reproduces the temporal cloud evolution in Fig. 22.The spread gives an idea about the cloud thickness due to the fact that if a cloud overpasses the lidar beam, it is very likely that cloud water near the cloud top is detected at the beginning and/or at the end of the overpass.However, in some cases, it indicates some single clouds that start at an altitude of 1600 m.At noon, the cloud base gets higher and the clouds are generally thinner, which can be seen the in range of variation in the lidar data, which is very low at this period of time.Overall, there is less variability in the LES data due to the spatial averaging over the downwind region.

Conclusions
We have presented a numerical study for investigating islandinduced effects on boundary layer structure, cloud generation and vertical mixing of aerosol layers at the easternmost Caribbean island of Barbados.The simulations were performed with the model ASAM on large eddy scale, where the horizontal resolution of x = y = 200 m is adequate to still resolve the lower wavenumber end of the inertial subrange of the turbulence spectrum.In order to generate inflow turbulence consistent with the upstream marine boundary layer forcing, the cell perturbation method based on finite amplitude perturbations was used.This method has been successfully adapted to moist boundary layer simulations with open lateral boundary conditions.Spectral analysis and examinations of vertical profiles of boundary layer quantities were used to determine the optimal Eckert number for the simulations.It was found that a value of Ec = 0.4 is most suitable for the island simulations, guaranteeing rapid development of turbulence that is very close to results from periodic BC simulations.This perturbation Eckert number differs from the optimum Ec = 0.2 derived for neutral stability (Muñoz-Esparza et al., 2015).These weaker-amplitude temperature perturbations appear to be due to a combination of the use of open boundary conditions and near production range grid resolution.
After the model was set up, several simulations were performed to analyze island effects.The atmospheric state is described by a horizontally homogeneous profile obtained via nighttime radiosonde launches.On the one hand, the used profile cannot be seen as a representative state, which could have been achieved by averaging multiple profiles.On the other hand, these quasi-idealized simulations en-abled the possibility to compare the LES output data with DWL measurements.Large-scale forcings are not applied during the simulation time, i.e., a time-invariant background state is used, which makes it easier to analyze daytime changes in boundary layer and cloud characteristics as a result of convective island activity.The radiosonde profile from 22 June 2013 served as a reference case for a crucial sensitivity study, which revealed the following: -Disadvantages of neglecting a turbulent inflow and modeling within the "terra incognita" or "gray zone" become apparent (Wyngaard, 2004).If the turbulent inflow generation is turned off, i.e., no marine boundary layer develops, cloud properties are drastically changed over and in the lee of the island.There, cloud cover, liquid water path and cloud base height have significantly higher values compared to the reference case.Also, the cloud morphology is different, having horizontally aligned cloud bands instead of scattered cumulus clouds.
-Using a generally coarser horizontal grid spacing of x = y = 400 m shows the same issues regarding cloud properties.Additionally, boundary layer turbulence is not well resolved but vertical energy transport is enhanced overall, which leads to the formation of a cloud street along the updraft band past Barbados.This has not been seen in the other cases and can thus be considered an artificial effect.
-For the two considered cases (22 and 27 June 2013), analysis of the daytime convective boundary layer shows generally similar vertical profiles for both, although slight differences were detected due to wind shear effects in particular.Also, for the case with the strong trade inversion, the boundary layer grew approximately 150 m deeper.Significant differences were again visible in cloud properties.In the latter case, the trade inversion inhibited cloud growth beyond 1800 m altitude, whereas local overshoots through the weaker inversion occurred in the 22 June case.
Vertical mixing of aerosols, in our case Saharan dust layers, has been qualitatively and quantitatively analyzed by in-  cluding passive tracers in the LES model.They are initialized at the same heights as they are seen in lidar data during the considered days.The model data suggested that a continuous subsidence velocity within these layers was present, which led to a mean sinking of 400 m toward the surface.Layers of turbulent downward mixing have also been detected between 1200 and 1700 m altitude.It can be concluded that if the dust layer reaches this height range, turbulent downward mixing of aerosol takes place, which is separated from large-scale subsidence effects.The exact position of these layers and the strength of downward turbulent mixing and transport are mainly controlled by atmospheric stability and wind shear.It became apparent that for the INV case the net direction of vertical transport is upwards for heights around the Saharan dust layer base.A larger number of simulation cases would provide further insight into this effect and the influence of trade inversion strength.
Qualitative comparisons with DWL measurements were carried out to validate the LES model results.A lot of similarities were found despite comparing different days, which, however, were very similar in terms of atmospheric stability, mean wind and temperature profile.First of all, large-scale changes in dynamics occur, which expectedly cannot be directly captured by the model.Furthermore, the vertical wind field shows a similar structure for both measurements and LES, highlighting the convective activity over the island and trapped gravity waves around the strong trade inversion.The amplitude of these waves was a bit higher in the measurements, which is also seen in the vertical wind variance profile.Taking these comparison results into account and con-necting them to the boundary layer and tracer analysis, it is possible that turbulent (downward) mixing of aerosol layers could be even more enhanced than the LES model results suggest.
From the stationary Raman lidar BERTHA, which was deployed near the west coast of Barbados, cloud base heights and thickness were estimated and compared with the LES data.With these two techniques combined, a consistent picture of the diurnal convective activity and cloud generation over the island was gained for the most part.
Possible future model development could focus on the direct and indirect aerosol effect as well as shadowing effects.This would lead to a better understanding of the effects of dust particles in connection with low-level clouds on their radiative feedback, e.g., as shown in Ge et al. (2014).
Furthermore, finer horizontal grid spacings than 200 m would be desirable to either confirm the robustness of the obtained results or reveal additional effects of finer resolved turbulent structures.This would lead to a better simulation of stably stratified areas, e.g., at the trade inversion.Also, it can be expected that, with a higher effective model grid resolution, the computed velocity variances would be closer to the measurement results of the Doppler lidar system.However, this requires a great deal of computational effort, especially when performing sensitivity tests with a larger number of simulations.
As a general conclusion of this work, the simulations performed provided a detailed image of downwind boundary layer structure, cloud and vertical mixing processes, which agree well with lidar measurements.The model data can also help to better interpret the ground-based observations gained during the SALTRACE campaign at the Barbados west coast (CIMH field site).

Figure 1 .
Figure 1.Left panel: part of the Caribbean including the Lesser Antilles.Right panel: topographical map of Barbados.The topographical data are obtained from the Consortium for Spatial Information (CGIAR-CSI) Shuttle Radar Topography Mission (SRTM) data set at 200 m resolution.The white star denotes the location of the CIMH, which is close to the measurement field site.

Figure 3 .
Figure 3. Parameterized (LES) and measured wind profiles from radiosondes on 22 June 2013 (left) and 27 June 2013 (right).Names indicate date and time in UTC.
is expected.Due to the vertical and temporal variability in the www.atmos-chem-phys.net/16

Figure 4 .
Figure 4. Parameterized diurnal variation in sensible (SHF) and latent heat fluxes (LHF) over island and ocean areas.

Figure 6 .
Figure 6.Spatial evolution of time-averaged energy spectra of u (left), v (middle) and w component (right) in the y direction at z/z i = 0.5 for three perturbation Eckert numbers and the NOCP case.Color lines indicate distance from the east lateral boundary in steps of 0.8 km, 51.2 : −0.8 : 1.6 km, from blue to red.The dashed black line corresponds to the reference spectrum from the periodic computation, additionally averaged in the streamwise direction.

Figure 7 .
Figure 7. Vertical profiles of horizontal wind speed U (top left panel), potential temperature θ (top middle), turbulent kinetic energy TKE (top right), momentum flux τ m (bottom left), sensible heat flux ρ 0 c pd w θ (bottom middle) and latent heat flux ρ 0 L 0 w q v (bottom right) at a downstream distance of 40 km from the east boundary for different Ec numbers and the NOCP case (averaged along the y direction and in time).The solid black line corresponds to the reference profile from the periodic computation, additionally averaged in the streamwise direction.The profiles are valid for the 22 June 2013 case.

Figure 8 .
Figure 8. Definitions of subdomains for spatial averaging to cover different boundary layer characteristics: upwind marine regime east of Barbados (right) and downwind island regime over the west coast area (left).

Figure 9 .
Figure 9. Perspective view of surface temperature and specific humidity field at the western boundary for the REF case.Clouds are visualized by 0.1 g kg −1 isosurface in white coloring.Red coloring depicts an isosurface of 2 m s −1 updrafts.Snapshot taken at 12:00 LT.A model area of 80 km × 60 km × 5 km is shown.
baum and Fairman, 2015).In the REF (and NOCP) case, the relatively low moisture load (RH = 80 % near the surface, decreasing below 60 % at z ≈ 1300 m) and a weak tradewind inversion leads to a suppression of the development of a cloud trail.Both moisture and stratification are increased in the INV case but the stronger mean trade winds (almost 12 m s −1 ) are the suppressing factor here(Kirshbaum and Fairman, 2015).Due to the absence of a turbulent inflow velocity field, the cumulus clouds over the island are horizontally aligned to the mean wind direction in the NOCP case.In the REF and INV cases, more realistic scattered cumulus cloud fields over the island area and downwind are modeled.

Figure 12 .
Figure 12.Time series of boundary layer height, maximum vertical velocity, cloud cover, cloud base height and liquid water path for the downwind domain.Spatial averaging as indicated in Fig. 8.

Figure 13 .
Figure 13.Vertical profiles of liquid water content LWC , resolved turbulent kinetic energy TKE , sensible and latent heat flux ρ 0 c pd w T and ρ 0 L 0 w q v for all four considered cases (REF, NOCP, DX400 and INV from top to bottom).Spatial averaging as indicated in Fig. 8. Black dashed lines represent the marine upwind area and are temporally averaged between 08:00 and 20:00 LT.Colored solid lines represent for hourly averages during different times of the day for the downwind area.

Figure 14 .
Figure 14.Meridional cut planes (y-z profiles, 3 km off the west coast of Barbados) of the relative passive tracer concentrations and relative humidity for the REF (top panel) and the INV case (bottom panel).Red/orange contour coloring represents the Saharan dust tracer concentration φ SDT and blue contour coloring represents the boundary layer tracer concentration φ BLT on a logarithmic scale.White circles denote the location of the northern and southern island edges.The snapshots are taken at 12:00 LT.

Figure 15 .Figure 16 .
Figure 15.Vertical profiles of the passive Saharan dust tracer φ SDT (left panel), the total vertical tracer flux wφ SDT (middle panel) and the resolved turbulent tracer flux w φ SDT (right panel) for the REF case (top row) and the INV case (bottom row).Solid line colors and spatial as well as temporal averaging as in Fig. 13.The horizontal brown dashed line represents the initial location of the Saharan dust tracer between 1800 m ≤ z ≤ 2800 m for REV and between 1700 m ≤ z ≤ 2900 m for INV.

Figure 17 .Figure 18 .
Figure17.Locations of the five considered vertical planes for the comparison between the LES results and the DWL measurements within a domain of 60 km (y) and 80 km (x) up to an altitude of 2.5 km (z).The vertical wind (in m s −1 ) along the planes is indicated in red (updrafts) and blue (downdrafts).Plane 2 is similar to the Falcon flight track.Planes 1 and 5 have no intersection with the island area, which is visualized by surface cut cells of the computational grid.Ocean area is in blue.

Figure 19 .
Figure 19.LES model output of vertical wind fields around the trade-wind inversion height at z = 1975 m at 10:00, 14:00 and 18:00 LT on 27 June 2013.The gray dashed line indicates the Falcon flight track (similar to cut plane 2 in Fig. 17).

Figure 20 .
Figure 19.LES model output of vertical wind fields around the trade-wind inversion height at z = 1975 m at 10:00, 14:00 and 18:00 LT on 27 June 2013.The gray dashed line indicates the Falcon flight track (similar to cut plane 2 in Fig. 17).

Figure 21 .Figure 22 .
Figure 21.Top panel: BERTHA lidar measurements at CIMH and three radiosonde profiles of 22 June 2013.The range-corrected signal of the 532 nm cross-polarized channel is shown.The aerosol layer (in yellow) reaches up to 3.7 km.The signal peaks (in red) are liquid clouds through which the lidar beam cannot penetrate, resulting in the lack of signal above the cloud.At 14:04 LT the temporal resolution was changed from 15 s to 3 s.Bottom panel: LES model results of relative humidity at the CIMH grid point for the REF case, corresponding to 22 June 2013.

Table 1 .
LES model configuration for the simulations performed in Sect.3.

Table 2 .
Parameter values of the cases to be examined: 22 June and 27 June 2013.

Table 3 .
Parameter choices for the sensitivity simulations performed.