Great Sand Dunes National Park Research Articles

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes.

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Earthquake swarm near Great Sand Dunes, Colorado, investigated with temportary seismic network and machine learning sesmic phase analysis

In November of 2018, residents living in the Zapata Subdivision south of Great Sand Dunes National Park and Preserve reported hearing and feeling multiple small earthquakes. Reports of additional earthquakes continued, escalating in late February of 2019, when the USGS recorded over 27 magnitude 0.9 and larger earthquakes over a two-day period. Subdivision residents became concerned that these could be foreshocks to a future, larger earthquake. To further study these earthquakes, we installed a temporary network of seismometers in the area during 2019 and used a convolution neural network seismic phase picker along with the GLASS3 associator to detect over 700 earthquakes in a 3.5-month period during the earthquake swarm. The earthquakes were located using a regional velocity model and a double-difference algorithm. The Northern Sangre de Cristo Fault (NSCF) cuts through the subdivision at the base of the Sangre de Cristo Mountains. Based on geologic evidence, it is one of the most active faults in Colorado but has been nearly aseismic historically. Initially, minor movement on the NSCF was suspected of being the geologic source of the earthquakes. However, nearly all recorded epicenters lie east of the trace of west-dipping fault and are not located on it. Instead, the earthquake epicenters define a narrow, linear, east-west-trending zone that projects eastward across the entire Northern Sangre de Cristo Range and into the headwaters of the Huerfano River Valley. We propose several possible geologic sources for the earthquakes including several mapped, but unnamed faults. Available evidence for any particular source in this geologically complex area is not conclusive. Additional geologic and geophysical investigations are needed to better understand the geology of the earthquake swarm and its implications for seismic hazards.

Dune Initiation in a Bimodal Wind Regime

AbstractEarly‐stage bedforms develop into mature dunes through complex interactions between wind, sand transport, and surface topography. Depending on varying environmental and wind conditions, the mechanisms driving dune formation and, ultimately, the shape of nascent dunes may differ markedly. In cases where sand availability is plentiful, the emergence and growth of dunes can be studied with a linear stability analysis of coupled transport and hydrodynamic equations. Until now, this analysis has only been applied using field evidence in unidirectional winds. However, in many areas of the world and on other planets, wind regimes are more often bimodal or multimodal. Here, we investigate field evidence of protodune formation under a bimodal wind regime by applying linear stability analysis to a developing protodune field. Employing recent development of the linear stability theory and experimental research, combined with in situ wind, sediment transport, and topographic measurements during a monthlong field campaign at Great Sand Dunes National Park, Colorado, USA, we predict the spatial characteristics (orientation and wavelength) and temporal evolution (growth rate and migration velocity) of a protodune field. We find that the theoretical predictions compare well with measured dunefield attributes as characterized by high‐resolution Digital Elevation Models measured using repeat terrestrial laser scanning. Our findings suggest that linear stability analysis is a quantitative predictor of protodune development on sandy surfaces with a bimodal wind regime. This result is significant as it offers critical validation of the linear stability analysis for explaining the initiation and development of dunes toward maturity in a complex natural environment.

Water balance as an indicator of natural resource condition: Case studies from Great Sand Dunes National Park and Preserve

Managing climate impacts to natural resources in protected areas can be hampered by lack of monitoring data, poor understanding of natural resource responses to climate, or lack of timely condition assessments that can inform management actions. Here we demonstrate the utility of water balance as a tool for understanding natural resource responses to climate by developing case studies focused on stream flow, vegetation production, and wildfire ignition at Great Sand Dunes National Park and Preserve (GSDNP), U.S.A. The efficacy of water balance to predict these responses stems from the explicit integration of climate with site conditions that modify the effects of climate. This in turn results in estimates of water availability, water use, and water need that are proximal drivers of aquatic and terrestrial natural resource conditions. The water balance model successfully forecasted stream flow (r2 = 0.69, P < 0.001); determined the critical water needs for maintaining annual vegetation production in different vegetation types spanning a large environmental gradient (r2 = 0.18–0.71); and predicted proportion of historic wildfire ignitions in forest (r2 = 0.96–0.99) and non-forest (r2 = 0.96–0.97) vegetation types. Collectively, these case studies demonstrate practical approaches to translate climate data into assessments of natural resource condition that inform long-term planning and near-term strategic actions needed for conservation of protected areas.

