Yung-Chieh Wang

Remote Sensing Data with the Conditional Latin Hypercube Sampling and Geostatistical Approach to Delineate Landscape Changes Induced by Large Chronological Physical Disturbances

This study applies variogram analyses of normalized difference vegetation index (NDVI) images derived from SPOT HRV images obtained before and after the ChiChi earthquake in the Chenyulan watershed, Taiwan, as well as images after four large typhoons, to delineate the spatial patterns, spatial structures and spatial variability of landscapes caused by these large disturbances. The conditional Latin hypercube sampling approach was applied to select samples from multiple NDVI images. Kriging and sequential Gaussian simulation with sufficient samples were then used to generate maps of NDVI images. The variography of NDVI image results demonstrate that spatial patterns of disturbed landscapes were successfully delineated by variogram analysis in study areas. The high-magnitude Chi-Chi earthquake created spatial landscape variations in the study area. After the earthquake, the cumulative impacts of typhoons on landscape patterns depended on the magnitudes and paths of typhoons, but were not always evident in the spatiotemporal variability of landscapes in the study area. The statistics and spatial structures of multiple NDVI images were captured by 3,000 samples from 62,500 grids in the NDVI images. Kriging and sequential Gaussian simulation with the 3,000 samples effectively reproduced spatial patterns of NDVI images. However, the proposed approach, which integrates the conditional Latin hypercube sampling approach, variogram, kriging and sequential Gaussian simulation in remotely sensed images, efficiently monitors, samples and maps the effects of large chronological disturbances on spatial characteristics of landscape changes including spatial variability and heterogeneity.

A System Dynamic Model and Sensitivity Analysis for Simulating Domestic Pollution Removal in a Free-Water Surface Constructed Wetland

This work develops a system dynamic simulation model for free-water surface constructed wetlands, as well as provides appropriate values for the parameters of constructed wetland management. The system dynamic model is calibrated and validated by using data from a 1-year study of a constructed wetland in Tainan of southern Taiwan. Additionally, the major parameters that affect the simulation output are obtained via sensitivity analysis by using generalized likelihood uncertainty estimation (GLUE). A high R2 and Nash–Sutcliffe coefficient of efficiency between the simulated and measured outflow values indicate that in addition to reproducing the changing trends of dissolved oxygen (DO), 5-day biological oxygen demand (BOD5), total nitrogen (TN), total suspended sediment (TSS), and total phosphorous (TP) concentrations, the model can simulate the variations of DO, BOD5, and TSS. Taken into account the interactions among parameters, the GLUE method successfully obtained the model sensitive parameters from the Monte Carlo parameter sets. Sensitivity analysis results indicate that the parameters of microorganisms are sensitive factors that affect DO, BOD5, and TN, while sediment diameter largely influences TP and TSS. Further elucidating environmental microorganisms would increase the model accuracy and provide a valuable reference for constructed wetland management and design.

A decision-making approach for delineating sites which are potentially contaminated by heavy metals via joint simulation.

This work develops a new approach for delineating sites that are contaminated by multiple soil heavy metals and applies it to a case study. First a number of contaminant sample data are transformed into multiple spatially un-correlated factors using Uniformly Weighted Exhaustive Diagonalization with Gauss iterations (U-WEDGE). Sequential Gaussian simulation (sGs) is then used to generate sets of realizations of each resultant factor. These are then transformed into sets of sGs contaminant distribution realizations, which are then used to analyze the local and spatial (global) uncertainties in the distribution and concentration of contaminants via joint simulation. Finally, Info-Gap Decision Theory (IGDT) is used to consider different monitoring and or remediation regimes based on the analysis of contaminant realization spatial uncertainty. In our case study each heavy metal contaminant was considered individually and together with all other heavy metals; as the number of heavy metals considered increased, higher critical proportion values of local probability were chosen to obtain a low global uncertainty (to provide high reliability). Info-Gap Decision Theory (IGDT) yielded the most appropriate critical proportion values which minimized information loss in terms of specific goals. When the false negative rate is set to zero, meaning that it is necessary to monitor all potentially polluted areas, the corresponding false positive rates are at least 63%, 65%, 66%, 68%, 70%, and 78% to yield robustness levels of 0.50, 0.60, 0.70, 0.80, 0.90, and 1.00 respectively. However, when the false negative rate tolerance threshold is raised to 50%, the false positive rate tolerance which yields robustness levels of 0.50, 0.60, 0.70, 0.80, 0.90 and 1.00 drop to 12%, 14%, 15%, 18%, 20%, and 39%. The case study demonstrates the effectiveness of the developed approach at making robust decisions concerning the delineation of sites contaminated by multiple heavy metals.

Assessing and mapping spatial associations among oral cancer mortality rates, concentrations of heavy metals in soil, and land use types based on multiple scale data.

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.

