EXPERIMENTAL STUDY OF THE SURFACE MARINE ATMOSPHERIC BOUNDARY LAYER OVER AEGEAN SEA, GREECE

Similar documents
Spectral characteristics of the wind components in the surface Atmospheric Boundary Layer

Surface Fluxes and Wind-Wave Interactions in Weak Wind Conditions

Conditions for Offshore Wind Energy Use

ON THE PHYSICAL PROCESSES THAT INFLUENCE THE DEVELOPMENT OF THE MARINE LOW-LEVEL JET

Gravity waves in stable atmospheric boundary layers

Study of the impact of a new Drag Coefficient parameterization on the forecasting skill of an Ocean Model in the Aegean and Levantine Seas

Observed Roughness Lengths for Momentum and Temperature on a Melting Glacier Surface

Flow modelling hills complex terrain and other issues

Dances with waves. Dances with waves: the knmi wave follower. by Cor Jacobs, Wiebe Oost, Cor van Oort and Ed Worrell

Air-Sea Interaction Spar Buoy Systems

RECTIFICATION OF THE MADDEN-JULIAN OSCILLATION INTO THE ENSO CYCLE

ANALYSIS OF TURBULENCE STRUCTURE IN THE URBAN BOUNDARY LAYER. Hitoshi Kono and Kae Koyabu University of Hyogo, Japan

An Investigation of Atmospheric Stability and Its Impact on Scatterometer Winds Across the Gulf Stream

2.5 SHIPBOARD TURBULENCE MEASUREMENTS OF THE MARINE ATMOSPHERIC BOUNDARY LAYER FROM HIRES EXPERIMENT

1 INTRODUCTION. Figure 2: Synoptical situation at the beginning of the simulation: 5th January 1999 at 12UTC.

E. Agu, M. Kasperski Ruhr-University Bochum Department of Civil and Environmental Engineering Sciences

The Influence of Ocean Surface Waves on Offshore Wind Turbine Aerodynamics. Ali Al Sam

Are Advanced Wind Flow Models More Accurate? A Test of Four Models

Super-parameterization of boundary layer roll vortices in tropical cyclone models

On the temporal evolution of the vertical momentum fluxes within the Marine Atmospheric Boundary Layer

Thorsten Mauritsen *, Gunilla Svensson Stockholm University, Stockholm, Sweden

2.4. Applications of Boundary Layer Meteorology

On the Interpretation of Scatterometer Winds near Sea Surface Temperature Fronts

Wind and Wind Stress Measurements in HiRes

Investigation of Vertical Wind Shear Characteristics Using 50m Meteorological Tower Data

CHANGE OF THE BRIGHTNESS TEMPERATURE IN THE MICROWAVE REGION DUE TO THE RELATIVE WIND DIRECTION

Refined Source Terms in WAVEWATCH III with Wave Breaking and Sea Spray Forecasts

10.6 The Dynamics of Drainage Flows Developed on a Low Angle Slope in a Large Valley Sharon Zhong 1 and C. David Whiteman 2

Wave-Phase-Resolved Air-Sea Interaction

Determination of the Spatial Variation of the Atmosphere and Ocean Wave Fields in Extremely Light Wind Regimes

Development of a Regional Coupled Ocean-Atmosphere Model

Real Life Turbulence and Model Simplifications. Jørgen Højstrup Wind Solutions/Højstrup Wind Energy VindKraftNet 28 May 2015

A Comparison of the UK Offshore Wind Resource from the Marine Data Exchange. P. Argyle, S. J. Watson CREST, Loughborough University, UK

UNDERSTANDING NEAR-SURFACE AND IN-CLOUD TURBULENT FLUXES IN THE COASTAL STRATOCUMULUS-TOPPED BOUNDARY LAYERS

MODELING OF CLIMATE CHANGE IMPACTS ON COASTAL STRUCTURES - CONTRIBUTION TO THEIR RE-DESIGN

RESOURCE DECREASE BY LARGE SCALE WIND FARMING

LONG TERM OCEAN WAVE FORECASTING ALONG INDIAN COAST

Wind resource assessment over a complex terrain covered by forest using CFD simulations of neutral atmospheric boundary layer with OpenFOAM

Wind Flow Modeling: Are computationally intensive models more accurate?