Ten unique and charismatic new species of Microgastrinae wasps (Hymenoptera, Braconidae) from North America.

Ten new species within four genera of Microgastrinae parasitoid wasps (Hymenoptera: Braconidae) are described from Canada and United States: Diolcogaster ichiroi, Diolcogaster miamensis, Glyptapanteles pseudotsugae, Microgaster archboldensis, Microgaster syntopic, Microplitis altissimus, Microplitis jorgeluisi, Microplitis juanmanueli, Microplitis julioalbertoi, and Microplitis mariamargaritae. The new taxa are significant because they represent the first North American records of a tropical group (species of the basimacula group in Diolcogaster), exemplify interesting ecological cases (niche-based host selection in Glyptapanteles, syntopic species in Microgaster), and showcase unique morphological features and/or altitudinal records (Microplitis). Most of the new species were collected in protected areas or areas with strong research programs (Archbold Biological Station and hammock forests near Miami, Florida; Great Sand Dunes National Park and Preserve, and Mount Evans Wilderness Area, Colorado; Sapelo Island, Georgia; Tonto National Forest, Arizona), and thus are also of value and interest for conservation and research efforts.

Determining mineralogical variations of aeolian deposits using thermal infrared emissivity and linear deconvolution methods

Abstract We apply linear deconvolution methods to derive mineral and glass proportions for eight field sample training sites at seven dune fields: (1) Algodones, California; (2) Big Dune, Nevada; (3) Bruneau, Idaho; (4) Great Kobuk Sand Dunes, Alaska; (5) Great Sand Dunes National Park and Preserve, Colorado; (6) Sunset Crater, Arizona; and (7) White Sands National Monument, New Mexico. These dune fields were chosen because they represent a wide range of mineral grain mixtures and allow us to gauge a better understanding of both compositional and sorting effects within terrestrial and extraterrestrial dune systems. We also use actual ASTER TIR emissivity imagery to map the spatial distribution of these minerals throughout the seven dune fields and evaluate the effects of degraded spectral resolution on the accuracy of mineral abundances retrieved. Our results show that hyperspectral data convolutions of our laboratory emissivity spectra outperformed multispectral data convolutions of the same data with respect to the mineral, glass and lithic abundances derived. Both the number and wavelength position of spectral bands greatly impacts the accuracy of linear deconvolution retrieval of feldspar proportions (e.g. K-feldspar vs. plagioclase) especially, as well as the detection of certain mafic and carbonate minerals. In particular, ASTER mapping results show that several of the dune sites display patterns such that less dense minerals typically have higher abundances near the center of the active and most evolved dunes in the field, while more dense minerals and glasses appear to be more abundant along the margins of the active dune fields.

Geophysical expression of buried range-front embayment structure: Great Sand Dunes National Park, Rio Grande rift, Colorado

Great Sand Dunes National Park and Preserve (GRSA, Colorado) lies along the eastern margin of the San Luis Basin and the tectonically active Sangre de Cristo fault system that are part of the northern Rio Grande rift. GRSA lies within a prominent embayment in the range front where two separate sections of the Sangre de Cristo fault system intersect. Fault scarps are observed along both intersecting fault zones within older basin-fill alluvium, but have been obscured by the actively migrating dunefield. The dune sand is also strongly magnetic, locally limiting the usefulness of aeromagnetic methods for mapping concealed structure. This study uses airborne geophysical methods, primarily airborne gravity gradient data, along with constraints from geologic mapping and limited subsurface data and groundwater modeling, to interpret the subsurface basin geometry and range-front structure of the embayment. Using forward modeling of the gravity gradient data and locations of faults inferred from gravity gradient and aeromagnetic lineaments, several previously unrecognized tectonic elements are interpreted adjacent to the range front. Some of the largest rift-­related fault offsets are demonstrated to be basinward of the normal fault zones mapped at the surface along the range front of the Sangre de Cristo Mountains, along faults concealed under the dunefield and subparallel to the two fault sections. A fault-bounded structural bench, likely composed of Proterozoic rocks, underlies most of the high dunefield at depths of 500 m to 1 km. The bench is truncated on its southwest margin by a northwest-trending, southwest-dipping normal fault. A northeast-trending, northwest-dipping normal fault with ∼600 m of estimated relief lies under the southern margin of the dunefield and bounds a structurally higher bench of Proterozoic rocks concealed at <400 m depth near the range front. The northwest- and northeast-trending geophysical lineaments generally correspond well with the trends of faults mapped at the surface, and with both pre- and syn-rift structures in the Sangre de Cristo Mountains. Aeromagnetic anomalies are explained by variations in the magnetization of pre-rift rocks, and the strongly magnetic dune sand.

