Browsing by master's degree program "Magisterprogrammet i atmosfärsvetenskaper"

In order to reliably quantify the impact of cities to the climate change, carbon dioxide (CO2) fluxes from urban areas need to be accurately estimated. Currently, eddy covariance (EC) measurements of net ecosystem exchange (NEE) are among the most used methods for CO2 balance assessment at ecosystem scale. They cannot, however, directly separate gross primary production (GPP) from different sources of CO2 , but instead different tracers and methods need to be used for the separation. In this thesis, urban carbonyl sulfide (COS) fluxes are reported and used as a tracer for biogenic CO2 uptake. For estimating GPP, a leaf-relative uptake (LRU) of CO2 and COS during photo- synthesis has to be calculated. Three different methods for ecosystem scale LRU estimations are compared in order to find which performs the best at assessing urban GPP from EC measurements. LRUNEE method is based on averaging LRU obtained directly from an equation using EC measurements of CO2 and COS. LRUP AR uses parameterisation of photosynthetically available radiation (PAR) to estimate LRU, and LRUCAP uses additional environmental parameters and information of soil water content (SWC), water vapour pressure deficit (VPD) and PAR. EC measurements were conducted at SMEAR III station (ICOS ecosystem associate site) in Helsinki during summer 2022 using a quantum cascade laser gas analyser. Mole fractions of COS, CO2, CO, and water vapour were collected with 10 Hz frequency, together with three-dimensional wind speed measurements. 30-min fluxes were calculated and used for estimating LRU and GPP with the three methods above. Source area of the station was divided into urban, street, and vegetation sectors. To compare different LRU methods, EC measured GPP from vegetation sector was used because it had the smallest anthropogenic emissions. CO and COS fluxes were compared to estimate co-emissions from fuel combustion. Median COS flux of −25.5 pmol m−2 s−1 was measured, indicating a biogenic CO2 uptake in the study area. For simplicity, soil related COS fluxes were neglected. LRUPAR performed the best at estimating GPP, giving an average CO2 uptake of 13.9 μmol m−2 s−1 for the vegetation sector. Furthermore, peak median anthropogenic emissions of 9.1 μmol m−2 s−1 were deduced with LRUPAR, which agrees with earlier research. However, due to problems with instrumentation the data collection period was so short that distinction between LRUPAR and LRUCAP was not unambiguous. Clear is that setting a constant value for LRU leads to significant underestimation of GPP, and a misunderstanding of its behaviour in heterogeneous urban landscape. Further development of COS based GPP estimation should include a more careful instrumentation and revision of the parameterisation of LRUCAP method.

In atmospheric sciences, measurements provided by remote-sensing instruments are crucial in observing the state of atmosphere. The associated uncertainties are important in nearly all data analyses. Random uncertainties reported by satellite instruments are typically estimated by inversion algorithms (ex-ante). They can be incomplete due to simplified or incomplete modelling of atmospheric processes used in the retrievals, and thus validating random uncertainties is important. However, such validation of uncertainties (or their estimates from statistical analysis afterwards, i.e. ex-post) is not a trivial task, because atmospheric measurements are obtained from the ever-changing atmosphere. This Thesis aims to explore the structure function method – an important approach in spatial statistics – and apply it to total ozone column measurements provided by the nadir-viewing satellite instrument TROPOMI. This method allows us to simultaneously perform validation of reported ex-ante random uncertainties and to explore of local-scale natural variability of atmospheric parameters. Two-dimensional structure functions of total ozone column have been evaluated based on spatial separations in latitudinal and longitudinal directions over selected months and latitude bands. Our results have indicated that the ex-post random uncertainties estimated agree considerably well with the reported ex-ante random uncertainties, which are within 1-2 DU. Discrepancies between them are very small in general. The morphology of ozone natural variability has also been illustrated: ozone variability is minimal in the tropics throughout the year, whereas in middle latitudes and polar regions they attain maxima in local spring and winter. In every scenario, the ozone structure functions are anisotropic with a stronger variability in the latitudinal direction, except at specific seasons in polar regions where isotropic behaviour is observed. Our analysis has demonstrated that the structure function method is a remarkable and promising tool for validating random uncertainties and exploring natural variability. It has a high potential for applications in other remote sensing measurements and atmospheric model data.

