Measured And Simulated Latent And Sensible Heat Fluxes At Two Marine Sites In The Baltic Sea

In this study, turbulent heat flux data from two sites within the Baltic Sea are compared with estimates from two models. The main focus is on the latent heat flux. The measuring sites are located on small islands close to the islands of Bornholm and Gotland. Both sites have a wide wind direction sector with undisturbed over-water fetch. Mean parameters and direct fluxes were measured on masts during May to December 1998.The two models used in this study are the regional-scale atmospheric model HIRLAM and the ocean model PROBE-Baltic. It is shown that both models overestimate the sensible and latent heat fluxes. The overestimation can, to a large extent, be explained by errors in the air-water temperature and humidity differences. From comparing observed and modelled data, the estimated 8-month mean errors in temperature and humidity are up to 1 °C and 1 g kg^-1, respectively. The mean errors in the sensible and latent heat fluxes for the same period are approximately 15 and 30 W m^-2, respectively.Bulk transfer coefficients used for calculating heat and humidity fluxes at the surface were shown to agree rather well with the measurements, at least for the unstable data. For stable stratification, the scatter in data is generally large, and it appears that the bulk formulation chosen overestimates turbulent heat fluxes.     

Area-Averaged Surface Fluxes Over the Litfass Region Based on Eddy-Covariance Measurements

Micrometeorological measurements (including eddy-covariance measurements of the surface fluxes of sensible and latent heat) were performed during the LITFASS-2003 experiment at 13 field sites over different types of land use (forest, lake, grassland, various agricultural crops) in a 20 × 20 km² area around the Meteorological Observatory Lindenberg (MOL) of the German Meteorological Service (Deutscher Wetterdienst, DWD). Significant differences in the energy fluxes could be found between the major land surface types (forest, farmland, water), but also between the different agricultural crops (cereals, rape, maize). Flux ratios between the different surfaces changed during the course of the experiment as a result of increased water temperature of the lake, changing soil moisture, and of the vegetation development at the farmland sites. The measurements over grass performed at the boundary-layer field site Falkenberg of the MOL were shown to be quite representative for the farmland part of the area. Measurements from the 13 sites were composed into a time series of the area-averaged surface flux by taking into account the data quality of the single flux values from the different sites and the relative occurrence of each surface type in the area. Such composite fluxes could be determined for about 80% of the whole measurement time during the LITFASS-2003 experiment. Comparison of these aggregated surface fluxes with area-averaged fluxes from long-range scintillometer measurements and from airborne measurements showed good agreement.     

Increased heat fluxes near a forest edge

Summary ?Observations of sensible and latent heat flux above forest downwind of a forest edge show these fluxes to be larger than the available energy over the forest. The enhancement averages to 56 W m⁻², or 16% of the net radiation, at fetches less than 400 m, equivalent to fetch to height ratios less than 15. The enhancement of turbulent energy fluxes is explained by advection and increases with the difference in temperature and humidity of the air over the upwind area as compared to the forest. The relatively high temperature and humidity of the upwind air are not caused by high surface heat fluxes, but are explained by the relatively low aerodynamic roughness of the upwind surface. Although the heat fluxes over forest are enhanced, the momentum fluxes are almost adjusted to the underlying forest. The different behaviour of heat and momentum fluxes is explained by absorption of momentum by pressure gradients near the forest edge. It is concluded that fetch requirements to obtain accurate surface fluxes from atmospheric observations need to be more stringent for scalar fluxes as compared to momentum fluxes. Received November 23, 2001; accepted May 13, 2002     

On the parameterisation of evaporation and sensible heat exchange for shallow lakes

Summary This paper demonstrates that typical air–sea interaction models can be transferred to shallow water conditions and lakes including a depth dependent function. As depth decreases, it is assumed that exchange conditions increase. From our results it follows that the model without the shallow water correction underestimates the latent and sensible heat fluxes of a lake by about 20%. This was shown through the comparison of the model with eddy covariance data for a shallow lake during the LITFASS-98 experiment in Germany. Furthermore, data selection according to the fetch and footprint conditions is presented. These models are also useful for gap filling procedures at measuring stations over lakes with limited fetch conditions in selected directions.     

On gas fluxes from small lakes and ponds

