• Nie Znaleziono Wyników

Field observations of tidal flow separation at a mega-scale beach nourishment

N/A
N/A
Protected

Academic year: 2021

Share "Field observations of tidal flow separation at a mega-scale beach nourishment"

Copied!
14
0
0

Pełen tekst

(1)

FIELD OBSERVATIONS OF TIDAL FLOW SEPARATION AT A MEGA-SCALE BEACH NOURISHMENT

MAX RADERMACHER1, WILMAR ZEELENBERG1, MATTHIEU DE SCHIPPER1, AD RENIERS1,2

1. Faculty of Civil Engineering and Geosciences, Delft University of Technology, Stevinweg 1, 2628CN Delft, The Netherlands. M.Radermacher@tudelft.nl.

2. Deltares, Boussinesqweg 1, 2629HV Delft, The Netherlands.

Abstract:

Large beach nourishments are often considered an efficient way to deal with coastal zone management issues, notably coastal erosion. Such large nourishments can act as a geometric perturbation, which might cause the tidal current to detach from the coastline. The Sand Motor is a mega-scale beach nourishment in the Netherlands. Field observations of the tidal flow field north of the Sand Motor are presented, showing evidence of flow separation and eddy formation during the flood phase of the tide. GPS-tracked drifters are transported offshore upon deployment near the tip of the nourishment. Flow convergence measured with nearshore current profilers is indicative of flow reversal at the lee side. Images of the temporal development of coherent patterns in time-averaged X-band radar backscatter align with in-situ observations of the flow field. Tidal flow separation at such large nourishments might have an impact on morphological development of the nourishment, as well as on swimmer safety.

Introduction

The impact of sand nourishments on nearshore hydrodynamics is a highly relevant, yet complex topic. The efficiency of a nourishment project depends on the observed and desired spatio-temporal spreading of the newly placed sand, which is governed by nearshore hydrodynamics (Hamm et al., 2002). If the nourished beach represents important recreational values as well, the safety of bathers and swimmers in relation to nearshore waves and currents is usually closely monitored by coastal managers, engineers and the general public alike (De Zeeuw et al., 2012).

Presently, possibilities for upscaling of nourishment volumes are being explored with the Sand Motor mega-nourishment project (Stive et al., 2013; also referred to as Sand Engine, see Figure 1) near The Hague in the Netherlands. A ~20 Mm3 hook-shaped peninsula of roughly 2 km long (alongshore) and 1 km wide (cross-shore) was constructed in 2011. The Sand Motor is expected to nourish the adjacent coastline throughout the coming decades without further human interference. Initially, the nourishment had a rather blunt, non-streamlined shape. As a result of large alongshore sediment transport gradients, by 2014 the shape of the nourishment has become more bell-shaped and less pronounced.

(2)

Fig. 1. Aerial pictures of the Sand Motor just after completion in July 2011 (left) and in a later stage of morphological development in September 2014 (right). In the left picture, suspended sediment

patterns around the tip of the mega-nourishment suggest tidal flow separation. Photography by Rijkswaterstaat / Joop van Houdt.

The presence of this protrusive coastline feature interacts with the tidal currents. The Delfland coast is a micro-tidal environment, with the tidal range varying between approximately 1.4 m and 1.8 m during a spring-neap cycle. The tide in the southern North Sea is a counter-clockwise rotating Kelvin wave, driving tidal currents in alongshore direction. Depth-averaged tidal currents are in the order of 0.6 m/s, with a phase difference of approximately 0° between the horizontal and vertical tide (Wijnberg, 2002). In the nearshore, the tidal currents contract around the Sand Motor. This paper aims to determine the presence of tidal flow separation at the mega-nourishment. If the restoring forces are not strong enough to keep the contracting streamlines attached to the coastline, the tidal current separates from the shoreline, creating a downstream region with large horizontal velocity shear discussed next.