The oxidation mechanisms of atmospheric organic compounds are an important puzzle piece for many atmospherically relevant topics, including but not limited to air quality and climate change. One poorly understood step in this oxidation process is peroxy radical recombination, in some conditions the most important sink reaction for peroxy radicals, which are formed in abundance due to gas phase reactions in the lower troposphere. After a few initial steps, the peroxy radical recombination reaction results in the ejection of O_2 leaving behind a pair of alkoxy radicals in close proximity. This reactive complex has three known reaction pathways: Hydrogen shift forming an alcohol and a carbonyl compound, radical recombination forming a ROOR dimer, and diffusive break-up forming two free alkoxy radicals. In this thesis, alkoxy bond scission followed by radical recombination resulting in the formation of a ROR is proposed as a fourth reaction pathway. To test the hypothesis, computational chemistry was used to determine alkoxy bond scission rates for radicals of atmospheric significance, and gas-phase oxidation experiments were realized on three peroxy radical precursor molecules to look for signs of ROR formation. More precisely, the Eyring equation was used to calculate the rate of alkoxy bond scission on a potential energy surface determined using density functionals, with corrections to electronic energy using coupled-cluster calculations. In the experiments, liquid phase alkenes were vaporized, and oxidized by O_3 in the gas phase, resulting in peroxy radical formation, after which the possible dimers were detected using a NO_3^- -atmospheric pressure chemical ionization time-of-flight mass spectrometer. A highly oxidized radical reaction partner was present in the chamber to improve the detectability of the formed dimers. The combined results of these two approaches suggest that the reaction pathway is possible in standard atmospheric conditions and may thus be important for a number of peroxy radicals.

Tässä tutkielmassa tarkastellaan horisontaalisten gravitaatiovaihteluiden vaikutusta ilmakehän perusyhtälöihin sekä yksinkertaisen ilmakehämallin tuloksiin erilaisissa simulaatioissa. Työn motivointina oli tutkia putoamiskiihtyvyyden vaikutusta mallinnustarkkuuteen, koska se on yksi monista säänennustus- ja ilmastosimulaatioihin liittyvistä epätarkkuustekijöistä. Ilmakehän perusyhtälöt johdettiin aluksi uudelleen huomioimalla gravitaation vaihtelu vaakasuunnassa. Tämän jälkeen vastaavat yhtälömuutokset tehtiin SPEEDY-mallin lähdekoodiin, ja mallin avulla tehtiin simulaatioita gravitaatiovaihteluiden vaikutusten selvittämiseksi. Jotta tulosten analysointi olisi mahdollisimman helppoa, käytettiin simulaatioissa paljon yksinkertaistuksia. Näistä merkittävin oli mallimaapallon korvaaminen vesiplaneetalla. Yhtälömuutosten oikeellisuus mallissa verifioitiin yhden aika-askeleen kokeilla, minkä jälkeen muokatuille perusyhtälöille tehtiin suuruusluokka-analyysi. Analyysin perusteella gravitaatiovaihteluista aiheutuvat lisätermit olivat pääosin yhdestä kahteen kertaluokkaa yhtälöiden muita termejä pienempiä. Lopuksi tehtiin kymmenen vuoden simulaatioita, joissa tarkasteltiin niin sanotun normaaligravitaatiojakauman vaikutuksia mallin tuloksiin. Näissä kokeissa havaittiin, että meteorologisten suureiden anomaliat olivat pääosin maltillisia, mutta eivät merkityksettömän pieniä. Esimerkiksi tuulikentässä havaitut muutokset olivat suurimmillaan noin 2 m/s, kun taas lämpötila-anomaliat jäivät globaalisti alle puoleen asteeseen. Meridionaalisen kiertoliikkeen anomaliassa havaittiin puolestaan selkeä antisymmetria pallonpuoliskojen välillä: intertrooppinen konvergenssivyöhyke siirtyi päiväntasaajalta leveyspiirin 10°S tienoille, kun taas leveyspiirillä 5°N nousuliike heikkeni. Lisäksi länsituulet hidastuivat pohjoisen pallonpuoliskon keskileveysasteilla, mutta voimistuivat eteläisellä pallonpuoliskolla. Tulosten perusteella aiheen tutkimista kannattaa jatkaa myös tulevaisuudessa.

