USGS Science⌕ Search

SEARCH · USGS Science

Results for “Forests”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83Linked to original sources
USGS1993

Snaring to control feral pigs sus scrofa in a remote Hawaiian rain forest

Feral pig Sus scrofa control in Kipahulu Valley, a remote rain forest in Haleakala National Park, Maui, Hawaiian Islands, has been achieved with snares over a 45-month period. Initial pig densities in fenced management units of 6·2 km 2 and 7·8 km 2 were estimated at 6 animals/km 2 and 14·3 animals/km 2 for the two units, based on population reconstruction from animals killed and aged. During the 45 months of the study, 1978 snares were set, and 1·6 million snare nights were logged. Snare density reached 96/km 2 and 200/km 2 for the two management units by the end of the study. A mean effort of 43 worker hours/pig was used to remove 53 pigs from the upper management unit, and a mean of 7 worker hours/pig to remove 175 animals from the more densely populated lower unit. Pig activity monitoring along transects provided a good measure of control effectiveness until densities of about 1 pig/km 2 were achieved, after which transects became less useful than scouting for determining pig activity.

Hawai'i↗
USGS1995

Denudation rates determined from the accumulation of in situ-produced 10Be in the luquillo experimental forest, Puerto Rico

We present a simple method for estimation of long-term mean denudation rates using in situ-produced cosmogenic 10Be in fluvial sediments. Procedures are discussed to account for the effects of soil bioturbation, mass wasting and attenuation of cosmic rays by biomass and by local topography. Our analyses of 10Be in quartz from bedrock outcrops, soils, mass-wasting sites and riverine sediment from the Icacos River basin in the Luquillo Experimental Forest, Puerto Rico, are used to characterize denudation for major landform elements in that basin. The 10Be concentration of a discharge-weighted average of size classes of river sediment corresponds to a long-term average denudation of ≈ 43 m Ma −1, consistent with mass balance results.

Earth and Planetary Science Letters↗
USGS1985

Chemistry and transport of soluble humic substances in forested watersheds of the Adirondack Park, New York

Studies were conducted in conjunction with the Integrated Lake-Watershed Acidification Study (ILWAS) to examine the chemistry and leaching patterns of soluble humic substances in forested watersheds of the Adirondack region. During the summer growing season, mean dissolved organic carbon (DOC) concentrations in the ILWAS watersheds ranged from 21–32 mg C l −1 in O/A horizon leachates, from 5–7 mg C l −1 in B horizon leachates, from 2–4 mg C l −1 in groundwater solutions, from 6–8 mg C l −1 in first order streams, from 3–8 mg C l −1 in lake inlets, and from 2–7 mg C l −1 in lake outlets. During the winter, mean DOC concentrations dropped significantly in the upper soil profile. Soil solutions from mixed and coniferous stands contained as much as twice the DOC concentration of lysimeter samples from hardwood stands. Results of DOC fractionation analysis showed that hydrophobia and hydrophilic acids dominate the organic solute composition of natural waters in these watersheds. Charge balance and titration results indicated that the general acid-base characteristics of the dissolved humic mixture in these natural waters can be accounted for by a model organic acid having an average p K a of 3.85, an average charge density of 4–5 μeq mg −1 C at ambient pH, and a total of 6–7 meq COOH per gram carbon.

New York↗
USGS1990

Effects of experimental overgrowth on survival and change in the turf assemblage of a giant kelp forest

Crustose coralline algae were the prevalent cover among sessile organisms that paved or grew near the substratum, and also the most commonly overgrown species in a giant kelp Macrocystis pyrifera (L.) C.A. Agardh forest located off San Nicolas Island, California. Giant kelp was the largest and most conspicuous species that overgrew large patches of the substrata; overgrowth among turf organisms also appeared common. To determine the effects of giant kelp holdfasts on crustose coralline algae and other turf organisms, “artificial holdfasts” were placed on 0.125-m 2 plots for 5, 8 and 12 months. In these treatments, 50–57% of the crustose coralline algae survived. Because these algae also recruited while covered, the total cover (survivorship plus recruitment) differed by only 7–26% from that sampled at the start of the study. The decline of these algae in control plots was similar to that in the treatment plots mostly because of overgrowth by sessile invertebrates. Bryozoans increased markedly on the control plots, whereas 0–12% survived in the treatment plots. Bryozoans and sponges also recruited under the artificial holdfasts. Some arborescent turf algae survived in the 5- and 8-month treatments; articulated coralline algae survived better than did foliose algae. High survival recruitment of crustose coralline algae while overgrown contributed to their prevalence in benthic communities.