Flow separation occurs when a boundary layer detaches from its associated wall and has significant consequences for the free-stream flow (Simpson, 1989). Generally, detachment of the boundary layer is caused by an adverse pressure gradient, which disturbs the balance of forces that constitutes a boundary layer. The flow near the wall changes sign at the separation point and forms a so-called return flow. In the case of contracting flow around a geometric disturbance in the wall, an adverse pressure gradient can be caused by flow expansion and deceleration in the lee of the disturbance. The context of tidal currents in coastal waters adds an important aspect to the general case of flow separation. The periodicity of tidal currents makes these flows unsteady. Unsteadiness does not influence the character of the separating flow if the timescale of the unsteady, organized motion is much larger than the timescale of energy-containing turbulent structures (Simpson, 1989). For the case of tidal flow around a large shoreline perturbation, the length scale of the tide (the tidal excursion) can be of the same order as the length scale at which turbulent energy is initialized (the size of the perturbation). Hence unsteadiness should be taken into account.

(3)

Signell and Geyer (1991) have put these processes in the right framework with their study of eddy formation at coastal headlands. The interplay between advection, friction and local acceleration, together with the headland geometry, determine the character of the tidal flow around the headland. When advection is not the dominant forcing mechanism, i.e. in case of large friction or in case of a wide headland compared to the tidal excursion, flow separation does not occur. In all other cases (i.e. with dominant advection), a region of strong vorticity develops in the lee of the headland, giving rise to the formation of a large-scale eddy. The periodicity of tidal flow and the influence of friction in shallow coastal waters avoid the formation of multiple, subsequent eddies that are released from the shoreline (Signell and Geyer, 1991), as was also observed in the field by Geyer and Signell (1990) at a headland in Vineyard Sound, Massachusetts.

At the Sand Motor, another factor might play a role, which is not treated in the framework of Signell and Geyer. Approximately 10 km South of the mega-nourishment, the fresh water runoff of the river Rhine is discharged into the North Sea, see Figure 2. This adds significant density differences to the system (De Boer, 2009). In the vicinity of the river mouth, fresh water is discharged in a pulsating fashion, with large outflow during low water (LW) and limited outflow during high water (HW). The near-surface freshwater discharge forms a buoyant lens, that is advected with the tidal currents. During flood, when the tidal flow is northward, this lens regularly approaches the coastline near the Sand Motor as a freshwater front (Souza and Simpson, 1997).

The Sand Motor is a smaller and more streamlined coastline perturbation than headlands and other geometries studied in literature (Geyer and Signell, 1990; Wolanski et al., 1984). This has raised questions about the presence of flow separation and tidal eddies at the Sand Motor, as flow separation might influence the morphological development of the nourished sand and could potentially affect the safety of swimmers and bathers. This research presents field observations of tidal flow separation at the Sand Motor. Measurements were conducted during the 6-week Mega Perturbation Experiment (MegaPEX) in fall 2014. Using data of one particular tidal cycle, an assessment is made of the measured tidal flow patterns and the occurrence of tidal flow separation is determined. The next section provides more information on the methods used in the field. It is followed by an overview of the results and a discussion on the complex interaction between the tidal currents and buoyancy forcing.

(4)

Fig. 2. Overview of the measurement location. The Sand Motor is located 10 km north of the mouth of the Rhine river. The red dot in the upper left panel indicates the location of the X-band radar station. Measurements were collected using a nearshore grid of current profilers, out of which two are used in the research presented here. Colors in the upper right panel represent bed level in m with

respect to the Dutch datum (~ Mean Sea Level).

Field methods

The MegaPEX field campaign took place at the Sand Motor between September 15th and October 26th of 2014. Its general aim was to acquire a multi-thematic dataset of all relevant processes affecting the Sand Motor’s morphological, ecological and recreational development. In this study, use is made of data obtained with GPS-tracked drifters, a fixed nearshore array of current profilers and an X-band radar station. Figure 2 gives a spatial overview of the measurement locations.

A set of 13 GPS-tracked drifters was used to perform Lagrangian current measurements. The drifters, modified after Schmidt et al. (2003) and MacMahan et al. (2009), consisted of a PVC tube for buoyancy with a circular plate and a

(5)

weight attached to the bottom to avoid surfing on the bores of broken waves. The GPS location of the drifters was tracked with a Locosys GT-31 handheld GPS data logger. Drifters were deployed such that their tracks would indicate whether flow separation was present (main flow detaching from the shore, eddy formation in the lee) or not (streamlines stay parallel to the coastline). Altogether, 15 separate drifter sessions were performed, 8 of which took place during the flood period of the tide and 7 during the ebb period. Per session, all available drifters were deployed one or several times depending on the visually observed drifter tracks. Sessions were started around the HW or LW times, as the process of flow separation, if present, would be fully developed by then (Signell and Geyer, 1991).