The legislation of the Paris Agreement obliges Finland to pursue actions that keep the global average temperature rise below 2°C and aim to limit the average temperature rise to 1.5°C. The current Finnish government has aligned the national goal of carbon neutrality by 2035. The role of municipalities in promoting or compensating carbon sinks has not yet been defined, although municipalities play an important role as a platform for climate work at local and regional levels. However, it is already known that the Finnish National Climate Act, which is being reformed at this moment, will be subject to an obligation to produce their own climate programs at municipal, regional or provincial level. Environmental competence and environmental development have been important in Lahti for several decades already. The City of Lahti has set its target for carbon neutrality for 2025 and it includes targets for reducing, compensating, and increasing carbon sinks. This work focused on the examination of carbon sequestration and sinks in an urban environment in Lahti, in the example area of approximately 82 hectares, through which a wider understanding of the city's potential to grow coal stocks and sinks in a tight urban structure within different land use classes and different ground cover between them. Based on the Finnish Environment Agency's CORINE land cover classification, the current potential of carbon sequestration for urban land use classes were calculated in this work and the actions to increase carbon sequestration capacity were identified. The work examined the availability of the finished spatial data and to supplement incomplete information, existing literature on the topic was used, as well as other existing spatial records of the city of Lahti and previously made surveys. The largest carbon sink was observed in forest areas, of which in mixed forests representing the largest forest type in the area. Through the calculations and literature carbon sinks and stocks in residential areas were also found to be significant in terms of vegetation, as well as in terms of soil based on the literature review. In planting street and park trees for the purpose of increasing the carbon sink, the most important thing was found to be the long lifetime of trees and securing it. Growing of carbon sinks is most effective in areas where carbon sequestration is already at a high level but increasing vegetation cover in all urban land covers will increase the carbon sink in the long run. One major conclusion of the work was that Lahti's current method of determining carbon sinks and stocks has been inadequate at least for the determining them in built areas, and future measures to maintain, preserve and increase carbon stocks and sinks would not be seen by the same calculation method in the computing. In general, the research data and methods are still largely based on observations and results from the operational processes of natural ecosystems, and these are utilized in urban planning, construction, and maintenance of urban green areas. An incomplete knowledge of the ecological processes in urban areas is a problem that produced challenges in this work as well. More research data is needed on carbon sinks in urban land use classes to gain a more secure understanding of carbon sinks and stocks, although the common importance of vegetation in urban areas is already clear. Although the work focused on carbon in an urban environment, it is necessary to remember the diversity of the urban environment and the other ecosystem services it produces. Land use planning, as well as the management of green spaces in the urban environment, can enhance both the size of carbon storages and sinks and biodiversity and they do not have to be entirely separate from each other.

Low-level wind shear is a significant aviation hazard. A sudden reduction in the headwind along an aircraft flight path can induce a loss of lift, from which an aircraft may not be able to recover when it is close to the ground. Airports therefore use low-level wind shear alert systems to monitor wind velocities within the airport terminal area and alert of any detected hazardous wind shear. There exist three ground-based sensor systems capable of independently observing low-level wind shear: a Doppler weather radar-based, a Doppler wind lidar-based, and an anemometer-based system. However, as no single sensor system is capable of all-weather wind shear observations, multiple alert systems are used simultaneously, and observations from each system are integrated to produce one set of integrated wind shear alerts. Algorithms for integrating Doppler weather radar and anemometer wind shear observations were originally developed in the early 1990s. However, the addition of the Doppler wind lidar-based alert system in more recent years warrants updates to the existing radar/anemometer integration algorithms. This thesis presents four different replacement candidates for the original radar/anemometer integration algorithms. A grid-based integration approach, where observations from different sensor systems are mapped onto a common grid and integrated, is found to best accommodate central integration considerations, and is recommended as the replacement to the original radar/anemometer algorithms in operational use. The grid-based approach is discussed in further detail, and a first possible implementation of the algorithm is presented. In addition, ways of validating the algorithm and adopting it for operational use are outlined.

