Optical Studies of the Antarctic Glacio-Oceanic System

The Antarctic system comprises of the continent itself, Antarctica, and the ocean surrounding it, the Southern Ocean. The system has an important part in the global climate due to its size, its high latitude location and the negative radiation balance of its large ice sheets. Antarctica has also been in focus for several decades due to increased ultraviolet (UV) levels caused by stratospheric ozone depletion, and the disintegration of its ice shelves. In this study, measurements were made during three Austral summers to study the optical properties of the Antarctic system and to produce radiation information for additional modeling studies. These are related to specific phenomena found in the system. During the summer of 1997-1998, measurements of beam absorption and beam attenuation coefficients, and downwelling and upwelling irradiance were made in the Southern Ocean along a S-N transect at 6°E. The attenuation of photosynthetically active radiation (PAR) was calculated and used together with hydrographic measurements to judge whether the phytoplankton in the investigated areas of the Southern Ocean are light limited. By using the Kirk formula the diffuse attenuation coefficient was linked to the absorption and scattering coefficients. The diffuse attenuation coefficients (Kpar) for PAR were found to vary between 0.03 and 0.09 1/m. Using the values for KPAR and the definition of the Sverdrup critical depth, the studied Southern Ocean plankton systems were found not to be light limited. Variabilities in the spectral and total albedo of snow were studied in the Queen Maud Land region of Antarctica during the summers of 1999-2000 and 2000-2001. The measurement areas were the vicinity of the South African Antarctic research station SANAE 4, and a traverse near the Finnish Antarctic research station Aboa. The midday mean total albedos for snow were between 0.83, for clear skies, and 0.86, for overcast skies, at Aboa and between 0.81 and 0.83 for SANAE 4. The mean spectral albedo levels at Aboa and SANAE 4 were very close to each other. The variations in the spectral albedos were due more to differences in ambient conditions than variations in snow properties. A Monte-Carlo model was developed to study the spectral albedo and to develop a novel nondestructive method to measure the diffuse attenuation coefficient of snow. The method was based on the decay of upwelling radiation moving horizontally away from a source of downwelling light. This was assumed to have a relation to the diffuse attenuation coefficient. In the model, the attenuation coefficient obtained from the upwelling irradiance was higher than that obtained using vertical profiles of downwelling irradiance. The model results were compared to field measurements made on dry snow in Finnish Lapland and they correlated reasonably well. Low-elevation (below 1000 m) blue-ice areas may experience substantial melt-freeze cycles due to absorbed solar radiation and the small heat conductivity in the ice. A two-dimensional (x-z) model has been developed to simulate the formation and water circulation in the subsurface ponds. The model results show that for a physically reasonable parameter set the formation of liquid water within the ice can be reproduced. The results however are sensitive to the chosen parameter values, and their exact values are not well known. Vertical convection and a weak overturning circulation is generated stratifying the fluid and transporting warmer water downward, thereby causing additional melting at the base of the pond. In a 50-year integration, a global warming scenario mimicked by a decadal scale increase of 3 degrees per 100 years in air temperature, leads to a general increase in subsurface water volume. The ice did not disintegrate due to the air temperature increase after the 50 year integration.

Grönlannin mannerjäätikön tulevaisuus: merenpinnan nousu maailmanlaajuisesti ja Itämerellä

Sea level rise is among the most worrying consequences of climate change, and the biggest uncertainty of sea level predictions lies in the future behaviour of the ice sheets of Greenland and Antarctica. In this work, a literature review is made concerning the future of the Greenland ice sheet and the effect of its melting on Baltic Sea level. The relation between sea level and ice sheets is also considered more generally from a theoretical and historical point of view. Lately, surprisingly rapid changes in the amount of ice discharging into the sea have been observed along the coastal areas of the ice sheets, and the mass deficit of Greenland and West Antarctic ice sheets which are considered vulnerable to warming has been increasing from the 1990s. The changes are probably related to atmospheric or oceanic temperature variations which affect the flow speed of ice either via meltwater penetrating to the bottom of the ice sheet or via changes in the flow resistance generated by the floating parts of an ice stream. These phenomena are assumed to increase the mass deficit of the ice sheets in the warming climate; however, there is no comprehensive theory to explain and model them. Thus, it is not yet possible to make reliable predictions of the ice sheet contribution to sea level rise. On the grounds of the historical evidence it appears that sea level can rise rather rapidly, 1 2 metres per century, even during warm climate periods. Sea level rise projections of similar magnitude have been made with so-called semiempirical methods that are based on modelling the link between sea level and global mean temperature. Such a rapid rise would require considerable acceleration of the ice sheet flow. Stronger rise appears rather unlikely, among other things because the mountainous coastline restricts ice discharge from Greenland. The upper limit of sea level rise from Greenland alone has been estimated at half a metre by the end of this century. Due to changes in the Earth s gravity field, the sea level rise caused by melting ice is not spatially uniform. Near the melting ice sheet the sea level rise is considerably smaller than the global average, whereas farther away it is slightly greater than the average. Because of this phenomenon, the effect of the Greenland ice sheet on Baltic Sea level will probably be rather small during this century, 15 cm at most. Melting of the Antarctic ice sheet is clearly more dangerous for the Baltic Sea, but also very uncertain. It is likely that the sea level predictions will become more accurate in the near future as the ice sheet models develop.

