Assessing the combined effect of land management and wildfire on runoff and soil erosion in north central Portugal

Strong and sometimes extreme responses in runoff and soil erosion following wildfires have been reported worldwide. However, in the case of North-Central Portugal, little research had been carried out regarding the hydrologic and erosive impacts of several land management activities in recently burnt areas (such as ground preparation, post-fire logging or post-fire mitigation treatments). This study aims to assess post-fire runoff and soil erosion response on Eucalypt and Maritime pine plantations during the first, second and third years following wildfires. The effect of several pre-fire ground preparation operations (ploughed down-slope, contour ploughed and inclined terraces), post-fire logging activities (on both the eucalypt and pine plantations), as well as the application of hydromulch (a post-fire emergency treatment) on overland flow and soil erosion were compared to burnt but undisturbed and untreated areas. The intensive monitoring of runoff, soil erosion and selected soil properties served to determine the main factors involved in post-fire runoff and soil erosion and their spatial and temporal variation. Soil water repellency deserved special attention, due to its supposed important role for overland flow generation. Repeated rainfall simulation experiments (RSE’s), micro-scale runoff plots and bounded sediment fences were carried out and/or installed immediately after the wildfire on seven burnt slopes. Micro-scale runoff plots results under natural rainfall conditions were also compared to the RSE’s results, which was useful for assessing the representativeness of the data obtained with artificial rainfall. The results showed comparable runoff coefficient (20-60%) but lower sediment losses (125-1000 g m-2) than prior studies in Portugal, but especially outside Portugal. Lower sediment losses were related with the historic intensive land use in the area. In evaluating these losses, however, the shallowness and stoniness of the soils, as well as the high organic matter fraction of the eroded sediments (50%) must not be overlooked. Sediment limited erosion was measured in all the ploughed sites, probably due to the time since ploughing (several years). The disturbance of the soil surface cover due to post-fire logging and wood extraction substantially increased sediment losses at both the pine and eucalypt sites. Hydromulch effectiveness in reducing the runoff (70%) and sediment losses (83%) was attributed to the protective high coverage provided by hydromulch. The hydromulch significantly affected the soil cover and other soil properties and these changes also reduced the soil erosion risk. The rainfall amount was the main factor explaining the variance in runoff. However, a shift from rainfall amount to rainfall intensity was detected when either the surface cover or the infiltration capacity (hydrophilic conditions) increased. Sediment losses were controlled by rainfall intensity and surface cover. The role of soil water repellency on runoff generation was not consistent; the overall repellency levels alone were not enough to assess its hydrological impact. Soil water repellency explained runoff generation in the specific-sites model better than in the overall model. Additionally, soil moisture content was a better predictor for soil water repellency than antecedent rainfall. The natural rainfall results confirmed that RSE’s were able to capture the specific sediment losses and its organic matter content as well as the differences between the ploughed and unploughed sites. Repeated RSE’s also captured the seasonal variations in runoff and sediment losses attributed to soil water repellency. These results have implications for post-fire soil erosion modelling and soil conservation practices in the region, or areas with the same land use, climate and soil characteristics. The measured sediment loss, as well as the increasing frequency of ploughing in recently burnt and unburnt eucalypt stands, suggests ploughing is not an effective as a soil conservation measure. Logging activities with less impact are recommended in order to maintain the forest litter protecting the soil surface. Due to its high effectiveness in reducing runoff and soil erosion, hydromulch is recommended for highly sensitive and vulnerable areas.

Understanding and modelling hydrological and soil erosion processes in burnt forest catchments