Journal of Experimental Marine Biology and Ecology↗
USGS2017

Dated wood from Alaska and the Yukon: Implications for forest refugia in Beringia

Postulations on the existence of forest refugia in parts of Beringia during the last glacial have been, in large part, based on ambiguous evidence. Existing data on radiocarbon-dated and identified fossil wood and macrofossils from Alaska and northwest Canada are synthesized here and are augmented by results of palynological studies in an effort to show the persistence of some, and total extinction of other, tree and large shrub species. Possible dispersal routes taken by species that reinvaded Beringia in postglacial times are also reconstructed from the fossil record. Macrofossil and pollen evidence, when combined with climatic factors, makes cottonwood a good candidate for survival during the last glacial. Larch and aspen are also candidates, though the evidence for them is less positive. Pollen and macrofossils of alder are very scarce in deposits of the last glacial age, and if it survived at all, it was probably in very isolated vegetatively reproducing clones. Shrub birch may have been present in Beringia, but tree birch probably was reintroduced during the Holocene. Spruce also appears to have been absent in Alaska from about 30,000 to 11,500 yr ago and probably reinvaded Beringia from a refugium south of the Laurentide ice sheet.

Alaska↗
USGS1989

Characterization of alkyl carbon in forest soils by CPMAS 13C NMR spectroscopy and dipolar dephasing

Samples obtained from forest soils at different stages of decomposition were treated sequentially with chloroform/methanol (extraction of lipids), sulfuric acid (hydrolysis), and sodium chlorite (delignification) to enrich them in refractory alkyl carbon. As revealed by NMR spectroscopy, this treatment yielded residues with high contents of alkyl carbon. In the NMR spectra of residues obtained from litter samples, resonances for carbohydrates are also present, indicating that these carbohydrates are tightly bound to the alkyl carbon structures. During decomposition in the soils this resistant carbohydrate fraction is lost almost completely. In the litter samples the alkyl carbon shows a dipolar dephasing behavior indicative of two structural components, a rigid and a more mobile component. As depth and decomposition increase, only the rigid component is observed. This fact could be due to selective degradation of the mobile component or to changes in molecular mobility during decomposition, e.g., because of an increase in cross linking or contact with the mineral matter of the soil.

Science of the Total Environment↗
USGS1988

Preliminary observations of streamflow generation during storms in a forested Piedmont watershed using temperature as a tracer

Variations in streamwater temperature at the outlet of a 41-ha forested watershed at Panola Mountain in the Georgia Piedmont indicate that the initial rapid hydrologic response is caused by a combination of groundwater discharge and channel interception of rainwater. A storm in May 1986 caused a rapid increase in discharge that was accompanied by a decrease in streamwater temperature and a rise in the water table level adjacent to the stream. The higher water table provided the hydraulic gradient necessary to increase the discharge of colder groundwater to the stream. Storms that occurred under very dry antecedent conditions in July 1986 and June 1987 caused a rapid hydrologic response but no change in water table level, indicating the response was caused by channel interception of rainwater. This conclusion was supported by increases in streamwater temperature in the June storm and by chemical changes in the July storm. When rainfall is sufficient, flow in the ephemeral part of the stream in the catchment headwaters generates a second and larger discharge peak that reflects the chemistry and temperature of runoff from a 3-ha granite outcrop in the headwaters; sulfate concentration and temperature increase and alkalinity decreases relative to prestorm conditions. The initial response, however, results from channel interception and groundwater discharge. Rapid rises in the water table level during some storms suggest that macropore flow may play a major role in the hydrologic response of the watershed to rainstorms.

Journal of Contaminant Hydrology↗
USGS2015

Stable-isotope and solute-chemistry approaches to flow characterization in a forested tropical watershed, Luquillo Mountains, Puerto Rico

