Suites of Best Management Practices (BMPs) are usually selected to be economically and environmentally efficient in reducing nonpoint source (NPS) pollutants from agricultural areas in a watershed. The objective of this research was to compare the selection and placement of BMPs in a pasture-dominated watershed using multiobjective optimization and targeting methods. Two objective functions were used in the optimization process, which minimize pollutant losses and the BMP placement areas. The optimization tool was an integration of a multi-objective genetic algorithm (GA) and a watershed model (Soil and Water Assessment Tool—SWAT). For the targeting method, an optimum BMP option was implemented in critical areas in the watershed that contribute the greatest pollutant losses. A total of 171 BMP combinations, which consist of grazing management, vegetated filter strips (VFS), and poultry litter applications were considered. The results showed that the optimization is less effective when vegetated filter strips (VFS) are not considered, and it requires much longer computation times than the targeting method to search for optimum BMPs. Although the targeting method is effective in selecting and placing an optimum BMP, larger areas are needed for BMP implementation to achieve the same pollutant reductions as the optimization method.

In Taiwan, the steep landscape and highly vulnerable geology make it difficult to predict soil erosion and sediment transportation via variable transport conditions. In this study, we integrated the Taiwan universal soil loss equation (TUSLE) and slope stability conditions in the soil and water assessment tool (SWAT) as the SWAT-Twn model to improve sediment simulation and assess the sediment transport functions in the Chenyulan watershed, a small mountainous catchment. The results showed that the simulation of streamflow was satisfactory for calibration and validation. Before model calibration and validation for sediment, SWAT-Twn with default sediment transport method performed better in sediment simulation than the official SWAT model (version 664). The SWAT-Twn model coupled with the simplified Bagnold equation could estimate sediment export more accurately and significantly reduce the overestimated sediment yield by 65.7%, especially in highly steep areas. Furthermore, five different sediment transport methods (simplified Bagnold equation with/without routing by particle size, Kodoatie equation, Molinas and Wu equation, and Yang sand and gravel equation) were evaluated. It is suggested that modelers who conduct sediment studies in the mountainous watersheds with extreme rainfall conditions should adjust the modified universal soil loss equation (MUSLE) factors and carefully evaluate the sediment transportation equations in SWAT.

In this study, a deconvolution procedure was used to create a variogram of oral cancer (OC) rates. Based on the variogram, area-to-point (ATP) Poisson kriging and p-field simulation were used to downscale and simulate, respectively, the OC rate data for Taiwan from the district scale to a 1 km × 1 km grid scale. Local cluster analysis (LCA) of OC mortality rates was then performed to identify OC mortality rate hot spots based on the downscaled and the p-field-simulated OC mortality maps. The relationship between OC mortality and land use was studied by overlapping the maps of the downscaled OC mortality, the LCA results, and the land uses. One thousand simulations were performed to quantify local and spatial uncertainties in the LCA to identify OC mortality hot spots. The scatter plots and Spearman’s rank correlation yielded the relationship between OC mortality and concentrations of the seven metals in the 1 km cell grid. The correlation analysis results for the 1 km scale revealed a weak correlation between OC mortality rate and concentrations of the seven studied heavy metals in soil. Accordingly, the heavy metal concentrations in soil are not major determinants of OC mortality rates at the 1 km scale at which soils were sampled. The LCA statistical results for local indicator of spatial association (LISA) revealed that the sites with high probability of high-high (high value surrounded by high values) OC mortality at the 1 km grid scale were clustered in southern, eastern, and mid-western Taiwan. The number of such sites was also significantly higher on agricultural land and in urban regions than on land with other uses. The proposed approach can be used to downscale and evaluate uncertainty in mortality data from a coarse scale to a fine scale at which useful additional information can be obtained for assessing and managing land use and risk.