Snow cover models in operational watershed forecasting

Yhteenveto: Lumimallit vesistöjen ennustemalleissa

Lumen sulamisesta ja sulamiskauden vesitaseesta Lammin Pääjärvellä 1970-72

English summary: On snowmelt and water balance during snowmelt period in Pääjärvi representative basin in 1970-72

Metsämaan vesitaseen muutokset Suomessa ilmaston muuttuessa

Muuttuva ilmasto, erityisesti kohoava lämpötila ja hiilidioksidipitoisuus, ei voi olla vaikuttamatta metsien kasvuun Suomessa. Metsätalouden sopeutumistoimilla voidaan pyrkiä lisääntyvän tuotospotentiaalin hyödyntämiseen ja metsätuhoriskien pienentämiseen. Sopeutumistoimien suunnittelemiseksi tarvitaan kuitenkin ensin tietoa siitä, miten ilmasto muuttuu ja mitä vaikutuksia ilmastonmuutoksella on. Työni oli osa vuonna 2005 käynnistettyä Ilmastonmuutoksen sopeutumistutkimusohjelmaa (ISTO). ISTO, ja näin myös oma tutkimukseni, tuottaa siis tietoa, jonka avulla voidaan mukauttaa metsänhoitoa vastaamaan muuttuvia ilmastooloja. Työssä tarkasteltiin maan vesipitoisuutta menneessä, nykyisessä ja muuttuvassa ilmastossa. Tavoitteena oli ennustaa, miten kuivuuden esiintyminen muuttuu Suomessa. Lisäksi menneiden kuivien vuosien kasvuvaikutuksia tutkittiin lustomittausten avulla. Vesitasemuutoksia lähdettiin tutkimaan mallilaskelmilla. Laskelmat perustuivat yksinkertaiseen open bucket – tyyppiseen vesimalliin. Mallissa maaperä ajatellaan tilana, johon sadanta ja lumen sulaminen tuovat vettä. Haihdunta ja valunta päinvastaisesti vähentävät sitä. Mallia sovellettiin Metlan yhdeksällä provenienssikoealueella, joihin lukeutui sekä kuusikoita että männiköitä eri kasvupaikoilta. Mallin yksinkertaisuuden takia koealueilta tarvittiin ainaostaan päivittäiset tiedot lämpötilasta, sadannasta, säteilystä ja VPD:stä. Lisäksi tuli selvittää kasvupaikkojen savi- ja hiekkapitoisuudet maaperäparametrien laskemiseksi. Vesimallilla tehtiin kahdet laskelmat. Ensimmäisissä laskettiin päivittäistä vesitasetta vuosina 1961-2008 mitattua sääaineistoa käyttäen. Toisissa laskelmissa käytettiin mitatusta säädatasta modifioitua sääennustetta. Lämpötilaa ja sadantaa oli kasvatettu vastaamaan vuosisadan lopulle tehtyä ennustetta. Näiden kahden laskelman tuloksia vertailtiin keskenään kuivuuspäivien lukumäärän osalta. Malliin oli siis ohjelmoitu kuivuuspäiväindeksi. Indeksin perusteella kuivuutta katsottiin olleen päivänä, jolloin kasveille käyttökelpoisen veden määrä laski alle 4 prosentin ja haihdunta leikkaantui 90 prosentilla normaalitasoon nähden. Kuivuuspäivät lisääntyivät kaikilla provenienssikoealueilla. Suurin muutos olisi tulosten perusteella odotettavissa etelärannikolla ja männyn kasvupaikoilla, jotka jo entuudestaan olivat alueista kuivimpia. Vähäisin muutos olisi taas odotettavissa kuusen kasvupaikoilla, etenkin jos sijainti ei ole aivan eteläisimmässä Suomessa. Kasvuvaikutuksia koskevien tulosten perusteella maan vesipitoisuus ei ole ollut merkittävä lustojen kasvua rajoittava tekijä provenienssikoealueilla menneinä vuosikymmeninä. Kuivuuspäivien määriin ei tule suhtautua absoluuttisina totuuksina, sillä ne ovat määritelmänsä mukaisesti vain suuntaa antavia. Niiden avulla voidaan tehdä ainoastaan päätelmiä metsiköiden ja vuosien välisistä suhteellisista eroista ja kehityssuunnista. Suoria johtopäätöksiä puuston kasvusta tai elinvoimaisuudesta ei pystytä tekemään. Eikä myöskään pystytä arvioimaan, reagoivatko mänty ja kuusi ennustettuihin muutoksiin eri tavoin. Voidaan kuitenkin päätellä, että keskimääräisten vuotuisten kuivuuspäivien lukumäärän kasvu johtaa brutto- ja nettoprimäärituotoksen vähenemiseen. Tätä kautta kuivuuden lisääntyminen alentaa myös kasvua verrattuna tilanteeseen, jossa kuivuuspäivien määrä pysyisi muuttumattomana.