11.4 WIND STRESS AND WIND WAVE OBSERVATIONS IN THE PRESENCE OF SWELL 2. METHODS

Evaluation of four numerical wind flow models

Wind Project Siting & Resource Assessment

Validation Study of the Lufft Ventus Wind Sensor

Review of Equivalent Neutral Winds and Stress

DIRECTION DEPENDENCY OF OFFSHORE TURBULENCE INTENSITY IN THE GERMAN BIGHT

Chapter 2. Turbulence and the Planetary Boundary Layer

Low level coastal jet

3.3 USING A SIMPLE PARCEL MODEL TO INVESTIGATE THE HAINES INDEX

Figure 1 Lake Ontario Offshore Study Area near East Toronto

Wind Stress Working Group 2015 IOVWST Meeting Portland, OR

Undertow - Zonation of Flow in Broken Wave Bores

WATER TRANSPORT VARIABILITY IN THE AEGEAN SEA AND ITS CONNECTION WITH NORTH SEA CASPIAN PATTERN (NCP)

The Air-Sea Interaction. Masanori Konda Kyoto University

Short-term Spectral and Statistical Analysis of Sea Surface Elevation Data from buoys located in the Greek Seas

Importance of thermal effects and sea surface roughness for offshore wind resource assessment

Abstract. 1 Introduction

ABNORMALLY HIGH STORM WAVES OBSERVED ON THE EAST COAST OF KOREA

Goal: Develop quantitative understanding of ENSO genesis, evolution, and impacts

Mesoscale air-sea interaction and feedback in the western Arabian Sea

An experimental study of internal wave generation through evanescent regions

SIO 210 MIDTERM, 26 October 2009.

Effects of very harsh 2012 winter conditions on the deep layer of the Southern Adriatic Sea

ASSESSMENT OF SEA BREEZE CHARACTERISTICS FROM SODAR ECHOGRAMS

Validation of Measurements from a ZephIR Lidar

Next Generation Modeling for Deep Water Wave Breaking and Langmuir Circulation

Meteorological Measurements OWEZ

Surface Wind Speed Distributions: Implications for Climate and Wind Power

Wind shear and its effect on wind turbine noise assessment Report by David McLaughlin MIOA, of SgurrEnergy

Undertow - Zonation of Flow in Broken Wave Bores

FLUXES ESTIMATION AND THE DERIVATION OF THE ATMOSPHERIC STABILITY AT THE OFFSHORE MAST FINO1.

Meteorological Measurements OWEZ

Observations and Modeling of Coupled Ocean-Atmosphere Interaction over the California Current System

Stefan Emeis

Observed directional characteristics of the wind, wind stress, and surface waves on the open ocean

Atmospheric Forcing and the Structure and Evolution of the Upper Ocean in the Bay of Bengal

Rokjin J. Park, Jaein I. Jeong, Minjoong Kim

Surface Waves NOAA Tech Refresh 20 Jan 2012 Kipp Shearman, OSU

Observations of Velocity Fields Under Moderately Forced Wind Waves

Atmospheric Rossby Waves in Fall 2011: Analysis of Zonal Wind Speed and 500hPa Heights in the Northern and Southern Hemispheres

3D Turbulence at the Offshore Wind Farm Egmond aan Zee J.W. Wagenaar P.J. Eecen

The Use of Bulk and Profile Methods for Determining Surface Heat Fluxes in the Presence of Glacier Winds

1. Atmospheric Diffusion of Stack Gases

MET Lecture 8 Atmospheric Stability

Mass coral mortality under local amplification of 2 C ocean warming

Profiles of Wind and Turbulence in the Coastal Atmospheric Boundary Layer of Lake Erie

8.4 COASTAL WIND ANOMALIES AND THEIR IMPACT ON SURFACE FLUXES AND PROCESSES OVER THE EASTERN PACIFIC DURING SUMMER

Importance of thermal effects and sea surface roughness for offshore wind resource assessment

Experience with the representation of stable conditions in the ECMWF model

AIS data analysis for vessel behavior during strong currents and during encounters in the Botlek area in the Port of Rotterdam

Modelling atmospheric stability with CFD: The importance of tall profiles

from a decade of CCD temperature data

ValidatingWindProfileEquationsduringTropicalStormDebbyin2012

Effect of Orography on Land and Ocean Surface Temperature