A semi-empirical model for estimating gas fluxes at the air-water interface in a neutrally stratified atmosphere is presented. Gas transfer velocity is assumed to depend on local friction velocity, whose variation with over-water fetch is estimated from Walmsleyet al.'s (1989) Guidelines. As guidance for limnologists estimating CO₂ and CH₄ fluxes from lakes and ponds in the Hudson Bay Lowlands region, we present model results of the variation of gas transfer velocity with over-water fetch,k(x), for comparison with the customary constant valuek ₀. Gas fluxes are estimated over well mixed, circular and elliptical lakes and ponds. The results of accounting for fetch dependence, wind direction, pond shapes and sizes lead to correction factors to fluxes obtained withk ₀. Recent studies of gas transfer velocity for gases of low solubility from lakes are reviewed, and the problem of time averaging on quantities determined from experimental measurements is also discussed.     

Simulation of surface energy fluxes using high-resolution non-hydrostatic simulations and comparisons with measurements for the LITFASS-2003 experiment

Land-surface heterogeneity effects on the subgrid scale of regional climate and numerical weather prediction models are of vital interest for the energy and mass exchange between the surface and the atmospheric boundary layer. High-resolution numerical model simulations can be used to quantify these effects, and are a tool used to obtain area-averaged surface fluxes over heterogeneous land surfaces. We present high-resolution model simulations for the LITFASS area near Berlin during the LITFASS-2003 experiment, which were carried out using the non-hydrostatic model FOOT3DK of the University of Köln with horizontal resolutions of 1 km and 250 m. The LITFASS-2003 experimental dataset is used for comparison. The screen level quantities show good quality for the simulated pressure, temperature, humidity and wind speed and direction. Averaged over the four week experimental period, simulated surface energy fluxes at land stations show a small bias for the turbulent heat fluxes and an underestimation of the net radiation caused by excessive cloudiness in the simulations. For eight selected days with low cloud amounts, the net radiation bias is close to zero, but the sensible heat flux shows a strong positive bias. Large differences are found for latent heat fluxes over a lake, which are partly due to local effects on the measurements, but an additional problem seems to be the overestimation of the turbulent exchange under stable conditions in the daytime internal boundary layer over the lake. In the area average over the LITFASS area of 20 ×  20 km², again a strong positive bias of 70 W m^?2 for the sensible heat is present. For the low soil moisture conditions during June 2003, the simulation of the turbulent heat fluxes is sensitive to variations in the soil type and its hydrological properties. Under these conditions, the supply of ground water to the lowest soil layer should be accounted for. Different area-averaging methods are tested. The experimental set-up of the LITFASS-2003 experiment is found to be well suited for the computation of area-averaged turbulent heat fluxes.     

Airborne measurements of turbulent fluxes during LITFASS-98: Comparison with ground measurements and remote sensing in a case study

Summary ?Simultaneous flight measurements with the research aircraft Do 128 and the helicopter-borne turbulence probe Helipod were performed on 18 June 1998 during the LITFASS-98 field experiment. The area-averaged turbulent vertical fluxes of momentum, sensible, and latent heat were determined on a 15 km × 15 km and a 10 km × 10 km flight pattern, respectively. The flights were carried out over heterogeneous terrain at different altitudes within a moderately convective boundary layer with Cumulus clouds. Co-spectra-analysis demonstrated that the small scale turbulent transport was completely sampled, while the comparatively small flight patterns were possibly of critical size regarding the large-scale turbulence. The phygoide of the airplane was identified as a significant peak in some co-spectra. The turbulent fluxes of momentum and sensible heat at 80 m above the ground showed systematic dependence on the location of the flight legs above the heterogeneous terrain. This was not observed for the latent heat flux, probably due to the vertical distribution of humidity in the boundary layer. Statistical error analysis of the fluxes F showed that the systematic statistical error ΔF was one order of magnitude smaller than the standard deviation σ _F . The difference between area-averaged fluxes derived from simultaneous Helipod and Do 128 measurements was much smaller than σ _F , indicating that the systematic statistical error was possibly over-estimated by the usual method. In the upper half of the boundary layer the airborne-measured sensible heat flux agreed well with windprofiler/RASS data. A linear fit was the best approximation for the height dependence of all three fluxes. The linear extrapolations of the latent and sensible heat fluxes to the ground were in good agreement with tower, scintillometer, and averaged ground-station measurements on various surface types. Systematic discrepancies between airborne and ground-based measurements were not found. Received June 18, 2001; revised December 21, 2001; accepted June 3, 2002     