Minimizing economic impacts from post-fire debris flows in the western United States

For individual burned drainage basins, existing hazard models and readily available data can be combined in a geographic information system to rapidly estimate debris-flow-related damages following a wildfire. The results can then be integrated into an optimization model, whose output guides allocation of emergency management funds and selection of cost-optimized debris-flow management strategies for burned areas consisting of multiple drainage basins. This paper describes methods to identify and value elements-at-risk from a range of possible post-fire debris-flow scenarios, methods to integrate these results with common debris-flow mitigation techniques and best management practices, and methods to apply this information to optimize the mitigation decisions for burned areas. Despite the potential to transform the way hazard managers approach debris-flow mitigation decisions following wildfires, natural hazard and social science management models have not previously been linked in the literature. Results from Santa Barbara (California), Great Sand Dunes National Park (Colorado), and Colfax/Las Animas Counties (Colorado, New Mexico) study sites indicate that optimization modeling can be used to select natural hazard management methods whose benefit for mitigation of post-fire debris flows can easily outweigh the cost of implementation.

Plant Community Changes Following Closure of Artesian Wells in Great Sand Dunes National Park, Colorado

Artificial artesian wells have existed in the San Luis Valley of south central Colorado for over 100 years, and they are an important source of water for livestock and wildlife. When Great Sand Dunes National Park and Preserve (GRSA) expanded its boundaries in 2000, ten of these wells were within the new park boundary. Because the habitat surrounding the wells was so severely disturbed by cattle and wildlife trampling, the National Park Service capped the wells in 2010 to promote restoration of the areas to a more natural state. To study changes after well-capping, we compared the plant communities in 2011 and 2012. We measured and compared plant cover, species frequency, and species diversity at recently capped wells on GRSA, at flowing wells (i.e., with water) on adjacent private lands of The Nature Conservancy, and at reference sites. In general, there was little difference in percentage of plant cover and species diversity between well types. For percent cover, annual variation and distance from wellheads were the best predictors for the native and exotic plant species, respectively. Plant species composition varied by well type, with reference sites having the greatest frequency of native plants per site and capped well sites having the greatest frequency of exotic plants per site. For native plants, reference sites had the highest species diversity, while capped wells had the highest community evenness. Overall, our results suggest that in the short term, native plant populations have changed little after well-capping, but the frequency of exotic plants has increased at capped well sites relative to reference sites.

Debris-flow hazard assessment and validation following the Medano Fire, Great Sand Dunes National Park and Preserve, Colorado

A debris-flow hazard assessment was conducted for the Medano Creek drainage basin, at Great Sand Dunes National Park and Preserve, following the 2,400 ha Medano Fire in 2010. Debris-flow probability and volume predictions were made using empirical regression models. Model parameters include burn severity, rainfall intensity, topographic characteristics, and soil properties. Model results provided park resource managers with information on potential sub-basin-specific hazards to park visitors, roads, and campsites. Monitoring of the first significant rainfall events following the fire and the resulting debris-flow responses provided means for model validation. Of the three probability models utilized, two predicted high probability of debris-flow occurrence for all sub-basins that produced debris flows, with numerous false positives, while the third failed to predict high probability in any of the sub-basins. The volume model predicted volumes within approximately one order of magnitude higher than those measured.

DIVERGENCE IS FOCUSED ON FEW GENOMIC REGIONS EARLY IN SPECIATION: INCIPIENT SPECIATION OF SUNFLOWER ECOTYPES

Early in speciation, as populations undergo the transition from local adaptation to incipient species, is when a number of transient, but potentially important, processes appear to be most easily detected. These include signatures of selective sweeps that can point to asymmetry in selection between habitats, divergence hitchhiking, and associations of adaptive genes with environments. In a genomic comparison of ecotypes of the prairie sunflower, Helianthus petiolaris, occurring at Great Sand Dunes National Park and Preserve (Colorado), we found that selective sweeps were mainly restricted to the dune ecotype and that there was variation across the genome in whether proximity to the nondune population constrained or promoted divergence. The major regions of divergence were few and large between ecotypes, in contrast with an interspecific comparison between H. petiolaris and a sympatric congener, Helianthus annuus. In general, the large regions of divergence observed in the ecotypic comparison swamped locus-specific associations with environmental variables. In both comparisons, regions of high divergence occurred in portions of the genetic map with high marker density, probably reflecting regions of low recombination. The difference in genomic distributions of highly divergent regions between ecotypic and interspecific comparisons highlights the value of studies spanning the spectrum of speciation in related taxa.