Forest fires implications in overland flow and soil erosion have been researched for several years. Therefore, is widely known that fires enhance hydrological and geomorphological activity worldwide as also in Mediterranean areas. Soil burn severity has been widely used to describe the impacts of fire on soils, and has being recognized as a decisive factor controlling post-fire erosion rates. However, there is no unique definition of the term and the relationship between soil burn severity and post-fire hydrological and erosion response has not yet been fully established. Few studies have assessed post-fire erosion over multiple years, and the authors are aware of none which assess runoff. Small amount of studies concerning pre-fire management practices were also found. In the case of soil erosion models, the Revised Universal Soil Loss Equation (RUSLE) and the revised Morgan–Morgan–Finney (MMF) are well-known models, but not much information is available as regards their suitability in predicting post-fire soil erosion in forest soils. The lack of information is even more pronounced as regards post-fire rehabilitation treatments. The aim of the thesis was to perform an extensive research under the post fire hydrologic and erosive response subject. By understanding the effect of burn severity in ecosystems and its implications regarding post fire hydrological and erosive responses worldwide. Test the effect of different pre-fire land management practices (unplowed, downslope plowed and contour plowed) and time-since-fire, in the post fire hydrological and erosive response, between the two most common land uses in Portugal (pine and eucalypt). Assess the performance of two widely-known erosion models (RUSLE and Revised MMF), to predict soil erosion rates during first year following two wildfires of distinctive burn severity. Furthermore, to apply these two models considering different post-fire rehabilitation treatments in an area severely affected by fire. Improve model estimations of post-fire runoff and erosion rates in two different land uses (pine and eucalypt) using the revised MMF. To assess these improvements by comparing estimations and measurements of runoff and erosion, in two recently burned sites, as also with their post fire rehabilitation treatments. Model modifications involved: (1) focusing on intra-annual changes in parameters to incorporate seasonal differences in runoff and erosion; and (2) inclusion of soil water repellency in runoff predictions. Additionally, validate these improvements with the application of the model to other pine and eucalypt sites in Central Portugal. The review and meta-analysis showed that fire occurrence had a significant effect on the hydrological and erosive response. However, this effect was only significantly higher with increasing soil burn severity for inter-rill erosion, and not for runoff. This study furthermore highlighted the incoherencies between existing burn severity classifications, and proposed an unambiguous classification. In the case of the erosion plots with natural rainfall, land use factor affected annual runoff while land management affected both annual runoff and erosion amounts significantly. Time-since-fire had an important effect in erosion amounts among unplowed sites, while for eucalypt sites time affected both annual runoff and erosion amounts. At all studied sites runoff coefficients increase over the four years of monitoring. In the other hand, sediment concentration in the runoff, recorded a decrease during the same period. Reasons for divergence from the classic post-fire recovery model were also explored. Short fire recurrence intervals and forest management practices are viewed as the main reasons for the observed severe and continuing soil degradation. The revised MMF model presented reasonable accuracy in the predictions while the RUSLE clearly overestimated the observed erosion rates. After improvements: the revised model was able to predict first-year post-fire plot-scale runoff and erosion rates for both forest types, these predictions were improved both by the seasonal changes in the model parameters; and by considering the effect of soil water repellency on the runoff, individual seasonal predictions were considered accurate, and the inclusion of the soil water repellency in the model also improved the model at this base. The revised MMF model proved capable of providing a simple set of criteria for management decisions about runoff and erosion mitigation measures in burned areas. The erosion predictions at the validation sites attested both to the robustness of the model and of the calibration parameters, suggesting a potential wider application.

Biogeochemical response of Tagus Estuary to climate change : a modelling study

Estuaries are highly dynamic systems which may be modified in a climate change context. These changes can affect the biogeochemical cycles. Among the major impacts of climate change, the increasing rainfall events and sea level rise can be considered. This study aims to research the impact of those events in biogeochemical dynamics in the Tagus Estuary, which is the largest and most important estuary along the Portuguese coast. In this context a 2D biophysical model (MOHID) was implemented, validated and explored, through comparison with in-situ data. In order to study the impact of extreme rainfall events, which can be characterized by an high increase in freshwater inflow, three scenarios were set by changing the inputs from the main tributaries, Tagus and Sorraia Rivers. A realistic scenario considering one day of Tagus and Sorraia River extreme discharge, a scenario considering one day of single extreme discharge of the Tagus River and finally one considering the extreme runoff just from Sorraia River. For the mean sea level rise, two scenarios were also established. The first with the actual mean sea level value and the second considering an increase of 0.42 m. For the extreme rainfall events simulations, the results suggest that the biogeochemical characteristics of the Tagus Estuary are mainly influenced by Tagus River discharge. For sea level rise scenario, the results suggest a dilution in nutrient concentrations and an increase in Chl-a in specific areas.For both scenarios, the suggested increase in Chl-a concentration for specific estuarine areas, under the tested scenarios, can lead to events that promote an abnormal growth of phytoplankton (blooms) causing the water quality to drop and the estuary to face severe quality issues risking all the activities that depend on it.

