Multispectral remote sensing data acquired by the Landsat 8 Operational Land Imager (OLI) sensor were analyzed using a new, automated technique to generate a map of exposed mineral and vegetation groups in the western San Juan Mountains, Colo., and the Four Corners Region of the United States. Band ratio results were combined into displayed mineral and vegetation groups using Boolean algebra. New analysis logic has been implemented to exploit the coastal aerosol band in Landsat 8 OLI data and identify concentrations of iron sulfate minerals. These may indicate the presence of near-surface pyrite, which can be a potential nonpoint source of acid rock drainage. Zoned occurrences of iron sulfate minerals in areas surrounding and down gradient of known sources of pyrite have been mapped using this technique. They show high correlation with occurrences of jarosite-bearing mineral assemblages, as mapped using airborne imaging spectrometer data and supporting field verification surveys. Mapping the occurrence of iron sulfate minerals produced by the weathering and oxidation of pyrite in exposed hydrothermally altered rocks can also provide an important indicator of the genetic environment of alteration and the associated mineral deposit type. The automated analysis methodology is being employed to rapidly and cost-effectively generate maps of large regions of the United States in support of U.S. Geological Survey mineral resource and mineral-environmental assessments. This map, which includes the ERDAS IMAGINE (.img) thematic raster format in the data release, has been attributed by pixel value with mineral and vegetation group classification data, which can be queried in most image processing and GIS software packages.