Nearshore currents were measured using an array of 10 Acoustic Doppler Current Profilers (ADCP’s), see Figure 2. The array consisted of two cross-shore transects (transect A at the tip of the Sand Motor and transect B along the northern edge) and one alongshore transect in between (transect R). The array was positioned such that a potential separation point was likely to be located in between transects A and B, which would enable identification of flow separation from different flow characteristics in both transects. Raw ADCP data were sampled at 1 Hz. All ADCP data presented in this study have been filtered with a uniform 15-minute window as to remove the signature of surface (infra-)gravity waves and were resampled at 1-minute intervals.

As a part of the ongoing Sand Motor monitoring campaign, an X-band radar station was deployed in the dunes near the village of Kijkduin, see Figure 2. The radar produced images of measured backscatter in space and time, which, for the marine half of the Sand Motor radar domain, represents so-called sea clutter (Bell, 1999). Through Bragg scattering of the radar signal on the sea surface, the radar effectively provides a measure of sea surface roughness (Dankert et al., 2003). Time averaging is commonly used to remove noise and reveal spatial (flow) patterns in radar images (e.g. Dankert et al., 2003; Haller et al., 2013). At the Sand Motor, 5-minute time-averaged X-band radar images revealed signatures of the incoming fresh water lens, as well as spatial patterns associated with flow separation. These qualitative measures were combined with the in-situ measurements in order to gain confidence in the radar signal and to create a synoptic picture of the flow patterns at tidal scales around the Sand Motor.

Results

Here, a subset of the full dataset is presented, concerning measurements on October 1st, 2014. Detailed analysis of processes happening during a single tidal cycle is necessary before any general trend in the full dataset can be pinpointed. On that day, a 5 m/s wind was blowing from SSW. Offshore measurements

(6)

show a significant wave height of 0.8 m, coming from SW. These mild wind and wave conditions are expected not to drive strong currents in the region of interest and will therefore not have a significant influence on the data presented in this study.

Fig. 3. Drifter tracks as recorded by the GPS loggers. The part of the track that was recorded after

arrival of the fresh water front is shown in red. Circles mark the end point of every drifter track. Background colors represent bed level in m with respect to chart datum.

First, the drifter data are presented. All 13 drifters were deployed in between the A5 and B3 ADCP stations (Figure 2) approximately 1 hour after the HW peak. The recorded drifter tracks are presented in Figure 3. Upon deployment, most drifters move obliquely to the general flood direction due to a steady offshore flow component. Recorded velocities range between 0.4 and 0.5 m/s during that phase. Four drifters that were released closest to station B3 are trailing behind, of which two even remain stagnant in the nearshore. Three drifters that were deployed at the same depth contour as the stagnant ones, but closer to station A5, are moving along with the main cluster of drifters. This implies the presence

(7)

of strong spatial gradients in the flow field. Approximately 40 minutes after deployment, the heading of all drifters changes to become more aligned with the flood direction. The velocity magnitude of the offshore cluster of drifters increases to 0.8 m/s. Between 60 and 75 minutes after deployment, most drifters are withdrawn from the water using a personal watercraft. Only the two drifters in the offshore cluster that covered the longest distance remained afloat, together with the two nearshore drifters.

The offshore velocity components that appear from the drifter data differ from an ordinary, alongshore, non-separating flow field. As this pattern cannot be explained by wind or wave forcing, separation of the tidal flow is likely to be the cause. Eulerian velocity measurements obtained from ADCP stations A5 and B3 are used now to complement that image. Both stations are located around the 4 m depth contour, which corresponds to the seaward slope of the sub-tidal bar. Filtered flow velocities are presented in Figure 4. The velocity signal has been decomposed into a cross-shore (u, first panel) and an alongshore (v, second panel) component, based on the local coastline orientation at every station. The shaded area marks the period of the drifter deployment. Alongshore flow velocities (second panel) at both stations contain a dominant signal at the tidal timescale. During the flood period, alongshore velocities at A5 and B3 rise simultaneously until 1h20 before HW. After that time, v at B3 steadily decreases, whereas at A5 it reaches its peak value 35 minutes before HW. For both stations, cross-shore u is increasing in that period and eventually turns offshore. From 1 hour after HW onwards, alongshore velocities at B3 are negative, which indicates reversal of the flow at that location. This suggests the presence of flow conversion and possibly a flow separation point in between both stations, which vanishes 2h30 after HW. During this period, u peaks in offshore direction at A5. After slack tide, which occurs approximately 3h30 after HW, v reverses at both stations and the ebb phase commences. Flow patterns during ebb are rather straightforward, although a distinct offshore component is present at A5.