Modelling the Minho river plume intrusion into the rias baixas

The Minho River, situated 30 km south of the Rias Baixas is the most important freshwater source flowing into the Western Galician Coast (NW of the Iberian Peninsula). This discharge is important to determine the hydrological patterns adjacent to its mouth, particularly close to the Galician coastal region. The buoyancy generated by the Minho plume can flood the Rias Baixas for long periods, reversing the normal estuarine density gradients. Thus, it becomes important to analyse its dynamics as well as the thermohaline patterns of the areas affected by the freshwater spreading. Thus, the main aim of this work was to study the propagation of the Minho estuarine plume to the Rias Baixas, establishing the conditions in which this plume affects the circulation and hydrographic features of these coastal systems, through the development and application of the numerical model MOHID. For this purpose, the hydrographic features of the Rias Baixas mouths were studied. It was observed that at the northern mouths, due to their shallowness, the heat fluxes between the atmosphere and ocean are the major forcing, influencing the water temperature, while at the southern mouths the influence of the upwelling events and the Minho River discharge were more frequent. The salinity increases from south to north, revealing that the observed low values may be caused by the Minho River freshwater discharge. An assessment of wind data along the Galician coast was carried out, in order to evaluate the applicability of the study to the dispersal of the Minho estuarine plume. Firstly, a comparative analysis between winds obtained from land meteorological stations and offshore QuikSCAT satellite were performed. This comparison revealed that satellite data constitute a good approach to study wind induced coastal phenomena. However, since the numerical model MOHID requires wind data with high spatial and temporal resolution close to the coast, results of the forecasted model WRF were added to the previous study. The analyses revealed that the WRF model data is a consistent tool to obtain representative wind data near the coast, showing good results when comparing with in situ wind observations from oceanographic buoys. To study the influence of the Minho buoyant discharge influence on the Rias Baixas, a set of three one-way nested models was developed and implemented, using the numerical model MOHID. The first model domain is a barotropic model and includes the whole Iberian Peninsula coast. The second and third domains are baroclinic models, where the second domain is a coarse representation of the Rias Baixas and adjacent coastal area, while the third includes the same area with a higher resolution. A bi-dimensional model was also implemented in the Minho estuary, in order to quantify the flow (and its properties) that the estuary injects into the ocean. The chosen period for the Minho estuarine plume propagation validation was the spring of 1998, since a high Minho River discharge was reported, as well as favourable wind patterns to advect the estuarine plume towards the Rias Baixas, and there was field data available to compare with the model predictions. The obtained results show that the adopted nesting methodology was successful implemented. Model predictions reproduce accurately the hydrodynamics and thermohaline patterns on the Minho estuary and Rias Baixas. The importance of the Minho river discharge and the wind forcing in the event of May 1998 was also studied. The model results showed that a continuous moderate Minho River discharge combined with southerly winds is enough to reverse the Rias Baixas circulation pattern, reducing the importance of the occurrence of specific events of high runoff values. The conditions in which the estuarine plume Minho affects circulation and hydrography of the Rias Baixas were evaluated. The numerical results revealed that the Minho estuarine plume responds rapidly to wind variations and is also influenced by the bathymetry and morphology of the coastline. Without wind forcing, the plume expands offshore, creating a bulge in front of the river mouth. When the wind blows southwards, the main feature is the offshore extension of the plume. Otherwise, northward wind spreads the river plume towards the Rias Baixas. The plume is confined close to the coast, reaching the Rias Baixas after 1.5 days. However, for Minho River discharges higher than 800 m3 s-1, the Minho estuarine plume reverses the circulation patterns in the Rias Baixas. It was also observed that the wind stress and Minho River discharge are the most important factors influencing the size and shape of the Minho estuarine plume. Under the same conditions, the water exchange between Rias Baixas was analysed following the trajectories particles released close to the Minho River mouth. Over 5 days, under Minho River discharges higher than 2100 m3 s-1 combined with southerly winds of 6 m s-1, an intense water exchange between Rias was observed. However, only 20% of the particles found in Ria de Pontevedra come directly from the Minho River. In summary, the model application developed in this study contributed to the characterization and understanding of the influence of the Minho River on the Rias Baixas circulation and hydrography, highlighting that this methodology can be replicated to other coastal systems.