HYDROLOGIC MODELING OF LAND PROCESSES IN PUERTO RICO USING REMOTELY SENSED DATA1

ABSTRACT: An integrated, multi-disciplinary effort to model land processes affecting Mayaguez Bay in western Puerto Rico is described. A modeling strategy was developed to take advantage of remotely sensed data. The spatial, interannual, and seasonal variability of sediment discharges to the bay were also evaluated. Classified images of remotely sensed data revealed the spatial distribution and quantities of land use classes in the region and aided in the discretization of the watershed into homogeneous regions. These regions were modeled using a geomorphic modeling technique based upon spatially averaged parameters. Simulation results from the modeling effort compared favorably with observations at two locations within the watershed. Results showed that runoff and sediment loads from the area exhibit a marked seasonal trend and that deforested areas located in the foothill regions of the watershed contribute a disproportionate share of the sediment load to the bay. In years when rainfall distributions are uniformly distributed over the area, the sediment yields may be up to 100 percent higher than years when the rainfall is concentrated in the heavily forested mountainous regions.     

HYDROLOGIC MODELING WITH REMOTELY SENSED DATABASES1

ABSTRACT: An integrated remotely sensed database was used as the basis for a hydrologic and sediment transport modeling effort for an agricultural area of western Puerto Rico. Classified spectral images of airborne radiance data provided ground cover information and were used in conjunction with topographic and soils data to guide model construction and provide input to the water balance and sediment yield simulations. Runoff and sediment discharge from hydrologically homogeneous regions were routed through the drainage network and combined at the basin outlet. The model was used to simulate four years of observed sediment discharge from the basin. Relative contributions to the total sediment yield of forested and agricultural areas were determined and compared.     

Selenium distribution and fluxes in intertidal wetlands, San Francisco Bay, California

Selenium (Se) concentrations exceeding ecological guidelines for sediments and suspended particulate matter (SPM) have been observed in the northern reach of the San Francisco Bay estuary. Longterm availability of elevated Se in wetland sediments depends in part on the fluxes controlling Se distribution. The relative contribution of sedimentary vs. post-depositional Se fluxes in two San Francisco Bay intertidal wetlands was estimated. Selenium concentrations on surface wetland sediments were compared with levels on SPM, and with previously established background levels in San Francisco Bay sediments. Sediment Se fluxes to the wetlands were measured directly using sediment traps. Although dissolved Se concentrations are higher than particulate Se concentrations in San Francisco Bay water, sediment input into the system provides the major flux of Se. Strong correlation between Se and C on SPM (r2 = 0.81) indicates the importance of organic particulate deposition. Dependence on sediment texture was qualitatively established by measuring Se on particle-size separates. Normalization to Al showed that 65% of Se spatial variability is related to sediment texture. Selenium is further enriched in the marsh via post-depositional inputs, probably due to in situ adsorption from overlying water and chemical reduction. According to sediment flux measurements, enrichment in the marsh is equivalent to 20 to 25% of the particulate Se flux, thereby defining the marsh as a Se sink. These findings highlight the need for more intensive monitoring of SPM as the major source of Se to intertidal wetlands.     

ESTIMATING THE SUSPENDED SEDIMENT LOAD IN RESERVOIRS1

ABSTRACT: The total suspended sediment loads of four north Mississippi reservoirs were determined from measurements of concentrations of suspended sediment in a vertical profile at several locations on each reservoir made during the year. These data were combined with the stage-height and known stage-volume relationships for each reservoir in a numerical integration to determine the total suspended sediment in the water body. Total suspended sediments were estimated using the product of the suspended sediment concentration in the surface water by the appropriate reservoir volume. The averaged ratios of the estimated to measured suspended sediment loads for each reservoir exceeded 0.90. Since the concentration of suspended sediments in surface waters of north Mississippi reservoirs has been shown as highly correlated with spectral reflectance, estimating the total suspended sediment of these reservoirs using remotely sensed spectral reflectance data is possible.     

Remotely Sensed Data for Ecosystem Analyses: Combining Hierarchy Theory and Scene Models