The measured flow patterns that were described above show signs of tidal flow separation, but are too localized in space to reveal the full tidal flow pattern in the lee of the Sand Motor. Furthermore, the influence of the fresh water runoff from the Rhine river remains unclear. These aspects can be addressed using X-band radar data, albeit in a more qualitative sense. Figure 5 shows the 5-minute averaged backscatter as received by the radar at 9 half hour intervals between 6h30 (1h20 before HW, marked with the first black, vertical line in Figure 4) and 10h30 (2h40 after HW). The Sand Motor is visible as a dark patch in the lower left corner of the images. Numerous bright spots, both on land and in the water, correspond with signposts, buildings, demarcation buoys, measurement buoys and watercraft. From 6h30 onwards, a white, curved streak develops. It

(8)

Fig. 4. Flow velocity and pressure as measured at ADCP stations A5 and B3.Currents in offshore and flood direction have a positive sign. Time is given in UTC+2. The shaded area marks the drifter deployment. The left black vertical line marks the initiation of flow separation in the x-band images. The right black vertical line marks the arrival of the fresh water front near the tip of the Sand Motor.

(9)

starts near the most seaward point of the Sand Motor, gradually grows longer and moves slightly in offshore direction. Shoreward of this line, coherent clockwise rotating patterns can be distinguished every now and then (only visible in a movie with 5 minute intervals). It is thought to be related to tidal flow separation: the main flow is located offshore of the bright streak, whereas the in-shore region is dominated by the development of a large-scale eddy. This supposition aligns with the facts that (1) the ADCP stations are situated on both sides of the separation line and have alongshore velocity components with opposite sign from 8h30 onwards and (2) the separation line starts to develop around 6h30, which is the same moment when the alongshore flow velocity measured in station B3 starts to decay.

At 9h00, a second white streak propagates shoreward into the radar domain. This streak marks the freshwater front associated with the outflow of the Rhine river, establishing a stratified regime with a freshwater layer on top of a saline layer. As soon as the front propagates past the separation streak, the latter disappears. The second black, vertical line at 9h25 in the time series of Figure 4 marks the moment at which the front is in the vicinity of station A5. From this moment onwards, the decaying alongshore velocities in both stations start to increase again in flood direction and reach a second local maximum before the ebb flow phase finally commences. The parts of the drifter tracks in Figure 3 that are affected by the front passage (> 9h25) have been given a red color. Especially the nearshore drifters show an abrupt change of heading and an increase in speed. The change of heading of the offshore drifters is more gradual. Visualization of the drifter tracks in space and time together with the radar images (not shown here) confirms that all drifters are first transported offshore along the flow separation line and subsequently alongshore and onshore with the front line. After entering the recirculation zone, the front gradually deforms and curls up, which is visible in the last radar image.

Discussion

The results presented in the previous section show traces of tidal flow separation and eddy creation downstream of the Sand Motor during the flood phase. The X-band radar images play a crucial role in showing that the ADCP and drifter measurements are part of a system with larger length scales than the in-situ measurement domain. The ADCP stations are positioned such that they capture flow reversal, which is a very characteristic feature of the separating flow field. Radar images illustrate that part of the time the separation point is positioned in between both ADCP’s, resulting in opposing alongshore velocities at both stations. The recirculation zone gradually expands further offshore and closer to the most seaward point of the Sand Motor. This is in line with the tidal flow patterns observed around a headland by Geyer and Signell (1990). By the time

(10)

Fig. 5. Time-averaged (5 min.) images of X-band radar. The two ADCP stations used in this research are indicated with red markers.

of front arrival, the cross-shore extent of the recirculation zone is approximately twice as wide as the cross-shore dimension of the Sand Motor.

It remains unclear what actually causes the bright flow separation streak in the time-averaged radar backscatter. For the freshwater front it is known from visual observations that it is associated with a slightly rougher water surface and floating foam around the interface. Which properties of the streamline that divides the mean flow and the recirculating flow in the lee of the Sand Motor lead to such high backscatter values in the X-band signal is presently unknown.

