Observations And Modelling Of The On-Ice And Off-Ice Air Flow Over The Northern Baltic Sea

Two cases of on-ice and off-ice air flow characterizing the opposite weather situations over the ice-edge zone in the northern Baltic Sea are analysed on the basis of aircraft observations, and modelled using atwo-dimensional mesoscale model. The stable boundary layer (SBL) during theon-ice flow exhibited little thermal modification, but a low-level jet (LLJ) was generated at the 250-m high top of the SBL. In the model, the LLJ was associated with inertial oscillations in space, while the baroclinicity explained the shape of the wind profile well above the SBL. Although the observed LLJ was most pronounced over the ice, the modelling suggests that it was not generated by the ice edge but by the coastline some 400 km upwind of the ice edge, where a much more drastic change in the thermal stratification and surface roughness took place. The generation, maintenance, and strength of the LLJ were very sensitive to the parameterization of turbulent mixing in the SBL. In the case of the off-ice flow, the modification of the air mass and the development of a convective boundary layer (CBL) both over the ice and open sea were reasonably well modelled. Sensitivity runs suggested that it was essential to take into account the effects of subgrid-scale leads, a forest in the archipelago (which was crossed by the air flow), and water vapour condensationinto ice crystals. The heat flux from leads was particularly important for the heatbudget of the CBL, and the observed growth of the CBL was partly due to theeffective mixing over the rough and relatively warm forest.     

The KABEG’97 field experiment: An aircraft-based study of katabatic wind dynamics over the Greenland ice sheet

The aircraft-based experiment KABEG97 (Katabatic wind and boundary-layer front experiment around Greenland) was performed in April/May 1997. During the experiment, surface stations were installed at five positions on the ice sheet and in the tundra near Kangerlussuaq, West Greenland. A total of nine katabatic wind flights were performed during quite different synoptic situations and surface conditions, and low-level jets with wind speeds up to 25m s^-1 were measured under strong synoptic forcing of the katabatic wind system. The KABEG data represent a unique data set for the investigation of katabatic winds. For the first time, high-resolution and accurate aircraft measurements can be used to investigate the three-dimensional structure of the katabatic wind system for a variety of synoptic situations.Surface station data show that a pronounced daily cycle of the near-surface wind is present for almost all days due to the nighttime development of the katabatic wind. In a detailed case study the stably-stratified boundary layer over the ice and the complex boundary-layer structure in the transition zone ice/tundra are investigated. The katabatic wind system is found to extend about 10 km over the tundra area and is associated with strong wind convergence and gravity waves. The investigation of the boundary-layer dynamics using the concept of a two-layer katabatic wind model yields the results that the katabatic flow is always a shooting flow and that the pure katabatic force is the main driving mechanism for the flow regime, although a considerable influence of the large-scale synoptic forcing is found as well.     

The Coupling State of an Idealized Stable Boundary Layer

The coupling state between the surface and the top of the stable boundary layer (SBL) is investigated using four different schemes to represent the turbulent exchange. An idealized SBL is assumed, with fixed wind speed and temperature at its top. At the surface, two cases are considered, first a constant temperature, 20 K lower than the SBL top, and later a constant 2 K h⁻¹ cooling rate is assumed for 10 h after a neutral initial condition. The idealized conditions have been chosen to isolate the influence of the turbulence formulations on the coupling state, and the intense stratification has the purpose of enhancing such a response. The formulations compared are those that solve a prognostic equation for turbulent kinetic energy (TKE) and those that directly prescribe turbulence intensity as a function of atmospheric stability. Two TKE formulations are considered, with and without a dependence of the exchange coefficients on stability, while short and long tail stability functions (SFs) are also compared. In each case, the dependence on the wind speed at the SBL top is considered and it is shown that, for all formulations, the SBL experiences a transition from a decoupled state to a coupled state at an intermediate value of mechanical forcing. The vertical profiles of potential temperature, wind speed and turbulence intensity are shown as a function of the wind speed at the SBL top, both for the decoupled and coupled states. The formulation influence on the coupling state is analyzed and it is concluded that, in general, the simple TKE formulation has a better response, although it also tends to overestimate turbulent mixing. The consequences are discussed.     