Expression of terrain and surface geology in high-resolution helicopter-borne gravity gradient (AGG) data: Examples from Great Sand Dunes National Park, Rio Grande Rift, Colorado

Airborne gravity gradient (AGG) data are rapidly becoming standard components of geophysical mapping programs because of their advantages in cost, access, and resolution over measurements of the gravity field on the ground. Unlike conventional techniques that measure the gravity field, AGG methods measure derivatives of the gravity field. This means that effects of terrain and near-surface geology are amplified in AGG data, and that proper terrain corrections are critically important for AGG data processing. However, terrain corrections require reasonable estimates of density for the rocks and sediments that make up the terrain. A recommended philosophical approach is to use the terrain and surface geology, with their strong expression in AGG data, to the interpreter's advantage. An example of such an approach is presented here for an area with difficult ground access and little ground gravity data. Nettleton-style profiling is used with AGG data to estimate the densities of the sand dunefield and adjacent Precambrian rocks from the area of Great Sand Dunes National Park in southern Colorado. Processing of the AGG data using the density estimate for the dunefield allows buried structures, including a hypothesized buried basement bench, to be mapped beneath the sand dunes.

Consequences of flight height and line spacing on airborne (helicopter) gravity gradient resolution in the Great Sand Dunes National Park and Preserve, Colorado

Line spacing and flight height are critical parameters in airborne gravity gradient surveys; the optimal trade-off between survey costs and desired resolution, however, is different for every situation. This article investigates the additional benefit of reducing the flight height and line spacing though a study of a survey conducted over the Great Sand Dunes National Park and Preserve, which is the highest-resolution public-domain airborne gravity gradient data set available, with overlapping high- and lower-resolution surveys. By using Fourier analysis and matched filtering, it is shown that while the lower-resolution survey delineates the target body, reducing the flight height from 80 m to 40 m and the line spacing from 100 m to 50 m improves the recoverable resolution even at basement depths.

Recent nonhybrid origin of sunflower ecotypes in a novel habitat

The genomics of local adaptation is an increasingly active field, providing insights into the forces driving ecological speciation and the repeatability of evolution. Demography and gene flow play an important role in determining the paths by which parallel evolution occurs and the genomic signatures of adaptation. In the annual sunflowers, hybridization between species has repeatedly led to the colonization of extreme habitats, such as sand dunes. In a new case of adaptation to sand dunes that occurs in populations of H. petiolaris growing at Great Sand Dunes National Park and Preserve (Colorado), we wished to determine the age and long-term migration patterns of the system, as well as its ancestry. We addressed these questions with restriction-associated DNA (RAD) sequence data, aligned to a reference transcriptome. In an isolation with migration model using RAD sequences, coalescent analysis showed that the dune ecotype originated since the last ice age, which is very recent compared with the hybrid dune species, H. anomalus. Large effective population sizes and substantial numbers of gene migrants per generation between dune and nondune ecotypes explained the highly heterogeneous divergence observed among loci. Analysis of RAD-derived SNPs identified heterogeneous divergence between the dune and nondune ecotypes, as well as identifying its nearest relative. Our results did not support the hypothesis that the dune ecotype has hybrid ancestry, suggesting that adaptation of sunflowers to dunes has occurred by multiple mechanisms. The ancestry and long-term history of gene flow between incipient sunflower species provides valuable context for our understanding of ecological speciation and parallel adaptation.

Variable wind ripple migration at Great Sand Dunes National Park and Preserve, observed by timelapse imaging

Granule ripples at Great Sand Dunes National Park and Preserve (GSDNPP) were observed using inexpensive digital timelapse cameras over a 70-day period in winter 2010–2011. The ripples migrated during a handful of discrete events — visible ripple movement occurred on only 11 days during the observation period. The movement conditions are documented with hourly and 15-minute records from two nearby weather stations, and by a cup anemometer at the site itself. During the most prominent movement episode, when local winds averaged ~ 10 m/s, ripples of several sizes were observed simultaneously and a reciprocal relationship of ripple size and propagation speed was seen, with small (~ 10 cm) ripples moving at ~ 1.4 cm/min, and larger (~ 80 cm) ripples at ~ 0.15 cm/min. Ripple sizes and morphologies evolve throughout the observation period.

