Geology ReportsSearch

SEARCH · Geology Reports

Results for “SOIL”

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 595 records · Page 33Linked to original sources

An improved method for quantifying soil macroporosity

Quantitative information on macroporosity is needed to predict water flow and solute transport in field soils. A method was developed for determining the number, shape, and size distribution of soil macropores. Horizontal serial sections sawed from paraffin-impregnated soil cores were photographed under ultraviolet (UV) light. Anthracene, mixed with the paraffin, fluoresces a bright bluish white under UV light and provides a sharp contrast between the soil matrix and the paraffin-filled pore space. Section photographs were converted to 256 level, grey-scale digital images using a flat-bed scanner. Image processing was used to classify each pixel in a digital image as pore space or soil matrix, to group the pore space pixels into pores, and to measure the area and perimeter of each pore. The method was able to measure pores with an equivalent radius ≥85 µm. Macroporosity in soil cores sampled form a tillage path and from an adjacent, undisturbed (notill) region was quantified. Tillage sections contained, on the average, 9.4 macropores/cm 2 with an equivalent macroporosity of 8%, while no-till sections contained 0.8 macropores/cm 2 with an equivalent macroporosity of 0.3%. Computed intrinsic permeabilities for tillage sections that included macropore information were significantly larger than values computed using micropore information alone, suggesting that macropores must be included in permeability calculations when the number of macropores is large. The developed method should be useful for quantifying macroporosity in nonskeletal soils.

Soil Science Society of America Journal

Pedogenic silica accumulation in chronosequence soils, southern California

Chronosequential analysis of soil properties has proven to be a valuable approach for estimating ages of geomorphic surfaces where no independent age control exists. In this study we examined pedogenic silica as an indicator of relative ages of soils and geomorphic surfaces, and assessed potential sources of the silica. Pedogenic opaline silica was quantified by tiron (4,5-dihydroxy-1,3-benzene-disulfonic acid [disodium salt], C 6 H 4 Na 2 O 8 S 2 ) extraction for pedons in two different chronosequences in southern California, one in the San Timoteo Badlands and one in Cajon Pass. The soils of both of these chronosequences are developed in arkosic sediments and span 11.5 to 500 ka. The amount of pedogenic silica increases with increasing duration of pedogenesis, and the depth of the maximum silica accumulation generally coincides with the maximum expression of the argillic horizon. Pedogenic silica has accumulated in all of the soils, ranging from 1.2% tiron-extractable Si (Si tn ) in the youngest soil to 4.6% in the oldest. Primary Si decreases with increasing duration of weathering, particularly in the upper horizons, where weathering conditions are most intense. The loss of Si coincides with the loss of Na and K, implicating the weathering of feldspars as the likely source of Si loss. The quantity of Si lost in the upper horizons is adequate to account for the pedogenic silica accumulation in the subsoil. Pedogenic silica was equally effective as pedogenic Fe oxides as an indicator of relative soil age in these soils.

California

Soil organic carbon stocks in Alaska estimated with spatial and pedon data

Temperatures in high-latitude ecosystems are increasing faster than the average rate of global warming, which may lead to a positive feedback for climate change by increasing the respiration rates of soil organic C. If a positive feedback is confirmed, soil C will represent a source of greenhouse gases that is not currently considered in international protocols to regulate C emissions. We present new estimates of the stocks of soil organic C in Alaska, calculated by linking spatial and field data developed by the USDA NRCS. The spatial data are from the State Soil Geographic database (STATSGO), and the field and laboratory data are from the National Soil Characterization Database, also known as the pedon database. The new estimates range from 32 to 53 Pg of soil organic C for Alaska, formed by linking the spatial and field data using the attributes of Soil Taxonomy. For modelers, we recommend an estimation method based on taxonomic subgroups with interpolation for missing areas, which yields an estimate of 48 Pg. This is a substantial increase over a magnitude of 13 Pg estimated from only the STATSGO data as originally distributed in 1994, but the increase reflects different estimation methods and is not a measure of the change in C on the landscape. Pedon samples were collected between 1952 and 2002, so the results do not represent a single point in time. The linked databases provide an improved basis for modeling the impacts of climate change on net ecosystem exchange.

Soil Science Society of America Journal

