Development of stochastic modeling systems using deterministic models and GIS: Principles and a case study in the Atlantic Zone of Costa Rica

Shuguang Liu; William A. Reiners; Kenneth G. Gerow; David S Schimel; Michael Keller

USGS · 70263458

Development of stochastic modeling systems using deterministic models and GIS: Principles and a case study in the Atlantic Zone of Costa Rica

Shuguang Liu·William A. Reiners·Kenneth G. Gerow·David S Schimel·Michael Keller

Abstract

The most important requirements for large-area environmental modeling are a tight integration between models and data, and a close match of the spatial scale at which the model is developed with the scale at which the model is to be applied. To better match the scale of data with that of the model, we propose a set of principles for the development of stochastic modeling systems based on linkage of deterministic models with GIS data. For modeling purposes, a region is usually rasterized into cells and the environmental conditions of those cells are specified by ranges or classes using GIS data layers. It is not necessary to simulate each and every GIS cell in the study area because many cells may have similar environmental conditions and can be grouped together to form cohorts. We define a cohort as the assembly of the cells sharing a unique combination of environmental conditions within the study region. Multiple model simulations can be performed for any given cohort. For each simulation, some of the parameter values can be randomly generated within the specified environmental conditions of the cohort according to a certain statistical distribution which, in turn, can be specified by GIS data layers. By this method the variance and covariance of environmental variables in space and time are integrated into the simulation processes with these modeling systems to make full use of the available data and to assess the uncertainties of the simulated results. An integrated simulation system between CENTURY model and GIS was developed to demonstrate the value of the concepts imbedded in stochastic simulation systems for large area studies.

Explore related subjects

Conference Paper

Source-reported bounding extent: 9.499325548625436° to 11.16801392273345° latitude; -84.16707378749192° to -82.44230009071387° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Shuguang Liu, William A. Reiners, Kenneth G. Gerow, David S Schimel, Michael Keller. Development of stochastic modeling systems using deterministic models and GIS: Principles and a case study in the Atlantic Zone of Costa Rica. https://pubs.usgs.gov/publication/70263458

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Authoritative knowledge infrastructure for GeoAI: Implications for scientific spatial data providers

Geospatial Artificial Intelligence (GeoAI) increasingly relies on automated analysis of authoritative geospatial data, exposing limitations in data infrastructures originally designed for simple dissemination rather than for knowledge-driven use. From the perspective of scientific spatial data providers, the challenge is enabling GeoAI responsibly while retaining semantic clarity, provenance, reproducibility, and governance. This short paper (i) examines how GeoAI exposed limitations within existing geospatial data infrastructures, and (ii) identifies provider-side requirements for an authoritative knowledge infrastructure. Drawing on GIScience, scientific workflow, and operational AI literature, the paper outlines implications for national mapping organizations that seek to support advanced analytics without jeopardizing scientific rigor or stewardship responsibilities.

Conference Paper↗

USGS2026

Springs of the Catskill Mountains: A conduit flow system in predominantly siliciclastic strata

The Catskill Mountains of New York state are composed of thick sequences of interbedded mudstone, sandstone, and siltstone that are intercalated with thin beds of paludal limestone and calcareous mudstone. Sandstones often contain rip-up clasts of calcareous material at their bases and contain minor amounts of carbonate cement throughout their mass. This minor amount of calcareous material weathers out to produce thin zones of solutional porosity and permeability integrated within perched aquifers that feed perennial hillslope springs. Solutional features are evident within thin carbonate lenses and along fractures in sandstones, occasionally resulting in explorable cave-sized conduits. This study focuses on the behavior of bedrock springs and seeps in the headwaters of the Neversink River basin above the Neversink Reservoir. Bedrock springs are numerous throughout the Catskill Mountains and serve to provide baseflow to the headwater streams that tends to be more persistent during dry times than other baseflow sources. We observed dissolution-influenced conduits within the bedrock, and suffosional conduits within unconsolidated materials. Continuous water temperature at six springs was monitored at hourly intervals from May to October 2022 to investigate the thermal response of groundwater to daily and seasonal changes in air temperature along with event-based meteorological events. Temperature monitoring was augmented at some springs with time-lapse photography to evaluate persistence of flow. The thermal responses of springs to precipitation events were highly variable across sites, with some springs exhibiting flashy responses indicating a strong connection to the surface, while other springs showed muted or no thermal response to events. The seasonal thermal response also varied across sites having lag times of 23-48 days between peak air temperature and peak spring water temperature. The water temperatures and geochemical compositions of these bedrock springs support a conceptual model of multiple residence times, showing both event and seasonal groundwater storage within the conduit flow system that feeds the springs of the Catskill Mountains.

New York↗

USGS2026

Stability changes after weak layer collapse: Insights from an unexpected dataset

While conducting avalanche mitigation or traveling through the backcountry, avalanche professionals sometimes experience a "whumpf": a collapse-mode crack accompanied by surface cracking that does not produce an avalanche. These whumpfs are a clear warning sign of instability in the surrounding terrain. But assessing the collapsed slope itself, if it is steep enough to avalanche, carries much more uncertainty. Can slopes that whumpf and are steep enough to slide still avalanche? The answer to this question supports avalanche mitigators and ski guides in assessing whether to open such slopes. Previous research on collapsed slopes is limited to one study conducted before tests that estimate crack propagation propensity, such as the Propagation Saw Test (PST), were developed. This work showed that weak layer shear strength decreased immediately after collapsing and then increased over time. However, no research has examined crack propagation propensity after a whumpf. Studying propagation propensity is key for assessing how whumpfs immediately affect slope stability. To improve our understanding of changes in propagation propensity following weak layer collapse, we conducted PSTs on a test slope in the hours and days after an unexpected collapse occurred on a layer of surface hoar. Two PSTs at our study site three days prior to this collapse provided an estimate of weak layer propagation propensity. PST cut lengths decreased immediately following the collapse, suggesting that the collapse increased the propagation propensity. However, over the next 4-5 hours, PST cut lengths steadily increased at our site, likely due to ongoing sintering and strengthening in the collapsed weak layer. By the end of day two, we were no longer able to initiate propagation, suggesting substantial stabilization of the weak layer. Our results support prior findings and observations suggesting that slopes may remain capable of producing avalanches immediately after a collapse but likely stabilize over the next day or two. Although this research does not provide answers about precisely when a slope becomes stable after a collapse, it does highlight uncertainties in our understanding and interpretation of the stability of slopes that have whumped.

Montana↗