Environmental effects of thermal and radioactive discharges from nuclear power plants in the boreal brackish-water conditions of the northern Baltic Sea

During recent decades, thermal and radioactive discharges from nuclear power plants into the aquatic environment have become the subject of lively debate as an ecological concern. The target of this thesis was to summarize the large quantity of results obtained in extensive monitoring programmes and studies carried out in recipient sea areas off the Finnish nuclear power plants at Loviisa and Olkiluoto during more than four decades. The Loviisa NPP is located on the coast of the Gulf of Finland and Olkiluoto NPP on that of the Bothnian Sea. The state of the Gulf of Finland is clearly more eutrophic; the nutrient concentrations in the surface water are about 1½ 2 times higher at Loviisa than at Olkiluoto, and the total phosphorus concentrations still increased in both areas (even doubled at Loviisa) between the early 1970s and 2000. Thus, it is a challenge to distinguish the local effects of thermal discharges from the general eutrophication process of the Gulf of Finland. The salinity is generally low in the brackish-water conditions of the northern Baltic Sea, being however about 1 higher at Olkiluoto than at Loviisa (the salinity of surface water varying at the latter from near to 0 in early spring to 4 6 in late autumn). Thus, many marine and fresh-water organisms live in the Loviisa area close to their limit of existence, which makes the biota sensitive to any additional stress. The characteristics of the discharge areas of the two sites differ from each other in many respects: the discharge area at Loviisa is a semi-enclosed bay in the inner archipelago, where the exchange of water is limited, while the discharge area at Olkiluoto is more open, and the exchange of water with the open Bothnian Sea is more effective. The effects of the cooling water discharged from the power plants on the temperatures in the sea were most obvious in winter. The formation of a permanent ice cover in the discharge areas has been delayed in early winter, and the break-up of the ice occurs earlier in spring. The prolonging of the growing season and the disturbance of the overwintering time, in conditions where the biota has adjusted to a distinct rest period in winter, have been the most significant biological effects of the thermal pollution. The soft-bottom macrofauna at Loviisa has deteriorated to the point of almost total extinction at many sampling stations during the past 40 years. A similar decline has been reported for the whole eastern Gulf of Finland. However, the local eutrophication process seems to have contributed into the decline of the zoobenthos in the discharge area at Loviisa. Thermal discharges have increased the production of organic matter, which again has led to more organic bottom deposits. These have in turn increased the tendency of the isolated deeps to a depletion of oxygen, and this has further caused strong remobilization of phosphorus from the bottom sediments. Phytoplankton primary production and primary production capacity doubled in the whole area between the late 1960s and the late 1990s, but started to decrease a little at the beginning of this century. The focus of the production shifted from spring to mid- and late summer. The general rise in the level of primary production was mainly due to the increase in nutrient concentrations over the whole Gulf of Finland, but the thermal discharge contributed to a stronger increase of production in the discharge area compared to that in the intake area. The eutrophication of littoral vegetation in the discharge area has been the most obvious, unambiguous and significant biological effect of the heated water. Myriophyllum spicatum, Potamogeton perfoliatus and Potamogeton pectinatus, and vigorous growths of numerous filamentous algae as their epiphytes have strongly increased in the vicinity of the cooling water outlet, where they have formed dense populations in the littoral zone in late summer. However, the strongest increase of phytobenthos has extended only to a distance of about 1 km from the outlet, i.e., the changes in vegetation have been largest in those areas that remain ice-free in winter. Similar trends were also discernible at Olkiluoto, but to a clearly smaller extent, which was due to the definitely weaker level of background eutrophy and nutrient concentrations in the Bothnian Sea, and the differing local hydrographical and biological factors prevailing in the Olkiluoto area. The level of primary production has also increased at Olkiluoto, but has remained at a clearly lower level than at Loviisa. In spite of the analogous changes observed in the macrozoobenthos, the benthic fauna has remained strong and diversified in the Olkiluoto area. Small amounts of local discharge nuclides were regularly detected in environmental samples taken from the discharge areas: tritium in seawater samples, and activation products, such as 60Co, 58Co, 54Mn, 110mAg, 51Cr, in suspended particulate matter, bottom sediments and in several indicator organisms (e.g., periphyton and Fucus vesiculosus) that effectively accumulate radioactive substances from the medium. The tritium discharges and the consequent detection frequency and concentrations of tritium in seawater were higher at Loviisa, but the concentrations of the activation products were higher at Olkiluoto, where traces of local discharge nuclides were also observed over a clearly wider area, due to the better exchange of water than at Loviisa, where local discharge nuclides were only detected outside Hästholmsfjärden Bay quite rarely and in smaller amounts. At the farthest, an insignificant trace amount (0.2 Bq kg-1 d.w.) of 60Co originating from Olkiluoto was detected in Fucus at a distance of 137 km from the power plant. Discharge nuclides from the local nuclear power plants were almost exclusively detected at the lower trophic levels of the ecosystems. Traces of local discharge nuclides were very seldom detected in fish, and even then only in very low quantities. As a consequence of the reduced discharges, the concentrations of local discharge nuclides in the environment have decreased noticeably in recent years at both Loviisa and Olkiluoto. Although the concentrations in environmental samples, and above all, the discharge data, are presented as seemingly large numbers, the radiation doses caused by them to the population and to the biota are very low, practically insignificant. The effects of the thermal discharges have been more significant, at least to the wildlife in the discharge areas of the cooling water, although the area of impact has been relatively small. The results show that the nutrient level and the exchange of water in the discharge area of a nuclear power plant are of crucial importance.