(11)

Possibly this is related to density differences between the fresher nearshore and more saline offshore water masses due to the presence of the ROFI (e.g. De Boer, 2009).

In the present study, only the flow field north of the Sand Motor has been discussed. It shows that tidal flow separation can develop at the lee side of a large beach nourishment. This might influence the spreading of the nourished sand, as suggested for headleands by Pingree (1978) and Geyer and Signell (1990). Observed grain size distributions around the Sand Motor show two patches of relatively fine sediment (Huisman et al., these proceedings) that coincide with the observed locations of the tidal eddies. Furthermore, offshore directed flow velocities around the tip of the nourishment, although limited to about 0.5 m/s, might pose a risk to swimmers. The results presented in the previous section raise questions about (1) the frequency of occurrence of flow separation at the northern side of the Sand Motor, (2) the presence of tidal flow separation south of the Sand Motor during ebb flow and (3) the presence of tidal flow separation at other, smaller and less pronounced beach nourishments. These questions will be subject to future research, partly based on the full dataset obtained during the MegaPEX field campaign. The initial topographic design of a (mega-)nourishment is of key importance to the occurrence of tidal flow separation. The results and considerations presented in this study can be taken into account during the design phase of a future nourishment project.

Conclusions

Large beach nourishments can act as a coastline perturbation, giving rise to separation of the tidal flow. The tidal flow field around the Sand Motor mega-nourishment in the Netherlands was measured using GPS-tracked drifters, a nearshore grid of current profilers and X-band radar images. Tidal flow separation was observed at the lee side of the nourishment. During the flood phase of the tide, drifters and nearshore current profilers recorded flow convergence and offshore flow near the tip of the nourishment. This is indicative of flow reversal and the formation of a large-scale tidal eddy. Spatial and temporal development of coherent structures in five-minute averaged X-band radar backscatter images were shown to align with in-situ measurements of the flow field.

The onset of flow separation was observed well before the high water peak. The size of the associated tidal eddy grows steadily over time and gradually extends further seaward, as far as two times the cross-shore width of the Sand Motor. The freshwater lens discharged from the Rhine river plays an important role in the flow field at the tidal timescale. During the second half of the flood phase,

(12)

the lens propagates shoreward. It deforms through interaction with coherent structures in the recirculation zone.

Tidal flow patterns measured at the lee side of the Sand Motor agree with observations around coastal headlands as published in literature. Potentially, the presence of flow separation entails a risk for swimmers and bathers. Furthermore, it might influence the morphological development of the nourishment. An assessment of the tidal flow field over more tidal cycles, south of the Sand Motor and around other beach nourishments could yield important information for future (mega-)nourishment design.

Acknowledgements

MR was supported by STW Grant 12686: Nature-driven Nourishments of Coastal Systems (NatureCoast), S1: Coastal safety. MdS was supported by the ERC-Advanced Grant 291206 - Nearshore Monitoring and Modeling (NEMO). We express our gratitude to the colleagues who helped us out in the field: Sierd, Bonny, Meagan, Martijn, Bas, Marion, Duong, Nick, Teun, Lisa, Jeroen, Emma, Pim, Marijke, Patrick and Carolin. Furthermore, we thank Jamie MacMahan for his valuable advice regarding design and realization of our field setup. The Naval Postgraduate School, Rijkswaterstaat and the University of Miami are greatly acknowledged for the use of their instrumentation.

References

Bell, P. S. (1999). “Shallow water bathymetry derived from an analysis of X-band marine radar images of waves,” Coastal Engineering, 37, 513-527. Dankert, H., Horstmann, J. and Rosenthal, W. (2003). “Ocean wind fields

retrieved from radar-image sequences,” Journal of Geophysical Research, 108(C11), 3352.

De Boer, G. J. (2009). “On the interaction between tides and stratification in the Rhine region of fresh water influence,” PhD-thesis, Delft University of Technology, Delft, The Netherlands.

De Zeeuw, R. C., De Schipper, M. A., Roelvink, D., De Vries, S., and Stive, M. J. F. (2012). “Impact of nourishments on nearshore currents and swimmer safety on the Dutch coast,” Proceedings of the 33rd International