The prospect of changing climate has led to uncertainty about the resilience of forested mountain watersheds in the tropics. In watersheds where frequent, high rainfall provides ample runoff, we often lack understanding of how the system will respond under conditions of decreased rainfall or drought. Factors that govern water supply, such as recharge rates and groundwater storage capacity, may be poorly quantified. This paper describes 8-year data sets of water stable isotope composition ( δ 2 H and δ 18 O) of precipitation (4 sites) and a stream (1 site), and four contemporaneous stream sample sets of solute chemistry and isotopes, used to investigate watershed response to precipitation inputs in the 1780-ha Río Mameyes basin in the Luquillo Mountains of northeastern Puerto Rico. Extreme δ 2 H and δ 18 O values from low-pressure storm systems and the deuterium excess ( d -excess) were useful tracers of watershed response in this tropical system. A hydrograph separation experiment performed in June 2011 yielded different but complementary information from stable isotope and solute chemistry data. The hydrograph separation results indicated that 36% of the storm rain that reached the soil surface left the watershed in a very short time as runoff. Weathering-derived solutes indicated near-stream groundwater was displaced into the stream at the beginning of the event, followed by significant dilution. The more biologically active solutes exhibited a net flushing behavior. The d- excess analysis suggested that streamflow typically has a recent rainfall component (∼25%) with transit time less than the sampling resolution of 7 days, and a more well-mixed groundwater component (∼75%). The contemporaneous stream sample sets showed an overall increase in dissolved solute concentrations with decreasing elevation that may be related to groundwater inputs, different geology, and slope position. A considerable amount of water from rain events runs off as quickflow and bypasses subsurface watershed flowpaths, and better understanding of shallow hillslope and deeper groundwater processes in the watershed will require sub-weekly data and detailed transit time modeling. A combined isotopic and solute chemistry approach can guide further studies to a more comprehensive model of the hydrology, and inform decisions for managing water supply with future changes in climate and land use.

Applied Geochemistry↗
USGS2008

Advances and limitations of individual-based models to analyze and predict dynamics of mangrove forests: A review

Mangrove ecosystems are considered vulnerable to climate change as coastal development limits the ecosystem services and adaptations important to their survival. Although they appear rather simple in terms of species diversity, their ecology is complex due to interacting geophysical forces of tides, surface runoff, river and groundwater discharge, waves, and constituents of sediment, nutrients and saltwater. These interactions limit developing a comprehensive framework for science-based sustainable management practices. A suite of models have been developed independently by various academic and government institutions worldwide to understand the dynamics of mangrove ecosystems and to provide ecological forecasting capabilities under different management scenarios and natural disturbance regimes. The models have progressed from statistical tables representing growth and yield to more sophisticated models describing various system components and processes. Among these models are three individual-based models (IBMs) (FORMAN, KIWI, and MANGRO). A comparison of models’ designs reveal differences in the details of process description, particularly, regarding neighbor competition among trees. Each model has thus its specific range of applications. Whereas FORMAN and KIWI are most suitable to address mangrove forest dynamics of stands, MANGRO focuses on landscape dynamics on larger spatial scale. A comparison of the models and a comparison of the models with empirical knowledge further reveal the general needs for further field and validation studies to advance our ecological understanding and management of mangrove wetlands.

Aquatic Botany↗
USGS2016

Atmospheric inputs of organic matter to a forested watershed: Variations from storm to storm over the seasons

The objectives of this study were to determine the quantity and chemical composition of precipitation inputs of dissolved organic carbon (DOC) to a forested watershed; and to characterize the associated temporal variability. We sampled most precipitation that occurred from May 2012 through August 2013 at the Susquehanna Shale Hills Critical Zone Observatory (Pennsylvania, USA). Sub-event precipitation samples (159) were collected sequentially during 90 events; covering various types of synoptic meteorological conditions in all climatic seasons. Precipitation DOC concentrations and rates of wet atmospheric DOC deposition were highly variable from storm to storm, ranging from 0.3 to 5.6 mg C L −1 and from 0.5 to 32.8 mg C m −2 h −1 , respectively. Seasonally, storms in spring and summer had higher concentrations of DOC and more optically active organic matter than in winter. Higher DOC concentrations resulted from weather types that favor air advection, where cold frontal systems, on average, delivered more than warm/stationary fronts and northeasters. A mixed modeling statistical approach revealed that factors related to storm properties, emission sources, and to the chemical composition of the atmosphere could explain more than 60% of the storm to storm variability in DOC concentrations. This study provided observations on changes in dissolved organic matter that can be useful in modeling of atmospheric oxidative chemistry, exploring relationships between organics and other elements of precipitation chemistry, and in considering temporal changes in ecosystem nutrient balances and microbial activity.

