The second-largest wildfire in the history of South Korea occurred in 2022 due to strong winds and dry climates. Quantitative evaluation of soil erosion is necessary to prevent subsequent sediment disasters in the wildfire areas. The erosion rates in two watersheds affected by the wildfires were assessed using the revised universal soil loss equation (RUSLE), a globally popular model, and the soil erosion model for mountain areas (SEMMA) developed in South Korea. The GIS-based models required the integration of maps of the erosivity factor, erodibility factor, length and slope factors, and cover and practice factors. The rainfall erosivity factor considering the 50-year and 80-year probability of rainfall increased from coastal to mountainous areas. For the LS factors, the traditional version (TV) was initially used, and the flow accumulation version (FAV) was additionally considered. The cover factor of the RUSLE and the vegetation index of the SEMMA were calculated using the normalized difference vegetation index (NDVI) extracted from Sentinel-2 images acquired before and after the wildfire. After one year following the wildfire, the NDVI increased compared to during the year of the wildfire. Although the RUSLE considered a low value of the P factor (0.28) for post-fire watersheds, it overestimated the erosion rate by from 3 to 15 times compared to the SEMMA. The erosion risk with the SEMMA simulation decreased with the elapsed time via the vegetation recovery and stabilization of topsoil. While the FAV of RUSLE oversimulated by 1.65~2.31 times compared to the TV, the FAV of SEMMA only increased by 1.03~1.19 times compared to the TV. The heavy rainfall of the 50-year probability due to Typhoon Khanun in 2023 generated rill and gully erosions, landslides, and sediment damage in the post-fire watershed on forest roads for transmission tower construction or logging. Both the RUSLE and SEMMA for the TV and FAV predicted high erosion risks for disturbed hillslopes; however, their accuracy varied in terms of the intensity and extent. According to a comparative analysis of the simulation results of the two models and the actual erosion situations caused by heavy rain, the FAV of SEMMA was found to simulate spatial heterogeneity and a reasonable erosion rate.

A large wildfire occurred due to strong winds and dry climates in the Gangwon province of South Korea. Thereafter, floods and sediment damage were caused by Typhoon Mitag in the burned areas. This study was an attempt to quantitatively evaluate the risk of soil erosion in wildfire areas using the Soil Erosion Model for Mountain Areas (SEMMA) based on GIS, which was developed in South Korea. The model required the integration of maps of the main factors involved, i.e., rainfall erosivity, vegetation index, soil erodibility, and slope length and steepness. According to the model simulation results, high erosion rates of over 100 t/ha were concentrated within the wildfire areas. Sediment yields from the study watershed, including the wildfire areas, were estimated to be 40.33 t/ha for the 30-year frequency of rainfall, which is similar to those of the typhoon. The high risk of erosion was predominantly observed in the upper mountains, which are characterized by steep slopes, silt loam, and shallow soil depths within the wildfire areas. Urgent and excessive logging of burned trees further increased the risk of erosion. However, various treatment strategies were implemented to control soil erosion and sediment transport from the post-fire watershed. This study confirmed that temporal and spatial BMPs should be selected and enforced to reduce sediment disasters in wildfire areas.

Hydraulic characteristics of rill incipience on steep slopes with sandy soils are different from those of gentle slopes in the mountainous watersheds. In other to better understand the processes of the incision and expansion of rill, rainfall simulations were conducted in a laboratory soil box with an inflow device installed at the top of the steep slope. The location and shape of the rill incision were tested with respect to different slopes (20, 24, and 28°), rainfall intensity (80 and 100 mm/h), and segment distances with various inflow rates (0.0, 8.8, 17.5, and 26.3 ml/s). The result showed that rills with heatcut were mostly incised on the 2.4 m segments and progressively moved upslope on the 3.2 m segments. The steps of rill evolution with increase of inflow are identified as four stages; (1) sheet flow with impact of raindrops, (2) overland flow with flowpaths, (3) microrill with headcut incision, and (4) rill enhanced by headcut expansion. It appears that incision and headcut of rill is prevalent when the Froude number (Fr) changes from subcritical flow to supercritical flow and the maximum velocity of rill reaches over 0.065 m/s. Unit stream power using the maximum velocity of rill was the best parameter for estimating rill erosion as rill evolution greatly depends on maximum velocity of the local flow enhanced by rill incision and expansion As a result, it was found that critical conditions for rill incipient of sandy soils on steep slope are observed greater than those of loamy soils or gentle slopes.

Rill erosion is an important source of sediment on steep disturbed hillslopes, although the interaction among runoff accumulation, rill micromorphological evolution, and sediment response remains insufficiently understood. This study investigated rill development on sandy slopes of 15° and 20° under controlled rainfall-only and rainfall-inflow conditions using direct measurements and LiDAR-derived three-dimensional topography. Under rainfall-only conditions, rills initiated approximately 3 min after the onset of rainfalland propagated upslope, whereas additional inflow induced rapid incision within approximately 30 s, thereby promoting downslope extension and network connectivity. Additional inflow substantially enhanced rill development across both slopes; however, the dominant mode varied, with the 15° slopes exhibiting a pronounced 2.9-fold increase in network length and density, while lateral enlargement was enhanced on the 20° slopes, resulting in a 67% increase in mean rill width. Manually estimated and LiDAR-derived rill volumes were strongly associated. The ratios of measured rill volume to the volume-equivalent of exported sediment were 0.62 and 0.73 at 15° and 20°, respectively. Sediment yield was strongly correlated with mean rill flow velocity, as well asgeometrically estimated rill volume, and LiDAR-derived rill volume, while runoff discharge and individual rill dimensions showed weaker relationships. These findings indicate that sediment yield was more closely associated with concentrated-flow hydraulics and integrated rill development than with runoff magnitude alone. These 3D morphological changes imply that rill volume can serve as a key structural indicator of sediment response, capturing the cumulative effects of concentrated flow on steep disturbed slopes.