Remotely sensed data have been used extensively for environmental monitoring and modeling at a number of spatial scales; however, a limited range of satellite imaging systems often constrained the scales of these analyses. A wider variety of data sets is now available, allowing image data to be selected to match the scale of environmental structure(s) or process(es) being examined. A framework is presented for use by environmental scientists and managers, enabling their spatial data collection needs to be linked to a suitable form of remotely sensed data. A six-step approach is used, combining image spatial analysis and scaling tools, within the context of hierarchy theory. The main steps involved are: (1) identification of information requirements for the monitoring or management problem; (2) development of ideal image dimensions (scene model), (3) exploratory analysis of existing remotely sensed data using scaling techniques, (4) selection and evaluation of suitable remotely sensed data based on the scene model, (5) selection of suitable spatial analytic techniques to meet information requirements, and (6) cost–benefit analysis. Results from a case study show that the framework provided an objective mechanism to identify relevant aspects of the monitoring problem and environmental characteristics for selecting remotely sensed data and analysis techniques.     

Measurements and simulation of forest leaf area index and net primary productivity in Northern China

Large scale process-based modeling is a useful approach to estimate distributions of global net primary productivity (NPP). In this paper, in order to validate an existing NPP model with observed data at site level, field experiments were conducted at three sites in northern China. One site is located in Qilian Mountain in Gansu Province, and the other two sites are in Changbaishan Natural Reserve and Dunhua County in Jilin Province. Detailed field experiments are discussed and field data are used to validate the simulated NPP. Remotely sensed images including Landsat Enhanced Thematic Mapper plus (ETM+, 30 m spatial resolution in visible and near infrared bands) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER, 15m spatial resolution in visible and near infrared bands) are used to derive maps of land cover, leaf area index, and biomass. Based on these maps, field measured data, soil texture and daily meteorological data, NPP of these sites are simulated for year 2001 with the boreal ecosystem productivity simulator (BEPS). The NPP in these sites ranges from 80 to 800 gCm(-2)a(-1). The observed NPP agrees well with the modeled NPP. This study suggests that BEPS can be used to estimate NPP in northern China if remotely sensed images of high spatial resolution are available.     

Analyzing the Variability of Remote Sensing and Hydrologic Model Evapotranspiration Products in a Watershed in Michigan

The goal of this paper was to statistically explore the spatiotemporal performance of remotely sensed actual evapotranspiration (ETa) datasets and a remotely sensed ensemble in a region that lacks observed data. The remotely sensed datasets were further compared with ETa results from a physically based hydrologic model (Soil and Water Assessment Tool) to examine the differences and determine the level of agreement between the ETa datasets and the model outputs. ETa datasets were compared on temporal (i.e., monthly and seasonal basis) and spatial (i.e., landuse) scales at both watershed and subbasin levels. The results showed a lack of consistent similarities and differences among the datasets when evaluating the monthly ETa variations; however, the seasonal aggregated data presented more consistent similarities and differences during the spring and summer compared to the fall and winter. Meanwhile, spatial analysis of the datasets showed the MOD16A2 500 m ETa product was the most versatile of the tested datasets, being able to differentiate between landuses during all seasons. Furthermore, the use of an averaging ensemble was able to improve overall ETa performance in the study area. This study showed that the remotely sensed ETa products are not similar throughout the year, but the appropriate application periods for different ETa products were identified. Finally, spatial variabilities of the ETa products are more in tune with landuse and climate characteristics.     

MULTIVARIATE ANALYSIS OF WATER QUALITY AND PHYSICAL CHARACTERISTICS OF SELECTED WATERSHEDS IN PUERTO RICO1

ABSTRACT: Multivariate analyses were used to develop equations that could predict certain water quality (WQ) conditions for unmonitored watersheds in Puerto Rico based on their physical characteristics. Long term WQ data were used to represent the WQ of 15 watersheds in Puerto Rico. A factor analysis (FA) was performed to reduce the number of chemical constituents. Cluster analysis (CA) was used to group watersheds with similar WQ characteristics. Finally, a discriminant analysis (DA) was performed to relate the WQ clusters to different physical parameters and generate predicting equations. The FA identified six factors (77 percent of variation explained): nutrients, dissolved ions, sodium and chloride, silicacious geology, red ox conditions, and discharge. From the FA, specific conductance, sodium, phosphorous, silica, and dissolved oxygen were selected to represent the WQ characteristics in the CA. The CA determined five groups of watersheds (forested, urban polluted, mixed urban/rural, forested plutonic, and limestone) with similar WQ properties. From the five WQ clusters, two categories can be observed: forested and urban watersheds. The DA found that changes in forest cover, percent of limestone, mean annual rainfall, and watershed shape factor were the most important physical features affecting the WQ of watersheds in Puerto Rico.     