Atmospheric Environment↗
USGS2007

Patch-occupancy models indicate human activity as major determinant of forest elephant Loxodonta cyclotis seasonal distribution in an industrial corridor in Gabon

The importance of human activity and ecological features in influencing African forest elephant ranging behaviour was investigated in the Rabi-Ndogo corridor of the Gamba Complex of Protected Areas in southwest Gabon. Locations in a wide geographical area with a range of environmental variables were selected for patch-occupancy surveys using elephant dung to assess seasonal presence and absence of elephants. Patch-occupancy procedures allowed for covariate modelling evaluating hypotheses for both occupancy in relation to human activity and ecological features, and detection probability in relation to vegetation density. The best fitting models for old and fresh dung data sets indicate that (1) detection probability for elephant dung is negatively related to the relative density of the vegetation, and (2) human activity, such as presence and infrastructure, are more closely associated with elephant distribution patterns than are ecological features, such as the presence of wetlands and preferred fresh fruit. Our findings emphasize the sensitivity of elephants to human disturbance, in this case infrastructure development associated with gas and oil production. Patch-occupancy methodology offers a viable alternative to current transect protocols for monitoring programs with multiple covariates.

Biological Conservation↗
USGS2013

Response of palila and other subalpine Hawaiian forest bird species to prolonged drought and habitat degradation by feral ungulates

Extinction has claimed half of all historically-known Hawaiian passerines, and today many extant species are increasingly threatened due to the combined effects of invasive species and climate change. Habitat disturbance has affected populations of feeding specialists most profoundly, and our results indicate that specialists continue to be most vulnerable, although even some abundant, introduced, generalist species also may be affected. Surveys of passerines during 1998–2011 in subalpine woodland habitat on Mauna Kea Volcano, Island of Hawai′i, revealed that the abundance of the critically endangered palila ( Loxioides bailleui ), a seed specialist, declined by 79% after 2003. The ′akiapōlā′au ( Hemignathus munroi ), an endangered specialist insectivore, was not detected in the survey area after 1998. The Hawai′i ′amakihi ( Hemignathus virens virens ), a generalist feeder and the most abundant species on Mauna Kea, was the only native species to maintain a stable population. The Japanese white-eye ( Zosterops japonicus ), a well-entrenched generalist and one of the three most common introduced species, declined. Drought prevailed in 74% of months during 2000–2011, and dry conditions contributed to the recent decline of the palila by reducing the annual māmane ( Sophora chrysophylla ) seed pod crop, which influences palila breeding and survival. Sustained browsing by introduced ungulates also lowered habitat carrying capacity, and their elimination should reduce the effects of drought and promote forest restoration. Our results illustrate how the feeding ecology of a species can influence its response to interacting environmental perturbations, and they underscore the value of long-term monitoring to detect population trends of sensitive species.

Hawai'i↗
USGS2015

Comparing spatial capture–recapture modeling and nest count methods to estimate orangutan densities in the Wehea Forest, East Kalimantan, Indonesia

Accurate information on the density and abundance of animal populations is essential for understanding species' ecology and for conservation planning, but is difficult to obtain. The endangered orangutan ( Pongo spp.) is an example; due to its elusive behavior and low densities, researchers have relied on methods that convert nest counts to orangutan densities and require substantial effort for reliable results. Camera trapping and spatial capture–recapture (SCR) models could provide an alternative but have not been used for primates. We compared density estimates calculated using the two methods for orangutans in the Wehea Forest, East Kalimantan, Indonesia. Camera trapping/SCR modeling produced a density estimate of 0.16 ± 0.09–0.29 indiv/km 2 , and nest counts produced a density estimate of 1.05 ± 0.18–6.01 indiv/km 2 . The large confidence interval of the nest count estimate is probably due to high variance in nest encounter rates, indicating the need for larger sample size and the substantial effort required to produce reliable results using this method. The SCR estimate produced a very low density estimate and had a narrower but still fairly wide confidence interval. This was likely due to unmodeled heterogeneity and small sample size, specifically a low number of individual captures and recaptures. We propose methodological fixes that could address these issues and improve precision. A comparison of the overall costs and benefits of the two methods suggests that camera trapping/SCR modeling can potentially be a useful tool for assessing the densities of orangutans and other elusive primates, and warrant further investigation to determine broad applicability and methodological adjustments needed.

