Assessing the effects of subsurface drainage on hydrology and nitrogen transport in Nordic fields

Salaojitus on pääasiallinen maatalouden vesien hallintamenetelmä pohjoisilla alueilla. Uusien salaojitusten asentaminen sekä vanhojen salaojitusten täydennysojittaminen muuttavat veden virtausreittejä maanpinnalla ja maaperässä sekä vaikuttavat pellon vesitaseeseen. Nämä muutokset vaikuttavat pellolla tapahtuviin typen prosesseihin sekä typpikuormaan. Salaojitusmenetelmän vaikutuksien arviointia monimutkaistavat kuivatukseen vaikuttavat sääolot, maaperän hydrauliset ominaisuudet sekä peltoa ympäröivät alueet. Tutkimuksen tavoitteena oli selvittää kuivatusjärjestelmien merkitystä veden virtaukseen ja typen kulkeutumiseen kerättyjen mittausaineistojen, tilastollisten menetelmien sekä matemaattisen mallintamisen avulla.

Kahden salaojitusmenetelmän toimintaa Sievin koepellolla tutkittiin tilastollisilla menetelmillä. Erot pohjavedenpinnan tasoissa johtuivat salaojitusmenetelmästä sekä maalajista salaojasyvyydessä. Absoluuttiset erot eri ojitusmenetelmillä käsiteltyjen lohkojen välillä olivat keskimäärin vain 0,1 m. Maan ominaisuuksien ja kuivatusjärjestelmän vaikutuksia veden virtausreitteihin tutkittiin tarkemmin prosessipohjaisella hydrologisella mallilla. Salaojitusvaihtoehtojen vaikutuksia simuloitiin 1D, 2D ja 3D mallisovelluksissa, joissa käytettiin mittausaineistoa eteläsuomalaiselta savipellolta (Nummela). Kuivatusjärjestelmä voitiin kuvata ojaväli- ja peltolohkomittakaavan malleilla. 3D ojavälisimulaatioista saatiin tietoa, kuinka maanäytteiden mittaustietojen suora käyttö mallin parametrisoinnissa voi toimia mallin kalibroinnin vaihtoehtona. Ojavälimittakaavan parametrisointi skaalattiin peltolohkomittakaavaan pitkän ajan 2D simulaatioita varten. Lyhyen ajan 3D simulaatioiden mukaan veden virtausta dominoi mallissa kuvatut makrohuokoset verrattuna kuivatusjärjestelmän vaikutuksiin. Pitkäaikaisten 2D simulaatiotulosten vertailu paljasti, että kuivatuksen tehostaminen muutti pellon valuntareittien suhteita, laskien myös pohjavesivalunnan määrään. Työssä kehitettiin yleinen aineiden kulkeutumista kuvaava matemaattinen malli ja sillä tutkittiin typen kierron prosesseja ja typen kulkeutumista. Typpimallia testattiin 1D sovelluksissa lyhyille syysjaksoille, jolloin huuhtoumaa kontrolloi pääosin kasvukauden jälkeiset maan typpivarastot ja sateiden ajoittuminen. Maan kosteus selitti kaasumaisen typpivirran suuruutta.

Pitkäaikaiset mittausaikasarjat, tilastolliset analyysimenetelmät ja prosessi-pohjainen mallinnus osoittivat, että salaojituksen vaikutukset ovat kohdekohtaisia. 1D, 2D ja 3D mallisovelluksilla pystyttiin selvittämään, mitkä tekijät vaikuttavat salaojavalunnan määrään, ja mistä alkulähteistä kuivatusvalunta muodostuu. Kuivatusjärjestelmää kuvattiin näissä sovelluksissa eri tavoille, mikä vaikutti veden virtausreitteihin. Peltoviljelyn ravinnekuormituksen vähentäminen vaatii pellon vesitaseen sekä kuormitusmekanismien kokonaisvaltaista ymmärtämistä. Tutkimustulosten mukaan salaojituksen tehostaminen siirtää typpikuormitusta pohjavesivalunnasta salaojavaluntaan eikä näin ollen oleellisesti kasvata typpikuorman määrää.

———————–

Subsurface drainage is the primary water management approach in field cultivation in Nordic areas. Installing and improving a subsurface drainage system change water flow dynamics and routes in the soil, which affect nitrogen (N) load from the field. The role of soil properties, drainage system improvements and surrounding areas in the formation of water flow routes is not fully understood. The objective of this study was to quantify the effects of subsurface drainage on water flow and N transport using field monitoring data, statistical analysis and mathematical modeling.

The performance of two drainage installation methods applied in the Sievi experimental field were investigated with statistical analysis. Differences in groundwater level occurred due to drainage installation method and soil type at the drain depth, but the absolute differences were small (0.1 m). A state-of-the-art process-based hydrological model was applied to investigate the effects of soil properties and drainage systems on water flow routes. Field subsurface drainage schemes were simulated with 3D, 2D and 1D model applications using data from a clayey field in southern Finland (Nummela). Model applications showed how field drainage can be described with models of different dimensions and scale (from drain spacing to field section scale). The 3D drain spacing simulations demonstrated the benefits of using detailed soil data in model parameterization as an alternative to model calibration. In the 2D long-term simulations, the 3D soil parameterization was up-scaled to field section scale. The short-term 3D model simulations showed the dominant nature of soil macropores over the drainage system description. Comparison of the long-term 2D model simulations revealed that the improved drainage installations changed the shares of all the water flow routes, including groundwater outflow. A generic solute transport component was developed and tailored to describe N cycle transport and processes in 1D model simulations. Autumn period simulations of a poorly and well drained field sections showed that nitrate N loading was mainly controlled by the initial soil N storages after harvest and the timing of the precipitation events, while the soil moisture content differences explained the magnitudes of gaseous N losses.

Long-term monitoring data series, statistical analysis and process-based modeling showed that the practical effects of subsurface drainage are site specific and comprehensive view on the local water and nutrient management is needed when controlling the environmental impacts of field cultivation. The 1D, 2D and 3D model applications could all be used to replicate the measured drain discharge data, even though the drainage system description differed between the cases due to the differences in water flow directions and the boundaries of the simulated domain. The finding of the research suggests that field water management moves the N load from one path to another rather than affecting to the total amount of the water volume or N loading.