Gulf of Bothnia was simulated with NEMO sea model and LIM3 sea ice model. The results were used to count ice area, ice thickness, ice season length and the dates for freezing ang thawing. Results contained years from 1975 to 2059 but only years from 2006 were mostly used. Before 2006 the model was forced with a statistical history run. It was compared to thickness observations and the ice model was considered reliable enough. Atmospheric forcings became from three earth system models and two emission scenarios. With smaller emissions results about ice conditions and its trends vary between model runs with different focings a lot more than with higher emissions. In many case a propability distribution was used on results unlike on earlier similar researches which have mostly used only medians of ice parameters. Propability distribution for ice area was made using cumulative Gumbel probability distribution which enabled counting time periods longer than the time series. It also increased the reliability of extreme results. According to the results the northen part of the Bay of Bothnia freezes in every winter at 2050s. Ice thickness is 80cm at the most and median is 50cm. In almost every result southern part of the Bothnian Sea has over 50% chance for ice to occur and the median of thickness is more than 10cm. Length of the ice season close to 2050s can exeed 150 days often at the northen Bay of Bothnia but elsewhere it is closer to 100 days and varies a lot. The counted propability distributions for ice area let us suspect that a part of the Bay of Bothnia would remain unfreezed at least once every 30 years or more propably every 20 years in years from 2006 to 2059. In contradiction, not much can be said about large areas, because results vary a lot about wheter the whole Gulf of Bothnia can be expected to freeze close to 2050s.

The parameterization of deep convection is simulated poorly over the Central and East Pacific. This could lead to issues in predicting the annual total precipitation in the tropics, such as the existence of a double intertropical convergence zone over the equatorial Pacific. Resolving tropical deep convection instead of parameterization leads to the presence of more linear systems. Observations over Atlantic indicate that shear-perpendicular lines (squall lines) propagate faster than shear- parallel lines, mainly due to their connection with the low-level vertical wind shear (VWS). The study examines the different movement speeds of mesoscale convective systems (MCSs) over the East- and West Pacific to determine whether this could explain the reason why climate models have problems with predicting deep convection. A higher proportion of fast moving MCSs (squall lines) could contribute to the prediction problems in the tropics. The MCS motion is determined by the sum of the mean wind and propagation speed. In squall lines, the MCS motion is mainly influenced by the propagation, which is associated with the low-level VWS. Therefore, the effect of the low-level VWS on the fast- and slow moving MCSs is also investigated. The Global High-Resolution Mesoscale Convective System Database is used, which provides infor- mation about the time, location and movement of MCSs. Additionally, ERA5 wind data is used to obtain the mean wind and VWS. Two specific areas over the northern equatorial Pacific are chosen to compare the different types of MCSs. These areas are over the East Pacific (120°W - 140°W and 5°N - 12°N) and West Pacific (140°E - 160°E and 0°N - 7°N). The movement speed is used to categorize the MCSs into three groups: slow moving MCSs (< 3 m/s), moderate moving MCSs (3 m s−1 − 7 m s−1) and fast moving MCSs (> 7 m/s). The study reveals that the share of the fast moving MCSs is 9.8% over the East Pacific and 13.8% over the West Pacific. This is only a 4 percent point difference between the two areas. Therefore, it is not shown that the fast moving MCSs contribute to the existing issues that models have in predicting the annual total precipitation over the East Pacific. Moreover, approximately 85-90% are categorized as slow- to moderate moving MCSs. Hence, the influence of fast moving MCSs is relatively small when compared to the other types. A difference is seen in the mean wind and VWS over the East Pacific, but do not explain the MCS motion vector. Therefore, the difference between fast- and slow moving MCSs cannot be explained by only the monthly averaged mean wind and low-level VWS over the East Pacific. Over the West Pacific, the mean wind direction and VWS are about the same in direction and speed. Therefore, the difference between fast- and slow moving MCSs is not explained by the low-level VWS over the West Pacific.

Tutkielmassa on selvitetty lumensyvyyden muutoksia ERA5-Land-uudelleenanalyysin antamille tuloksille ajanjaksolla 1950-2021. Datan analyysi ja käsittely on toteutettu Pythonilla. ERA5-Land:n etuihin lukeutuu muun muassa parempi erotuskyky kuin ERA5-uudelleenanalyysiin, mikä parantaa aineiston tarkkuutta huomattavasti. ERA5-Land ei suoraan käytä havaintoarvoja vaan lumensyvyys lasketaan muitten sääsuureitten, kuten lämpötilan ja sademäärän, aikasarjoja hyödyntäen. Näin ollen tutkielmassa käsiteltyihin suureisiin on sisällytetty muitakin suureita kuin vain lumensyvyys jotta syitä lumensyvyyden muutoksiin olisi helpompaa hahmottaa. Tutkielmaan valikoitiin kolme tutkimusaluetta; Suomi kokonaisuudessaan, Fennoskandian alueelta Norjan, Ruotsin sekä Suomen yhteinen pinta-ala, sekä Japanista Hokkaidon ja Honshun saarten muodostama maa-alue. Aineiston pohjalta voidaan todeta varsin yksiselitteisesti lumensyvyyden kuukausikeskiarvojen olevan, varsinkin lumensyvyyden huippuarvokuukausina, vertailukaudelle 1991-2020 pienempiä kuin vertailukaudelle 1951-1980. Sama trendi näkyy myös vuotuisten keskiarvojen kehityksessä. Muutosten suuruus on jonkin verran sidoksissa alueellisiin erityispiirteisiin mutta päätrendi on kaikille alueille jokseenkin samansuuntainen; nykyisentyyppinen ilmastollinen kehitystrendi laskee lumensyvyyttä pitkällä aikatähtäimellä. Toisin sanoen ilmastonmuutos vähentää lumimäärää pohjoisella pallonpuoliskolla.