Antarctic ice shelf melting and its impact on the global sea ice-ocean system

Earth s ice shelves are mainly located in Antarctica. They cover about 44% of the Antarctic coastline and are a salient feature of the continent. Antarctic ice shelf melting (AISM) removes heat from and inputs freshwater into the adjacent Southern Ocean. Although playing an important role in the global climate, AISM is one of the most important components currently absent in the IPCC climate model. In this study, AISM is introduced into a global sea ice-ocean climate model ORCA2-LIM, following the approach of Beckmann and Goosse (2003; BG03) for the thermodynamic interaction between the ice shelf and ocean. This forms the model ORCA2-LIM-ISP (ISP: ice shelf parameterization), in which not only all the major Antarctic ice shelves but also a number of minor ice shelves are included. Using these two models, ORCA2-LIM and ORCA2-LIM-ISP, the impact of addition of AISM and increasing AISM have been investigated. Using the ORCA2-LIM model, numerical experiments are performed to investigate the sensitivity of the polar sea ice cover and the Antarctic Circumpolar Current (ACC) transport through Drake Passage (DP) to the variations of three sea ice parameters, namely the thickness of newly formed ice in leads (h0), the compressive strength of ice (P*), and the turning angle in the oceanic boundary layer beneath sea ice (θ). It is found that the magnitudes of h0 and P* have little impact on the seasonal sea ice extent, but lead to large changes in the seasonal sea ice volume. The variation in turning angle has little impact on the sea ice extent and volume in the Arctic but tends to reduce them in the Antarctica when ignored. The magnitude of P* has the least impact on the DP transport, while the other two parameters have much larger influences. Numerical results from ORCA2-LIM and ORCA2-LIM-ISP are analyzed to investigate how the inclusion of AISM affects the representation of the Southern Ocean hydrography. Comparisons with data from the World Ocean Circulation Experiment (WOCE) show that the addition of AISM significantly improves the simulated hydrography. It not only warms and freshens the originally too cold and too saline bottom water (AABW), but also warms and enriches the salinity of the originally too cold and too fresh warm deep water (WDW). Addition of AISM also improves the simulated stratification. The close agreement between the simulation with AISM and the observations suggests that the applied parameterization is an adequate way to include the effect of AISM in a global sea ice-ocean climate model. We also investigate the models capability to represent the sea ice-ocean system in the North Atlantic Ocean and the Arctic regions. Our study shows both models (with and without AISM) can successfully reproduce the main features of the sea ice-ocean system. However, both tend to overestimate the ice flux through the Nares Strait, produce a lower temperature and salinity in the Hudson Bay, Baffin Bay and Davis Strait, and miss the deep convection in the Labrador Sea. These deficiencies are mainly attributed to the artificial enlargement of the Nares Strait in the model. In this study, the impact of increasing AISM on the global sea ice-ocean system is thoroughly investigated. This provides a first idea regarding changes induced by increasing AISM. It is shown that the impact of increasing AISM is global and most significant in the Southern Ocean. There, increasing AISM tends to freshen the surface water, to warm the intermediate and deep waters, and to freshen and warm the bottom water. In addition, increasing AISM also leads to changes in the mixed layer depths (MLD) in the deep convection sites in the Southern Ocean, deepening in the Antarctic continental shelf while shoaling in the ACC region. Furthermore, increasing AISM influences the current system in the Southern Ocean. It tends to weaken the ACC, and strengthen the Antarctic coastal current (ACoC) as well as the Weddell Gyre and the Ross Gyre. In addition to the ocean system, increasing AISM also has a notable impact on the Antarctic sea ice cover. Due to the cooling of seawater, sea ice concentration and thickness generally become higher. In austral winter, noticeable increases in sea ice concentration mainly take place near the ice edge. In regards with sea ice thickness, large increases are mainly found along the coast of the Weddell Sea, the Bellingshausen and Amundsen Seas, and the Ross Sea. The overall thickening of sea ice leads to a larger volume of sea ice in Antarctica. In the North Atlantic, increasing AISM leads to remarkable changes in temperature, salinity and density. The water generally becomes warmer, more saline and denser. The most significant warming occurs in the subsurface layer. In contrast, the maximum salinity increase is found at the surface. In addition, the MLD becomes larger along the Greenland-Scotland-Iceland ridge. Global teleconnections due to AISM are studied. The AISM signal is transported with the surface current: the additional freshwater from AISM tends to enhance the northward spreading of the surface water. As a result, more warm and saline water is transported from the tropical region to the North Atlantic Ocean, resulting in warming and salt enrichment there. It would take about 30 40 years to establish a systematic noticeable change in temperature, salinity and MLD in the North Atlantic Ocean according to this study. The changes in hydrography due to increasing AISM are compared with observations. Consistency suggests that increasing AISM is highly likely a major contributor to the recent observed changes in the Southern Ocean. In addition, the AISM might contribute to the salinity contrast between the North Atlantic and North Pacific, which is important for the global thermohaline circulation.