Hormone therapy in elderly women; a comparative study with alendronate of the effects on bone, cardiovascular risk factors, periodontal conditions and quality of life

Thirty percent of 70-year-old women have osteoporosis; after age of 80 its prevalence is up to 70%. Postmenopausal women with osteoporosis seem to be at an increased risk for cardiovascular events, and deterioration of oral health, as shown by attachment loss of teeth, which is proportional to the severity of osteoporosis. Osteoporosis can be treated with many different medication, e.g. estrogen and alendronate. We randomized 90 elderly osteoporotic women (65-80 years of age) to receive hormone therapy (HT)(2mg E2+NETA), 10mg alendronate, and their combination for two years and compared their effects on bone mineral density (BMD) and turnover, two surrogate markers of the risk of cardiovascular diseases, C-reactive protein (CRP) and E-selectin, as well as oral health. The effect of HT on health-related quality of life (HRQoL) was studied in the population-based cohort of 1663 postmenopausal women (mean age 68 yr) (585 estrogen users and 1078 non-users). BMD was measured with dual-energy X-ray absorptiometry (DXA) at 0, 12 and 24 months. Urinary N-telopeptide (NTX) of type I collagen, a marker of bone resorption, and serum aminoterminal propeptide of human type I procollagen (PINP), a marker of bone formation, were measured every six months of treatment. Serum CRP and E-selectin, were measured at 0, 6, and 12 months. Dental, and periodontal conditions, and gingival crevicular fluid (GCF) matrix metalloproteinase (MMP)-8 levels were studied to evaluate the oral health status and for the mouth symptoms a structured questionnaire was used. The HRQoL was measured with 15D questionnaire. Lumbar spine BMD increased similarly in all treatment groups (6.8-8.4% and 9.1-11.2%). Only HT increased femoral neck BMD at both 12 (4.9%) and 24 months (5.8%), at the latter time point the HT group differed significantly from the other groups. HT reduced bone marker levels of NTX and PINP significantly less than other two groups.Oral HT significantly increased serum CRP level by 76.5% at 6 and by 47.1% (NS) at 12 months, and decreased serum E-selectin level by 24.3% and 30.0%. Alendronate had no effect on these surrogate markers. Alendronate caused a decrease in the resting salivary flow rate and tended to increase GCF MMP-8 levels. Otherwise, there was no effect on the parameters of oral health. HT improved the HRQoL of elderly women significantly on the dimensions of usual activities, vitality and sexual activity, but the overall improvement in HRQoL was neither statistically significant nor clinically important. In conclusion, bisphosphonates might be the first option to start the treatment of postmenopausal osteoporosis in the old age.

