LimnologicalReview4(2004)17–24 SecularfluctuationsoficephenomenainUpperRaduniaLake, KashubianLakeland DariuszBorowiak,JacekBarańczuk UniversityofGdańsk,DepartmentofLimnology, Dmowskiego16a,80–264Gdańsk Abstract: ThepaperisanattemptatreconstructingmultiannualfluctuationsofthedurationoficephenomenaonUpperRadunia Lake,usingtheirclosedependenceontheindexofwinterNorthAtlanticOscillation(NAO).Againstthesecularchangesthere wereanalysedthefluctuationsofparticularphasesofthedevelopmentoficephenomenainthelast40years,whichprovedtobethe periodofthegreatesttransformationsoftheiceregimeofthelake.Intheyears19612004therewasobservedatendencyofashortening oftheicecoverduration(9.8days/10yrs)andthedurationoficephenomena(9.0days/10yrs).Thepercentageshareoftheicecover durationwithrespecttothetotaldurationoficephenomenadecreasedatarateof6.3%/10ys.Therearealsoobservedshiftsofdates ofthebeginningandendofindividualphasesoftheiceregimemanifestedingradientsofarangeof1.7–4.9days/10years. Keywords: icephenomena,icecover,NorthAtlanticOscillation. Introduction research was performed on the basis of short measurement series covering data sets of periods Theoccurrenceoficephenomena,widespreadin no longer than 25 years. Only the recent years lakesoftemperateclimaticzone,isoneofthebasic broughtstudieswhichuseddataseriesoftheyears characteristicsofthethermalregimeofeachlake. 1961–2000 (Skowron, 1997, 2003; Sziwa, 2002), Moreover,theintensityoftheseprocessesaswell significantlyenrichingtheknowledgeofspatialand astheirdurationaresignificantfactorsshapingthe temporaldeversityoficeregimesofPolishlakes. ecologicalconditionsoflakesandfactorswhichmay Alhough they present a new approach to the havesocioeconomicconsequences. discussedproblem,concentratingbasicallyonthe TheinvestigationoficeregimesofPolishlakes analysis of multiannual changes of chosen wasstartedinthelate1930’sbyMatusewicz(1939) parameters of ice regime of lakes, they still reveal andformanyyearsafterthattheproblemwasonly basic interpretation limitations caused by restricting marginally dealt with. Its main directions were the measurement series to a period of maximum 40 basicallydeterminedby(1)regionaldescriptionof years (Skowron, 2003). Hence, there arises the the diversity of the course of ice phenomena question if the tendencies of changes observed at (Gołek, 1957; Pasławski, 1982; Skowron and present in Polish lakes are a persistent process or Szczepanik, 1988), (2) forecasting parameters of onlyapeculiardisturbanceinthesecularcycle. ice regime (Lityńska, 1966, 1969; Pasławski, Theaimofthisstudyisdeterminedbythesearch 1982)and(3)determinationofthecourseandnatural foranswerstotheabovequestionsandcoversan conditions of ice phenomena in chosen lakes attempt at the reconstruction of multiannual (Chojnowski, 1964; Grześ, 1974). Regional fluctuations of the duration of ice phenomena in 18 DariuszBorowiak,JacekBarańczuk Upper Radunia Lake from the beginning of the the highest situated within the European Lowland 19th century, atthedeterminationofcharacteristic (WieŜyca 329 m above sea level), mean air periodsintheircourseandcomparingtheobtained temperaturesareherelower(about1.5°C)thanin results with the tendencies in the course of this theadjacentareas.UpperRaduniaLakesituatedin phenomenon observed presently in other lakes of a postglacial channel is subject to characteristic PolishLowland. wind operation. Winds occur mainly along the longitudinal axix of the lake from SW or NE differinginspeedandfrequencyofoccurrencein Areaofstudy comparison with winds on the morainic plateau (Okulanis, 1981). Winds blowing along the UpperRaduniaLakeissituatedinthecatchmentof channelaxisconstituteintotalasmuchas73.3% theupperriverRadunia(KashubianLakeland)at of the sum of winds. Also wind speeds show a analtidudeof162mandcoversanareaof387.2ha. similartendency.Thehighestspeedsareofsouthern Its maximum depth is 43 m, with mean depth and southwestern winds: 4.5 and 4.4 m/s reaching15.5mandcapacityis60.2mlnm 3.The respectively(BarańczukandMarchlewicz,2003). morphologyofthelakebottomischaracterisedby Measurements of watertemperatureand ice theoccurrenceofoverdeepenings,bottombarsand phenomenahavebeenperformedintheLimnological islands dividing the basin into three clearly StationoftheUniversityofGdańsk,situatedinthe distinctive parts. The lake is a flowthrough northwesternBorucnoBasin,continuouslysince reservoir with the intensity of horizontal water 1961. exchange of FR = 0.467 and retention time RT = 2.14 (Borowiak, 2000), which indicates a passive hydrologicaltype.Acharacteristicfeatureofthe Results lake is its special structure of water balance in whichthedominatingroleonthesideofincomeis Iceconditionsbetween1961–2004 playedbyundergourndinflow,constitutingover50% ofthetotalsupplyvalue. Theanalysisofarchivaldataresultsintheconclusion KashubianLakelandisacharacteristicregion thatthecourseoficephenomenainUpperRadunia inyoungglaciallandscape.Intermsofclimateit Lakerevealsconsiderableyeartoyearfluctuations belongs to the province of Lakeland manifested in large oscillations in dates of distinguished by a big number of frosty days, on formation,disappearanceanddurationofvarious theaverige130(Kwiecień,1979).Astheregionis formsofice(Fig.1). SecularfluctuationsoficephenomenainUpperRaduniaLake,KashubianLakeland 19