Eddy Covariance Tilt Corrections over a Coastal Mountain Area in South-east China: Significance for Near-Surface Turbulence Characteristics

Turbulence characteristics of an atmospheric surface layer over a coastal mountain area were investigated under different coordinate frames.Performances of three methods of coordinate rotation:double rotation(DR),triple rotation(TR),and classic planar-fit rotation(PF) were examined in terms of correction of eddy covariance flux.Using the commonly used DR and TR methods,unreasonable rotation angles are encountered at low wind speeds and cause significant run-to-run errors of some turbulence characteristics.The PF method rotates the coordinate system to an ensemble-averaged plane,and shows large tilt error due to an inaccurate fit plane over variable terrain slopes.In this paper,we propose another coordinate rotation scheme.The observational data were separated into two groups according to wind direction.The PF method was adapted to find an ensemble-averaged streamline plane for each group of hourly runs with wind speed exceeding 1.0 m s 1.Then,the coordinate systems were rotated to their respective bestfit planes for all available hourly observations.We call this the PF10 method.The implications of tilt corrections for the turbulence characteristics are discussed with a focus on integral turbulence characteristics,the spectra of wind-velocity components,and sensible heat and momentum fluxes under various atmospheric stabilities.Our results show that the adapted application of PF provides greatly improved estimates of integral turbulence characteristics in complex terrain and maintains data quality.The comparisons of the sensible heat fluxes for four coordinate rotation methods to fluxes before correction indicate that the PF10 scheme is the best to preserve consistency between fluxes.     

Reynolds-Number Dependence of Gas Dispersion Over a Wavy Wall

Direct numerical simulations of an Ekman layer are performed to study flow evolution during the response of an initially neutral boundary layer to stable stratification. The Obukhov length, L, is varied among cases by imposing a range of stable buoyancy fluxes at the surface to mimic ground cooling. The imposition of constant surface buoyancy flux , i.e. constant-flux stability, leads to a buoyancy difference between the ground and background that tends to increase with time, unlike the constant-temperature stability case where a constant surface temperature is imposed. The initial collapse of turbulence in the surface layer owing to surface cooling that occurs over a time scale proportional to \(L/u_*\), where \(u_*\) is the friction velocity, is followed by turbulence recovery. The flow accelerates, and a “low-level jet” (LLJ) with inertial oscillations forms during the turbulence collapse. Turbulence statistics and budgets are examined to understand the recovery of turbulence. Vertical turbulence exchange, primarily by pressure transport, is found to initiate fluctuations in the surface layer and there is rebirth of turbulence through enhanced turbulence production as the LLJ shear increases. The turbulence recovery is not monotonic and exhibits temporal intermittency with several collapse/rebirth episodes. The boundary layer adjusts to an increase in the surface buoyancy flux by increased super-geostrophic velocity and surface stress such that the Obukhov length becomes similar among the cases and sufficiently large to allow fluctuations with sustained momentum and heat fluxes. The eventual state of fluctuations, achieved after about two inertial periods (\(ft \approx 4\pi \)), corresponds to global intermittency with turbulent patches in an otherwise quiescent background. Our simplified configuration is sufficient to identify turbulence collapse and rebirth, global and temporal intermittency, as well as formation of low-level jets, as in observations of the stratified atmospheric boundary layer.     

Large‐eddy simulation of small‐scale surface effects on the convective boundary‐layer structure

Abstract

The effects of small‐scale surface inhomogeneities on the turbulence structure in the convective boundary layer are investigated using a high‐resolution large‐eddy simulation model. Surface heat flux variations are sinusoidal and two‐dimensional, dividing the total domain into a checkerboard‐like pattern of surface hot spots with a 500‐m wavelength in the x and y directions, or 1/4 of the domain size. The selected wind speeds were 1 and 4 m s^‐l, respectively. As a comparison, a simulation of the turbulence structure was performed over a homogeneous surface.

When the wind speed is light, surface heat flux variations influence the horizontally averaged turbulence statistics, including the higher moments despite the small characteristic length of the surface perturbation. Stronger mean wind speeds weaken the effects of inhomogeneous surface conditions on the turbulence structure in the convective boundary layer.