Lake area changes in the middle Yangtze region of China over the 20th century

The Jianghan Plain and the Dongting lake area, located in the middle reaches of the Yangtze River are famous for their abundant freshwater lakes. The lakes have undergone changes in size and number over thousands of years due to natural causes and human activities. The 20th century particularly, witnessed dramatic changes in the freshwater resources of this region. This paper traces and analyzes lake evolution in the middle reaches of the Yangtze River over the past century. Lakes greater than 0.1 km(2) in size are mapped using Geographic Information System. Data is acquired and integrated from drainage network maps, relief maps, historical maps and remotely sensed images for different time periods. The results indicate that while there has been little change in the number of lakes over the past century, the lake area has experienced a dramatic decrease of 58.06%. The paper also examines the natural processes and human activities that may have contributed to the decrease in lake area. The results show that the decrease in total lake area appears to coincide with periods of rapid land reclamation in the middle reaches of the Yangtze River. Moreover, uncontrolled land reclamation activities can create an increase in sediment deposition in lakes thereby further reducing the lake size. Reduction of the lake area directly affects flood-control and has a negative ecological impact on the environment and on human life and property.     

A Model for Predicting Daily Peak Visitation and Implications for Recreation Management and Water Quality: Evidence from Two Rivers in Puerto Rico

Visitor use surveys and water quality data indicates that high visitor use levels of two rivers in Puerto Rico does not appear to adversely affect several water quality parameters. Optimum visitor use to maximize visitor defined satisfaction is a more constraining limit on visitor use than water quality. Our multiple regression analysis suggests that visitor use of about 150 visitors per day yields the highest level of visitor reported satisfaction, a level that does not appear to affect turbidity of the river. This high level of visitor use may be related to the gregarious nature of Puerto Ricans and their tolerance for crowding on this densely populated island. The daily peak visitation model indicates that regulating the number of parking spaces may be the most effective way to keep visitor use within the social carrying capacity.     

Road Erosion and Delivery Index (READI): A Model for Evaluating Unpaved Road Erosion and Stream Sediment Delivery

The Road Erosion and Delivery Index (READI) is a new geographic information system–based model to assess erosion and delivery of water and sediment from unpaved road networks to streams. READI quantifies the effectiveness of existing road surfacing and drain placements in reducing road sediment delivery and guides upgrades to optimize future reductions. Roads are draped on a digital elevation model and parsed into hydrologically distinct segments. Segments are further divided by engineered drainage structures. For each segment, a kinematic wave approximation generates runoff hydrographs for specified storms, with discharge directly to streams at road–stream crossings and onto overland‐flow plumes at other discharge points. Plumes are attenuated by soil infiltration, which limits their length, with delivery occurring if plumes intersect streams. Sediment production and sediment delivery can be calculated as a relative dimensionless index. READI predicts only a small proportion of new drains and new surfacing results in the majority of sediment delivery reductions. The model illustrates how the spatial relationships between road and stream networks, controlled by topography and network geometries, influence patterns of road–stream connectivity. READI was applied in seven northern California basins. The model was also applied in a recent burn area to examine how reduced hillslope infiltration can result in increased hydrologic connectivity and sediment delivery.     

PREDICTABILITY OF SURFACE WATER POLLUTION LOADING IN PENNSYLVANIA USING WATERSHED‐BASED LANDSCAPE MEASUREMENTS1