Strategies for soil quality assessment using VNIR gyperspectral spectroscopy in a western Kenya Chronosequence

Visible and near-infrared reflectance spectroscopy (VNIRS) is a rapid and nondestructive method that can predict multiple soil properties simultaneously, but its application in multidimensional soil quality (SQ) assessment in the tropics still needs to be further assessed. In this study, VNIRS (350–2500 nm) was employed to analyze 227 air-dried soil samples of Ultisols from a soil chronosequence in western Kenya and assess 16 SQ indicators. Partial least squares regression (PLSR) was validated using the full-site cross-validation method by grouping samples from each farm or forest site. Most suitable models successfully predicted SQ indicators (R 2 ≥ 0.80; ratio of performance to deviation [RPD] ≥ 2.00) including soil organic matter (OM LOI ), active C, Ca, cation exchange capacity (CEC), and clay. Moderately-well predicted indicators (0.50 ≤ R 2 < 0.80; 1.40 ≤ RPD < 2.00) were water stable aggregation (WSA), Cu, silt, Mg, pH, sand, water content at permanent wilting point (Θ pwp ), and field capacity (Θ fc ). Poorly predicted indicators (R 2 < 0.50; RPD < 1.40) were EC, S, P, available water capacity (AWC), K, Zn, and penetration resistance. Combining VNIRS with selected field- and laboratory-measured SQ indicator values increased predictability. Furthermore, VNIRS showed moderate to substantial agreement in predicting interpretive SQ scores and a composite soil quality index (CSQI) especially when combined with directly measured SQ indicator values. In conclusion, VNIRS has good potential for low cost, rapid assessment of physical and biological SQ indicators but conventional soil chemical tests may need to be retained to provide comprehensive SQ assessments.

Western Kenya

Can the impact of canopy trees on soil and understory be altered using litter additions?

Trees can have large effects on soil nutrients in ways that alter succession, particularly in the case of nitrogen-(N)-fixing trees. In Hawaiʻi, forest restoration relies heavily on use of a native N-fixing tree, Acacia koa (koa), but this species increases soil-available N and likely facilitates competitive dominance of exotic pasture grasses. In contrast, Metrosideros polymorpha (‘ōhi‘a), the dominant native tree in Hawaiʻi, is less often planted because it is slow growing; yet it is typically associated with lower soil N and grass biomass, and greater native understory recruitment. We experimentally tested whether it is possible to reverse high soil N under koa by adding ‘ōhi‘a litter, using additions of koa litter or no litter as controls, over 2.5 yr. We then quantified natural litterfall and decomposition rates of ‘ōhi‘a and koa litter to place litter additions in perspective. Finally, we quantified whether litter additions altered grass biomass and if this had effects on native outplants. Adding ‘ōhi‘a litter increased soil carbon, but increased rather than decreased inorganic soil N pools. Contrary to expectations, koa litter decomposed more slowly than ‘ōhi‘a, although it released more N per unit of litter. We saw no reduction in grass biomass due to ‘ōhi‘a litter addition, and no change in native outplanted understory survival or growth. We conclude that the high N soil conditions under koa are difficult to reverse. However, we also found that outplanted native woody species were able to decrease exotic grass biomass over time, regardless of the litter environment, making this a better strategy for lowering exotic species impacts.

Ecological Applications

Fire modifies plant–soil feedbacks

Although plant–soil feedbacks (interactions between plants and soils, often mediated by soil microbes, abbreviated as PSFs) are widely known to influence patterns of plant diversity at local and landscape scales, these interactions are rarely examined in the context of important environmental factors. Resolving the roles of environmental factors is important because the environmental context may alter PSF patterns by modifying the strength or even direction of PSFs for certain species. One important environmental factor that is increasing in scale and frequency with climate change is fire, though the influence of fire on PSFs remains essentially unexamined. By changing microbial community composition, fire may alter the microbes available to colonize the roots of plants and thus seedling growth post-fire. This has potential to change the strength and/or direction of PSFs, depending on how such changes in microbial community composition occur and the plant species with which the microbes interact. We examined how a recent fire altered PSFs of two leguminous, nitrogen-fixing tree species in Hawaiʻi. For both species, growing in conspecific soil resulted in higher plant performance (as measured by biomass production) than growing in heterospecific soil. This pattern was mediated by nodule formation, an important process for growth for legume species. Fire weakened PSFs for these species and therefore pairwise PSFs, which were significant in unburned soils, but were nonsignificant in burned soils. Theory suggests that positive PSFs such as those found in unburned sites would reinforce the dominance of species where they are locally dominant. The change in pairwise PSFs with burn status shows PSF-mediated dominance might diminish after fire. Our results demonstrate that fire can modify PSFs by weakening the legume-rhizobia symbiosis, which may alter local competitive dynamics between two canopy dominant tree species. These findings illustrate the importance of considering environmental context when evaluating the role of PSFs for plants.