Results from conditional sampling show that when the mean wind speed is small, weak mean circulations occur, with updraft branches above the high heat flux regions and down‐draft branches above the low heat flux regions. The inhomogeneous surface induces significant differences in the turbulence statistics between the high and low heat flux regions. However, the effect of the surface perturbations weaken rapidly when the mean wind speed increases. This research has implications in the explanation of the large‐scale variability commonly encountered in aircraft observations of atmospheric turbulence.     

Characteristics of intermittent turbulence in the upper stable boundary layer over Greenland

The experiment IGLOS (Investigation of the Greenland Boundary Layer Over Summit) was conducted in June and July 2002 in the central plateau of the Greenland inland ice. The German research aircraft Polar2, equipped with the turbulence measurement system Meteopod, was used to investigate turbulence and radiation flux profiles near research station “Summit Camp”. Aircraft measurements are combined with measurements of radiation fluxes and turbulent quantities made from a 50 m tower at Summit Camp operated by Eidgenössische Technische Hochschule (ETH) Zürich. During all six flight missions, well-developed stable boundary layers were found. Even in high-wind conditions, the surface inversion thickness did not exceed roughly 100 m. The turbulent height of the stable boundary layer (SBL) was found to be much smaller than the surface inversion thickness. Above the surface layer, significant turbulent fluxes occurred only intermittently in intervals on the order of a few kilometres. Turbulent event fraction in the upper SBL shows the same dependence on gradient Richardson number as reported for near-surface measurements. Clear-air longwave radiation divergence was always found to contribute significantly to the SBL heat budget. In low-wind cases, radiative cooling even turned out to be dominant.     

A study of the planetary boundary layer over the polynya downwind of St. Lawrence Island in the Bering Sea using aircraft data

Turbulence data collected with the gust probe system on the NOAA P-3 aircraft over the polynya downwind of St. Lawrence Island in the Bering Sea are used to study the fluxes of heat, momentum, and moisture from the polynya. The data also allow study of the effect of the topography of St. Lawrence Island on the atmospheric boundary-layer flow over the polynya and ultimately on ice production in the polynya. Two cases are studied: one (Feb. 15, 1982) where the topographic effects are minimal and the other (Feb. 18, 1983) where the topographic effects are dominant. Calculation of the surface drag coefficient, C _D, for the Feb. 15, 1982 case over young grey/white ice gave a value of 1.2 × 10^-3, which is in close agreement with previous results. The value of the drag coefficient for the grey/white ice regime on Feb. 18, 1983, where the upstream topography on St. Lawrence Island had an important influence on the flow over the polynya, was 3.2 × 10^-3. It was determined that this higher value was related to the more efficient mixing of momentum downward by turbulent eddies generated by flow over and around the topography. The area-averaged heat transfer coefficient, C _H, over the polynya was on the order of 1.1 × 10^-3 for both days, but there were large variations in heat flux across the polynya due to variations in the flow caused by the topography. Conditional sampling techniques applied to the turbulence data showed that the fractional areas occupied by updrafts and downdrafts were 28% and 36%, respectively, and that these results were within the range of values found in previous studies for over-land and over-ocean conditions.     

Nocturnal CO2 exchange over a tall forest canopy associated with intermittent low-level jet activity

Summary Nocturnal eddy-covariance carbon dioxide fluxes have uncertainties arising from non-stationary atmospheric processes. Low-level jets (LLJ) are one of the prominent nocturnal boundary-layer phenomena observed over non-mountainous terrain, and are capable of generating shear and turbulence close to the ground. The influence of intermittent LLJ activity on nocturnal carbon dioxide exchange measurements is investigated using wind profile observations and eddy-covariance flux measurements over a tall forest canopy. Results suggest that the buildup and venting of CO₂ are closely associated with LLJ activity during the night. Of significance in quantifying nocturnal fluxes, this paper demonstrates how low-level jet activity introduces sporadic coupling between the canopy and the atmosphere.     

Impact of the Nocturnal Low-Level Jet and Orographic Waves on Turbulent Motions and Energy Fluxes in the Lower Atmospheric Boundary Layer