Mars-planeetan kaasukehä koostuu enimmäkseen hiilidioksidista, kun taas vesihöyryä on hyvin vähän. Kaasukehän lämpötila vaihtelee noin +10 ja -130 Celsius-asteen välillä ja pintapaine on vain noin sadasosa Maan ilmakehän paineesta. Marsin kaasukehässä on usein paljon hienojakoista pölyä, joka absorboi tehokkaasti auringonsäteilyä ja täten vaikuttaa kaasukehän toimintaan. Marsin pinnan reagoidessa erittäin nopeasti auringonsäteilyn määrän muutoksiin sekä kaasukehässä olevan pölyn vuoksi rajakerroksen mallinnuksessa käytettävissä malleissa säteilyn parametrisaatioiden täytyy olla mahdollisimman hyviä. Helsingin yliopisto ja Ilmatieteen laitos ovat kehittäneet Marsin kaasukehän tutkimukseen tarkoitetun 1-ulotteisen pylväsmallin. Malli on erittäin nopea ja helposti muokattavissa, joten sillä voidaan testata uusia ilmakehäfysiikan lainalaisuuksia ja algoritmeja, joita voidaan mahdollisesti lisätä kolmiulotteisiin Marsin kaasukehän malleihin. Tämä työ tehtiin osana Ilmatieteen laitoksen Marsin tutkimusryhmää ja työssä tutustutaan Marsin kaasukehän rajakerrokseen sekä pylväsmalliin. Lisäksi mallin antamia tuloksia esitellään ja verrataan Curiosity mönkijän (toiselta nimeltään Mars Science Laboratory, MSL) havaintoihin sekä tutkitaan mallin herkkyyttä sen alustusparametreihin. Mallin ennustamia lämpötilan, vesihöyryn tilavuuden sekoitussuhteen ja suhteellisen kosteuden vuorokausisyklejä verrattiin MSL:n havaintoihin eri vuodenaikoina. MSL laskeutui vuonna 2012 lähelle Marsin päiväntasaajaa Gale-kraatterin pohjalle ja se sisältää Ilmatieteen laitoksen suunnittelemat ja rakentamat mittalaitteet paineelle ja suhteelliselle kosteudelle. Mallin ennustamat vuorokausisyklit vastasivat hyvin mönkijän havaintoja ja tuloksista nähtiin myös lämpötilan suuri vuorokausivaihtelu kaasukehän reagoidessa nopeasti auringonsäteilyn muutoksiin. MSL:n paineen mittauksista (yli 3000 Marsin vuorokautta) nähtiin selvästi hiilidioksidin vuodenaikaiskierto etelänavalta pohjoisnavalle ja päinvastoin. Lisäksi vuoden 2018 globaali pölymyrsky näkyi monissa eri mittaustuloksissa. Mallin herkkyyttä tutkittiin muuttamalla neljää eri alustusparametria: pinnan lämpötilaa ja painetta, ilmapylvään vesisisältöä (PWC) sekä pölyn optista paksuutta (tau). Näiden testien perusteella mallin ennustamiin vuorokauden lämpötilaprofiileihin eniten vaikuttivat pinnan lämpötilan ja pölyn optisen paksuuden alustus, kun taas kosteusprofiileihin eniten vaikuttivat PWC:n ja pölyn optisen paksuuden alustus. Näistä parametreista pinnan paineen alustuksella oli vähiten vaikutusta mallin ennustamiin profiileihin.