Ecology

Wetland and hydric soils

Soil and the inherent biogeochemical processes in wetlands contrast starkly with those in upland forests and rangelands. The differences stem from extended periods of anoxia, or the lack of oxygen in the soil, that characterize wetland soils; in contrast, upland soils are nearly always oxic. As a result, wetland soil biogeochemistry is characterized by anaerobic processes, and wetland vegetation exhibits specific adaptations to grow under these conditions. However, many wetlands may also have periods during the year where the soils are unsaturated and aerated. This fluctuation between aerated and nonaerated soil conditions, along with the specialized vegetation, gives rise to a wide variety of highly valued ecosystem services.

Book chapter

Responses of amphibian populations to water and soil factors in experimentally-treated aquatic macrocosms

Survival of anuran embryos and tadpoles is reduced in acidic (pH<5.0) waters under laboratory conditions. However, field data on the presence-absence of amphibian species and acidity are equivocal. This study attempts to reconcile some of this discrepancy by using macrocosms to examine the interaction of soil type and water acidification on free-ranging tadpole populations. Tadpoles were caught with activity traps in 24 aquatic macrocosms experimentally treated with H 2 SO 4 and Al 2 (SO 4 ) 3 and lined with either comparatively high metal, low organic matter clay soils or lower metal, higher organic matter loams. Northern cricket frog ( Acris crepitans ) tadpole abundance was less in acidified macrocosms than in circumneutral ones (p<0.05) and less in those with loam soils than in macrocosms with clay soils (p<0.04). Gray treefrog ( Hyla versicolor ) abundance was affected by an interaction between soil and acidification (p<0.07) in that treatment effects were only observed in macrocosms with clay soils (p<0.01). No differences were observed among treatments for green frog ( Rana clamitans ) or southern leopard frog ( R. utricularia ) tadpoles. The study shows that soil type may interact with water conditions to affect amphibian populations in acidified waters.

Archives of Environmental Contamination and Toxico

Survey and evaluation of contaminants in earthworms and in soils derived from dredged material at confined disposal facilities in the Great Lakes region

Soils derived from dredged material were collected, together with earthworms from nine confined disposal facilities located in the Great Lakes Region. These samples were analyzed for 18 elements, 11 organochlorine pesticides, PCBs, and 24 polycyclic aromatic hydrocarbons. The concentrations detected in earthworms were evaluated in terms of their potential hazard to wildlife, which for the sake of the evaluation were assumed to prey entirely either on earthworms or on other soil invertebrates having similar concentrations. The soil concentrations (dry wt.) of the contaminants of greatest concern were <1.9 to 32 ppm Cd, <0.053 to 0.94 ppm Hg, 4.6 to 550 ppm Pb, and <0.1 to 1.0 ppm PCBs. The concentrations in earthworms (dry wt., ingested soil included) were as high as 91 ppm Cd, 1.6 ppm Hg, 200 ppm Pb, and 1.8 ppm PCBs. Based on laboratory toxicity studies of relatively sensitive species, and on concentration factors calculated from the earthworm and soil data, we estimated that lethal or serious sublethal effects on wildlife might be expected at concentrations of 10 ppm Cd, 3 ppm Hg, 670 ppm Pb, and 1.7 ppm PCBs in alkaline surface soils derived from dredged material. Concentrations of polycyclic aromatic hydrocarbons in earthworms were well below those in soil.

Environmental Monitoring and Assessment

Regional investigations of soil and overburden analysis and plant uptake of metals

