A model for the height of the internal boundary layer over an area with an irregular coastline

A model for the time and space variation of the internal boundary-layer height over a land area with an irregular coastline is presented. It is based on the analytical model of the boundary-layer height proposed by Gryning and Batchvarova (1990) and Batchvarova and Gryning (1991), The model accounts for the temperature jump and the mean vertical air motion at the top of the internal boundary-layer. Four cases from experiments in Nanticoke and Vancouver are used for model validation. The agreement between the calculated and measured internal boundary layer height at the observational sites is fairly good. The input information for the model consist of wind speed and direction, friction velocity and kinematic heat flux in time and space for the area, and the potential temperature gradient and the mean vertical air motion above the internal boundary layer. For the experiments used in the validation the effect of subsidence is relatively important in the afternoon under low wind speed high pressure conditions, lowering the height of the internal boundary layer by up to 10%, and it is negligible in the morning hours. The effect of the mixing height over the sea is found to be negligible.     

Marine Boundary Layer And Turbulent Fluxes Over The Baltic Sea: Measurements And Modelling

Two weeks of measurements of the boundary-layer height over a small island (Christiansø) in the Baltic Sea are discussed. The meteorological conditions are characterised by positive heat flux over the sea. The boundary-layer height was simulated with two models, a simple applied high-resolution (2 km × 2 km) model, and the operational numerical weather prediction model HIRLAM (grid resolution of 22.5 km × 22.5 km). For southwesterly winds it was foundthat a relatively large island (Bornholm) lying 20-km upwind of the measuring site influences the boundary-layer height. In this situation the high-resolution simple applied model reproduces the characteristics of the boundary-layer height over the measuring site. Richardson-number based methods using data from simulations with the HIRLAM model fail, most likely because the island and the water fetch to the measuring site are about the size of the grid resolution of the HIRLAM model and therefore poorly resolved. For northerly winds, the water fetch to the measuring site is about 100 km. Both models reproduce the characteristics of the height of the marine boundary layer. This suggests that the HIRLAM model adequately resolves a water fetch of 100 km with respect to predictions of the height of the marine boundary layer.     

Estimation Of The Lagrangian Structure Function Constant C0 From Surface-Layer Wind Data

Eulerian turbulence observations, madein the surface layer under unstable conditions (z/L > 0),by a sonic anemometer were used to estimatethe Lagrangian structure function constant _O. Twomethods were considered. The first one makes use of arelationship, widely used in the Lagrangian stochasticdispersion models, relating _O to the turbulent kineticenergy dissipation rate , wind velocity variance andLagrangian decorrelation time. The second one employsa novel equation, connecting _O to the constant of thesecond-order Eulerian structure function. Beforeestimating _O, the measurements were processed in orderto discard non-stationary cases at least to a firstapproximation and cases in which local isotropy couldnot be assumed. The dissipation was estimated eitherfrom the best fit of the energy spectrum in theinertial subrange or from the best fit of the third-orderlongitudinal Eulerian structure function. Thefirst method was preferred and applied to the subsequentpart of the analysis. Both methods predict thepartitioning of _O in different spatial components as aconsequence of the directional dependence of theEulerian correlation functions due to the isotropy.The _O values computed by both methods are presented anddiscussed. In conclusion, both methods providerealistic estimates of _O that compare well withprevious estimations reported in the literature, evenif a preference is to be attributed to the second method.     

Precipitation and evaporation budgets over the Baltic Proper: observations and modelling

Precipitation and evaporation budgets over the Baltic Sea were studied in a concerted project called PEP in BALTEX (Pilot study of Evaporation and Precipitation in the Baltic Sea), combining extensive field measurements and modelling efforts. Eddy-correlation-measurements of turbulent heat flux were made on a semi-continuous basis for a 12 month period at four well-exposed coastal sites in the Baltic Proper (the main basin of the Baltic Sea). Precipitation was measured at land-based sites with standard gauges and on four merchant ships travelling between Germany and Finland with the aid of specially designed ship rain gauges (SRGs). The evaporation and precipitation regime of the Baltic Sea was modelled for a 12 month period by applying a wide range of numerical models: the operational atmospheric High Resolution Limited Area Model (HIRLAM, Swedish and Finnish versions), the German atmospheric REgional-scale MOdel, REMO, the operational German Europe Model (only precipitation), the oceanographic model PROBE-Baltic, and two models that use interpolation of ground-based data, the Swedish MESAN model of SMHI and a German model of IFM-GEOMAR Kiel. Modelled precipitation was compared with SRG measurements on board the ships. A reasonable correlation was obtained, but the regional-scale models and MESAN gave some 20% higher precipitation over the sea than is measured. Bulk parameterisation schemes for evaporation were evaluated against measurements. A constant value of C_HN and C_EN with wind speed, underestimated large fluxes of both sensible and latent heat flux. The limited area models do not resolve the influence of the height of the marine boundary layer in coastal zones and the entrainment (on the surface fluxes), which may explain the observed low correlations between modelled and measured latent heat fluxes. Estimates of evaporation, E, and precipitation, P, for the entire Baltic Proper were made with several models for a 12 month period. While the annual variation was well represented by all predictions, there are still important differences in the annual means. Evaporation ranges from 509 to 625 mm year^-1 and precipitation between 624 and 805 mm year^-1 for this particular 12 month period. Taking the results of model verification from the present study into account, the best estimate of P-E is about 100 ± 50 mm for this particular 12 month period. But the annual mean of P-E varies considerably from year to year. This is reflected in simulations with the PROBE-Baltic model for an 18 year period, which gave 95 mm year^-1 for the 12 month period studied here and 32 mm year^-1 as an average for 18 years.     