ABSTRACT: We formally evaluated the relationship between landscape characteristics and surface water quality in the state of Pennsylvania (USA) by regressing two different types of pollutant responses on landscape variables that were measured for whole watersheds. One response was the monthly exported mass of nitrogen estimated from field measurements, while the other response was a GIS‐modeled pollution potential index. Regression models were built by the stepwise selection protocol, choosing an optimal set of landscape predictors. After factoring out the effect of physiography, the dominant predictors were the proportion of “annual herbaceous” land and “total herbaceous” land for the nitrogen loading and pollution potential index, respectively. The strength of these single predictors is encouraging because the marginal land cover proportions are the simplest landscape measurements to obtain once a land cover map is in hand; however, the optimal set of predictors also included several measurements of spatial pattern. Thus, for watersheds at this general hierarchical scale, gross landscape pattern may be an important influence on instream pollution loading. Overall, there is strong evidence that using landscape measurements alone, obtained solely from remotely sensed data, can explain most of the water quality variability (R²= approx. 0.75) within these watersheds.     

Satellite-based indices in the analysis of land cover for municipalities in the province of Siena, Italy

Territorial indicators based on the satellite measured reflected or emitted energy can provide valuable information on the spatial evolution of a territory. Information regarding land use, biomass coverage and radiant temperature were obtained from remotely obtained measurements for each municipality in the Province of Siena. Remotely sensed data were combined into different indices to compare characteristics of land cover between territories and create an information base for continued study. A classification was performed on a municipal level and land use classes were grouped together and a general index of use pressure was created. A vegetation index was used to compare biomass densities. A radiant temperature index was calculated using measure thermal infrared emissions. The results of these analyses allowed researchers to examine the intra and inter-municipal spatial heterogeneity within the Province. A comparison between the results of each index was made and the problems related to the use of such indices were examined.     

CHANGING NEAR‐STREAM LAND USE AND RWER CHANNEL MORPHOLOGY IN THE VENEZUELAN ANDES1

ABSTRACT: The shape of a river channel is linked to surrounding land use through interacting hydrologic and geologic processes. This study analyzes the relationship between the change in near‐stream land use and the shape of the adjacent river channel over time. Three watersheds in the foothills of the Venezuelan Andes that have experienced differing degrees of development were studied to determine river channel width, sinuosity, and position relative to surrounding land use. Change in land use over time was obtained from multiple‐date aerial photographs (1946 and 1980) referenced to 1996 Landsat Thematic Mapper (TM) satellite imagery, and verified by field inspection. Measurements of land‐use type and amount and river channel morphology from the two dates were made using geographic information system (GIS) methods. The three watersheds differed in the extent of deforestation, the location of remaining forested land, and how much land‐use change had already occurred by 1946. Change in river channel morphology was greatest at the most deforested sites. Valley shape and channel constraint also had a discernible effect on change in channel morphology. This study introduces a method for analyzing change in coupled terrestrial‐aquatic systems based on multiple‐date, remotely sensed data and GIS analysis of spatial properties. The results document human impacts on river channels through a comparison of multiple watersheds over a 35‐year time interval.     

Recreation benefits of natural area characteristics at the El Yunque National Forest

The paper applies the Contingent Valuation Method using a visitor survey to determine trip values to three rivers in northeastern Puerto Rico. Sample data were obtained from 984 interviews conducted in 13 sites during the summer months. The trip value is statistically decomposed into components representing natural features and recreation activities. The median net economic value of a river visit without waterfalls or foot trails is $96. Recreation at rivers with both waterfalls and foot trails has a median value of $138 per trip. The presence of waterfalls accounts for $23 and foot trails $19 per trip. This information can be pertinent for managers in selecting recreation site locations and protecting them from competing uses.     

Manatee Mortality in Puerto Rico

Trichechus manatus ) in Puerto Rico is mortality due to human activities. We assessed 90 cases of manatee strandings in Puerto Rico based on historical data and a coordinated carcass salvage effort from 1990 through 1995. We determined patterns of mortality, including type of event, condition of carcasses, spatial and temporal distribution, gender, size/age class, and the cause of death. The spatial distribution of stranding events was not uniform, with the north, northeast, and south coasts having the highest numbers. Six clusters representing the highest incidence included the areas of Fajardo and Ceiba, Bahía de Jobos, Toa Baja, Guayanilla, Cabo Rojo, and Rio Grande to Luquillo. The number of reported cases has increased at an average rate of 9.6%/yr since 1990. The seasonality of stranding events showed a bimodal pattern, from February through April and in August and September. Most identified causes of death were due to human interaction, especially captures and watercraft collisions. Natural causes usually involved dependent calves. From 1990 through 1995, most deaths were attributed to watercraft collisions. A reduction in anthropogenic mortality of this endangered species can be accomplished only through education and a proactive management and conservation plan that includes law enforcement, mortality assessment, scientific research, rescue and rehabilitation, and inter- and intraagency cooperation.     