Cardiac function before and after treatment for various types of loading conditions and in myocardial restriction : A prospective study on pediatric patients with two- and three-dimensional echocardiography and measurement of serum natriuretic peptides

Background: The incidence of all forms of congenital heart defects is 0.75%. For patients with congenital heart defects, life-expectancy has improved with new treatment modalities. Structural heart defects may require surgical or catheter treatment which may be corrective or palliative. Even those with corrective therapy need regular follow-up due to residual lesions, late sequelae, and possible complications after interventions. Aims: The aim of this thesis was to evaluate cardiac function before and after treatment for volume overload of the right ventricle (RV) caused by atrial septal defect (ASD), volume overload of the left ventricle (LV) caused by patent ductus arteriosus (PDA), and pressure overload of the LV caused by coarctation of the aorta (CoA), and to evaluate cardiac function in patients with Mulibrey nanism. Methods: In Study I, of the 24 children with ASD, 7 underwent surgical correction and 17 percutaneous occlusion of ASD. Study II had 33 patients with PDA undergoing percutaneous occlusion. In Study III, 28 patients with CoA underwent either surgical correction or percutaneous balloon dilatation of CoA. Study IV comprised 26 children with Mulibrey nanism. A total of 76 healthy voluntary children were examined as a control group. In each study, controls were matched to patients. All patients and controls underwent clinical cardiovascular examinations, two-dimensional (2D) and three-dimensional (3D) echocardiographic examinations, and blood sampling for measurement of natriuretic peptides prior to the intervention and twice or three times thereafter. Control children were examined once by 2D and 3D echocardiography. M-mode echocardiography was performed from the parasternal long axis view directed by 2D echocardiography. The left atrium-to-aorta (LA/Ao) ratio was calculated as an index of LA size. The end-diastolic and end-systolic dimensions of LV as well as the end-diastolic thicknesses of the interventricular septum and LV posterior wall were measured. LV volumes, and the fractional shortening (FS) and ejection fraction (EF) as indices of contractility were then calculated, and the z scores of LV dimensions determined. Diastolic function of LV was estimated from the mitral inflow signal obtained by Doppler echocardiography. In three-dimensional echocardiography, time-volume curves were used to determine end-diastolic and end-systolic volumes, stroke volume, and EF. Diastolic and systolic function of LV was estimated from the calculated first derivatives of these curves. Results: (I): In all children with ASD, during the one-year follow-up, the z score of the RV end-diastolic diameter decreased and that of LV increased. However, dilatation of RV did not resolve entirely during the follow-up in either treatment group. In addition, the size of LV increased more slowly in the surgical subgroup but reached control levels in both groups. Concentrations of natriuretic peptides in patients treated percutaneously increased during the first month after ASD closure and normalized thereafter, but in patients treated surgically, they remained higher than in controls. (II): In the PDA group, at baseline, the end-diastolic diameter of LV measured over 2SD in 5 of 33 patients. The median N-terminal pro-brain natriuretic peptide (proBNP) concentration before closure measured 72 ng/l in the control group and 141 ng/l in the PDA group (P = 0.001) and 6 months after closure measured 78.5 ng/l (P = NS). Patients differed from control subjects in indices of LV diastolic and systolic function at baseline, but by the end of follow-up, all these differences had disappeared. Even in the subgroup of patients with normal-sized LV at baseline, the LV end-diastolic volume decreased significantly during follow-up. (III): Before repair, the size and wall thickness of LV were higher in patients with CoA than in controls. Systolic blood pressure measured a median 123 mm Hg in patients before repair (P < 0.001) and 103 mm Hg one year thereafter, and 101 mm Hg in controls. The diameter of the coarctation segment measured a median 3.0 mm at baseline, and 7.9 at the 12-month (P = 0.006) follow-up. Thicknesses of the interventricular septum and posterior wall of the LV decreased after repair but increased to the initial level one year thereafter. The velocity time integrals of mitral inflow increased, but no changes were evident in LV dimensions or contractility. During follow-up, serum levels of natriuretic peptides decreased correlating with diastolic and systolic indices of LV function in 2D and 3D echocardiography. (IV): In 2D echocardiography, the interventricular septum and LV posterior wall were thicker, and velocity time integrals of mitral inflow shorter in patients with Mulibrey nanism than in controls. In 3D echocardiography, LV end-diastolic volume measured a median 51.9 (range 33.3 to 73.4) ml/m² in patients and 59.7 (range 37.6 to 87.6) ml/m² in controls (P = 0.040), and serum levels of ANPN and proBNP a median 0.54 (range 0.04 to 4.7) nmol/l and 289 (range 18 to 9170) ng/l, in patients and 0.28 (range 0.09 to 0.72) nmol/l (P < 0.001) and 54 (range 26 to 139) ng/l (P < 0.001) in controls. They correlated with several indices of diastolic LV function. Conclusions (I): During the one-year follow-up after the ASD closure, RV size decreased but did not normalize in all patients. The size of the LV normalized after ASD closure but the increase in LV size was slower in patients treated surgically than in those treated with the percutaneous technique. Serum levels of ANPN and proBNP were elevated prior to ASD closure but decreased thereafter to control levels in patients treated with the percutaneous technique but not in those treated surgically. (II): Changes in LV volume and function caused by PDA disappeared by 6 months after percutaneous closure. Even the children with normal-sized LV benefited from the procedure. (III): After repair of CoA, the RV size and the velocity time integrals of mitral inflow increased, and serum levels of natriuretic peptides decreased. Patients need close follow-up, despite cessation of LV pressure overload, since LV hypertrophy persisted even in normotensive patients with normal growth of the coarctation segment. (IV): In children with Mulibrey nanism, the LV wall was hypertrophied, with myocardial restriction and impairment of LV function. Significant correlations appeared between indices of LV function, size of the left atrium, and levels of natriuretic peptides, indicating that measurement of serum levels of natriuretic peptides can be used in the clinical follow-up of this patient group despite its dependence on loading conditions.