Atmospheric Waves James Cayer, Wesley Rondinelli, Kayla Schuster. Abstract

6.28 PREDICTION OF FOG EPISODES AT THE AIRPORT OF MADRID- BARAJAS USING DIFFERENT MODELING APPROACHES

Are Hurricanes Becoming More Furious Under Global Warming?

CROSS-SHORE SEDIMENT PROCESSES

Impact of fine-scale wind stress curl structures on coastal upwelling dynamics : The Benguela system as a case of study.

JCOMM Technical Workshop on Wave Measurements from Buoys

Anemometry. Anemometry. Wind Conventions and Characteristics. Anemometry. Wind Variability. Anemometry. Function of an anemometer:

Transcription:

CEST213 Athens, Greece Ref no: XXX EXPERIMENTAL STUDY OF THE SURFACE MARINE ATMOSPHERIC BOUNDARY LAYER OVER AEGEAN SEA, GREECE KOSTOPOULOS* V. E. and HELMIS* C. G. *Department of Environmental Physics and Meteorology, Faculty of Physics, University of Athens, University Campus, Zografou, 15784, Athens, Greece, vkostopoulos@phys.uoa.gr EXTENDED ABSTRACT Micrometeorological measurements within the surface Marine Atmospheric Boundary Layer (MABL) took place at the shoreline of island of Karpathos, at the South East end of the Aegean Sea (Karpathio Pelagos), during summer 212. Α meteorological mast was installed at a distance less than 3m from the shoreline, instrumented with a sonic anemometer and a fast krypton hygrometer at 14.5m height with a sampling frequency of 1Hz, in order to measure the turbulent transport of mass and energy of the surface MABL. At the same mast slow response sensors (1Hz) were also installed providing vertical profiles of wind speed and direction, temperature and humidity at three levels (3, 8 and 12 m). Satellite sea surface temperature images were also obtained. The spacial and temporal scales of the Aegean Pelagos have not been studied in the past while it is expected to present differences compared with those measured over oceans or open seas. It is worth to mention that this study is a follow-up of a previous study which took place during summer 211 at Skyros island, northern Aegean provided similar measurements but with different fetch. During Karpathos campaign, more than of 12 hours stationary surface MABL observations of near neutral wind flows were recorded. The vertical profiles of temperature reveal the presence of a coastal internal boundary layer, mostly below 8m height. Eddy correlation analysis is used to study the turbulent fluxes of mass and energy. According to first results on the stable MABL, momentum fluxes depend on stability and give increasing values with the wind speed. The estimated values are higher, almost by a factor of two, compared to typical values measured over the ocean. These increased values could be attributed to a developed sea state since wind stress is greater over young and developing wave fields. Sensible heat flux is found to be increased with stability and rather independent from wind speed while moisture flux seems to not be dependent neither to stability or wind speed. Some part of the records present upward heat fluxes under near neutral conditions and mostly for positive sea-air temperature differences, possibly associated with sea spray or mesoscale heterogeneity of the sea surface temperature. Ongoing investigation of the turbulent exchanges is expected to provide a more thorough insight on the coupling processes between the surface MABL and Aegean sea. Keywords: surface marine ABL, turbulent fluxes, Aegean Please do not use page numbers

