On the consideration of mesoscale transports in climate modelling

Summary The dynamical effect of land surface heterogeneity on heat fluxes in the atmospheric boundary layer (ABL) is investigated using numerical simulations with a non-hydrostatic model over a wide range of grid resolutions. It is commonly assumed that mesoscale or dynamical fluxes associated with mesoscale and convective circulations simulated by a high-resolution model (subgrid (SG) model) on the subgrid scale of a climate model (large-scale (LS) model) represent additional processes in the ABL, which are not considered by the turbulence scheme of the LS-model, and which can be parameterized using the SG-model. The present study investigates the usefulness of this methodology for small-scale and large-scale idealized heterogeneities using a SG-model resolving mesoscale or even microscale circulations to compute the mesoscale fluxes on the scale of the LS-model. It is shown that the dynamical transports as derived from the SG-model should not be used to correct the parameterized turbulent fluxes of the LS-model. The reason is that the subgrid circulations simulated by the SG-model interact with the fields of wind and scalars in the ABL, which results in reduced turbulent fluxes in the ABL. Thus the methodology of previous studies to use mesoscale/dynamical fluxes for the correction of flux profiles simulated by climate models seems to be questionable.     

Validation of the “Lokal-Modell” over heterogeneous land surfaces using aircraft-based measurements of the REEEFA experiment and comparison with micro-scale simulations

Summary An aircraft-based experimental investigation of the atmospheric boundary layer (ABL) structure and of the energy exchange processes over heterogeneous land surfaces is presented. The measurements are used for the validation of the mesoscale atmospheric model “Lokal-Modell” (LM) of the German Weather Service with 2.8 km resolution. In addition, high-resolution simulations using the non-hydrostatic model FOOT3DK with 250 m resolution are performed in order to resolve detailed surface heterogeneities. Two special observation periods in May 1999 show comparable convective boundary layer (CBL) conditions. For one case study vertical profiles and area averages of meteorological quantities and energy fluxes are investigated in detail. The measured net radiation is highly dependent on surface albedo, and the latent heat flux exhibits a strong temporal variability in the investigation area. A reduction of this variability is possible by aggregation of multiple flight patterns. To calculate surface fluxes from aircraft measurements at low altitude, turbulent energy fluxes were extrapolated to the ground by the budget method, which turned out to be well applicable for the sensible heat flux, but not for the latent flux. The development of the ABL is well captured by the LM simulation. The comparison of spatiotemporal averages shows an underestimation of the observed net radiation, which is mainly caused by thin low-level clouds in the LM compared to observed scattered CBL clouds. The sensible heat flux is reproduced very well, while the latent flux is highly overestimated especially above forests. The realistic representation of surface heterogeneities in the investigation area in the FOOT3DK simulations leads to improvements for the energy fluxes, but an overestimation of the latent heat flux still persists. A study of upscaling effects yields more structures than the LM fields when averaged to the same scale, which are partly caused by the non-linear effects of parameter aggregation on the LM scale.     

Impact of coherent eddies on airborne measurements of vertical turbulent fluxes