Effects of climate and morphology on temperature conditions of lakes

To a large extent, lakes can be described with a one-dimensional approach, as their main features can be characterized by the vertical temperature profile of the water. The development of the profiles during the year follows the seasonal climate variations. Depending on conditions, lakes become stratified during the warm summer. After cooling, overturn occurs, water cools and an ice cover forms. Typically, water is inversely stratified under the ice, and another overturn occurs in spring after the ice has melted. Features of this circulation have been used in studies to distinguish between lakes in different areas, as basis for observation systems and even as climate indicators. Numerical models can be used to calculate temperature in the lake, on the basis of the meteorological input at the surface. The simple form is to solve the surface temperature. The depth of the lake affects heat transfer, together with other morphological features, the shape and size of the lake. Also the surrounding landscape affects the formation of the meteorological fields over the lake and the energy input. For small lakes the shading by the shores affects both over the lake and inside the water body bringing limitations for the one-dimensional approach. A two-layer model gives an approximation for the basic stratification in the lake. A turbulence model can simulate vertical temperature profile in a more detailed way. If the shape of the temperature profile is very abrupt, vertical transfer is hindered, having many important consequences for lake biology. One-dimensional modelling approach was successfully studied comparing a one-layer model, a two-layer model and a turbulence model. The turbulence model was applied to lakes with different sizes, shapes and locations. Lake models need data from the lakes for model adjustment. The use of the meteorological input data on different scales was analysed, ranging from momentary turbulent changes over the lake to the use of the synoptical data with three hour intervals. Data over about 100 past years were used on the mesoscale at the range of about 100 km and climate change scenarios for future changes. Increasing air temperature typically increases water temperature in epilimnion and decreases ice cover. Lake ice data were used for modelling different kinds of lakes. They were also analyzed statistically in global context. The results were also compared with results of a hydrological watershed model and data from very small lakes for seasonal development.