Regional studies on the bioavailability of metals at native and disturbed sites were conducted over the past seven years by the USGS. The work was concentrated in the Fort Union, Powder River, and Green River coal resource regions where measures of extractable metals in soils were found to have limited use in predicting metal levels in plants. Correlations between Cu, Fe, and Zn in plants and extractable (DTPA, EDTA, and oxalate) or total levels in native A- and C-horizons of soil were occasionally significant. A simple linear model is generally not adequate, however, in estimating element uptake by plants. Prediction capabilities were improved when a number of soil chemical and physical parameters were included as independent variables in a stepwise linear multiple regression analysis; however, never more than 54% of the total variability in the data was explained by the equations for these metals. Soil pH was the most important variable relating soil chemistry to plant chemistry. This relation was always positive and apparently a response to soil levels of metal carbonates and not Fe and Mn oxides. Studies that compared the metal uptake by rehabilitation species to extractable (DTPA) metal levels in mice soils produced similar results. ?? 1984 Science and Technology Letters.

Minerals and the Environment

Soil base saturation combines with Beech Bark Disease to influence composition and structure of Sugar Maple-Beech forests in an acid rain-impacted region

Sugar maple, an abundant and highly valued tree species in eastern North America, has experienced decline from soil calcium (Ca) depletion by acidic deposition, while beech, which often coexists with sugar maple, has been afflicted with beech bark disease (BBD) over the same period. To investigate how variations in soil base saturation combine with effects of BBD in influencing stand composition and structure, measurements of soils, canopy, subcanopy, and seedlings were taken in 21 watersheds in the Adirondack region of NY (USA), where sugar maple and beech were the predominant canopy species and base saturation of the upper B horizon ranged from 4.4 to 67%. The base saturation value corresponding to the threshold for Al mobilization (16.8%) helped to define the species composition of canopy trees and seedlings. Canopy vigor and diameter at breast height (DBH) were positively correlated (P < 0.05) with base saturation for sugar maple, but unrelated for beech. However, beech occupied lower canopy positions than sugar maple, and as base saturation increased, the average canopy position of beech decreased relative to sugar maple (P < 0.10). In low-base saturation soils, soil-Ca depletion and BBD may have created opportunities for gap-exploiting species such as red maple and black cherry, whereas in high-base saturation soils, sugar maple dominated the canopy. Where soils were beginning to recover from acidic deposition effects, sugar maple DBH and basal area increased progressively from 2000 to 2015, whereas for beech, average DBH did not change and basal area did not increase after 2010.

New York

Forest vegetation change and its impacts on soil water following 47 years of managed wildfire

Managed wildfire is an increasingly relevant management option to restore variability in vegetation structure within fire-suppressed montane forests in western North America. Managed wildfire often reduces tree cover and density, potentially leading to increases in soil moisture availability, water storage in soils and groundwater, and streamflow. However, the potential hydrologic impacts of managed wildfire in montane watersheds remain uncertain and are likely context-dependent. Here we characterize the response of vegetation and soil moisture to 47 years (1971-2018) of managed wildfire in Sugarloaf Creek Basin (SCB) in Sequoia-Kings Canyon National Park in the Sierra Nevada, California, USA, using repeat plot-measurements, remote-sensing of vegetation, and a combination of continuous in-situ and episodic spatially-distributed soil moisture measurements. We find that, by comparison to a nearby watershed with higher vegetation productivity and greater fire frequency, the managed wildfire regime at SCB caused relatively little change in dominant vegetation over the 47 year period, and relatively little response of soil moisture. Fire occurrence was limited to drier mixed-conifer sites; fire-caused overstory tree mortality patches were generally < 10 ha, and fires had little effect on removing mid- and lower strata trees. Few dense meadow areas were created by fire, with most forest conversion leading to sparse meadow and shrub areas, which had similar soil moisture profiles to nearby mixed-conifer vegetation. Future fires in SCB could be managed to encourage greater tree mortality adjacent to wetlands to increase soil moisture, although the potential hydrologic benefits of the program in drier basins such as this one may be limited.

California

Relationships between soil macroinvertebrates and nonnative feral pigs (Sus scrofa) in Hawaiian tropical montane wet forests