Physical properties of a boreal clay soil under differently managed perennial vegetation

The physical properties of surface soil horizons, essentially pore size, shape, continuity and affinity for water, regulate water entry into the soil. These properties are prone to changes caused by natural forces and human activity. The hydraulic properties of the surface soil greatly impact the generation of surface runoff and accompanied erosion, the major concern of agricultural water protection. The general target of this thesis was to improve our understanding of the structural and hydraulic properties of boreal clay soils. Physical properties of a clayey surface soil (0 - 10 cm, clay content 51%), with a micaceous/illitic mineralogy subjected to three different management practices of perennial vegetation, were studied. The study sites were vegetated buffer zones located side by side in SW Finland: 1) natural vegetation with no management, 2) harvested once a year, and 3) grazed by cattle. The soil structure, hydraulic properties, shrinkage properties and soil water repellency were determined at all sites. Two distinct flow domains were evident. The surface soil was characterized by subangular blocky, angular blocky and platy aggregates. Hence, large, partially accommodated, irregular elongated pores dominated the macropore domain at all sites. The intra-aggregate pore system was mostly comprised of pores smaller than 30 μm, which are responsible for water storage. Macropores at the grazed site, compacted by hoof pressure, were horizontally oriented and pore connectivity was poorest, which decreased water and air flux compared with other sites. Drying of the soil greatly altered its structure. The decrease in soil volume between wet and dry soil was 7 - 10%, most of which occurred in the moisture range of field conditions. Structural changes, including irreversible collapse of interaggregate pores, began at matric potentials around -6 kPa indicating, instability of soil structure against increasing hydraulic stress. Water saturation and several freezethaw cycles between autumn and spring likely weakened the soil structure. Soil water repellency was observed at all sites at the time of sampling and when soil was dryer than about 40 vol.%. (matric potential < -6 kPa). Therefore, water repellency contributes to water flow over a wide moisture range. Water repellency was also observed in soils with low organic carbon content (< 2%), which suggests that this phenomenon is common in agricultural soils of Finland due to their relatively high organic carbon content. Aggregate-related pedofeatures of dense infillings described as clay intrusions were found at all sites. The formation of these intrusions was attributed to clay dispersion and/or translocation during spring thaw and drying of the suspension in situ. These processes generate very new aggregates whose physical properties are most probably different from those of the bulk soil aggregates. Formation of the clay infillings suggested that prolonged wetness in autumn and spring impairs soil structure due to clay dispersion, while on the other hand it contributes to the pedogenesis of the soil. The results emphasize the dynamic nature of the physical properties of clay soils, essentially driven by their moisture state. In a dry soil, fast preferential flow is favoured by abundant macropores including shrinkage cracks and is further enhanced by water repellency. Increase in soil moisture reduces water repellency, and swelling of accommodated pores lowers the saturated hydraulic conductivity. Moisture- and temperature-related processes significantly alter soil structure over a time span of 1 yr. Thus, the pore characteristics as well as the hydraulic properties of soil are time-dependent.