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.     

Correction to: Spatio-temporal variability of atmospheric rivers and associated atmospheric parameters in the Euro-Atlantic region

Theoretical and Applied Climatology -     

Spatio-temporal variability of atmospheric rivers and associated atmospheric parameters in the Euro-Atlantic region

We study the spatio-temporal variability of Atmospheric Rivers (ARs) and associated integrated water vapor and atmospheric parameters over the Euro-Atlantic region using long-term reanalysis datasets. Winds, temperature, and specific humidity at different pressure levels during 1979–2018 are used to study the water vapor transport integrated between 1000 and 300 hPa (IVT300) in mapping ARs. The intensity of ARs in the North Atlantic has been increasing in recent times (2009–2018) with large decadal variability and poleward shift (~ 5° towards the North) in landfall during 1999–2018. Though different reanalysis datasets show similar spatial patterns of IVT300 in mapping ARs, bias in specific humidity and wind components led to IVT300 mean bias of 50 kg m⁻¹ s⁻¹ in different reanalysis products compared to ERA5. The magnitude of winds and specific humidity in the lower atmosphere (below 750 hPa) dominates the total column water vapor and intensity of ARs in the North Atlantic. Reanalysis datasets in the central North Atlantic show an IVT300 standard deviation of 200 kg m⁻¹ s⁻¹ which is around 33% of the ARs climatology (~ 600 kg m⁻¹ s⁻¹). Though ARs have a higher frequency of landfalling over Western Europe in winter half-year, the intensity of IVT300 in winter ARs is 3% lower than the annual mean. The lower frequency of ARs in the summer half-year shows 3% higher IVT300 than the annual mean. While ARs in the North Atlantic show a strong decadal change in frequency and path, the impact of the North Atlantic Oscillation (NAO) and Scandinavian blocking on the location of landfall of ARs are significant. Furthermore, there is a strong latitudinal dependence of the source of moisture flux in the open ocean, contributing to the formation and strengthening ARs.

     

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⁻².     

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.     

Methods for Estimating Air–Sea Fluxes of CO<Subscript>2</Subscript> Using High-Frequency Measurements

The most direct method for flux estimation uses eddy covariance, which is also the most commonly used method for land-based measurements of surface fluxes. Moving platforms are frequently used to make measurements over the sea, in which case motion can disturb the measurements. An alternative method for flux estimation should be considered if the effects of platform motion cannot be properly corrected for. Three methods for estimating CO₂ fluxes are studied here: the eddy-covariance, the inertial-dissipation, and the cospectral-peak methods. High-frequency measurements made at the land-based Östergarnsholm marine station in the Baltic Sea and measurements made from a ship during the Galathea 3 expedition are used. The Kolmogorov constant for CO₂, used in the inertial-dissipation method, is estimated to be 0.68 and is determined using direct flux measurements made at the Östergarnsholm site. The cospectral-peak method, originally developed for neutral stratification, is modified to be applicable in all stratifications. With these modifications, the CO₂ fluxes estimated using the three methods agree well. Using data from the Östergarnsholm site, the mean absolute error between the eddy-covariance and inertial-dissipation methods is 0.25 μmol  m^?2 s^?1. The corresponding mean absolute error between the eddy-covariance and cospectral-peak methods is 0.26 μmol m^?2 s^?1, while between the inertial-dissipation and cospectral-peak methods it is 0.14 μmol m^?2 s^?1.     

Spectra of CO<Subscript>2</Subscript> and Water Vapour in the Marine Atmospheric Surface Layer

Spectra of CO₂ and water vapour fluctuations from measurements made in the marine atmospheric surface layer have been analyzed. A normalization of spectra based on Monin–Obukhov similarity theory, originally developed for wind speed and temperature, has been successfully extended also to CO₂ and humidity spectra. The normalized CO₂ spectra were observed to have somewhat larger contributions from low frequencies compared to humidity spectra during unstable stratification. However, overall, the CO₂ and humidity spectra showed good agreement as did the cospectra of vertical velocity with water vapour and CO₂ respectively. During stable stratification the spectra and cospectra displayed a well-defined spectral gap separating the mesoscale and small-scale turbulent fluctuations. Two-dimensional turbulence was suggested as a possible source for the mesoscale fluctuations, which in combination with wave activity in the vertical wind is likely to explain the increase in the cospectral energy for the corresponding frequency range. Prior to the analysis the turbulence time series of the density measurements were converted to time series of mixing ratios relative to dry air. Some differences were observed when the spectra based on the original density measurements were compared to the spectra based on the mixing ratio time series. It is thus recommended to always convert the density time series to mixing ratio before performing spectral analysis.     

The effect of idealized water waves on the turbulence structure and kinetic energy budgets in the overlying airflow

The influence of an idealized moving wavy surface on the overlying airflow is investigated using direct numerical simulations (DNS). In the present simulations, the bulk Reynolds number is Re = 8000 (; where U₀ is the forcing velocity of the flow, h the height of the domain and v the kinematic viscosity) and the phase speed of the imposed waves relative to the friction velocity, i.e., the wave age varies from very slow to fast waves. The wave signal is clearly present in the airflow up to at least 0.15λ (where λ is the wave length) and is present up to higher levels for faster waves. In the kinetic energy budgets, pressure transport is mainly of importance for slow waves. For fast waves, viscous transport and turbulent transport dominate near the surface. Kinetic energy budgets for the wave and turbulent perturbations show a non-negligible transport of turbulent kinetic energy directed from turbulence to the wave perturbation in the airflow. The wave-turbulent energy transport depends on the size, tilt, and phase of the wave-induced part of the turbulent Reynolds stresses.According to the DNS data, slow waves are more efficient in generating isotropic turbulence than fast waves.Despite the differences in wave-shape as well as in Reynolds number between the idealized direct numerical simulations and the atmosphere, there are intriguing similarities in the turbulence structure. Important information about the turbulence above waves in the atmosphere can be obtained from DNS—the data must, however, be interpreted with care.     

Air–Sea Interaction Features in the Baltic Sea and at a Pacific Trade-Wind Site: An Inter-comparison Study

A systematic comparison of wind profiles and momentum exchange at a trade wind site outside Oahu, Hawaii and corresponding data from the Baltic Sea is presented. The trade wind data are to a very high degree swell dominated, whereas the Baltic Sea data include a more varied assortment of wave conditions, ranging from a pure growing sea to swell. In the trade wind region swell waves travel predominantly in the wind direction, while in the Baltic, significant cross-wind swells are also present. Showing the drag coefficient as a function of the 10-m wind speed demonstrates striking differences for unstable conditions with swell for the wind-speed range 2 m s^?1 < U ₁₀ < 7 m s^?1, where the trade-wind site drag values are significantly larger than the corresponding Baltic Sea values. In striking contrast to this disagreement, other features studied are surprisingly similar between the two sites. Thus, exactly as found previously in Baltic Sea studies during unstable conditions and swell, the wind profile in light winds (3 m s^?1) shows a wind maximum at around 7–8 m above the water, with close to constant wind speed above. Also, for slightly higher wind speeds (4 m s^?1 < U ₁₀ < 7 m s^?1), the similarity between wind profiles is striking, with a strong wind-speed increase below a height of about 7–8 m followed by a layer of virtually constant wind speed above. A consequence of these wind-profile features is that Monin–Obukhov similarity is no longer valid. At the trade-wind site this was observed to be the case even for wind speeds as high as 10 m s^?1. The turbulence kinetic energy budget was evaluated for four cases of 8–16 30- min periods at the trade-wind site, giving results that agree very well with corresponding figures from the Baltic Sea.     

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.