During the Hydrological-Atmospheric Pilot Experiment (HAPEX)-Sahel, which took place in Niger in the transitional period between the wet and dry seasons, two French aircraft probed the Sahelian boundary layer to measure sensible and latent heat fluxes. The measurements over the Niamey area often revealed organised structures of a few km scale that were associated with both thermals and dry intrusions. We study the impact of these coherent structures using a single day’s aircraft-measured fluxes and a numerical simulation of that day with a mesoscale model. The numerical simulation at high horizontal resolution (250 m) contains structures that evolve from streaks in the early morning to cells by noon. This simulation shows distribution, variance and skewness similar to the observations. In particular, the numerical simulation shows dry intrusions that can penetrate deeply into the atmospheric boundary layer (ABL), and even reach the surface in some cases, which is in accordance with the observed highly negatively skewed water vapour fluctuations. Dry intrusions and thermals organised at a few km scale give skewed flux statistics and can introduce large errors in measured fluxes. We use the numerical simulation to: (i) evaluate the contribution of the organised structures to the total flux, and (ii) estimate the impact of the organised structures on the systematic and random errors resulting from the 1D sampling of the aircraft as opposed to the 2D numerical simulation estimate. We find a significant contribution by the organised structures to the total resolved fluxes. When rolls occur, and for a leg length of about 30 times the ABL depth, the 1D sampled flux is shown to be sometimes 20% lower than the corresponding 2D flux when the 1D sampling direction is the same as the main axis of the rolls, whereas the systematic error is much lower when the direction of the leg is transverse to the rolls. In the case of cells, an underestimate of around 10% can still be observed with the 1D approach independent of direction, due to poor sampling of the energy-containing scales.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

On the Development of a Dynamic Non-linear Closure for Large-Eddy Simulation of the Atmospheric Boundary Layer

A dynamic procedure is developed to compute the model coefficients in the recently introduced modulated gradient models for both momentum and scalar fluxes. The magnitudes of the subgrid-scale (SGS) stress and the SGS flux are estimated using the local equilibrium hypothesis, and their structures (relative magnitude of each of the components) are given by the normalized gradient terms, which are derived from the Taylor expansion of the exact SGS stress/flux. Previously, the two model coefficients have been specified on the basis of theoretical arguments. Here, we develop a dynamic SGS procedure, wherein the model coefficients are computed dynamically according to the statistics of the resolved turbulence, rather than provided a priori or ad hoc. Results show that the two dynamically calculated coefficients have median values that are approximately constant throughout the turbulent atmospheric boundary layer (ABL), and their fluctuations follow a near log-normal distribution. These findings are consistent with the fact that, unlike eddy-viscosity/diffusivity models, modulated gradient models have been found to yield satisfactory results even with constant model coefficients. Results from large-eddy simulations of a neutral ABL and a stable ABL using the new closure show good agreement with reference results, including well-established theoretical predictions. For instance, the closure delivers the expected surface-layer similarity profiles and power-law scaling of the power spectra of velocity and scalar fluctuations. Further, the Lagrangian version of the model is tested in the neutral ABL case, and gives satisfactory results.     

Spatial and Temporal Characteristics of Horizontal Rolls and Cells in the Atmospheric Boundary Layer Based on Radar and in Situ Observations

This paper is the second of a series devoted to the observation and analysis of coherent structures in the cloudy Atmospheric Boundary Layer (ABL) such as horizontal rolls or thermal cells. In the first paper, the TRAC (Turbulence Radar Aircraft Cells) experiment which is the observational support of this investigation based on coupled radar-aircraft measurements, was presented along with an overview of the main results of this campaign held in June 1993 in France. Here the analysis is focused on the spatial characteristics (length-scale, orientation ... ) of the coherent structures, their temporal evolution and vertical distribution deduced from the radar reflectivity fields acquired in clear air at several levels for five different ABLs. For that, an original and efficient image processing method able to extract the major mode of the organisation was developed. These characteristics are examined in relation to the dynamic and thermodynamic state of the ABL using mean and turbulent information as observed by the in-situ aircraft. These experimental results are the basis of a future Large Eddy Simulation modeling of an organised ABL which is in progress and will be the concern of the third paper in the series.     

Large-Scale Turbulence Structures and Their Contributions to the Momentum Flux and Turbulence in the Near-Neutral Atmospheric Boundary Layer Observed from a 213-m Tall Meteorological Tower