1. INTRODUCTION Air-sea interaction is known to strongly influence the world s climate. The complex nature of the physical processes involved is still not fully understood although numerous studies addressed this problem over the past decades. Flux measurements of mass and energy within the surface Marine Atmospheric Boundary Layer (MABL) are of fundamental importance in order to understand the structure and the characteristics of the MABL and considered essential for weather and climate predictions since they provide modelers with parameterization schemes for surface turbulent fluxes of momentum, heat and water vapor. Up today, the vast majority of such measurements were conducted mostly over open seas and oceans. The semi-closed sea of Aegean Pelagos consisting of different spatial and temporal scales than the ocean (Sofianos et al. 22) is thought to have a different impact on the exchanges at the air-sea interface. Momentum transfer is accepted as an increasing function of the wind speed and the differences found between independent studies are a function of the sea state (Foreman et al. 21). Short fetched flows wind stress is documented greater over young and developing wave fields compared to older wave fields, which are in more equilibrium with the wind. Heat and moisture exchanges on the other hand were treated up to recent without including the contribution of the sea spray which is found to significant augment the interfacial transfers (Andreas 29), especially at high winds or of effects of mesoscale heterogeneity of the sea surface temperature (SST), mostly for weak wind cases (Mahrt et al. 212). Evidence of such contributions, which was found in our data, is presented. The purpose of this paper is to present the conducted experiment as well as results regarding the turbulent transfer over the Aegean Sea. This study is a follow-up of a previous study, which took place during summer 211 at Skyros Island, at the northern Aegean sea, which gave similar measurements but with different fetch (Kostopoulos et al. 212). It is worth to mention that there are no other published measurements over the Aegean. 2. EXPERIMENTAL SITE AND INSTRUMENTATION Micrometeorological measurements of the surface MABL, took place at the shoreline of island of Karpathos, at the South East of the Aegean Sea (Karpathio Pelagos), during summer 212. The experimental campaign lasted almost two weeks from the end of July to August, under mostly etesian flows. Karpathos Island is located at the southeastern Aegean approximately 15km from the islands of the center Aegean. Apart from the south end which is relatively flat and has mostly gentle slope shoreline (<1m), the island is characterized by intense orography with dense forests. Mountains extend to about 6m at the north side and over 1km height (1.22m) at the center of the island. The experimental site as well as the placement height of the fast response sensors was carefully chosen so that the measurements would not be affected from the developed Internal Boundary Layer (IBL) over land. The selection of the mast location was based on the fact that the developed IBL, at a distance of 3 m from the shoreline, will be characterized with heights much less than 1 m (Savelyev et al. 24). Α meteorological mast (figure 1) was installed at the south part of the island, at a distance less than 3m from the shoreline, instrumented with a sonic anemometer and a fast krypton hygrometer at 14.5m height with a sampling frequency of 1Hz, in order to measure the turbulent transport of mass and energy of the surface MABL. At the same mast slow response sensors (1Hz) were also installed providing vertical profiles of wind speed and direction, temperature and relative humidity at three levels (3, 8 and 12 m). The calculation of the turbulent fluxes is conducted using eddy correlation analysis with averaging time of

1min. Wind speed at 1m is produced using similarity theory wind profile equations (Large et al.198). More than 12 hours of observations of near neutral and stable marine wind flows were recorded. Data was initially treated with quality control procedures (Lee et al. 22) in order to find and exclude instrumentation errors. Less than 1-3 % of corrections over about 1-1 % of the total records were applied. Tilt corrections due to the deviation of the mast from true vertical where also made using the planar fit method (Wilczak et al. 21). Records were furthermore examined for stationarity (Mahrt et al. 1996). Finally additional data of daily SST images from the satellite INMARSAT-C were collected. Figure 1. The meteorological mast (facing North) at the experimental site (left) located at the south-west end of Karpathos island (right). 3. MEASUREMENTS OF TURBULENT FLUXES OVER THE AEGEAN SEA Moderate to strong etesian winds were recorded during the experimental period, ranging up to 11m/s while wind direction varied from west to mostly north-west (NW). Only few records with directions blowing within the golf (>345 degrees) were excluded. Temperature fluctuated from 25ºC to 3ºC and relative humidity from 6% to 85%. During the whole experimental period the sky was cloudless. The sign of the temperature gradients with height as well as the stability parameter values recorded from the sonic anemometer at 14.5 m, revealed the presence of a coastal internal boundary layer, mostly below 8m height as expected (figure 2). According to this figure the temperature data at 12 m height are characterized by higher values than the ones at 8 m (due to the marine flow), while the stability parameter reveal a stable layer at 14.5 m which was expected for the same reason. Thus, flux measurements should be considered practically uninfected from the presence of the land and the developed IBL. When taking into account the difference between SST and the potential temperature of the air, several of the recorded heat fluxes are found to strongly correlate with temperature differences. Upward heat fluxes are observed during near neutral conditions and small air-sea temperature differences, mostly positive (figure 3), fluxes that increase with the wind of speed (not shown in this paper) and are possibly related to sea spray conditions or effects of mesoscale heterogeneity of the sea surface temperature. It is worth mentioning that Skyros measurements also revealed several upward heat flux conditions corresponding as well to positive or very small negative air-sea temperature