Atmospheric particle formation in spatially and temporally varying conditions

Atmospheric particles affect the radiation balance of the Earth and thus the climate. New particle formation from nucleation has been observed in diverse atmospheric conditions but the actual formation path is still unknown. The prevailing conditions can be exploited to evaluate proposed formation mechanisms. This study aims to improve our understanding of new particle formation from the view of atmospheric conditions. The role of atmospheric conditions on particle formation was studied by atmospheric measurements, theoretical model simulations and simulations based on observations. Two separate column models were further developed for aerosol and chemical simulations. Model simulations allowed us to expand the study from local conditions to varying conditions in the atmospheric boundary layer, while the long-term measurements described especially characteristic mean conditions associated with new particle formation. The observations show statistically significant difference in meteorological and back-ground aerosol conditions between observed event and non-event days. New particle formation above boreal forest is associated with strong convective activity, low humidity and low condensation sink. The probability of a particle formation event is predicted by an equation formulated for upper boundary layer conditions. The model simulations call into question if kinetic sulphuric acid induced nucleation is the primary particle formation mechanism in the presence of organic vapours. Simultaneously the simulations show that ignoring spatial and temporal variation in new particle formation studies may lead to faulty conclusions. On the other hand, the theoretical simulations indicate that short-scale variations in temperature and humidity unlikely have a significant effect on mean binary water sulphuric acid nucleation rate. The study emphasizes the significance of mixing and fluxes in particle formation studies, especially in the atmospheric boundary layer. The further developed models allow extensive aerosol physical and chemical studies in the future.

Water content variations of sandy soil under winter conditions

Yhteenveto: Maankosteusvaihtelut talvella hiekkamaassa

Snow cover models in operational watershed forecasting

Yhteenveto: Lumimallit vesistöjen ennustemalleissa

On the use of weather radar for mesoscale applications in northern conditions

Mesoscale weather phenomena, such as the sea breeze circulation or lake effect snow bands, are typically too large to be observed at one point, yet too small to be caught in a traditional network of weather stations. Hence, the weather radar is one of the best tools for observing, analyzing and understanding their behavior and development. A weather radar network is a complex system, which has many structural and technical features to be tuned, from the location of each radar to the number of pulses averaged in the signal processing. These design parameters have no universal optimal values, but their selection depends on the nature of the weather phenomena to be monitored as well as on the applications for which the data will be used. The priorities and critical values are different for forest fire forecasting, aviation weather service or the planning of snow ploughing, to name a few radar-based applications. The main objective of the work performed within this thesis has been to combine knowledge of technical properties of the radar systems and our understanding of weather conditions in order to produce better applications able to efficiently support decision making in service duties for modern society related to weather and safety in northern conditions. When a new application is developed, it must be tested against ground truth . Two new verification approaches for radar-based hail estimates are introduced in this thesis. For mesoscale applications, finding the representative reference can be challenging since these phenomena are by definition difficult to catch with surface observations. Hence, almost any valuable information, which can be distilled from unconventional data sources such as newspapers and holiday shots is welcome. However, as important as getting data is to obtain estimates of data quality, and to judge to what extent the two disparate information sources can be compared. The presented new applications do not rely on radar data alone, but ingest information from auxiliary sources such as temperature fields. The author concludes that in the future the radar will continue to be a key source of data and information especially when used together in an effective way with other meteorological data.