Large-scale turbulence structures in the near-neutral atmospheric boundary layer (ABL) are investigated on the basis of observations made from the 213-m tall meteorological tower at Tsukuba, Japan. Vertical profiles of wind speed and turbulent fluxes in the ABL were obtained with sonic anemometer-thermometers at six levels of the tower. From the archived data, 31 near-neutral cases are selected for the analysis of turbulence structures. For the typical case, event detection by the integral wavelet transform with a large time scale (180 s) from the streamwise velocity component (u) at the highest level (200 m) reveals a descending high-speed structure with a time scale of approximately 100 s (a spatial scale of 1 km at the 200-m height). By applying the wavelet transform to the u velocity component at each level, the intermittent appearance of large-scale high-speed structures extending also in the vertical is detected. These structures usually make a large contribution to the downward momentum transfer and induce the enhancement of turbulent kinetic energy. This behaviour is like that of “active” turbulent motions. From the analysis of the two-point space–time correlation of wavelet coefficients for the u velocity component, the vertical extent and the downward influence of large-scale structures are examined. Large fluctuations in the large-scale range (wavelet variance at the selected time scale) at the 200-m level tend to induce the large correlation between the higher and lower levels.     

Effects Of Changing Surface Heat Flux On Atmospheric Boundary-Layer Flow Over Flat Terrain

We examine the unsteady response of a neutral atmospheric boundary layer (ABL) of depth h and friction velocity u _* when a uniform surface heat flux is applied abruptly or decreased rapidly over a time scale t<inf>θ</inf> less than about h /(10u _*). Standard Monin–Obukhov (MO) relationships are used for the perturbed eddy viscosity profile in terms of the changes to the heat flux and mean shear. Analytical solutions for changes in temperature, mean wind and shear stress profile are obtained for the surface layer, when there are small changes in h /|L_MO| over the time scale t_MO~|L _MO|/(10u_*) (where L _MO and t _MO are the length and time scales, respectively). They show that a maximum in the wind speed profile occurs at the top of the thermal boundary layer for weak surface cooling, i.e. a wind jet, whereas there is a flattening of the profile and no marked maximum for weak surface heating. The modelled profiles are approximately the same as those obtained from the U.K. Met Office Unified Model when operating as a mesoscale model at 12-km horizontal resolution. The theoretical model is modified when strong surface heating is suddenly applied, resulting in a large change in h /|L _MO| (>>1), over the time scale t _MO. The eddy structure is predicted to change significantly and the addition of convective turbulence increases the shear turbulence at the ground. A low-level wind jet can form, with convective turbulence adding to the mean momentum of the flow. This was verified by our laboratory experiment and direct numerical simulations. Additionally, it is shown that the effects of Coriolis acceleration diminish (rather than as suggested in the literature, amplify) the formation of the wind jets in the situations considered here. Hence, only when the surface heat flux changes over time scales greater than 1/f (where f is the Coriolis parameter) does the ABL adjust monotonically between its equilibrium states. These results are also applicable to the ABL passing over spatially varying surface heat fluxes.     

Current problems in the stratocumulus-topped atmospheric boundary layer

Extended sheets of stratocumulus (Sc) in the upper part of the atmospheric boundary layer (ABL) often occur under appropriate meteorological conditions. These cloud decks are important both in climate studies and in weather forecasting. We review the current knowledge of the turbulent structure of the ABL capped by a cloud deck, in the light of recent observations and model studies. The most important physical processes determining this structure are longwave radiative cooling at cloud top, shortwave radiative wanning by absorption in the cloud, surface buoyancy flux, and wind shear in the ABL. As a result, turbulence can cause entrainment against the buoyancy jump at cloud top. In cases where only longwave radiative fluxes and surface buoyancy fluxes are important, the turbulent structure is relatively well understood. When shortwave radiative fluxes and/or wind shear are also important, the resulting turbulent structure may change considerably. A decoupling of the cloud from the sub-cloud layer or of the top of the cloud from the rest of the ABL is then regularly observed. In no cases are the details of the entrainment at cloud top understood well enough to derive a relatively simple formulation that is consistent with observations. Cloud-top entrainment instability may lead to the break-up of a cloud deck (but also to cloud deepening). The role of mesoscale circulations in determining fractional cloudiness is not yet well understood.     

