Amplified eruptive outbreaks of bark beetles as a consequence of climate change can cause tree mortality that significantly affects terrestrial water and carbon fluxes. However, the lack of field-scale observations of underlying physiological mechanisms currently hampers the expression of such ecosystem disturbances in predictive modelling. Based on a unique flux tower dataset from a subalpine forest located in the Rocky Mountains, mechanisms of stomatal response to an extensive bark beetle outbreak were investigated using various models and parametrizations. The datasets cover a decade, including the periods of pre-infestation, infestation, and post-infestation. Field measurements showed considerable decreases in evapotranspiration (ET), transpiration (T), and leaf area index (LAI) during the two-year infestation period compared to the pre-infestation period. Model interpretations of observed water and carbon fluxes indicated that the overall reductions in T were not solely due to decreased LAI, but also to changes in physiological behaviours. The summer season's canopy-scale stomatal conductance was significantly reduced during the infestation period, from 0.0018 to 0.0011 m s⁻¹. One primary reason for the observed variations is likely that the bark beetle infestation hampers the water transport in the xylem. The damage of xylem has important implications for water use efficiency (WUE), which also significantly influences the parameterization of stomatal conductance. When using stomatal conductance models to forecast ecosystem dynamics, it is crucial to recalibrate the model's parameters to ensure the accurate depiction of stomatal dynamics during various infestation periods. The neglect of the temporal variability of canopy-scale stomatal conductance under ecosystem disturbances (e.g., bark beetle infestations) in current earth system models, therefore, requires specific attention in assessments of large-scale water and carbon balances.

The process of evaporation interacts with the soil, which has various comprehensive mechanisms. Multiphase flow models solve air, vapour, water, and heat transport equations to simulate non-isothermal soil moisture transport of both liquid water and vapor flow, but are only applied in non-vegetated soils. For (sparsely) vegetated soils often energy balance models are used, however these lack the detailed information on non-isothermal soil moisture transport. In this study we coupled a multiphase flow model with a two-layer energy balance model to study the impact of non-isothermal soil moisture transport on evaporation fluxes (i.e., interception, transpiration, and soil evaporation) for vegetated soils. The proposed model was implemented at an experimental agricultural site in Florida, US, covering an entire maize-growing season (67 days). As the crops grew, transpiration and interception became gradually dominated, while the fraction of soil evaporation dropped from 100% to less than 20%. The mechanisms of soil evaporation vary depending on the soil moisture content. After precipitation the soil moisture content increased, exfiltration of the liquid water flow could transport sufficient water to sustain evaporation from soil, and the soil vapor transport was not significant. However, after a sufficient dry-down period, the soil moisture content significantly reduced, and the soil vapour flow significantly contributed to the upward moisture transport in topmost soil. A sensitivity analysis found that the simulations of moisture content and temperature at the soil surface varied substantially when including the advective (i.e., advection and mechanical dispersion) vapour transport in simulation, including the mechanism of advective vapour transport decreased soil evaporation rate under wet condition, while vice versa under dry condition. The results showed that the formulation of advective soil vapor transport in a soil-vegetation-atmosphere transfer continuum can affect the simulated evaporation fluxes, especially under dry condition.

Vegetation can affect slope hydrology and stability via plant transpiration and induced matric suction. Previous work suggested that the presence of plant roots would induce preferential flow, and its effects may be more significant when the planting density is high. However, there is a lack of numerical studies on how planting density affects soil pore-water pressure and shear strength during heavy rainfall. This study aims to investigate the impact of plant root-induced preferential flow on hydromechanical processes of vegetated soils under different planting densities. Two modelling approaches, namely single- and dual-permeability models, were integrated with an infinite slope stability approach to simulate pore-water pressure dynamics and slope stability. Laboratory tests on soils with two different planting densities for a plant species, Schefflera heptaphylla, were conducted for numerical simulations. The single-permeability model overestimated the pore-water pressure in shallow soil and underestimated the infiltration depth. The dual-permeability model, which is able to model the effects of preferential flow, can better capture the observations of rapid increase of pore-water pressure and deeper pressure response in the vegetated soil. However, caution should be taken on the choice of pore-water pressure when using the dual-permeability model to assess the factor of safety. The dual-permeability model using the pore-water pressure in the preferential flow domain and that in the matrix domain would result in a lower and higher factor of safety, respectively.