East Kalimantan↗
USGS2019

Satellite-detected forest disturbance forecasts American marten population decline: The case for supportive space-based monitoring

Limited monitoring resources often constrain rigorous monitoring practices to species or populations of conservation concern. Insufficient monitoring can induce a tautology as lack of monitoring resources makes it difficult to determine whether a species or population deserves additional monitoring resources. When in-situ monitoring resources are limited, remote habitat monitoring could be a useful supplementary tool, as linking parameterized species distribution models to spatially explicit time-series of environmental correlates allows iterative prediction of population change. Yet the performance of predictive forecasts or hindcasts has been difficult to evaluate. We paired contemporary field data, historic population estimates, and a remotely-sensed archive of landscape change to evaluate predictions of American marten ( Martes americana ) population decline owing to habitat loss in Maine, USA. We estimated contemporary spatial patterns in marten density relative to landscape disturbance with spatial capture-recapture models. We compared current density estimates to historical density calculations to evaluate population decline, and compared historical calculations to habitat-based model predictions to evaluate the efficacy of habitat monitoring as a proxy for direct monitoring. Marten density was negatively associated with the proportion of surrounding regenerating forest, and point estimates within focal townships were 50–80% lower than historical calculations. Habitat-based hindcasts of marten density across our entire focal area interest suggested a smaller population decline (roughly 50%) within our focal area. Thus, although habitat-based predictions underpredicted marten decline, they provided correct directional inference. Habitat monitoring and predictions from species distribution models may provide useful inference about population changes given trends in habitat at limited expense when in-situ information is lacking.

Maine↗
USGS2022

Rainforest carnivore ecology in a managed forest reserve: Differential seasonal correlates between habitat components and relative abundance

Studies of relationships between seasons and Neotropical carnivore distributions tend to focus on water and prey availability without considering other habitat components such as escape, foraging, and resting cover. Our goal was to evaluate habitat characteristics that may be important for predicting the seasonal (dry or rainy) relative abundance of four commonly captured Neotropical carnivores (i.e., jaguar [ Panthera onca ], puma [ Puma concolor ], ocelot [ Leopardus pardalis ], and grey fox [ Urocyon cinereoargenteus ]) in Chiquibul Forest Reserve in Belize, Central America. We used trail camera data and random-effect Poisson models to investigate how prey ratios (number of prey detections/total detections), cover (e.g., logs and stumps used for hiding cover from predators), vegetation structure, and environmental site characteristics (e.g., site harvest-history, slope, aspect) were related to carnivore relative abundance. Both prey ratios and vegetation structure appeared in supported models more frequently than other environmental site characteristics and were negatively correlated with carnivore relative abundance. Supported models differed for each season for all species except jaguars for which mammalian prey ratios and prey cover at sites was always negatively correlated with jaguar relative abundance. Carnivores appeared to avoid sites where vegetation created ideal escape and hiding cover for prey even though prey may be less abundant. Our data suggest that vegetation structure and composition can create conditions conducive to carnivore foraging and that these characteristics can differ by season in the tropics.

Biological Conservation↗
USGS2024

Climatic drought and trophic disruption in an endemic subalpine Hawaiian forest bird