Very-Large-Scale Motions in the Atmospheric Boundary Layer Educed by Snapshot Proper Orthogonal Decomposition

Large-eddy simulations of the atmospheric boundary layer (ABL) under a wide range of stabilities are conducted to educe very-large-scale motions and then to study their dynamics and how they are influenced by buoyancy. Preliminary flow visualizations suggest that smaller-scale motions that resemble hairpins are embedded in much larger scale streamwise meandering rolls. Using simulations that represent more than 150 h of physical time, many snapshots in the \(xy\) -, \(yz\) - and \(xz\) -planes are then collected to perform snapshot proper orthogonal decomposition and further investigate the large structures. These analyses confirm that large streamwise rolls that share several features with the very-large-scale motions observed in laboratory studies arise as the dominant modes under most stabilities, but the effect of the surface kinematic buoyancy flux on the energy content of these dominant modes is very significant. The first two modes in the \(yz\) -plane in the neutral case contain up to 3 % of the total turbulent kinetic energy; they also have a vertical tilt angle in the \(yz\) -plane of about 0 to 30 \(^\circ \) due to the turning effect associated with the Coriolis force. Unstable cases also feature streamwise rolls, but in the convective ABL they are strengthened by rising plumes in between them, with two to four rolls spanning the whole domain in the first few modes; the Coriolis effect is much weaker in the unstable ABL. These rolls are no longer the dominant modes under stable conditions where the first mode is observed to contain sheet-like motions with high turbulent kinetic energy. Using these proper orthogonal decomposition modes, we are also able to extract the vertical velocity fields corresponding to individual modes and then to correlate them with the horizontal velocity or temperature fields to obtain the momentum and heat flux carried by individual modes. Structurally, the fluxes are explained by the topology of their corresponding modes. However, the fraction of the fluxes produced by the modes is invariably smaller than the fraction of energy they contain, particularly under stable conditions where the first modes are found to perform weak counter-gradient fluxes.     

Observations of Coherent Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer

Turbulence structures of high Reynolds number flow in the near-neutral atmospheric boundary layer (ABL) are investigated based on observations at Shionomisaki and Shigaraki, Japan. A Doppler sodar measured the vertical profiles of winds in the ABL. Using the integral wavelet transform for the time series of surface wind data, the pattern of a descending high-speed structure with large vertical extent (from the surface to more than 200-m level) is depicted from the Doppler sodar data. Essentially this structure is a specific type of coherent structure that has been previously shown in experiments on turbulent boundary-layer flows. Large-scale high-speed structures in the ABL are extracted using a long time scale (240 s) for the wavelet transform. The non-dimensional interval of time between structures is evaluated as 3.0–6.2 in most cases. These structures make a large contribution to downward momentum transfer in the surface layer. Quadrant analyses of the turbulent motion measured by the sonic anemometer (20-m height) suggest that the sweep motion (high-speed downward motion) plays a substantial role in the downward momentum transfer. In general, the contribution of sweep motions to the momentum flux is nearly equal to that of ejection motions (low-speed upward motions). This contribution of sweep motions is related to the large-scale high-speed structures.     

Dynamical Initialization of Boundary-Layer Data

A procedure for the dynamic initialization of wind and temperature fields within the atmospheric boundary layer (ABL) is presented. The procedure uses transilient turbulence theory to assess the turbulent fluxes from observed wind and temperature data, which are used to estimate the eddy transfer coefficients, so enabling a local closure ABL model to be integrated forward in time to a new steady state. The method has been applied to initialize kytoon data taken at Kharagpur (22.3° N, 87.2° E) during the MONTBLEX field programme. Results of a case study for 17 June 1990 are discussed.     

Scaling Variances of Scalars in a Convective Boundary Layer Under Different Entrainment Regimes