Summary of the Second International Planetary Dunes Workshop: Planetary Analogs — Integrating Models, Remote Sensing, and Field Data, Alamosa, Colorado, USA, May 18–21, 2010

The Second International Planetary Dunes Workshop took place in Alamosa, Colorado, USA from May 18–21, 2010. The workshop brought together researchers from diverse backgrounds to foster discussion and collaboration regarding terrestrial and extra-terrestrial dunes and dune systems. Two and a half days were spent on five oral sessions and one poster session, a full-day field trip to Great Sand Dunes National Park, with a great deal of time purposefully left open for discussion. On the last day of the workshop, participants assembled a list of thirteen priorities for future research on planetary dune systems.

Water Rights for New Federal Land Conservation Programs: A Turn-of-The-Century Evaluation

This article addresses how to integrate the water requirements associated with new federal land conservation initiatives – such as new national parks, wildlife refuges, conservation areas, recreation areas, and wilderness any policy areas – with state water law systems. It looks at past controversies over federal reserved rights, and suggests a kind of third way – rights substantively defined by federal law and procedural perfected through state water rights systems, as a reasonable compromise between state and national interests. It touts the approach Congress used in designating the Great Sand Dunes National Park as a model.

Assessment of lake sensitivity to acidic deposition in national parks of the Rocky Mountains

The sensitivity of high-elevation lakes to acidic deposition was evaluated in five national parks of the Rocky Mountains based on statistical relations between lake acid-neutralizing capacity concentrations and basin characteristics. Acid-neutralizing capacity (ANC) of 151 lakes sampled during synoptic surveys and basin-characteristic information derived from geographic information system (GIS) data sets were used to calibrate the statistical models. The explanatory basin variables that were considered included topographic parameters, bedrock type, and vegetation type. A logistic regression model was developed, and modeling results were cross-validated through lake sampling during fall 2004 at 58 lakes. The model was applied to lake basins greater than 1 ha in area in Glacier National Park (n = 244 lakes), Grand Teton National Park (n = 106 lakes), Great Sand Dunes National Park and Preserve (n = 11 lakes), Rocky Mountain National Park (n = 114 lakes), and Yellowstone National Park (n = 294 lakes). Lakes that had a high probability of having an ANC concentration <100 microeq/L, and therefore sensitive to acidic deposition, are located in basins with elevations >3000 m, with <30% of the catchment having northeast aspect and with >80% of the catchment bedrock having low buffering capacity. The modeling results indicate that the most sensitive lakes are located in Rocky Mountain National Park and Grand Teton National Park. This technique for evaluating the lake sensitivity to acidic deposition is useful for designing long-term monitoring plans and is potentially transferable to other remote mountain areas of the United States and the world.

The rate of granule ripple movement on Earth and Mars

Abstract The rate of movement for 3- and 10-cm-high granule ripples was documented in September of 2006 at Great Sand Dunes National Park and Preserve during a particularly strong wind event. Impact creep induced by saltating sand caused ∼24 granules min −1 to cross each cm of crest length during wind that averaged ∼9 m s −1 (at a height well above 1 m), which is substantially larger than the threshold for saltation of sand. Extension of this documented granule movement rate to Mars suggests that a 25-cm-high granule ripple should require from hundreds to thousands of Earth-years to move 1 cm under present atmospheric conditions.

Estimating the Economic Value of National Parks with Count Data Models Using On-Site, Secondary Data: The Case of the Great Sand Dunes National Park and Preserve

We estimate an individual travel cost model for Great Sand Dunes National Park and Preserve (GSD) in Colorado using on-site, secondary data. The purpose of the on-site survey was to help the National Park Service better understand the visitors of GSD; it was not intended for a travel cost model. Variables such as travel cost and income were estimated based on respondents' Zip Codes. Following approaches found in the literature, a negative binomial model corrected for truncation and endogenous stratification fit the data the best. We estimate a recreational benefit of U.S. $89/visitor/year or U.S. $54/visitor/24-h recreational day (in 2002 U.S. $). Based on the approach presented here, there are other data sets for national parks, preserves, and battlefields where travel cost models could be estimated and used to support National Park Service management decisions.