Modelling sensible and latent heat fluxes over sea during unstable,very close to neutral conditions

During slightly unstable but still very close to neutral conditions new results from two previous investigations have shown a significant increase of sensible and latent heat fluxes over the sea. The vertical heat transport during these conditions is dominated by detached eddies originating at the top of the boundary layer, bringing relatively cold and dry air to the surface. This effect can be described in numerical models by either enhanced heat transfer coefficients for sensible and latent heat (Stanton and Dalton numbers respectively) or with an additional roughness length, added to the original roughness lengths for heat and humidity. Such new expressions are developed using turbulence measurements from the Baltic Sea valid for wind speeds up to 14 m s⁻¹. The effect of including the increased heat fluxes is investigated using two different numerical models: a regional three-dimensional climate model covering northern Europe, and a process-oriented ocean model for the Baltic Sea. During periods of several days, the latent heat flux can be increased by as much as 100 W m⁻². The increase in sensible heat flux is significantly smaller since the process is only of importance in the very near-neutral regime where the sensible heat flux is very small. The long-term average effect over the Baltic Sea is of the order of several W m⁻².     

Turbulent fluxes from Helipod flights above quasi-homogeneous patches within the LITFASS area

Turbulent fluxes of sensible and latent heat were measured with the helicopter-borne turbulence probe Helipod over a heterogeneous landscape around the Meteorological Observatory Lindenberg during the STINHO-2 and LITFASS-2003 field experiments. Besides the determination of area-averaged heat fluxes, the analysis focused on different aspects of the response of the turbulent structure of the convective boundary layer (CBL) on the surface heterogeneity. A special flight pattern was designed to study flux profiles both over quasi-homogeneous sub-areas of the study region (representing the major land use types—forest, farmland, water) and over a typical mixture of the different surfaces. Significant differences were found between the heat fluxes over the individual surfaces along flight legs at about 80 m above ground level, in agreement with large-aperture scintillometer measurements. This flux separation was still present during some flights at levels near the middle of the CBL. Different scales for the blending height and horizontal heterogeneity were calculated, but none of them could be identified as a reliable indicator of the mixing state of the lower CBL. With the exception of the flights over water, the latent heat flux measurements generally showed a larger statistical error when compared with the sensible heat flux. Correlation coefficients a nd integral length scales were used to characterise the interplay between the vertical transport of sensible and latent heat, which was found to vary between ‘fairly correlated’ and ‘decoupled’, also depending on the soil moisture conditions.     

Water vapour flux profiles in the convective boundary layer

Summary Water vapour flux profiles in the atmospheric boundary layer have been derived from measurements of water vapour density fluctuations by a ground-based Differential Absorption Lidar (DIAL) and of vertical wind fluctuations by a ground-based Doppler lidar. The data were collected during the field experiment LITFASS-2003 in May/June 2003 in the area of Lindenberg, Germany. The eddy-correlation method was applied, and error estimates of ±50 W/m² for latent heat flux were found. Since the sampling error dominates the overall measurement accuracy, time intervals between 60 and 120 min were required for a reliable flux calculation, depending on wind speed. Rather large errors may occur with low wind speed because the diurnal cycle restricts the useful interval length. In the lower height range, these measurements are compared with DIAL/radar-RASS fluxes. The agreement is good when comparing covariance and error values. The lidar flux profiles are well complemented by tower measurements at 50 and 90 m above ground and by area-averaged near surface fluxes from a network of micrometeorological stations. Water vapour flux profiles in the convective boundary layer exhibit different structures mainly depending on the magnitude of the entrainment flux. In situations with dry air above the boundary layer a positive entrainment flux is observed which can even exceed the surface flux. Flux profiles which linearly increase from the surface to the top of the boundary layer are observed as well as profiles which decrease in the lower part and increase in the upper part of the boundary layer. In situations with humid air above the boundary layer the entrainment flux is about zero in the upper part of the boundary layer and the profiles in most cases show a linear decrease.     

Heat Flux Apportionment to Heterogeneous Surfaces Using Flux Footprint Analysis