Detecting the Land-Cover Changes Induced by Large-Physical Disturbances Using Landscape Metrics, Spatial Sampling, Simulation and Spatial Analysis

The objectives of the study are to integrate the conditional Latin Hypercube Sampling (cLHS), sequential Gaussian simulation (SGS) and spatial analysis in remotely sensed images, to monitor the effects of large chronological disturbances on spatial characteristics of landscape changes including spatial heterogeneity and variability. The multiple NDVI images demonstrate that spatial patterns of disturbed landscapes were successfully delineated by spatial analysis such as variogram, Moran’I and landscape metrics in the study area. The hybrid method delineates the spatial patterns and spatial variability of landscapes caused by these large disturbances. The cLHS approach is applied to select samples from Normalized Difference Vegetation Index (NDVI) images from SPOT HRV images in the Chenyulan watershed of Taiwan, and then SGS with sufficient samples is used to generate maps of NDVI images. In final, the NDVI simulated maps are verified using indexes such as the correlation coefficient and mean absolute error (MAE). Therefore, the statistics and spatial structures of multiple NDVI images present a very robust behavior, which advocates the use of the index for the quantification of the landscape spatial patterns and land cover change. In addition, the results transferred by Open Geospatial techniques can be accessed from web-based and end-user applications of the watershed management.

Using CV-GLUE procedure in analysis of wetland model predictive uncertainty

This study develops a procedure that is related to Generalized Likelihood Uncertainty Estimation (GLUE), called the CV-GLUE procedure, for assessing the predictive uncertainty that is associated with different model structures with varying degrees of complexity. The proposed procedure comprises model calibration, validation, and predictive uncertainty estimation in terms of a characteristic coefficient of variation (characteristic CV). The procedure first performed two-stage Monte-Carlo simulations to ensure predictive accuracy by obtaining behavior parameter sets, and then the estimation of CV-values of the model outcomes, which represent the predictive uncertainties for a model structure of interest with its associated behavior parameter sets. Three commonly used wetland models (the first-order K-C model, the plug flow with dispersion model, and the Wetland Water Quality Model; WWQM) were compared based on data that were collected from a free water surface constructed wetland with paddy cultivation in Taipei, Taiwan. The results show that the first-order K-C model, which is simpler than the other two models, has greater predictive uncertainty. This finding shows that predictive uncertainty does not necessarily increase with the complexity of the model structure because in this case, the more simplistic representation (first-order K-C model) of reality results in a higher uncertainty in the prediction made by the model. The CV-GLUE procedure is suggested to be a useful tool not only for designing constructed wetlands but also for other aspects of environmental management.

Conservation planning to zone protected areas under optimal landscape management for bird conservation

This study proposes a two-stage conservation planning approach. Firstly, the Land-Use Pattern Optimization-library is used to maximize the suitability of habitats for target species by optimizing configuration based on the current landscape. Secondly, the systematic conservation planning tool, Marxan is used to identify protected areas based on the estimated species distributions from the optimal landscape configuration. We compared our conservation plan for three target bird species from a highland farm with the conservation plan produced using Marxan alone. Our comparison showed the effectiveness of our approach by selecting a reserve network with higher habitat suitability, better connection, and smaller size after relatively minor landscape modification. The proposed approach advances previous reserve site selection algorithms by considering optimal landscape configuration and potential species distributions for a reserve network design. Our approach yields priority maps to guide the design of a reserve network as well as identify landscape management for conservation.

Estimating and classifying spatial and temporal distributions of flow conditions for fish habitats by using geostatistical approaches with measured flow and fish data

This study investigated the relationship between the distribution of Sicyopterus japonicus, current velocity and water depth in four reaches along Datuan Stream in northern Taiwan during winter of 2007 and spring of 2008. The spatial distributions of current velocity and water depth were estimated by kriging from the stream mouth to the upstream section. The empirical rule method and the Froude number method with kriging estimated distributions of the current velocity and water depth were incorporated into a geographical information system (GIS) and used to classify flow conditions at the investigated reaches. Indicator kriging was used to estimate the probability of the presence of S. japonicus and superimposed on the estimated flow conditions at each reach. The field results showed that, in each investigated season, the average current velocity was low in the downstream and upstream reaches, but high in the middle stream reach. The flow conditions based on kriging estimated distribution of water velocity and water depth at the investigated reaches accurately reflect flow conditions at each reach along the stream. Geostatistical approaches, such as kriging and indicator kriging, can be used to estimate the flows and the appearance probability of fish accurately. Moreover, the overlapping maps of the flow conditions and the probabilities indicate that the preferences ofS. japonicus vary in different reaches and seasons. Based on the migration behavior of S. japonicus, the classification of the empirical rule flow classification method with geostatistical approaches in GIS can be used to estimate the preferences of S. japonicus in Datuan stream effectively.