For the presentation and analysis of atmospheric boundary-layer (ABL) data, scales are used to non-dimensionalise the observed quantities and independent variables. Usually, the ABL height, surface sensible heat flux and surface scalar flux are used. This works well, so long as the absolute values of the entrainment ratio for both the scalar and temperature are similar. The entrainment ratio for temperature naturally ranges from −0.4 to −0.1. However, the entrainment ratio for passive scalars can vary widely in magnitude and sign. Then the entrainment flux becomes relevant as well. The only customary scalar scale that takes into account both the surface flux and the entrainment flux is the bulk scalar scale, but this scale is not well-behaved for large negative entrainment ratios and for an entrainment ratio equal to −1. We derive a new scalar scale, using previously published large-eddy simulation results for the convective ABL. The scale is derived under the constraint that scaled scalar variance profiles are similar at those heights where the variance producing mechanisms are identical (i.e., either near the entrainment layer or near the surface). The new scale takes into account that scalar variance in the ABL is not only related to the surface flux of that scalar, but to the scalar entrainment flux as well. Furthermore, it takes into account that the production of variance by the entrainment flux is an order of magnitude larger than the production of variance by the surface flux (per unit flux). Other desirable features of the new scale are that it is always positive (which is relevant when scaling standard deviations) and that the scaled variances are always of order 1–10.     

Water and Heat Transport in the Desert Soil and Atmospheric Boundary Layer in Western China

In order to understand the exchange and transferprocesses of water and energy in the desert soil andthe atmospheric boundary layer (ABL), we have developeda coupled model, in which a desert soil modelincluding water movement of both liquid and vapourphase, and an ABL model based on a non-local transilientturbulence closure scheme, are coupled together. Withthis model, the evolution of potential temperature andspecific humidity, the distribution of net radiationamong sensible, latent and soil heat fluxes, and thewater and heat flux profiles both in the soil and ABLhave been simulated. The HEIFE (HEIhe River Basin FieldExperiment) observational data are used to calibrate calculation of the water and heat flux both in thesoil and the ABL. The sensible and latent heatfluxes warm and moisten the bottom grid box (100m) of theABL. In this way the ABL model and the desert soil model are coupled together.The simulated results show that when the flux of watervapour in the soil is neglected, the evaporation rateand the flux profiles of specific humidity in the ABLshow great changes, hence the importance of watervapour movement in the desert soil for the calculationof specific humidity in the ABL. In the upper 5cm of thesoil, which is called an active layer, water andheat transport are more effective than in the substrate(soil below 5 cm).     

Aircraft-Based Measurements Of Turbulence Structures In The Katabatic Flow Over Greenland

Turbulence structures in the katabatic flow in the stable boundary layer (SBL) over the ice sheet are studied for two case studies with high wind speeds during the aircraft-based experiment KABEG (Katabatic wind and boundary layer front experiment around Greenland) in the area of southern Greenland. The aircraft data allow the direct determination of turbulence structures in the katabatic flow. For the first time, this allows the study of the turbulence structure in the katabatic wind system over the whole boundary layer and over a horizontal scale of 80 km.The katabatic flow is associated with a low-level jet (LLJ), with maximum wind speeds up to 25 m s^-1. Turbulent kinetic energy (TKE) and the magnitude of the turbulent fluxes show a strong decrease below the LLJ. Sensible heat fluxes at the lowest level have values down to -25 W m^-2. Latent heat fluxes are small in general, but evaporation values of up to +13 W m^-2 are also measured. Turbulence spectra show a well-defined inertial subrange and a clear spectral gap around 250-m wavelength. While turbulence intensity decreases monotonously with height above the LLJ for the upper part of the slope, high spectral intensities are also present at upper levels close to the ice edge. Normalized fluxes and variances generally follow power-law profiles in the SBL.Terms of the TKE budget are computed from the aircraft data. The TKE destruction by the negative buoyancy is found to be very small, and the dissipation rate exceeds the dynamical production.     