The nocturnal low-level jet (LLJ) and orographic (gravity) waves play an important role in the generation of turbulence and pollutant dispersion and can affect the energy production by wind turbines. Additionally, gravity waves have an influence on the local mixing and turbulence within the surface layer and the vertical flux of mass into the lower atmosphere. On 25 September 2017, during a field campaign, a persistent easterly LLJ and gravity waves were observed simultaneously in a coastal area in the north of France. We explore the variability of the wind speed, turbulent eddies, and turbulence kinetic energy in the time–frequency and space domain using an ultrasonic anemometer and a scanning wind lidar. The results reveal a significant enhancement of the turbulence-kinetic-energy dissipation (by?50%) due to gravity waves in the LLJ shear layer (below the jet core) during the period of wave propagation. Large magnitudes of zonal and vertical components of the shear stress (approximately 0.4 and 1.5 m² s^?2, respectively) are found during that period. Large eddies (scales of 110 to 280 m) matching the high-wind-speed regime are found to propagate the momentum downwards, which enhances the mass transport from the LLJ shear layer to the roughness layer. Furthermore, these large-scale eddies are associated with the crests while comparatively small-scale eddies are associated with the troughs of the gravity wave.

     

Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

Boundary-layer structure near the cold front of a marine cyclone during “ERICA”

The boundary layer in the warm sector of a moderately deepening winter cyclone during the Experiment on Rapidly Intensifying Cyclones over the Atlantic (ERICA) is studied near the cold front. Data from the National Center for Atmospheric Research Electra research aircraft are used to examine mean and turbulence quantities. The aircraft data and supplemental data from ships, drifting buoys and moored buoys reveal an equivalent-barotropic pressure field. The area is found to be dominated by gradients in temperature and in turbulent fluxes, with changes occurring over 100 km horizontally being comparable to changes over 350 m vertically. The horizontal components of the gradients are found to be a maximum in a direction perpendicular to the front. Cross-sections perpendicular to the front are used to illustrate boundary-layer structure. Profiles of wind speed, stress, wind direction and stress direction are estimated from an Ekman model that is modified to take into account the equivalent-barotropic pressure field. Comparison of profiles from the model to the aircraft-measured data show reasonable agreement far from the front (100 km) when the model uses a constant eddy viscosity of approximately 6 kg m^–1 s^–1. Near the front there is less agreement with the model. Profiles of turbulent fluxes of momentum, heat and latent heat are divergent, with along-wind momentum flux negative and decreasing upward, cross-wind momentum flux positive and increasing upward, and heat flux and latent heat flux small, positive and decreasing upward. Far from the front, the turbulent kinetic energy budget shows that dissipation balances shear production. However, near-front behavior has an imbalance at low altitude, with shear production appearing as a TKE sink.     

Nocturnal Low-Level Jet Characteristics Over Kansas During Cases-99

Characteristics and evolution of the low-level jet (LLJ)over southeastern Kansas were investigated during the 1999 Cooperative Surface-AtmosphereExchange Study (CASES–99) field campaign with an instrument complement consisting of ahigh-resolution Doppler lidar (HRDL), a 60 m instrumented tower, and a triangle of Dopplermini-sodar/profiler combinations. Using this collection of instrumentation we determined thespeed U_X, height Z_X and direction D_X of the LLJ. We investigate here the frequencyof occurrence, the spatial distribution, and the evolution through the night, of these LLJcharacteristics. The jet of interest in this study was that which generates the shear and turbulencebelow the jet and near the surface. This was represented by the lowest wind maximum.We found that this wind maximum, which was most often between 7 and 10 m s¹,was often at or just below 100 m above ground level as measured by HRDL at the CASEScentral site. Over the 60 km profiler–sodararray, the topography varied by 100 m. The wind speed anddirection were relatively constant over this distance (with some tendency for strongerwinds at the highest site), but Z_X was more variable. Z_X was occasionally about equal at allthree sites, indicating that the jet was following the terrain, but more often it seemed to berelatively level, i.e., at about the same height above sea level. Z_X was also more variable thanU_X in the behaviour of the LLJ with time through the night, and on some nights $U_X wasremarkably steady. Examples of two nights with strong turbulence below jet level were furtherinvestigated using the 60 m tower at the main CASES–99 site. Evidence of TKE increasing withheight and downward turbulent transport of TKE indicates that turbulence was primarilygenerated aloft and mixed downward, supporting the upside–down boundary layer notion in thestable boundary layer.     