Nonnative feral pigs (Sus scrofa) are recognized throughout the New World as a highly significant introduced species in terms of ecosystem alteration. Similarly, nonnative soil macroinvertebrates (e.g. earthworms, ground beetles) invade and alter the structure and function of native habitats globally. However, the relationship between feral pigs and soil macroinvertebrates remains largely unknown. This study analyzed relationships between these taxa using nine sites located inside and outside of feral pig management units representing a ~ 25 year chronosequence of removal in tropical montane wet forests in Hawai‘i. Soil macroinvertebrates were sampled from plots categorized as: actively trampled by feral pigs, actively rooted by feral pigs, feral pigs present with no signs of recent activity, or feral pigs removed over time. In total, we found 13 families of primarily nonnative soil macroinvertebrates. Plots with active trampling correlated with lower total macroinvertebrate abundance, biomass, and family richness. Plots with active rooting were correlated with higher abundance of nonnative earthworms (Lumbricidae and Megascolicidae) and ground beetles (Carabidae). The abundance, biomass, and biodiversity of macroinvertebrates did not vary with time since feral pig removal. Collectively, these results indicate: (1) trampling by feral pigs negatively influences soil macroinvertebrates; (2) feral pigs either modify habitats while rooting thereby facilitating earthworm and ground beetle habitat use or selectively seek out target prey species of soil macroinvertebrates; and (3) removal of feral pigs has minimal impacts on soil macroinvertebrates over time. These results are important globally due to the broadly overlapping ranges of S. scrofa and nonnative macroinvertebrates.

Hawaii

The concurrent use of novel soil surface microclimate measurements to evaluate CO2 pulses in biocrusted interspaces in a cool desert ecosystem

Carbon cycling associated with biological soil crusts, which occupy interspaces between vascular plants in drylands globally, may be an important part of the coupled climate-carbon cycle of the Earth system. A major challenge to understanding CO 2 fluxes in these systems is that much of the biotic and biogeochemical activity occurs in the upper few mm of the soil surface layer (i.e., the ‘mantle of fertility’), which exhibits highly dynamic and difficult to measure temperature and moisture fluctuations. Here, we report a multi-sensor approach to simultaneously measuring temperature and moisture of this biocrust surface layer (0–2 mm), and the deeper soil profile, concurrent with automated measurement of surface soil CO 2 effluxes. Our results illuminate robust relationships between biocrust water content and field CO 2 pulses that have previously been difficult to detect and explain. All observed CO 2 pulses over the measurement period corresponded to surface wetting events, including when the wetting events did not penetrate into the soil below the biocrust layer (0–2 mm). The variability of temperature and moisture of the biocrust surface layer was much greater than even in the 0–5 cm layer of the soil beneath the biocrust, or deeper in the soil profile. We therefore suggest that coupling surface measurements of biocrust moisture and temperature to automated CO 2 flux measurements may greatly improve our understanding of the climatic sensitivity of carbon cycling in biocrusted interspaces in our study region, and that this method may be globally relevant and applicable.

Biogeochemistry

Elevated temperature and nutrients lead to increased N2O emissions from salt marsh soils from cold and warm climates

Salt marshes can attenuate nutrient pollution and store large amounts of ‘blue carbon’ in their soils, however, the value of sequestered carbon may be partially offset by nitrous oxide (N 2 O) emissions. Global climate and land use changes result in higher temperatures and inputs of reactive nitrogen (Nr) into coastal zones. Here, we investigated the combined effects of elevated temperature (ambient + 5℃) and Nr (double ambient concentrations) on nitrogen processing in marsh soils from two climatic regions (Quebec, Canada and Louisiana, U.S.) with two vegetation types, Sporobolus alterniflorus (= Spartina alterniflora ) and Sporobolus pumilus (= Spartina patens ), using 24-h laboratory incubation experiments. Potential N 2 O fluxes increased from minor sinks to major sources following elevated treatments across all four marsh sites. One day of potential N 2 O emissions under elevated treatments (representing either long-term sea surface warming or short-term ocean heatwaves effects on coastal marsh soil temperatures alongside pulses of N loading) offset 15–60% of the potential annual ambient N 2 O sink, depending on marsh site and vegetation type. Rates of potential denitrification were generally higher in high latitude than in low latitude marsh soils under ambient treatments, with low ratios of N 2 O:N 2 indicating complete denitrification in high latitude marsh soils. Under elevated temperature and Nr treatments, potential denitrification was lower in high latitude soil but higher in low latitude soil as compared to ambient conditions, with incomplete denitrification observed except in Louisiana S. pumilus . Overall, our findings suggest that a combined increase in temperature and Nr has the potential to reduce salt marsh greenhouse gas (GHG) sinks under future global change scenarios.