Observations Of The Turbulence Structure Within Two Stratocumulus-Topped, Marine Boundary Layers

Two situations observed during the second Aerosol Characterization Experiment (ACE-2) are analysed from aircraft measurements in the broken stratocumulus (Sc)-topped marine boundary layer. The first one (26 June 1997), characterized by a non-polluted, oceanic air mass, presents a decoupling between the Sc layer (1400–1520 m) and the turbulent mixed layer, this latter extending from the surface up to 580 m. In contrast, the second case (9 July 1997), during which continental air had been advected over the experimental area, presents a well-coupled layer extending from the surface up to the top of the Sc layer(910 m). This coupling, uncommon in this area in the middle of the day, isrelated to the relative shallowness of the boundary layer. For both situations,it is shown that the turbulent fluxes can be computed with reasonably goodaccuracy (better than 10 %), taking into account both the random and thesystematic errors involved in the eddy-correlation technique. Estimationof random error is based on the computation of the integral scale of thecovariance, and systematic error is estimated from the parameterizationof Mann and Lenschow. The fluxes show that the buoyancy, as a sourceof turbulence, is due to latent heat flux rather than sensible heat flux,with values comparable to previous experiments in the Azores-Canariesbasin. In addition, we propose a method to analyse, for coupled situations,the relationship between the fractional cloudiness and the organization ofthe turbulent field below the clouds. This method is based on a conditionalsampling technique. It is shown that this organization cannot be deducedfrom the analysis of the velocity signal, which is dominated by turbulence.However, when the signals are conditionally sampled according to thepresence or absence of clouds, a weak cloud-related organization can beshown, and the cloud-related transports quantified; the values found areof the order of 10 % of the total transfers, i.e. the same order of magnitude asthe errors on the total flux computation. The method developed is thereforepromising, provided that the uncertainties can be reduced by analyzing a highamount of data.     

A Case Study on the Effects of Heterogeneous Soil Moisture on Mesoscale Boundary-Layer Structure in the Southern Great Plains, U.S.A. Part I: Simple Prognostic Model

The atmospheric boundary-layer (ABL) depth was observed by airborne lidar and balloon soundings during the Southern Great Plains 1997 field study (SGP97). This paper is Part I of a two-part case study examining the relationship of surface heterogeneity to observed ABL structure. Part I focuses on observations. During two days (12–13 July 1997) following rain, midday convective ABL depth varied by as much as 1.5 km across 400 km, even with moderate winds. Variability in ABL depth was driven primarily by the spatial variation in surface buoyancy flux as measured from short towers and aircraft within the SGP97 domain. Strong correlation was found between time-integrated buoyancy flux and airborne remotely sensed surface soil moisture for the two case-study days, but only a weak correlation was found between surface energy fluxes and vegetation greenness as measured by satellite. A simple prognostic one-dimensional ABL model was applied to test to what extent the soil moisture spatial heterogeneity explained the variation in north–south ABL depth across the SGP97 domain. The model was able to better predict mean ABL depth and variations on horizontal scales of approximately 100 km using observed soil moisture instead of constant soil moisture. Subsidence, advection, convergence/divergence and spatial variability of temperature inversion strength also contributed to ABL depth variations. In Part II, assimilation of high-resolution soil moisture into a three-dimensional mesoscale model (MM5) is discussed and shown to improve predictions of ABL structure. These results have implications for ABL models and the influence of soil moisture on mesoscale meteorology     

Regional climate simulation with a high resolution GCM: surface radiative fluxes