The Chenyulan watershed, located in the central mountain area of Taiwan, has been suffering from earthquakes, typhoons, and heavy rainfalls in recent decades. These sequential natural disturbances have a cumulative impact on the watershed, leading to more fragile and fragmented land cover and loss of capacity of soil water conservation. In this study, the Soil and Water Assessment Tool (SWAT) and a landscape metrics tool (FRAGSTATS) were used to assess the direct impact (e.g., by annual rainfall) and indirect impact (e.g., by landscape configuration and composition) of natural disturbances on the ecohydrological processes of the Chenyulan watershed. Six SPOT satellite images from 2008 to 2013 were analyzed by using the nearest feature line embedding (NFLE) approach and reclassified into six land cover types: forest, cultivated land, grassland, river, landslide, and built-up. Forest was found to have the largest patch size, indicating that it is more resilient to disturbances, while agricultural land tended to expand from the river side toward the hill. Two land cover change scenarios were compared in the SWAT model. The results showed that there was no significant difference in simulated streamflow during 2004–2015 and sediment loading during 2004–2009; however, the model performed better for sediment loading during 2010–2015 with dynamic land cover change (coefficient of determination (R²) = 0.66, Nash-Sutcliffe efficiency coefficient (NSE) = 0.62, percent bias (PBIAS) = 10.5%, root mean square error observation standard deviation ratio (RSR) = 0.62) than with constant land cover (R² = 0.61, NSE = 0.54, PBIAS = −17.3%, RSR = 0.68), indicating that long-term land cover change should be considered in hydrologic modeling. Changes in landslides during 2008–2013 were found to significantly affect ecohydrological processes, especially after 2011. In general, annual precipitation plays a dominant role, and landscape composition had by far the strongest influence on water yield and sediment yield compared to landscape configuration. The results can be useful for understanding the effects of land cover change on ecohydrological processes in the Chenyulan watershed and the potential impact of ecohydrological changes on the environment and public health.

Soil erosion and landslide triggered by heavy rainfall are serious problems that have threatened water resources in Taiwan watersheds. This study investigated the relationship among streamflow, sediment load, sediment concentration and typhoon characteristics (path and rainfall amount) during 2000–2017 for nine gauging stations in five basins (Tamshui River basin, Zhuoshui River basin, Zengwen River basin, Gaoping River basin, and Hualien River basin) representing the diverse geomorphologic conditions in Taiwan. The results showed that streamflow and sediment load were positively correlated, and the correlation was improved when the sediment load data were grouped by sediment concentration. Among these basins, the Zhuoshui River basin has the highest unit-discharge sediment load and unit-area sediment load. The soil in the upstream was more erodible than the downstream soil during the normal discharge conditions, indicating its unique geological characteristics and how typhoons magnified sediment export. The spatiotemporal variation in sediment loads from different watersheds was further categorized by typhoons of different paths. Although typhoon path types matter, the Zhuoshui and Hualien River basin were usually impacted by typhoons of any path type. The results indicated that sediment concentration, the watershed soil characteristics, and typhoons paths were the key factors for sediment loads. This study can be useful for developing strategies of soil and water conservation implementation for sustainable watershed management.

Sustainable conservation aims to ensure the sustained conservation of landscape multi-functionality which in turn requires ensuring ecosystem service (ES) and habitat quality (HQ) sustainability with inclusive landscape-scale conservation planning. This study proposes a landscape conservation planning (LCP) framework for landscape-scale ES-HQ conservation and sustainability. Spatially explicit hotspots for five ESs and HQs are identified via InVEST and LISA software. Spatiotemporal changes in ES-HQ hotspots, in terms of stability and resilience, are delineated. The Zonation technique is applied to prioritize areas for conservation based on ES-HQ hotspot stability and resilience maps. High priority conservation areas are identified and are used as reserve area inputs for land use modeling with CLUE-S software to simulate future land use change under climate change scenarios. This study reports that varied rainfall and climate are major driving factors of ES-HQ sustainability disturbance in the study area. Furthermore, our proposed conservation Strategy 2 demonstrates that a larger extent of landscape multi-functionality can be sustained when the existing conservation area includes the total area of identified ES-HQ resilient hotspots. This study effectively identifies the stability and resiliency of ES-HQ hotspot areas affected by disturbances for high priority landscape conservation requirements to ensure ES-HQ sustainability and landscape multi-functionality in the study area.