Biogeochemistry

The influence of soil resources and plant traits on invasion and restoration in a subtropical woodland

It has been shown in some cases that nitrogen (N) addition to soil will increase abundance of plant invaders because many invaders have traits that promote rapid growth in response to high resource supply. Similarly, it has been suggested, and sometimes shown, that decreasing soil N via carbon (C) additions can facilitate native species recovery. Yet all species are unlikely to respond to resource supply in the same way. We asked how soil nutrients and competition affect native and exotic woody species in a restoration experiment where we added N or C, and crossed soil manipulation with the manipulation of dominant exotic grass abundance in a Hawaiian subtropical woodland. We related changes in survival and growth of outplanted individuals to native/exotic status and plant traits. As a group, N-fixers showed reduced survival compared to non-fixers in response to added N, with Morella faya (exotic) and Acacia koa (native) having dramatic negative responses. Among non-fixers, species with greater foliar %N had more positive survival responses to increasing soil N. Specific leaf area was not predictive of responses to nutrients or competition. In general, responses to carbon addition were weak, although reducing competition from existing exotic grasses was beneficial for all outplanted species, with N-fixers showing the most positive response. We conclude that commonly used restoration strategies to clear exotic species or lower soil resources with C addition will most greatly benefit N-fixing species, which themselves may be unwanted invaders. Thus statements about the influence of increased soil N on invasions should be carefully dissected by considering the traits (such as N-fixation status) of the regional species pool.

Hawai'i

Influence of soil microbiota on Taxodium distichum seedling performance during extreme flooding events

Plant associations with soil microbiota can modulate tree seedling growth and survival via mutualistic or antagonistic interactions. It is uncertain, however, whether soil microbiota influence seedling growth of coastal trees when exposed to extreme flooding regimes. We evaluated the role of soil microbes in promoting baldcypress ( Taxodium distichum ) seedling performance under different inundation scenarios and determined the influence of flooding on the colonization of in planta beneficial microbes. Seedlings reared in sterile and non-sterile soil were exposed to three different flooding regimes historically experienced in Louisiana swamps. Seedling growth was assessed, and the colonization by beneficial symbionts such as arbuscular mycorrizhal fungi (AMF), and dark septate endophytes (DSE) was evaluated in harvested roots. Seedlings grown in sterile soil had six times higher growth than seedlings reared in non-sterile soil. As a result, we evaluated pathogen load in the roots by assessing oomycete colonization. Flooding influenced the in planta colonization of DSE and oomycetes, but did not affect the colonization of mutualist AMF fungi. DSE and oomycetes were rarer in flooded conditions, while AMF remained abundant. Seedling biomass production was not correlated with in planta fungal colonization or pathogen load. Soil microbiota can negatively influence baldcypress seedling growth, and no growth benefit was evidenced from the root colonization of mutualist fungi. Flooding can modify baldcypress-fungal interactions by diminishing colonization of DSE. Overall, baldycpress seedlings were more sensitive to the presence of microbiota than flooding, and thus restoration efforts should focus on having a better understanding of plant–microbe interactions in swamps.

Louisiana

Moisture and temperature controls on nitrification differ among ammonia oxidizer communities from three alpine soil habitats

Climate change is altering the timing and magnitude of biogeochemical fluxes in many high elevation ecosystems. The consequent changes in alpine nitrification rates have the potential to influence ecosystem scale responses. In order to better understand how changing temperature and moisture conditions may influence ammonia oxidizers and nitrification activity, we conducted laboratory incubations on soils collected in a Colorado watershed from three alpine habitats (glacial outwash, talus, and meadow). We found that bacteria, not archaea, dominated all ammonia oxidizer communities. Nitrification increased with moisture in all soils and under all temperature treatments. However, temperature was not correlated with nitrification rates in all soils. Site-specific temperature trends suggest the development of generalist ammonia oxidizer communities in soils with greater in situ temperature fluctuations and specialists in soils with more steady temperature regimes. Rapidly increasing temperatures and changing soil moisture conditions could explain recent observations of increased nitrate production in some alpine soils.

Colorado