Effect of soil preparation of boreal spruce forest on air and soil temperature conditions in forest regeneration areas.

The effect of scarification, ploughing and cross-directional plouhing on temperature conditions in the soil and adjacent air layer have been studied during 11 consecutive growth periods by using an unprepared clear-cut area as a control site. The maximum and minimum temperatures were measured daily in the summer months, and other temperature observations were made at four-hour intervals by means of a Grant measuring instrument. The development of the seedling stand was also followed in order to determine its shading effect on the soil surface. Soil preparation decreased the daily temperature amplitude of the air at the height of 10 cm. The maximum temperatures on sunny days were lower in the tilts of the ploughed and in the humps of the cross-directional ploughed sites compared with the unprepared area. Correspondingly, the night temperatures were higher and so the soil preparation considerably reduced the risk of night frost. In the soil at the depth of 5 cm, soil preparation increased daytime temperatures and reduced night temperatures compared with unprepared area. The maximum increase in monthly mean temperatures was almost 5 °C, and the daily variation in the surface parts of the tilts and humps increased so that excessively high temperatures for the optimal growth of the root system were measured from time to time. The temperature also rose at the depths of 50 and 100 cm. Soil preparation also increased the cumulative temperature sum. The highest sums accumulated during the summer months were recorded at the depth of 5 cm in the humps of cross-directional ploughed area (1127 dd.) and in the tilts of the ploughed area (1106 dd.), while the corresponding figure in the unprepared soil was 718 dd. At the height of 10 cm the highest temperature sum was 1020 dd. in the hump, the corresponding figure in the unprepared area being 925 dd. The incidence of high temperature amplitudes and percentage of high temperatures at the depth of 5 cm decreased most rapidly in the humps of cross-directional ploughed area and in the ploughing tilts towards the end of the measurement period. The decrease was attributed principally to the compressing of tilts, the ground vegetation succession and the growth of seedlings. The mean summer temperature in the unprepared area was lower than in the prepared area and the difference did not diminish during the period studied. The increase in temperature brought about by soil preparation thus lasts at least more than 10 years.

Decomposition of Scots pine fine woody debris in boreal conditions

"Litter quality and environmental effects on Scots pine (Pinus sylvestris L.) fine woody debris (FWD) decomposition were examined in three forestry-drained peatlands representing different site types along a climatic gradient from the north boreal (Northern Finland) to south (Southern Finland) and hemiboreal (Central Estonia) conditions. Decomposition (percent mass loss) of FWD with diameter <= 10 mm (twigs) and FWD with diameter > 10 mm (branches) was measured using the litter bag method over 1-4-year periods. Overall, decomposition rates increased from north to south, the rate constants (k values) varying from 0.128 to 0.188 year(-1) and from 0.066 to 0.127 year(-1) for twigs and branches, respectively. On average, twigs had lost 34%, 19% and 19%, and branches 25%, 17% and 11% of their initial mass after 2 years of decomposition at the hemiboreal, south boreal and north boreal sites, respectively. After 4 years at the south boreal site the values were 48% for twigs and 42% for branches. Based on earlier studies, we suggest that the decomposition rates that we determined may be used for estimating Scots pine FWD decomposition in the boreal zone, also in upland forests. Explanatory models accounted for 50.4% and 71.2% of the total variation in FWD decomposition rates when the first two and all years were considered, respectively. The variables most related to FWD decomposition included the initial ash, water extractives and Klason lignin content of litter, and cumulative site precipitation minus potential evapotranspiration. Simulations of inputs and decomposition of Scots pine FWD and needle litter in south boreal conditions over a 60-year period showed that 72 g m(-2) of organic matter from FWD vs. 365 g m(-2) from needles accumulated in the forest floor. The annual inputs varied from 5.7 to 15.6 g m(-2) and from 92 to 152 g m(-2) for FWD and needles, respectively. Each thinning caused an increase in FWD inputs, Up to 510 g m(-2), while the needle inputs did not change dramatically. Because the annual FWD inputs were lowered following the thinnings, the overall effect of thinnings on C accumulation from FWD was slightly negative. The contribution of FWD to soil C accumulation, relative to needle litter, seems to be rather minor in boreal Scots pine forests. (C) 2008 Elsevier B.V. All rights reserved."