Fig.1.DurationoftheicephenomenaphasesinUpperRaduniaLakeintheyears1961–2004 Thebeginningoficephenomenaonthelake, days,whichis77.2%ofthewholeicephenomena formationofbankice,occursafter1stDecember period. (years1989,1994,1999),butnotlaterthanon15th Infiveyears(1961,1962,1967,1980,1981) February (1975). The mean date of bank ice theicecoverwasnotprecededbytheoccurrence appearance marks the 28th December (SD ofbankiceand/orfloatingice.Bankiceremained /standard deviation/ = 18.1 days). The thelongestin1974,for75days,whiletheshortest stabilisation of ice cover occurs only about two in1982–onlyoneday.Themeandurationofbank weeks later, on 10th January (SD = 20.6 days). icewas11days,whichis13.9%ofthetotalduration Extreme dates of occurrence of continuous ice oficephenomena.FloatingiceonUpperRadunia cover fall on the following extreme dates: 7th Lakeoccurredin42examinedperiods.Itdidnot December(1981)and3rdMarch(1988).However, occurin1961and1990.Thelargestnumberoftimes astableicecovermaynotformatall,asitwas itappearedin1973anditwasrecordedon23days thecasein1990,whilein1989itsoccurrencewas intotal.Themeandurationoffloatingiceonthelake limited to one day only, and ice thickness was was7days–8.9%ofthedurationoficephenomena. then just 1 cm. The disappearance of ice cover Intheperiod1961–2004therewereobserved takes place between 21st January (2002) and 21st tendenciesofchangesofindividualphasesoftheice April(1970),whereasthemeandateoftheendof regimeofthelake.Theyweremainifestedbasically ice cover falls on 18th March (SD = 21.5 days). inchangesinthedurationofparticularphasesas The termination of ice phenomena occurs only wellasashiftofthedatesoftheirbeginningand/or about 27th March (SD = 20.7 days), and threshold ending. The total duration of ice phenomena dates are 8th January (1990) and 30th April gradually decreased by 9.0 days/10 years on the (1969).Hence,theoretically,thelengthoftheperiod average. At the same time, there was observed an withicephenomenamayevenbefourmonths,i.e. earlier ending of ice phenomena with a 151daysintotal.Infact,thelongestdurationofice simultaneousearlierbeginning,by4.9days/10years phenomena,closetothetheoreticalone,wasrecorded and 1.7 days/10 years respectively. Whereas, the in1970anditwas141days,andastableicecover icecoverdurationdecreasedby9.8days/10years remainedfor123days.Theshortestdurationofice on the average. In this case, there was not phenomena was in 1990 (4 days in total), and a observedacleartendencyofashiftofthedayof continuousicecoverdidnotoccuratall.Themean icecoverappearance,thoughthedispersionfrom ice coverduration on Upper Radunia Lake is 61 themeanexpressedasstandarddeviationishere asmuchas20.6days.Ontheotherhand,thereis 20 DariuszBorowiak,JacekBarańczuk clearlymarkedchangeinthedatesoftheicecover ofheatresourcesismorejustifiedasdepthis,insuch disappearance.Theicecoverdisappearsearlier,and acase,aconstantvaluenotinfluencingtheresult theintensityofthesechangesismanifestedinthe of forecasting. Since, as it was shown earlier, ice gradientofabout4.9days/10years.Theresultant phenomena start on Upper Radunia Lake at the of the above observed tendencies is also the beginningofDecember,thesumofmeanmonthly changeinthestabilityofthedurationoficecover. temperatures of water surface of the period The percentage share of the ice cover duration October–Novemberwasadoptedastheindexof withregardtothetotaldurationoficephenomena heatresourcesofthelake. reducedatarateof6.3%/10years. Asaresultofcalculationstherewereobtained Themeanofthemaximumthicknessesofice,for equationsoflinearregressionforsuchelementsof the examined period, was 23.3 cm (SD = 11.9). iceregimeasthedurationoficephenomena,duration Conditionstheclosesttothemeanoccurredinthe oficecoverandmeanicethickness(Tab.1).Asthe years1967,1993,1997.Thehighestvalueoftheice degree of explanation of the variability of the ice coverwasobservedin1970whenitwasasmuch regime elements on the basis of equations with as 50 cm and occured for subsequent 44 days oneandtwopredictorsissimilar,andinthecaseof between25thFebruaryand9thApril.Extremewas determining the ice thickness, the degree of alsotheyear1990(lackoficecover)and1989with explaining of the dependent variable of one icecoverofonly1cmofthickness.Themeanice parameterequationisover50%,andduetoalack thicknessinthewholeperiodwas16.5cm(SD= of available observations of water temperatures, 8.6). The most typical, in this respect, were the which is a decisive reason, oneparameter years1966,1967,1977,1994and2003. equationswereusedtoreconstructiceconditions inthelakeintheyears1825–2004.Thisapproach isthemoreacceptableasthelevelofinformation Fluctuationsofelementsoftheiceregimeof lossinthecaseofreconstructingthedurationofice UpperRaduniaLakeintheyears1825–2004 phenomena (the biggest differences between both equations) is very small, and calculated values Intheresearchconcerningtheicingoflakes,air similar (Fig. 2). Finally, there was reconstructed temperatureandlakedepthwereindicatedasthe thechronologicalcourseoficephenomenafrom mainfactorsshapingtheconditionsoftheicecover 1825.Itindicatesahighinterannualdiversityof formation(Lityńska,1966,Pasławski,1982).Asair icephenomenainthediscussed180yearsperiod. temperatureisshapedbythedirectionsandintensity However,smoothingofthediagramandremoving of advection of air masses, it may be indirectly fluctuations of high frequency (short duration) determined by gradients of the pressure field or allows for distinguishing several distinctive indices describing them (KoŜuchowski and periods of the development of ice phenomena śmudzka2002). One of such indicators, strongly characterised by different values of ice regime conditioning the circulation of air masses over parameters, as well as different tendencies of northernEuropeandalsosignificantlyinfluencing occurringchanges(Fig.3). airtemperaturedistributionsin,especially in winter, is the index of winter North Atlantic Table1.Equationsdescribingselectedparametersoftheice Oscillation (NAO) (Hurrell, 1995; Marsz, 1999; regimeofUpperRaduniaLake(DIP–durationofice phenomena,DIC–durationoficecover,MIT–mean Marsz and Styszyńska, 2001). In this paper there ice thickness) as a function of the index of winter was used winter NAO index (mean value for NorthAtlanticOscillation(NAO)andthesumof period December – March) determined by Jones meanmonthlytemperaturesofwatersurfaceinOcto (Jones et al .,1997;CRC2004). Inturn,the mean berandNovember(WST).Calculationsfortheperiod depth of the lake indirectly characterises the 1961–2003 capacityofheatstorageintheperiodprecedingice EQUATION R2 phenomena and may be substituted for by a DIP=23,074*NAO–4,957*WST+174,300 0.667 parameter describing the quantity of the DIC=22,431*NAO–2,022*WST+105,140 0.576 MIT=5,258*NAO–0,453*WST+26,549 0.498 accumulated heat resources in the surface water DIP=21,796*NAO+87,553 0.613 layer. In the case of examining the course of ice DIC=21,904*NAO+69,854 0.567 phenomenainonlyonereservoir,adoptingtheindex MIT=5,282*NAO+19,196 0.527 SecularfluctuationsoficephenomenainUpperRaduniaLake,KashubianLakeland 21