(13)

Geyer, W. R. and Signell, R. (1990). “Measurements of tidal flow around a headland with a shipboard Acoustic Doppler Current Profiler,” Journal of

Geophysical Research, 95(C3), 3189-3197.

Haller, M. C., Honegger, D. and Catalan, P. A. (2013). “Rip current observations via marine radar,” Journal of Waterway, Port, Coastal, and

Ocean Engineering, 140(2), 115-124.

Hamm, L., Capobianco, M., Dette, H. H., Lechuga, A., Spanhoff, R., and Stive, M. J. F. (2002). “A summary of European experience with shore nourishment,” Coastal Engineering, 47, 237-264.

Huisman, B. J. A., Van der Zwaag, J. J., Luijendijk, A. P. and Ruessink, B. G. (2015). “Practical considerations on numerical modeling of sediment sorting at a large scale sand nourishment,” Proceedings of Coastal

Sediments 2015, San Diego, CA.

MacMahan, J., Brown, J. and Thornton, E. (2009). “Low-cost handheld global positioning system for measuring surf-zone currents,” Journal of Coastal

Research, 25(3), 744-754.

Pingree, R. D. (1978). “The formation of the Shambles and other banks by tidal stirring of the seas,” Journal of the Marine Biological Association of the

UK, 58, 211-226.

Schmidt, W. E., Woodward, B. T., Millikan, K. S., Guza, R. T., Raubenheimer, B. and Elgar, S. (2003). “A GPS-tracked surf zone drifter,” Journal of

Athmospheric and Oceanic Technology, 20, 1069-1075.

Signell, R. P. and Geyer, W. R. (1991). “Transient eddy formation around headlands,” Journal of Geophysical Research, 96(C2), 2561-2575.

Simpson, R. L. (1989). “Turbulent boundary-layer separation,” Annual Review

of Fluid Mechanics, 21, 205-234.

Souza, A. J. and Simpson, J. H. (1997). “Controls on stratification in the Rhine ROFI system,” Journal of Marine Systems, 12, 311-232.

(14)

Stive, M. J. F., De Schipper, M. A., Luijendijk, A. P., Aarninkhof, S. G. J. Van Gelder-Maas, C., Van Thiel de Vries, J. S. M., De Vries, S., Henriquez, M., Marx, S., and Ranasinghe, R. (2013). “A new alternative to saving our beaches from sea-level rise: The Sand Engine,” Journal of Coastal

Research, 29(5), 1001-1008.

Wijnberg, K. M. (2002). “Environmental controls on decadal morphologic behavior of the Holland coast,” Marine Geology, 189, 227-247.

Wolanski, E., Imberger, J. and Heron, M. L. (1984). “Island wakes in shallow coastal waters,” Journal of Geophysical Research, 89(C6), 10553-10569.

Cytaty

Powiązane dokumenty

Dotychczasowe spotkania wokół dyskursu autopromocyjnego przyniosły wiele interesujących rozważań i dyskusji, których pokłosiem są trzy publikacje: Dyskurs autopromocyjny i

Typowe dla wyżej scharakteryzowanego procesu zjawiska spoty­ kamy w Galii w VI — V II stuleciu, Widzimy, jak rwą się tam węzły łączące ludność

Szczególnie interesująca i do dziś dnia szeroko dyskutowana jest spra­ w a regenta spartańskiego Pauzaniasza, zwycięzcy, spod Platei, Cypru i Bizancjum.. A le

En ce qui concerne le matériel de guerre à transporter en Pologne, la France est disposée, dans les conditions prévues p ar les accords politiques et

Автор обращает внимание на определяющую роль Европы в этот период, на ее влияние на судьбы других континентов, на совпадение рубежей

Artykuł umieszczony jest w kolekcji cyfrowej bazhum.muzhp.pl, gromadzącej zawartość polskich czasopism humanistycznych i społecznych, tworzonej przez Muzeum Historii Polski

Korzystając z pracy (Zasady… 2009) , w której podano klasyfikację dotyczącą likwido- wanych terenów górniczych ze względu na ograniczenia w wykorzystaniu dla celów

Procesy wzbogacania grawitacyjnego węgla, także prowadzone metodami mokrymi, oprócz usuwania z nich niepożądanych z uwagi na parametry energetyczne zanieczyszczeń,