Lumettoman maan routaolojen mallintaminen ja ennustettavuus muuttuvassa ilmastossa

Raportissa on arvioitu ilmastonmuutoksen vaikutusta Suomen maaperän talviaikaiseen jäätymiseen lämpösummien perusteella. Laskelmat kuvaavat roudan paksuutta nimenomaisesti lumettomilla alueilla, esimerkiksi teillä, joilta satanut lumi aurataan pois. Luonnossa lämpöä eristävän lumipeitteen alla routaa on ohuemmin kuin tällaisilla lumettomilla alueilla. Toisaalta luonnollisessa ympäristössä paikalliset erot korostuvat johtuen mm. maalajeista ja kasvillisuudesta. Roudan paksuudet laskettiin ensin perusjakson 1971–2000 ilmasto-oloissa talviaikaisten säähavaintotietoihin pohjautuvien lämpötilojen perusteella. Sen jälkeen laskelmat toistettiin kolmelle tulevalle ajanjaksolle (2010–2039, 2040–2069 ja 2070–2099) kohottamalla lämpötiloja ilmastonmuutosmallien ennustamalla tavalla. Laskelman pohjana käytettiin 19 ilmastomallin A1B-skenaarioajojen keskimäärin simuloimaa lämpötilan muutosta. Tulosten herkkyyden arvioimiseksi joitakin laskelmia tehtiin myös tätä selvästi heikompaa ja voimakkaampaa lämpenemisarviota käyttäen. A1B-skenaarion mukaisen lämpötilan nousun toteutuessa nykyisiä mallituloksia vastaavasti routakerros ohenee sadan vuoden aikana Pohjois-Suomessa 30–40 %, suuressa osassa maan keski- ja eteläosissa 50–70 %. Jo lähivuosikymmeninä roudan ennustetaan ohentuvan 10–30 %, saaristossa enemmän. Mikäli lämpeneminen toteutuisi voimakkaimman tarkastellun vaihtoehdon mukaisesti, roudan syvyys pienenisi tätäkin enemmän. Roudan paksuuden vuosienvälistä vaihtelua ja sen muuttumista tulevaisuudessa pyrittiin myös arvioimaan. Leutoina talvina routa ohenee enemmän kuin normaaleina tai ankarina pakkastalvina. Päivittäistä sään vaihtelua simuloineen säägeneraattorin tuottamassa aineistoissa esiintyi kuitenkin liian vähän hyvin alhaisia ja hyvin korkeita lämpötiloja. Siksi näitten lämpötilatietojen pohjalta laskettu roudan paksuuskin ilmeisesti vaihtelee liian vähän vuodesta toiseen. Kelirikkotilanteita voi esiintyä myös kesken routakauden, jos useamman päivän suojasää ja samanaikainen runsas vesisade pääsevät sulattamaan maata. Tällaiset routakauden aikana sattuvat säätilat näyttävätkin yleistyvän lähivuosikymmeninä. Vuosisadan loppua kohti ne sen sijaan maan eteläosissa jälleen vähenevät, koska routakausi lyhenee oleellisesti. Tulevia vuosikymmeniä koskevien ilmastonmuutosennusteiden ohella routaa ja kelirikon esiintymistä on periaatteessa mahdollista ennustaa myös lähiaikojen sääennusteita hyödyntäen. Pitkät, viikkojen tai kuukausien mittaiset sääennusteet eivät tosin ole ainakaan vielä erityisen luotettavia, mutta myös lyhyemmistä ennusteista voisi olla hyötyä mm. tienpitoa suunniteltaessa.