Heat flux data collected from the Baiyangdian Heterogeneous Field Experiment were analyzed using the footprint method. High resolution (25 m) Landsat-5 satellite imaging was used to determine the land cover as one of four surface types: farmland, lake, wetland, or village. Data from two observation sites in September 2005 were used. One site (Wangjiazhai) was characterized by highly heterogeneous surfaces in the central area of the Baiyangdian: lake/wetland. The other site (Xiongxian) was on land with more uniform surface cover. An improved Eulerian analytical flux footprint model was used to determine "source areas" of the heat fluxes measured at towers located at each site from surrounding landscapes of mixed surface types.In relative terms results show that wetland and lake areas generally contributed most to the observed heat flux at Wangjiazhai, while farmland contributed most at Xiongxian. Given the areal distribution of surface type contributions, calculations were made to obtain the magnitudes of the heat flux from lake, wetland and farmland to the total observed flux and apportioned contributions of each surface type to the sensible and latent heat fluxes. Results show that on average the sensible heat flux from wetland and farmland were comparable over the diurnal cycle, while the latent heat flux from farmland was somewhat larger by about 30-50 W m-2 during daytime. The latent and sensible fluxes from the lake source in daytime were about 50 W m-2 and 100 W m-2 less, respectively, than from wetland and farmland. The results are judged reasonable and serve to demonstrate the potential for flux apportionment over heterogeneous surfaces.     

Heat flux apportionment to heterogeneous surfaces using flux footprint analysis

Heat flux data collected from the Baiyangdian Heterogeneous Field Experiment were analyzed using the footprint method. High resolution （25 m） Landsat-5 satellite imaging was used to determine the land cover as one of four surface types： farmland, lake, wetland, or village. Data from two observation sites in September 2005 were used. One site （Wangjiazhai） was characterized by highly heterogeneous surfaces in the central area of the Baiyangdian： lake/wetland. The other site （Xiongxian） was on land with more uniform surface cover. An improved Eulerian analytical flux footprint model was used to determine “source areas” of the heat fluxes measured at towers located at each site from surrounding landscapes of mixed surface types. In relative terms results show that wetland and lake areas generally contributed most to the observed heat flux at Wangjiazhai, while farmland contributed most at Xiongxian. Given the areal distribution of surface type contributions, calculations were made to obtain the magnitudes of the heat flux from lake, wetland and farmland to the total observed flux and apportioned contributions of each surface type to the sensible and latent heat fluxes. Results show that on average the sensible heat flux from wetland and farmland were comparable over the diurnal cycle, while the latent heat flux from farmland was somewhat larger by about 30-50 W m-2 during daytime. The latent and sensible fluxes from the lake source in daytime were about 50 W m-2 and 100 W m-2 less, respectively, than from wetland and farmland. The results are judged reasonable and serve to demonstrate the potential for flux apportionment over heterogeneous surfaces.     

Some evaluations of drag and bulk transfer coefficients over water bodies of different sizes

Three recent experiments allow evaluation of the bulk transfer coefficients for momentum, water vapour and sensible heat over water bodies of different sizes. As part of a study of evaporation rates from a swamp, measurements of latent and sensible heat fluxes were made over Lake Wyangan in southern N.S.W., Australia. This lake is of several kilometers diameter. In a later experiment, Reynolds stress and sensible heat transfer were measured from a natural-gas platform standing in Bass Strait, south of mainland Australia. The most recent experiment involved the direct measurement of each of these turbulent fluxes from a fixed tower erected in Lake Michigan, U.S.A.Perhaps the most important of the results is the finding that drag coefficients measured over Bass Strait are not significantly different from those over Lake Michigan, despite the obvious differences in depth, fetch, and hence surface wave structure. At both locations, drag coefficients are found to increase slightly with increasing wind speed, while at low wind speeds they are not significantly different from those corresponding to aerodynamic smoothness.Near-neutral bulk transfer coefficients for sensible heat and for water vapour are found to be similar. An average value of about 1.4 × 10^–3 is obtained.It is emphasized that stability effects should be considered in any discussion of drag coefficients or bulk transfer coefficients. Large errors can result if near-neutrality is incorrectly assumed.     

Flux Footprint Simulation Downwind of a Forest Edge