Fig.2.Observed(thicksolidline)andcalculatedvaluesofthedurationoficephenomena.SolidlinedenotesNAOmodelandbroken lineNAO–WST(watersurfacetemperature)model

Fig.3.IcephenomenainUpperRaduniaLakeintheyears18252004.Movingaverage(7yearsstep)ismarkedinthickline Discussion beginning and end. There is also a noticeable decreaseintheshareofthedurationoficecover TheresearchonicingoflakesofPolishLowland with respect to the length of the whole period of performedbySkowron(2003)unambiguouslyprove ice phenomena. Against this background, Upper theoccurrenceintheyears1961–2000ofatendency Radunia Lake does not differ from other Polish ofshorteningofthedurationoficecoverandice lakes.Thecharacterofthechangesobservedinthis phenomena.Therearealsomarkedclearchanges lakeisveryclosetothetendenciesoccurringinthe of dates of the beginning and end of individual majority of lakes of Pomerania and Masurian phasesofthedevelopmentoficephenomena.These Lakelands(Tab.2).Onlywithreferencetothedate changes are manifested through ever earlier of the beginning of the ice cover formation there 22 DariuszBorowiak,JacekBarańczuk were not observed significant changes. In the śnińskieDuŜe).Theshorteningoftheicecover quoted research also this parameter revealed a durationisalsoreflectedinagradualreduction of highly ambiguous character of changes. There can the ice cover thickness. The maximum ice belistednumerouslakesonwhichtheformation thicknessisdecreasingatarateof3.9cm/10years of ice cover was delayed (Lakes Osiek, Lubie, andfallwithintherangeobtainedinearlierresearch Charzykowskie,Bukowo)orrelativelystable(Lake (Skowron,2003). Table2.Meanvaluesofthedurationofparticularphasesoficephenomena(indays)andmaximumicethickness(incm)inUpper RaduniaLakeintheyears1961–2004dividedintodecades YEARS BANKICE ICECOVER FLOATINGICE ICEPHENOMENA ICETHICKNESS 1961–1970 6.9 84.4 8.7 100.0 32.8 1971–1980 13.9 61.4 6.9 82.2 24.0 1981–1990 9.6 55.8 5.6 71.0 19.9 1991–2000 11.1 43.9 9.7 64.7 19.1 2000–2004 12.0 51.8 4.3 67.8 17.3 However,thesetendencieshadavariedcourse Ascanbeseenfromthepresentedsurvey,the duringtheyears1825–2004(Fig.3).Tilltheendof evaluationofthetendenciesofchangesofelements the19thcentury,icephenomenashowedarelative oftheiceregimeofPolishlakesbasedondataof stability,andtheirfluctiationswereshortterm.In theyears1961–2000concernsaspecialperiodintheir thatperiodtheelementsoftheregimewereclose secular fluctuations. It covers the period of the tothemeanvalues( DIP =77 days, DIC=60 days, strongestintensificationoficephenomena(1960– 1970),followedaphaseofwinterwarmingstrongly MIT =17 cm)andtheywereDIP=78days(SD limitingtheirdevelopmentintheyears1981–2000. =23.5),DIC=60days(SD=23.6)andMIT= Thus,thisisaperiodoftheoccurrenceofthelargest 17.2cm(SD=5.6)respectively.Atthebeginning gradients of individual elements of the regime. ofthe20thcentury(1919–1930)thereoccursaclear Hence,thisperiodcanbetreatedasanidealone warmingcausedbytheintensificationofwestern fortheidentificationandevaluationoftheinfluence circulation,whichisreflected,amongothers,in of factors disturbing the development of ice shorteningofthedurationoficephenomenabyabout phenomenaonPolishlakes.Thereversalorstrong 10days(DIP=68±20.6days,DIC=50±20.7). reductionoftheobservedtrendsmustbearesult Simultaneously, the mean ice thickness is of various local conditions overlapping with the decreasing; its value in that period fell to 14.7 generalpictureofchangescausedbyextraregional ±4.9cm.Afterthatperiod,upto1980thereoccurs influences. an intensification of ice phenomena. Their duration increases to 86 ±23.8 days. The ice cover, in comparisonwiththevalueobtainedforthewhole References analysedperiod,occuredontheaverage7days longer(DIC=67±21.3).Themeanthicknessof Borowiak D., 2000, ReŜimy wodne i funkcje hydrolo icecoverincreasesto18.8±5.0cm.Theheightof gicznejeziorNiŜuPolskiego.Bad.Limnol.2,Wyd. the intensification of ice phenomena fell on the KLUG,Gdańsk. years 1961–1970 when the length of individual Barańczuk J., Marchlewicz R., 2003, Diversity of phasesoftheiceregimeincreasedbyover20%in development of ice phenomena in chosen lakes of comparisonwiththemeanvalues(DIP=99±31.0 KaszubskieLakelandinwinter2003,Limnological days, DIC=74 ±25.7 days), and the ice cover Review,vol.3,Kielce,3–8. thicknessincreasedby20.5%(MIT=20.5±6.0cm). ChojnowskiS.,(edit.),1964,Zbadańnadpokrywąlodo wąJezioraMikołajskiego,Biul.PIHM8. The end of the 20th century is distinguished by ClimaticResearchCenter,2004,NorthAtlanticOscillation, another increase in the intensity of western Retrieved June 24, 2004, from University of East circulation, and ice regime parameters reveal Anglia, Climatic Research Center Unit. Available similaritytotheonesrecordedatthebeginningof athttp://www.cru.uea.ac.uk/cru/data/nao.htm thecentury(DIP=68±31.8days,53±25.4days, GołekJ.,1957,Zjawiskalodowenarzekachpolskich, 15.5±6.0cm). PracePIHM48. SecularfluctuationsoficephenomenainUpperRaduniaLake,KashubianLakeland 23 GrześM.,1974,Badanianadtermikąizlodzeniemjeziora dejściedoomawianegozagadnienia,koncentrującsięna Gopło,Dok.Geogr.3,IGPAN. analiziewieloletnichzmianwybranychparametrówreŜi Hurrell,J.W.,1995,DecadaltrendsintheNorthAtlantic mulodowegojezior,wciąŜujawniajązasadniczeograni Oscillationandrelationshipstoregionaltemperature czenia interpretacyjne spowodowane redukcją serii po andprecipitation,Science269,676–679. miarowych do okresu maksimum 40letniego (Skow JonesP.D.,JonssonT.,WheelerD.,1997,Extensionto ron,2003).Pojawiasięzatempytanieczyobserwowa theNorthAtlanticOscillationusingearlyinstrumental newspółcześniewjeziorachpolskichtendencjezmian pressureobservationsfromGibraltarandSouthWest są procesem trwałym, czy stanowią jedynie swego ro Iceland,Int.J.Climatol.17,1433–1450. dzajuzaburzeniewcyklusekularnym?Celopracowa KoŜuchowskiK.,śmudzkaE.,2002,Cyrkulacjaatmosfe nia wyznacza zatem poszukiwanie odpowiedzi na po rycznaijejwpływnazmiennośćtemperaturypowie stawionepowyŜejpytanieiobejmujepróbęodtworze trzawPolsce,Prz.Geogr.74(4),591–604. nia wieloletnich,poczynającodpoczątkuXIX wieku, KwiecieńK.,1979,Warunkiklimatyczne,In:Augustowski wahań czasu trwania zjawisk lodowych na Jeziorze B.,(Edit),PojezierzeKaszubskie,GTN,Gdańsk. Raduńskim Górnym, wyznaczenie charakterystycznych LityńskaZ.,1966,ForecastofMazurianLakeDistrict okresówwichprzebieguorazodniesienieuzyskanych freeze,Hydrologyoflakesandreservoirs,IAHSPub. wynikówdoobserwowanychwspółcześnietendencji. No70(1),392–396. Zdokonanejanalizydanych archiwalnych wyni MarszA.,1999,OscylacjaPółnocnoatlantyckaareŜim ka,ŜeprzebiegzjawisklodowychwJeziorzeRaduńskim termicznyzimnaobszarzepółnocnozachodniejPolski Górnym wykazujeznacznefluktuacjemiędzyroczne, inapolskimwybrzeŜuBałtyku,PraceGeograficzne,71 przejawiające się duŜymi wahaniami w terminach po (3),225–245. wstania,zanikuiczasutrwaniaróŜnychformzlodzenia MarszA.,StyszyńskaA.,2001,OscylacjaPółnocnego (rys. 1). Początek zjawisk lodowych na jeziorze, for AtlantykuatemperaturapowietrzanadPolską,Wyd. mowaniesięlodubrzegowego,następujeniewcześniej WSM,Gdynia. aniŜeli1grudnia(lata1989,1994,1999)iniepóźniej MatusewiczJ.,1939,Badanianadzlodzeniemjezior niŜ15lutego(1975). Zakończeniezjawisk lodowych na wPolsce,Wiad.Sł.Hydrogr.2(1),49–59. stępujewJeziorzeRaduńskimdopierookoło27marca Okulanis E., 1981, Studium limnologiczne Jezior Raduń (SD=20,7dnia),adatyprogowewyznaczają:8stycz skoOstrzyckich,OssolineumGTN,WrocławGdańsk. nia(1990)oraz30kwietnia(1969).Wlatach1961–2004 PasławskiZ.,1982,ZlodzeniejeziorwPolsce,Prz. zaobserwowano tendencje zmian poszczególnych faz Geofiz.27(1–2):79–92. reŜimulodowego.Przejawiałysięonezasadniczozmia SkowronR.,1997,Tendencjezmiantemperaturywody namiwczasietrwaniaposzczególnychfazjakteŜprze powierzchniowejizjawisklodowychwjeziorachna sunięciemdatichpoczątkui/lubzaniku.Całkowityczas obszarachpojeziernychwPolsce,In:ChoińskiA., trwaniazjawisklodowychulegałstopniowemuskróce (edit.),Wpływantropopresjinajeziora,Wyd.Homini, niuśrednioo9,0dni/10lat.Jednocześniezaobserwo Poznań–Bydgoszcz,143–151. wano wcześniejsze zanikanie zjawisk lodowych przy Skowron R., 2003, Ice sheet in the lakes of the Polish wcześniejszymichpojawianiusię,odpowiednioo4,9 lowland, Distribution, differences and trends, dni/10latoraz1,7dni/10lat.Zkoleiczaszaleganiapo Limnological.Review,vol.3,Kielce,205–212. krywylodowejskróciłsięprzeciętnieo9,8dnia/10lat. SkowronR.,SzczepanikW.,1988,ZróŜnicowanieprze Wykorzystującdanezlodzeniajeziorazlat1961– strzenneitendencjezmianprzebieguzjawisklodo 2003orazzimowyindeksNAO(średniawartośćzaokres wychnajeziorachPolskipólnocnej,In:ChurskiZ., grudzień–marzec)opracowanomodelepozwalającena (edit.),Naturalneiantropogeniczneprzemianyjezior predykcjęwybranychcharakterystykreŜimulodowego imokradełwPolsce,Wyd.UMK,Toruń,125–136. jeziora(tab.1).PoniewaŜstopieńwyjaśnianiazmienności SkowronR.,SzczepanikW.,1993,ZróŜnicowanieprze elementówreŜimulodowegonapodstawierównańjedno strzenneitendencjezmiantemperaturywodyizlo idwuparametrowychjestzbliŜony(rys.2)dorekonstruk dzeniajezior,In:DynowskaI.,(edit.),Przemianysto cjiwarunkówlodowychpanującychwjeziorzewlatach sunków wodnych w Polsce w wyniku procesów 1825–2004 wykorzystano równania jednoparametro naturalnychiantropogenicznych,Wyd.UJ,Kraków, we. Analiza sekularnych wahań parametrów reŜimu lo 121–136. dowego pozwala na wydzielenie kilku odrębnych okre SziwaR.,2002,Maximumicecoverthicknessonlakes sów rozwoju zjawisk lodowych cechujących się od of the basin, Przymorze and the lower miennymiwartościamiśrednimiparametrów,atakŜe basin,LimnologicalReview,vol.2,Lublin,391–397. róŜnymitendencjamizachodzącychzmian(rys.3). Streszczenie UprzedniebadaniaSkowrona(2003)dowodząwy stępowaniawostatnim40leciutendencjiskracaniaczasu Ostatniepracedotyczącerozpoznaniaprzebieguzjawisk zalegania pokrywy lodowej oraz trwania zjawisk lodo lodowych najeziorachPolski choćprezentująnowepo wych.WskazująteŜnaczytelneprzesunięciadatpocząt 24 DariuszBorowiak,JacekBarańczuk kuikońcaposzczególnych fazrozwojuzjawisklodo geum przypada na lata 19611970 kiedy to długość wych.NatymtlejezioroRaduńskieGórnenieodbiega poszczególnych faz reŜimu lodowego ulega wydłuŜe odpozostałychjeziorPolski,acharakterzmianobser niuoponad20%wporównaniuzwartościamiśredni wowanych na tym jeziorze jest zbliŜony do tendencji mi.KoniecwiekuXXzaznaczasięponownymwzrostem obserwowanychnawiększościjeziorPojezierzaPomor natęŜeniacyrkulacjizachodniej,aparametryreŜimulodo skiegojakiMazurskiego(tab.2). wegoujawniająpodobieństwodonotowanychnapoczątku JednakŜe tendencje te róŜnie kształtowały się na stulecia(DIP=68dni,DIC=53dni,MIT=15,5cm). przestrzenilat1825–2004(rys.3).DokońcaXIXw. Ocena tendencji zmian elementów reŜimu lodo zjawiska lodowe wykazywały względną stałość, a ich wegojeziorpolskichopartanadanychzlat19612000 wahaniabyłykrótkookresowe.NapoczątkuwiekuXX odnosisięzatemdospecyficznegookresuwichseku (1991–1930)zaznaczasięociepleniepowodowaneinten larnychfluktuacjach.Obejmujebowiemokresonajsil syfikacjącyrkulacjizachodniej,comanifestujesięskra niejszejintensyfikacjizjawisklodowych(19601970), ceniemczasutrwaniazjawisklodowychookoło10dni poktórymzkoleinastępujefazazimowegoocieplenia (DIP=68dni,DIC=50dni).Jednocześnieredukcji,do silnieograniczającaichrozwójwlatach19812000.Tym 14,7cm,ulegaśredniagrubośćpokrywylodowej.Pofazie samymjesttookresnajwiększychgradientówposzcze tej,aŜdoroku1980,mamiejsceintensyfikacjazjawisk gólnych elementów reŜimu, moŜe więc być potrakto lodowych. Czas ichtrwania wydłuŜa siędo86dni.Po wanyjakoidealnydoidentyfikacjiiocenywpływuczyn krywalodowazalegaprzeciętnieo7dnidłuŜejaśred nikówzakłócającychrozwójzjawisklodowychnajezio niagrubośćpokrywylodowejwzrastado18,8cm.Apo rachpolskich.