differences. Such conditions are generally excluded from analysis since relatively small errors in the air-sea temperature difference may strongly contaminate the estimated transfer coefficient. On the other hand under such conditions although sensible heat flux is found to be small (less than 1W/m 2 ), may cover a large fraction of the sea and contribute significantly to the greater area s heat budget thus affecting not only weather but perhaps the local climatology. Analysis results of these records are excluded from this paper and will be addressed in a future study. In this study, records (~2) for negative sea-air temperature differences (<-1 C) are only presented. Figure 2. The vertical profile of temperature and the stability parameter (z/l). m K/s w'tv'.1.5 -.5 -.1 -.15 -.2-3 -2.5-2 -1.5-1 -.5.5 1 SST-Θair (K) Figure 3. Sensible heat flux as a function of SST-θ(z).

Momentum flux decreases with stability probably reflecting the insufficient generation of turbulence needed to overcome the absorption of energy by buoyant forces and gives increasing values with the wind speed as seen in figure 4. The estimated values are higher, almost by a factor of two, compared to typical values measured over the ocean or even when compared with Skyros different fetch results (Kostopoulos et al. 212). It is known that the fetch and the duration of an air flow over the sea plays crucial role to the development of the wave field on which momentum flux also depends on (Hasselmann et al.1973). In this case such values should be attributed to the special spatial and time scales found over Kapathio Pelagos. u'w' (m 2 /s 2 ) -.45 -.4 -.35 -.3 -.25 -.2 -.15 -.1 -.5 4 5 6 7 8 9 1 11 12 U 1 (m/s) u'w' (m 2 /s 2 ) -.45 -.4 -.35 -.3 -.25 -.2 -.15 -.1 -.5.5 1 1.5 2 2.5 3 Figure 4. Momentum flux towards wind speed and stability. Heat and moisture fluxes do not seem to follow the same consequence although appear enhanced compared to Skyros results corresponding to larger air-sea temperature differences. According to figure 5 sensible heat flux is found to be independent from the wind speed even though it increases with turbulence as expected (not shown). When plotted against negative air-sea temperature differences (figure 6), sensible heat flux is seen normally to increase with increasing difference. These values are close to the annual mean net heat flux which is estimated for the Aegean sea at 26 W/m 2 (Poulos et al., 1997) and which implies that the Aegean Sea is, on the average, losing heat through its surface. z/l

w't' (mk/s) -.7 -.6 -.5 -.4 -.3 -.2 -.1 4 5 6 7 8 9 1 11 12 U 1 (m/s) Figure 5. Sensible heat flux towards wind speed w't' (mk/s) -.6 -.5 -.4 -.3 -.2 -.1 w'q' (g h2o m / Kg air s).25.2.15.1.5-5 -4-3 -2-1 SST-Θ(z) (K) Figure 6. Sensible heat and specific humidity flux as a function of SST-θ(z) -.5 Moisture fluxes (figure 7) appear with rather typical positive values ranging from 1 1 to 1 2 W/m 2, independent from the wind speed, air-sea temperature difference or stability, except perhaps for increased stability conditions where it seems to slightly decrease possibly as a result of the reducing in turbulence. Further analysis is in progress in order to investigate the physical processes and the influences of various meteorological and sea state factors to the turbulent transport of momentum and energy in the surface MABL of the Aegean.

w'q' (g v m / Kg air s).3.2.1.5 1 1.5 2 2.5 3 -.1 z/l w'q' (g v m / Kg air s).3.2.1 4 5 6 7 8 9 1 11 12 -.1 4. CONCLUSIONS U 1 (m/s) Figure 7. Specific humidity flux towards wind speed and stability Onshore turbulent fluxes measurements of the surface MABL were conducted at the southern east Aegean sea during summer 212. This study is a follow-up of a recent study which took place during summer 211 at Skyros island in the northern Aegean. Upward heat fluxes were found on both campaigns measurements for near neutral conditions with positive sea-air temperature differences. Such processes which are possibly related to sea spray or mesoscale heterogeneity of the sea surface temperature are still not fully understood. Although these fluxes are measured and also documented small (mostly less than 1W/m 2 ) they could contribute significantly to the heat budget over the Aegean sea where the annual mean estimation of net heat flux is 26W/m 2. For negative sea-air temperature differences, momentum transfer presents enhanced relative to typical values measured over open sea. Such values could be attributed to the special spatial and time scales found over Kapathio Pelagos. Moisture fluxes are found positive