Surface fluxes, originating from forest patches, are commonly calculated from atmospheric flux measurements at some height above that patch using a correction for flux arising from upwind surfaces. Footprint models have been developed to calculate such a correction. These models commonly assume homogeneous turbulence, resulting in a simulated atmospheric flux equal to the average surface flux in the footprint area. However, atmospheric scalar fluxes downwind of a forest edge have been observed to exceed surface fluxes in the footprint area. Variations in atmospheric turbulence downwind of the forest edge, as simulated with an E – model, can explain enhanced atmospheric scalar fluxes. This E – model is used to calculate the footprint of atmospheric measurements downwind of a forest edge. Atmospheric fluxes appear mainly enhanced as a result of a stronger sensitivity to fluxes from the upwind surface. A sensitivity analysis shows that the fetch over forest, necessary to reach equilibrium between atmospheric fluxes and surface fluxes, tends to be longer for scalar fluxes as compared to momentum fluxes. With increasing forest density, atmospheric fluxes deviate even more strongly from surface fluxes, but over shorter fetches. It is concluded that scalar fluxes over forests are commonly affected by inhomogeneous turbulence over large fetches downwind of an edge. It is recommended to take horizontal variations in turbulence into account when the footprint is calculated for atmospheric flux measurements downwind of a forest edge. The spatially integrated footprint is recommended to describe the ratio between the atmospheric flux and the average surface flux in the footprint.     

A New Method for the Determination of Area-Averaged Turbulent Surface Fluxes from Low-Level Flights Using Inverse Models

The low-level flight method (LLF) has been combined with linear inverse models (IM) resulting in an LLF+IM method for the determination of area-averaged turbulent surface fluxes. With this combination, the vertical divergences of the turbulent latent and sensible heat fluxes were calculated from horizontal flights. The statistical errors of the derived turbulent surface fluxes were significantly reduced. The LLF+IM method was tested both in numerical and field experiments. Large-eddy simulations (LES) were performed to compare ‘true’ flux profiles with ‘measurements’ of simulated flights in an idealised convective boundary layer. Small differences between the ‘true’ and the ‘measured’ fluxes were found, but the vertical flux divergences were correctly calculated by the LLF+IM method. The LLF+IM method was then applied to data collected during two flights with the Helipod, a turbulence probe carried by a helicopter, and with the research aircraft Do 128 in the LITFASS-98 field campaign. The derived surface fluxes were compared with results from eddy-covariance surface stations and with large-aperture scintillometer data. The comparison showed that the LLF+IM method worked well for the sensible heat flux at 77 and 200 m flight levels, and also for the latent heat flux at the lowest level. The model quality control indicated failures for the latent heat flux at the 200 m level (and higher), which were probably due to large moisture fluctuations that could not be modelled using linear assumptions. Finally the LLF+IM method was applied to more than twenty low-level flights from the LITFASS-2003 experiment. Comparison with aggregated surface flux data revealed good agreement for the sensible heat flux but larger discrepancies and a higher statistical uncertainty for the latent heat flux     

Radiative surface temperature and energy balance of a wheat canopy

The sensible heat loss from a stand of winter wheat was calculated from radiometric measurements of crop surface temperature, measurements of air temperature, and an atmospheric resistance to momentum transfer; corresponding latent heat flux was obtained through the energy balance equation. These estimates of sensible and latent heat were compared with fluxes from the Bowen Ratio method. When radiative temperature was derived using a measured canopy emissivity of 0.98, calculations of sensible heat flux were systematically 50–100 W m^-2 less than Bowen Ratio values. The two techniques agreed more closely when an apparent emissivity of 0.96 was used with an apparent reflectivity of 0.03. The mean difference between the estimates of latent heat flux was then -16 ± 32 W m^-2.The surface temperature method showed less systematic error in comparison with the Bowen Ratio values than did estimates using the aerodynamic method.On leave from: University of Nottingham, School of Agriculture, Loughborough LE12 5RD.     

Local boundary-layer development over burnt and unburnt tropical savanna: an observational study

Fire scars have the ability to radically alter the surface energy budget within a tropical savanna by reducing surface albedo, increasing available energy for partitioning into sensible and latent heat fluxes and increasing substrate heat flux. These changes have the potential to alter boundary-layer conditions and ultimately feedback to local and regional climate. We measured radiative and energy fluxes over burnt and unburnt tropical savanna near Howard Springs, Darwin, Australia. At the burnt site a low to moderate intensity fire, ranging between 1,000 and 3,500 kW m⁻¹, initially affected the land surface by removing all understorey vegetation, charring and blackening the ground surface, scorching the overstorey canopy and reducing the albedo. A reduction in latent heat fluxes to almost zero was seen immediately after the fire when the canopy was scorched. This was then followed by an increase in the sensible heat flux and a large increase in the ground heat flux over the burnt surface. Tethered balloon measurements showed that, despite the presence of pre-monsoonal rain events occurring during the measurement period, the lower boundary layer over the burnt site was up to 2^°C warmer than that over the unburnt site. This increase in boundary-layer heating when applied to fire scars at the landscape scale can have the ability to form or alter local mesoscale circulations and ultimately create a feedback to regional heating and precipitation patterns that may affect larger-scale processes such as the Australian monsoon.     