Boundary-layer modification in wintertime cold-air outbreaks from the Arctic sea ice

During the field experiment ARKTIS 1993 ten cases of boundary-layer modification in wintertime cold-air outbreaks from the Arctic sea ice in the Spitsbergen region were observed by aircraft over a distance ranging from about 50 km over the ice to about 300 km over the water. The modification depends decisively on the initial conditions over the ice, the boundary conditions at the bottom and top of the boundary layer and on the conditions of the large-scale flow. The modification of the bulk boundary-layer characteristics in relation to these conditions is presented.Besides the air-sea temperature contrast, the most important role for the boundary-layer modification is played by the stability on top of the boundary layer and by the divergence of the large-scale flow. According to the high variability of these conditions the observed boundary-layer modifications were very variable ranging from 100 to 300 m thick boundary layers with air temperatures between -32 and -22 °C over the ice to thicknesses between 900 and 2200 m and air temperatures between -15 and -5 °C after 300 km fetch over the open water. In most cases the large-scale flow was anticyclonic and divergent over the ice and changed to cyclonic and convergent over the water and an ice-sea breeze was superimposed on it.The sensible and latent heat fluxes are the dominant terms in the surface energy budget over the open water and ranged between 200 and 700 W m^-2 whereas the net longwave radiation is the dominating term over the ice with the heat fluxes only about 10 W m^-2.     

Present and future near-surface wind climate of Greenland from high resolution regional climate modelling

The present and twenty-first century near-surface wind climate of Greenland is presented using output from the regional atmospheric climate model RACMO2. The modelled wind variability and wind distribution compare favourably to observations from three automatic weather stations in the ablation zone of southwest Greenland. The Weibull shape parameter is used to classify the wind climate. High values (κ > 4) are found in northern Greenland, indicative of uniform winds and a dominant katabatic forcing, while lower values (κ < 3) are found over the ocean and southern Greenland, where the synoptic forcing dominates. Very high values of the shape parameter are found over concave topography where confluence strengthens the katabatic circulation, while very low values are found in a narrow band along the coast due to barrier winds. To simulate the future (2081–2098) wind climate RACMO2 was forced with the HadGEM2-ES general circulation model using a scenario of mid-range radiative forcing of +4.5 W m^?2 by 2100. For the future simulated climate, the near-surface potential temperature deficit reduces in all seasons in regions where the surface temperature is below the freezing point, indicating a reduction in strength of the near-surface temperature inversion layer. This leads to a wind speed reduction over the central ice sheet where katabatic forcing dominates, and a wind speed increase over steep coastal topography due to counteracting effects of thermal and katabatic forcing. Thermally forced winds over the seasonally sea ice covered region of the Greenland Sea are reduced by up to 2.5 m s^?1.     

Structure and Dynamics of Katabatic Flow Jumps: Idealised Simulations

For the first time, results from a high-resolution numerical simulation (with horizontal grid spacing of 35m) were used to reveal the detailed structure near an atmospheric katabatic jump over an idealized slope. The simulation represents flow over the slopes of Coats Land, Antarctica for austral winter conditions. The katabatic jump is characterised by an updraft with vertical velocities of order 1ms⁻¹ and serves as a possible forcing mechanism for the gravity waves frequently observed over the ice shelves around the Antarctic. Results also indicate that strong turbulence is generally confined within a mixing zone near the top of the katabatic layer upstream of the jump and extends downstream through the top of the strong updraft associated with the jump. Detailed analyses of momentum and heat budgets across the katabatic jump indicate that, upstream of the jump, turbulent mixing is important in decelerating the upper part of the katabatic layer, while within the jump the upslope pressure gradient force associated with the pool of cold air plays a role in decelerating the flow near the surface. The heat budget near the jump reveals a simple two-term balance: the turbulent heat flux divergence is balanced by the advection. A comparison of model results with available theories indicates that mixing between layers of different potential temperature structure indeed plays some role in the development of katabatic flow jumps, especially for strong jumps. Theories used to study katabatic jumps should include this mixing process, of which the amount depends on the intensity of the jump. A conceptual model of a katabatic jump, including the main dynamical processes, is constructed from these detailed analyses.     