Overexploitation, habitat conversion, and introduced species have caused unprecedented extinctions and heavily degraded native bird populations in island ecosystems. In the Hawaiian Islands, stemming these losses has proven difficult as the highly specialized avifauna are often impacted – among other things – by poorly understood trophic disruptions as well as persistent climatic shifts. Here we investigate these dynamics by examining the trophic positions of 7 different taxa of producers and consumers across the last century in the subalpine māmane ( Sophora chrysophylla ) forest ecosystem on the island of Hawaiʻi. From museum collections and contemporary sampling, we analyzed the stable nitrogen (δ 15 N) values of producers and consumers to identify trophic and source amino acids and derive trophic enrichment constants specific to this food web. This enables us to reconstruct the diet of the palila ( Loxioides bailleui ), a critically endangered finch whose population recently declined 90 %. Our results show that from 1890 to 2006, the palila trophic position declined from 2.6 to 2.2, with cascading implications for its diet. Bayesian mixing model reconstructions indicate that palila trophic position changes likely arose from a 76 % decline (69.3 % to 16.6 % of diet) in the consumption of native moth caterpillars, and a 172 % increase (30.7 % to 83.4 %) in native plants. From the available ecosystem variables, exploratory Bayesian multiple regressions selected surface temperature changes, and the interactions of surface temperatures with drought and caterpillar parasitism as the primary drivers of these trophic changes. Despite the predicted increases of warming and drought, management interventions may build resiliency in this unique island ecosystem.

Hawaii↗
USGS2010

Effects of temperature on silicate weathering: Solute fluxes and chemical weathering in a temperate rain forest watershed, Jamieson Creek, British Columbia

Chemical weathering of silicate minerals has long been known as a sink for atmospheric CO 2 , and feedbacks between weathering and climate are believed to affect global climate. While warmer temperatures are believed to increase rates of weathering, weathering in cool climates can be accelerated by increased mineral exposure due to mechanical weathering by ice. In this study, chemical weathering of silicate minerals is investigated in a small temperate watershed. The Jamieson Creek watershed is covered by mature coniferous forest and receives high annual precipitation (4000 mm), mostly in the form of rainfall, and is underlain by quartz diorite bedrock and glacial till. Analysis of pore water concentration gradients indicates that weathering in hydraulically unsaturated ablation till is dominated by dissolution of plagioclase and hornblende. However, a watershed scale solute mass balance indicates high relative fluxes of K and Ca, indicating preferential leaching of these solutes possibly from the relatively unweathered lodgement till. Weathering rates for plagioclase and hornblende calculated from a watershed scale solute mass balance are similar in magnitude to rates determined using pore water concentration gradients. When compared to the Rio Icacos basin in Puerto Rico, a pristine tropical watershed with similar annual precipitation and bedrock, but with dissimilar regolith properties, fluxes of weathering products in stream discharge from the warmer site are 1.8 to 16.2-fold higher, respectively, and regolith profile-averaged plagioclase weathering rates are 3.8 to 9.0-fold higher. This suggests that the Arrhenius effect, which predicts a 3.5- to 9-fold increase in the dissolution rate of plagioclase as temperature is increased from 3.4° to 22 °C, may explain the greater weathering fluxes and rates at the Rio Icacos site. However, more modest differences in K and Ca fluxes between the two sites are attributed to accelerated leaching of those solutes from glacial till at Jamieson Creek. Our findings suggest that under conditions of high rainfall and favorable topography, weathering rates of silicate minerals in warm tropical systems will tend to be higher than in cool temperate systems, even if the temperate system is has been perturbed by an episode of glaciation that deposits regolith high in fresh mineral surface area.

British Columbia↗
USGS2015

Phenology and climate relationships in aspen (Populus tremuloides Michx.) forest and woodland communities of southwestern Colorado

Trembling aspen ( Populus tremuloides Michx.) occurs over wide geographical, latitudinal, elevational, and environmental gradients, making it a favorable candidate for a study of phenology and climate relationships. Aspen forests and woodlands provide numerous ecosystem services, such as high primary productivity and biodiversity, retention and storage of environmental variables (precipitation, temperature, snow–water equivalent) that affect the spring and fall phenology of the aspen woodland communities of southwestern Colorado. We assessed the land surface phenology of aspen woodlands using two phenology indices, start of season time (SOST) and end of season time (EOST), from the U.S. Geological Survey (USGS) database of conterminous U.S. phenological indicators over an 11-year time period (2001–2011). These indicators were developed with 250 m resolution remotely sensed data from the Moderate Resolution Imaging Spectroradiometer processed to highlight vegetation response. We compiled data on SOST, EOST, elevation, precipitation, air temperature, and snow water equivalent (SWE) for selected sites having more than 80% cover by aspen woodland communities. In the 11-year time frame of our study, EOST had significant positive correlation with minimum fall temperature and significant negative correlation with fall precipitation. SOST had a significant positive correlation with spring SWE and spring maximum temperature.

Colorado↗