Environmental mapping based on spatial variability

Environmental maps show the probable environmental states of different types of land use or development of landscape in a geographic context. Remotely sensed data are particularly efficient for environmental mapping in order to outline major environmental types. Multiple schemes of image classification used in environmental mapping are either traditionally statistical or heuristic. While the former methods do not take account of spatial variability in space and aerial data, the latter ones does not lend themselves to optimal solutions we present. Novel probabilistic models of piecewise-homogeneous images are used in environmental mapping to segment real images. The models consider both an image and a land cover map. Such a pair constitutes an example of a Markov random field specified by a joint Gibbs probability distribution of images and maps. Parameters of the model are estimated by using a stochastic approximation technique. Its convergence to the desired values is studied experimentally. Addition of spatial attributes appears to be necessary in most areas where the differences in spatial data between regions in the image occur. Experiments in generating the pairs of images and environmental maps and in segmenting the simulated as well as real images are discussed.     

RESERVOIR SEDIMENTATION WITH RANDOM DEPOSITS1

ABSTRACT: The accumulated volume of sediment in reservoris is investigated in this paper using a statistical method. Based on the fact that sediment load and river discharge are highly correlated, a relationship between them is established. With longer records of stream flows, this relationship will serve either as a data generation mechanism which produces a sediment S-sequence having the same length as the river discharge Q-sequence, or as a transformation of variables, by which the distribution of S is transferred from the distribution of Q. The mean and variance of the sediment volume accumulated in the design life of a reservoir are then derived, assuming the S- and Q-sequences follow the log-normal distribution. Finally, case studies are given for illustrating the proposed method.     

Operational Evapotranspiration Mapping Using Remote Sensing and Weather Datasets: A New Parameterization for the SSEB Approach

The increasing availability of multi‐scale remotely sensed data and global weather datasets is allowing the estimation of evapotranspiration (ET) at multiple scales. We present a simple but robust method that uses remotely sensed thermal data and model‐assimilated weather fields to produce ET for the contiguous United States (CONUS) at monthly and seasonal time scales. The method is based on the Simplified Surface Energy Balance (SSEB) model, which is now parameterized for operational applications, renamed as SSEBop. The innovative aspect of the SSEBop is that it uses predefined boundary conditions that are unique to each pixel for the “hot” and “cold” reference conditions. The SSEBop model was used for computing ET for 12 years (2000‐2011) using the MODIS and Global Data Assimilation System (GDAS) data streams. SSEBop ET results compared reasonably well with monthly eddy covariance ET data explaining 64% of the observed variability across diverse ecosystems in the CONUS during 2005. Twelve annual ET anomalies (2000‐2011) depicted the spatial extent and severity of the commonly known drought years in the CONUS. More research is required to improve the representation of the predefined boundary conditions in complex terrain at small spatial scales. SSEBop model was found to be a promising approach to conduct water use studies in the CONUS, with a similar opportunity in other parts of the world. The approach can also be applied with other thermal sensors such as Landsat.     

Baseline assessment for environmental services payments from satellite imagery: a case study from Costa Rica and Mexico

In this study we evaluate the accuracy of four global and regional forest cover assessments (MODIS, IGBP, GLC2000, PROARCA) as tools for baseline estimation. We conduct this research at the national scale for Costa Rica and for two tropical dry forest study sites in Costa Rica (Santa Rosa) and Mexico (Chamela-Cuixmala). We found that at the national level, the total forest cover accuracy of the four land cover maps was inflated due to an overestimation of forest in areas with an evergreen canopy. However, the four maps greatly underestimated the extent of the deciduous forest (dry forest); an ecosystem that faces high deforestation pressure and poses complications to the mapping of its extent from remotely sensed data. For the tropical dry forest sites, all maps have low forest cover accuracies (mean for Santa Rosa: 27%; mean for Chamela-Cuixmala: 56%). This has implications for policy implementation.