夏季青海湖局地环流及大气边界层特征的数值模拟

使用美国NCAR新版MM5V3.6非静力模式,采用两重嵌套方法,模拟了青海湖区域的局地环流及大气边界层特征,并且与无湖试验进行了比较。结果表明:白天由于青海湖的存在有很好的降温作用,夜晚则有保温效应,表现出明显的冷(暖)湖效应;青海湖对感热和潜热的影响有很强的日变化,白天湖面感热、潜热都小,夜间情况相反,这使得白天青海湖是冷干岛,夜间是暖湿岛;青海湖使得白天湖面边界层顶低,陆面边界层顶高,夜间相反。这样的边界层顶高度和温度、地面能量通量相配合,形成了一个很好的保护机制,对青海湖的水土保持和生态环境的维持有正效应;青海湖使得湖面上空大气下沉,陆面上空大气上升,从而产生了湖面上空大气冷干,陆面上空大气暖湿的边界层特征;青海湖边缘的陆面形成的较大的湿气柱围绕着湖面,起到了保护湖面的作用;青海湖低空白天有明显的湖面向四周的辐散气流,而夜间则为从北偏东方向来的陆风。     

Creating Surface Flux Maps from Airborne Measurements: Application to the Mackenzie Area GEWEX Study MAGS 1999

The objective of this study is to produce two-dimensional maps of the sensible and the latent heat fluxes from airborne measurements, based on the analysis of a flight pattern, called grid flights. A footprint model with along-wind and cross-wind components was used to project the measured fluxes onto the surface map. The method was applied to measurements over Arctic tundra during the Mackenzie Area GEWEX (Global Energy and Water Cycle Experiment) Study (MAGS) 1999. The resulting flux estimates were computed by integration of a wavelet transform, and corrected for long wavelength losses using information from 100 km long regional runs that were conducted close to the grid flights. The random flux error was estimated based on the flight length that is represented in each map element, and a map resolution of 3 × 3 km was chosen in order to keep the average random error of the latent heat flux below 25%. The random error of the sensible heat flux was smaller by a factor of 1.4 on average. An analysis of airborne flux measurements at different altitudes showed no significant increase of flux estimates for measurement heights below 90 m. Thus, the fluxes measured at heights between 48 and 64 m were not corrected for vertical flux divergence. The resulting flux maps provide quantitative two-dimensional estimates of the energy exchange between the surface and the atmosphere during the snow melt period in an Arctic environment, which are well-suited for calibration and validation of numerical models.     

Influence from Surrounding Land on the Turbulence Measurements Above a Lake

Turbulence measurements taken at a Swedish lake are analyzed. Although the measurements took place over a relatively large lake with several km of undisturbed fetch, the turbulence structure was found to be highly influenced by the surrounding land during daytime. Variance spectra of both horizontal velocity and scalars during both unstable and stable stratification displayed a low frequency peak. The energy at lower frequencies showed a daily variation, increasing in the morning and decreasing in the afternoon. This behaviour is explained by spectral lag, where the low frequency energy due to large eddies that originate from the convective boundary layer above the surrounding land. When the air is advected over the lake the small eddies rapidly equilibrate with the new surface forcing. However, the large eddies remain for an appreciable distance and influence the turbulence in the developing lake boundary layer. The variances of the horizontal velocity and scalars are increased by these large eddies, while the turbulent fluxes are mainly unaffected. The drag coefficient, Stanton number and Dalton number used to parametrize the momentum flux, heat flux and latent heat flux respectively all compare well with current parametrizations developed for open sea conditions. The diurnal cycle of the partial pressure of methane, $p\mathrm{CH}_{4}$ , observed at this site is closely related to the diurnal cycle of the lake-air methane flux. An idealized two-dimensional model simulation of the boundary layer at a lake site indicates that the strong response of $p\mathrm{CH}_{4}$ to the surface methane flux is due to the shallow internal boundary layer that develops above the lake, allowing methane to accumulate in a relatively small volume.