and on values ranging from 1 1 to 1 2 W/m 2. Sensible heat fluxes are downward and range from 1 to 1 1 W/m 2. Momentum and moisture fluxes decrease with stability while heat flux increases with augmenting scatter. Heat fluxes are seen to act as a function of the measured air-sea temperature difference. Ongoing investigation of the turbulent structure and the energy transfer in the recorded surface MABL is expected to provide further understanding on the relative processes over an Archipelago such as the Aegean sea. ACKNOWLEDGMENTS This research has been co-financed by the European Union (European Social Fund ESF) and Greek national funds through the Operational Program "Education and Lifelong Learning" of the National Strategic Reference Framework (NSRF) - Research Funding Program: Heracleitus II, Investing in knowledge society through the European Social Fund. REFERENCES 1. Kostopoulos V., Helmis C., Raptis I., Experimental Study of the Turbulent Structure of the Surface Marine Atmospheric Boundary Layer over the Aegean Pelagos Under Etesian Winds, 213, Advances in Meteorology, Climatology and Atmospheric Physics Springer Atmospheric Sciences, 213, pp 163-168 2. Mahrt L, Vickers D. Andreas E. Khelif D., Sensible Heat Flux in Near-Neutral Conditions over the Sea, 212, Journal of Physical Oceanography, vol.42, p 1134 1142. 3. Mahrt L. and D. Khelif, Heat fluxes over weak SST heterogeneity, 21, Journal of Geophysical Research, vol. 115, D1113, doi:1.129/29jd13161. 4. Andreas, Edgar L., 21: Spray-Mediated Enthalpy Flux to the Atmosphere and Salt Flux to the Ocean in High Winds. J. Phys. Oceanogr., 4, 68 619. doi: http://dx.doi.org/1.1175/29jpo4232.1 5. Foreman, Richard J., Stefan Emeis, 21: Revisiting the Definition of the Drag Coefficient in the Marine Atmospheric Boundary Layer. J. Phys. Oceanogr., 4, 2325 2332. doi: http://dx.doi.org/1.1175/21jpo442.1 6. Savelyev S, Taylor P (24) Internal Boundary Layers: I. Height Formulae for Neutral and Diabatic Flows. Bound Layer Meteorol 115:1-25. doi: 1.17/s1546-4-2122-z 7. Lee X et al (22) Handbook of Micrometeorology: A Guide for Surface Flux Measurement and Analysis, Kluwer, Dordrecht 8. Sofianos S, Johns W, Lascaratos A, Murray S, Olson D, Theocharis A (22) Draft Report of the Aegean Sea Workshop. http://www.oc.phys.uoa.gr/workshop/aegean_draft_report_f.htm 9. Wilczak JM, Oncley SP, Stage SA (21) Sonic anemometer tilt correction algorithms. Bound Layer Meteorol 99:127-15. doi:1.123/a:11896624465 1. Poulos, S. E., P. G. Drakopoulos, M. B. Collins, Seasonal variability in sea surface oceanographic conditions in the Aegean Sea (Eastern Mediterranean): an overview, J. Mar. Sys., 13, 225-244, 1997 11. Mahrt L, Vickers D, Howell J, Hojstrup J, Wilczak JM, Edson J, Hare J (1996) Sea surface drag coefficients in the Riso Air Sea Experiment. J Geophys 11:14327 14335. doi:1.129/96jc748 12. Smith SD, Fairall CW, Geernaert GL, Hasse L (1996) Air-sea fluxes: 25 years of progress. Bound Layer Meteorol 78:247-29. doi: 1.17/BF12938 13. Andreas E, Edson J., Monahan E, Rouault M., Smith D., The spray contribution to net evaporisation from the sea: a review of recent progress, 1995, Boundary-Layer Meteorology, Volume 72, pp 3-52. 14. Large WG, Pond S (198) Open sea momentum flux measurements in moderate to strong winds. J Phys Oceangr 11:324-336, doi: http://dx.doi.org/1.1175/152-485(1981)11<324:oomfmi>2..co;2 15. Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP), Hasselmann, K., and Coauthors, 1973: Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Dtsch. Hydrogr. Z., vol. 12, pp 1 95.