The ability of a high resolution (T106) version of the ECHAM3 general circulation model to simulate regional scale surface radiative fluxes has been assessed using observations from a new compilation of worldwide instrumentally-measured surface fluxes (Global Energy Balance Archive, GEBA). The focus is on the European region where the highest density of observations is found, and their use for the validation of global and regional climate models is demonstrated. The available data allow a separate assessment of the simulated fluxes of surface shortwave, longwave, and net radiation for this region. In summer, the incoming shortwave radiation calculated by the ECHAM3/T106 model is overestimated by 45 W m^–2 over most of Europe, which implies a largely unrealistic forcing on the model surface scheme and excessive surface temperatures. In winter, too little incoming shortwave radiation reaches the model surface. Similar tendencies are found over large areas of the mid-latitudes. These biases are consistent with deficiencies in the simulation of cloud amount, relative humidity and clear sky radiative transfer. The incoming longwave radiation is underestimated at the European GEBA stations predominantly in summer. This largely compensates for the excessive shortwave flux, leading to annual mean net radiation values over Europe close to observations due to error cancellation, a feature already noted in the simulated global mean values in an earlier study. Furthermore, the annual cycle of the simulated surface net radiation is strongly affected by the deficiencies in the simulated incoming shortwave radiation. The high horizontal resolution of the GCM allows an assessment of orographically induced flux gradients based on observations from the European Alps. Although the model-calculated and observed flux fields substantially differ in their absolute values, several aspects of their gradients are realistically captured. The deficiencies identified in the model fields are generally consistent at most stations, indicating a high degree of representativeness of the measurements for their larger scale setting.     

TRAC98: Detection of Coherent Structures in a Convective Boundary Layer using Airborne Measurements

The TRAC98 experimental campaign (Turbulence Radar Aircraft Cells) devoted to coherent structures analysis took place over the Beauce plain (France) during summer 1998. It allowed us to collect a large dataset of airborne measurements in addition to various ground measurements. This study aims at diagnosing the occurrence of coherent structures within the atmospheric boundary layer (ABL) through airborne measurements. The statistical analysis performed as a first step from turbulent parameters underlined the homogeneity of the ABL over the Beauce plain. However mixed-layer scaling failed at the top of the ABL, even when taking into account the entrainment rate. Coherent structures were detected through the analysis of ABL isotropy, using the opportunity of sampling with two perpendicular crossing planes, one of them being aligned with the wind. This approach allowed us to determine an organization scheme of the ABL for three of the five flights (ARAT30, MIV30 and MIV27). For the ARAT30 flight, the analysis was pursued by focusing on measurements of fluctuations in the inner flight legs. In this way, the low-level cloud cover has been investigated from the downward visible radiation (VISD). The results indicated an anisotropy of the horizontal cloud size. Secondly, the variations of some parameters were analysed through lagged correlation functions. This allowed us to infer relationships between the vertical velocity, VISD, mixing ratio and lifting condensation level. Length scales have also been extracted, and confirmed the ABL organization during the ARAT30 flight. Finally, the anisotropy observed in various flights has been investigated with respect to the underestimation of the latent heat fluxes revealed by the imbalance of the surface energy budget.     

Modelling of cloud formation due to air-sea interactions in an energy-active zone

A mesoscale 3D numerical model is described, with which detailed calculations have been made of turbulence and wind characteristics in the atmospheric boundary layer (ABL), as well as cloud particle size distribution, longwave and solar radiation fluxes and flux divergences, and atmosphere-ocean heat exchange. Based on numerical experiments simulating winter conditions of the Newfoundland energy-active zone of the ocean (EAZO), atmosphere-ocean energy exchange is investigated. It is shown that the basic mechanisms for the EAZO formation involve the following processes: (i) at the hydrological front between cold and warm ocean currents, the fluxes of sensible and latent heat grow significantly; (ii) at this front, in a particular synoptic situation, overcast low-level cloudiness forms, screening solar radiation so that in winter, the radiation budget at the front is reduced, and the radiative flux into the ocean is less than the energy release to the atmosphere; (iii) frequent occurrence of such synoptic situations with cloudiness decreases the oceanic enthalpy and creates negative SST anomalies. The transport of these anomalies by currents to the western coasts of the continents causes anomalies of weather and climate.