The Atmospheric Boundary Layer In An Arctic Wintertime On-Ice Air Flow

A warm on-ice air flow from the open water over the Arctic sea ice in the Fram Straitwas, for the first time, systematically measured on 12 March 1998 by aircraft in thelowest 3 km over a 300-km long distance. The air mass modification and the processesinvolved are discussed.Over the water, air temperature was lower than water temperature so that a convectiveboundary layer (CBL) was present as initial condition. As soon as the CBL passed theice edge, a shallow stable internal boundary layer (IBL) was formed. In the residual CBL, turbulence and pre-existing convective clouds dissolved within about 20 km. Within about the same distance, due to the transition from unstable to stable stratification, the influence of surface friction increased in the IBL and decreased above the IBL with consequent generation of a low-level jet at IBL top. The IBL was strongly stratified with respect to both temperature and wind. The wind shear was around 0.1 s^-1 so that the Richardson number in the IBL was subcritical and turbulence was generated. The IBL top grew to about 145 m over 230 km distance. The growth of the IBL was not monotonic and was influenced by (a) inhomogeneous ice surface temperatures causedby both different ice thickness and changes in the cloud conditions, and (b) leads in theice deck. At the front side of the on-ice flow, the air mass boundary between the warmair and the cold Arctic air was sharp (12 K over 10 km) at low levels and tilted withheight. Observations suggest that the stratified IBL was lifted as a slab on top of thecold air.     

Mesoscale modeling of katabatic winds over Greenland and comparisons with AWS and aircraft data

Summary Simulations of the katabatic wind system over the Greenland ice sheet for the two months April and May 1997 were performed using the Norwegian Limited Area Model (NORLAM) with a horizontal resolution of 25 km. The model results are intercompared and validated against observational data from automatic weather stations (AWS), global atmospheric analyses and instrumented aircraft observations of individual cases during that period. The NORLAM is able to simulate the synoptic developments and daily cycle of the katabatic wind system realistically. For most of the cases covered by aircraft observations, the model results agree very well with the measured developments and structures of the katabatic wind system in the lowest 400 m. Despite NORLAM’s general ability of reproducing the four-dimensional structure of the katabatic wind, problems occur in cases, when the synoptic background is not well captured by the analyses used as initial and boundary conditions for the model runs or where NORLAM fails to correctly predict the synoptic development. The katabatic wind intensity in the stable boundary layer is underestimated by the model in cases when the simulated synoptic forcing is too weak. An additional problem becomes obvious in cases when the model simulates clouds in contrast to the observations or when the simulated clouds are too thick compared to the observed cloud cover. In these cases, the excessive cloud amount prevents development of the katabatic wind in the model. Received September 22, 2000/Revised March 16, 2001     

Large-Eddy Simulation of the Daytime Boundary Layer in an Idealized Valley Using the Weather Research and Forecasting Numerical Model

A three-dimensional numerical meteorological model is used to perform large-eddy simulations of the upslope flow circulation over a periodic ridge-valley terrain. The subgrid-scale quantities are modelled using a prognostic turbulence kinetic energy (TKE) scheme, with a grid that has a constant horizontal resolution of 50 m and is stretched along the vertical direction. To account for the grid anisotropy, a modified subgrid length scale is used. To allow for the response of the surface fluxes to the valley-flow circulation, the soil surface temperature is imposed and the surface heat and momentum fluxes are computed based on Monin–Obukhov similarity theory. The model is designed with a symmetrical geometry using periodic boundary conditions in both the x and y directions. Two cases are simulated to study the influence of along-valley geostrophic wind forcing with different intensities. The presence of the orography introduces numerous complexities both in the mean properties of the flow and in the turbulent features, even for the idealized symmetric geometry. Classical definitions for the height of the planetary boundary layer (PBL) are revisited and redefined to capture the complex structure of the boundary layer. Analysis of first- and second-moment statistics, along with TKE budget, highlights the different structure of the PBL at different regions of the domain.