Mixing Height Over Water And Its Role On The Correlation Between Temperature And Humidity Fluctuations In The Unstable Surface Layer

Results from an experimental investigation of themixing height over inner Danish waters carriedout from September 1990 to October 1992, are discussed.The statistical analysis of the mixed-layer height (z_i)over the sea does not exhibit the dailyvariation that is characteristic of the mixed layerover land, but it is nearly constant over a24-hour cycle. During summer, the mixed layer ishigher than during winter. A second inversionwas often observed.A case study of the development of the mixed layerover the sea under near-neutral and unstableatmospheric conditions during six consecutivedays is presented. A zero-order mixed-layer heightmodel is applied. In addition to momentum and heatfluxes the effect of subsidence was found to be importantfor the evolution of the mixed layer over the sea. Themodelled evolution of z_i compared successfullywith measurements.We have investigated the influence of themixed-layer height on the correlation coefficient R_qTbetween temperature and humidity fluctuations usingthe values obtained with the model.We found that the evolution of R_qT follows theevolution of the mixing height. An empirical modellinking the surface values of R_qT to z_i and the Obukhov scaling length L has been suggested. The modelreproduces the experimental features.     

Modelling Internal Boundary-Layer Development in a Region with a Complex Coastline

The purpose of this paper is to test the ability of two quite different models to simulate the combined spatial and temporal variability of the internal boundary layer in an area of complex terrain and coastline during one day. The simple applied slab model of Gryning and Batchvarova, and the Colorado State University Regional Atmospheric Modelling System (CSU-RAMS) are tested by comparison with data gathered during a field study (called Pacific '93) of photochemical pollution in the Lower Fraser Valley of British Columbia, Canada. The data utilised here are drawn from tethered balloon flights, free flying balloon ascents, and downlooking lidar operated from an aircraft flown at roughly 3500 m above sea level. Both models are found to represent the temporal and spatial development of the internal boundary-layer depth over the Lower Fraser Valley very well, and reproduce many of the finer details revealed by the measurements.     

Humidity fluctuations in the marine boundary layer measured at a coastal site with an infrared humidity sensor

An extensive set of humidity turbulence data has been analyzed from 22-m height in the marine boundary layer. Fluctuations of humidity were measured by an OPHIR, an infrared humidity sensor with a 10 Hz scanning frequency and humidity spectra were produced. The shapes of the normalized spectra follow the established similarity functions. However the 10-min time averaged measurements underestimate the value of the absolute humidity. The importance of the humidity flux contribution in a marine environment in calculating the Obukhov stability length has been studied. Deviations from Monin-Obukhov similarity theory seem to be connected to a low correlation between humidity and temperature.     

On the extension of the wind profile over homogeneous terrain beyond the surface boundary layer

Analysis of profiles of meteorological measurements from a 160 m high mast at the National Test Site for wind turbines at Høvsøre (Denmark) and at a 250 m high TV tower at Hamburg (Germany) shows that the wind profile based on surface-layer theory and Monin-Obukhov scaling is valid up to a height of 50–80 m. At higher levels deviations from the measurements progressively occur. For applied use an extension to the wind profile in the surface layer is formulated for the entire boundary layer, with emphasis on the lowest 200–300 m and considering only wind speeds above 3 m s^?1 at 10 m height. The friction velocity is taken to decrease linearly through the boundary layer. The wind profile length scale is composed of three component length scales. In the surface layer the first length scale is taken to increase linearly with height with a stability correction following Monin-Obukhov similarity. Above the surface layer the second length scale (L _MBL ) becomes independent of height but not of stability, and at the top of the boundary layer the third length scale is assumed to be negligible. A simple model for the combined length scale that controls the wind profile and its stability dependence is formulated by inverse summation. Based on these assumptions the wind profile for the entire boundary layer is derived. A parameterization of L _MBL is formulated using the geostrophic drag law, which relates friction velocity and geostrophic wind. The empirical parameterization of the resistance law functions A and B in the geostrophic drag law is uncertain, making it impractical. Therefore an expression for the length scale, L _MBL , for applied use is suggested, based on measurements from the two sites.