USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Data”

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.

939 records · Page 29Linked to original sources
USGS2026

Nonstationary demographic state-space models using unreplicated counts for species undergoing environmental stressors

A fundamental task in ecological statistics is to estimate abundance and growth rate distributions from wildlife monitoring data to inform conservation management. Modeling time series of wildlife populations presents a number of challenges from both statistical and ecological perspectives, including discreteness; lack of replication; nonstationarity; and observation, demographic, and other phenomenological processes. Nonstationary dynamics are often exhibited by populations undergoing environmental stressors. Models must account for these characteristics to produce reliable estimates of abundance and trends, yet estimation can be challenging with unreplicated data. We propose nonstationary demographic state-space models using unreplicated counts for populations undergoing environmental stressors. A reduced growth rate model matches the complexity of the unreplicated count data, and a fecundity bound on growth rate distributions allows the separation of processes affecting growth rates like environmental stressors from those affecting abundance external to growth rates like migration. NDSSMs allow for the embedding of nonstationary model components, and we explore the use of changepoints, volatility clustering, and migration processes. We apply the proposed nonstationary models in case studies of herons affected by predator/competitor reestablishment and three bat species affected by a fungal pathogen causing white-nose syndrome. Nonstationary models outperform stationary models and generalized linear mixed effects models according to model scoring and visual inspection of predictions, and provide estimates more consistent with published values. Incorporating migration improves model fit universally, even with approximate one-way immigration, most likely because populations are extirpated, recolonized, and increase multiple-fold over the upper bound set by species fecundity. In addition, estimates of the timing and severity of the environmental stressor differed for models with migration. Including nonstationary and demographic components in a fecundity-bounded growth rate model improves inference and benefits interpretability of hyperparameters. In turn, this adjusts uncertainties in predictions of abundance and growth rates over time, providing the ingredients needed for informed conservation analysis and for directing future monitoring of at-risk species.

Journal of Agricultural, Biological and Environmen↗
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2025

Multi-scale geophysical mapping of the brine and bedrock surfaces along the Dolores River, Paradox Valley, Colorado, December 2023

Total dissolved solids derived from salt dome–sourced brine in the underlying alluvial aquifer substantially increase with distance in the reach of the Dolores River that passes through Paradox Valley in southwestern Colorado. The area has been the site of salinity control operations since the 1990s to reduce salt loading to the downstream Colorado River. Previous airborne and ground/water-based electromagnetic (EM) geophysical data have successfully mapped the top of the brine surface, albeit with relatively coarse near-surface resolution and limited spatial coverage. This present December 2023 study used ground-based high-resolution EM and passive seismic (horizontal-to-vertical spectral ratio, HVSR) tools to map in detail the depth and thickness of the brine zone in the alluvial aquifer (top of the brine down to bedrock contact) in areas immediately surrounding the Dolores River where previous airborne EM (AEM) results indicated brine within 10 m of land surface. Results indicate the deepest bedrock is generally associated with the shallowest brine and local depressions in the collapse breccia (caprock to the Paradox Formation salt) may facilitate vertical migration of brine into the alluvial aquifer. Additionally, the ground-based EM mapping corroborated general patterns in depth to brine that were observed in previous AEM results while also revealing additional detail, including suspected focused brine discharge zones to the Dolores River. A river-based survey wherein EM data, channel depth and river water–specific conductance information were collected augmented these findings and indicated specific areas on both the western and eastern sides of the river where focused brine discharge may occur. This study comprises a large-scale, ground- and water-based geophysical mapping effort, including hundreds of HVSR soundings and 100s of kilometres of EM data, which were successfully translated into digital brine and bedrock surfaces that could be incorporated into groundwater modelling efforts, future well siting or other decision-making.

Colorado↗
USGS2025

Airborne geophysical efforts for critical mineral systems mapping in the southern Midcontinent, USA

The increasing global demand for critical minerals to support energy and technological advancement has accelerated exploration and research efforts for these essential resources. Since 2019, the United States Geological Survey (USGS) Earth Mapping Resources Initiative (EMRI) has worked to modernize geologic mapping of the Nation to better understand its critical mineral resources. To further this initiative, the USGS has flown a series of high-resolution airborne magnetic and radiometric surveys over large areas of the southern Midcontinent. The surveys cover known critical mineral deposits and areas with the potential to host additional critical minerals based on the presence of one or more overlapping mineral systems. One aspect of EMRI emphasizes close collaboration between the USGS and the Association of American State Geologists, as well as other government and industry partners to leverage geophysical, geological, and geochemical expertise on both regional and local scales. The EMRI high-resolution airborne survey data provide new insights into the geophysical framework of the southern Midcontinent and its critical mineral endowment. Additionally, discoveries made from the data have directed new studies for critical mineral exploration.

southern Midcontinent↗
USGS2025

The influence of nodule versus crust morphology on the composition of seamount-hosted ferromanganese minerals

Ferromanganese mineral precipitation in the global oceans is ubiquitous, occurring in the form of both crusts and nodules at a broad range of depths and seafloor terrains. Although ferromanganese crusts and nodules are both composed of ferromanganese minerals, mineralogy and mean element concentrations compiled for regional crust versus nodule occurrences differ. Notably, most published compilations compare nodules from abyssal plain sediments to crusts forming on seamounts, and do not address the question as to whether mineral morphology may affect the composition of crusts and nodules in addition to environmental factors. To address this, we have compiled a dataset of co-located hydrogenetic (seawater sourced) open-ocean ferromanganese crusts and seamount-hosted nodules. The preliminary data presented here support our hypothesis that co-located crusts and nodules exhibit similar compositions, both in terms of their mineralogy and mean element concentrations, as well as overall variance in concentrations for any given element. These data indicate that local processes influence sample-to-sample variation in both crusts and nodules, in contrast to broader ferromanganese variation reported between open-ocean seamounts- hosted crusts and abyssal plain-hosted nodules.

Conference Paper↗
USGS2026

Multiple machine-learning estimation of groundwater levels and trends for the regional Mississippi River Valley alluvial aquifer

The Mississippi River Valley alluvial aquifer provides irrigation, public, and domestic water supplies across the south-central United States. Declining groundwater levels require improved characterization of changing conditions. Traditional potentiometric-surface mapping does not use all available water-level data or quantify uncertainty. To address these limitations, we developed a data-driven multiple machine-learning (MML) framework delivered through two open-source R packages. The covMRVAgen1 software assembles covariates to 155,960 monthly groundwater levels from 57,695 wells; the mmlMRVAgen1 software trains Cubist and Random Forest models, blends them, and makes 1-kilometer gridded predictions of monthly potentiometric surfaces for the period January 1980–December 2022. The MML approach provides a methodological foundation for region-scale spatiotemporal groundwater prediction and uncertainty quantification, generating 90-percent prediction limits with appropriate empirical coverage. Model performance is acceptable, with a root-mean-square error of about 4.2 feet, standard deviation of 24.82 feet, and a normalized Nash–Sutcliffe efficiency of 0.973.

Arkansas, Illinois, Louisiana, Mississippi, Missou↗
USGS2025

Rebuilding a volcano one lava flow at a time—Visualizing the lava dome-building eruption in the crater of Mount St. Helens, 1982–1986

Between 1980 and 1986, the U.S. Geological Survey made a series of 1:2,000-scale topographic contour maps from aerial photographic surveys to monitor the eruption. These maps were made for operational purposes and were not intended for publication. Since then, advances in technology made it possible to digitize the original, highly detailed hardcopy maps and derive new digital data elevation models of the surface of the lava dome. These digital elevation models allow for the visualization of the progression of the eruption and reveal the rubbly, chaotic surface of the lava flows and dome. Additionally, these new data help fill gaps in the long-term record of topographic changes that have occurred at the volcano since the May 18, 1980, eruption.

Washington↗
USGS2025

Regional characterization of coal resources in the U.S. Gulf Coast

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The U.S. Gulf Coast and the Williston basins are the two major lignite-bearing basins within the country. Recent REE and CM studies of the lignite in these basins have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of the U.S. Gulf Coast lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from 31,181 drill holes and other data points. The estimated total coal resource in the Gulf Coast is about 83 billion metric tons in the upper 90 m (~ 300 ft) of the subsurface. Texas accounted for 40 percent (32 billion metric tons) of the total resource, followed by Mississippi (24 %, 20 billion metric tons), Louisiana (14 %, 12 billion metric tons), Tennessee (10 %, 8.5 billion metric tons), and Arkansas (6 %, 5.1 billion metric tons). The remaining states each accounted for less than 5 percent of the total resource. Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known state-wide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast and Williston Basin resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses. Coal-power electric generation has sharply decreased in past decades and mining of these coals for CM and REE could provide additional co-products such as activated carbon and other uses such as fertilizer (soil enhancer).

EarthArXiv↗
USGS2026

Standardized method for logging drill core at the Idaho National Laboratory, Idaho

The U.S. Geological Survey’s (USGS) Lithologic Core Storage Library (CSL) at the Idaho National Laboratory stores more than 120,000 feet of drill core that is accessible to the public for research and sampling. To effectively convey the physical and descriptive properties of the drill core, USGS staff at the Idaho National Laboratory Project Office log the drill core and publish the lithologic logs as data releases. The logs provide essential data on the lithology, texture, mineralogy, alteration, and other physical properties of the core, which serve as valuable information for researchers to guide their research and sampling efforts. To ensure consistent, quality, and dependable lithologic logs, this document outlines the procedures and expectations for logging drill core at the CSL. This document describes the processes for storing, photographing, and logging core, and includes a variety of resources, reference materials, and appendixes designed to standardize and aid the logging process. Following the procedures outlined in this document will produce consistent, detailed logs that facilitate dependable observations and serve as an easy reference for researchers and other interested parties.

Idaho↗
USGS2025

2024 Surprise Inlet landslides: Insights from a prototype landslide‐triggered tsunami monitoring system in Prince William Sound, Alaska

Alaska's coastal communities face growing landslide hazards owing to glacier retreat and extreme weather intensified by the warming climate, yet hazard monitoring remains challenging. As part of ongoing experimental monitoring in Prince William Sound, we detected three large landslides (0.5–2.3 M m 3 ) at Surprise Inlet on 20 September 2024, within the span of an hour. These events were identified in near real-time through seismic data and later confirmed using satellite imagery, tidal records, and infrasound. The landslides generated a modest tsunami, and a 4 cm wave was recorded by a tide gauge 18 km away, marking the first recorded landslide to reach water since monitoring began in this region in 2021. Here, we examine the detection and interpretation of these landslides using multiple data sources and modeling. We demonstrate the effectiveness of this regional seismic monitoring system and show how complementary instrumentation, where available, can enhance detection capabilities.

Alaska↗
USGS2026

Emergent headwater types for watershed monitoring, research, and planning

Headwaters—that is headwater streams and the watersheds they drain—strongly influence water quality, ecosystem services, and hydrologic connectivity, yet remain poorly mapped, monitored, and understood at large spatial scales. Using available spatial data to classify headwater “types” could therefore enhance the utility of sparse monitoring data to advance research and watershed planning—and help update literature-based perceptual models of headwaters as steep, forested, and minimally disturbed watersheds. Here, we analyze headwaters across the contiguous United States to quantify how their attributes (climate, topography, land cover, soils) vary regionally and differ from their downstream watersheds. Our analysis (a) identified seven dominant headwater types with distinct combinations of biophysical and climatic attributes and (b) revealed systematic contrasts with downstream watersheds, with strongest differences in attributes representing potential flowpath depth. Our findings offer a transferable framework for classifying, modeling, and managing these critical yet understudied components of the hydrologic network.

conterminous United States↗
USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗
USGS2024

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)↗
USGS2025

Near-surface material and topography generate anomalous high-frequency ground motion amplification in Chugiak, Alaska

An ∼3 km long nodal array oriented approximately east–west was deployed in Chugiak, Alaska, by the U.S. Geological Survey during 2021. The array intersects with the permanent NetQuakes station NP.ARTY, where peak ground acceleration (PGA) value of 1.98 g was recorded during the 2018 M w 7.1 Anchorage, Alaska, earthquake, in sharp contrast to the PGA of ∼0.3 g at a site just 4 km to the west. Seismic data for M w 1.8–4.3 aftershocks from the M w 7.1 event recorded by the nodal array confirm the anomalously large ground motions obtained at NP.ARTY as well as similar amplifications at nodes within ∼1 km to the east. Here, we performed 0–10 Hz 3D finite‐difference simulations, including high‐resolution surface topography, to explore the cause of the unexpectedly large amplification. As expected, the simulations computed with a regional 3D tomography velocity model severely underpredict the 0–10 Hz acceleration records at almost all sites. Adding a near‐surface low‐velocity taper to 300 m depth amplifies the accelerations by up to a factor of 5 and enables a reasonable match between the nodal data and simulations at sites to the west of NP.ARTY. However, this model still underpredicts the spectral energy in the area covered by glacial sediments by up to an order of magnitude. The addition of a till layer using a depth‐dependent shear‐wave velocity ( ⁠⁠ V s ) profile along with a homogeneous, 8 m thick low‐velocity layer with V s = 250 m/s representing the kame terraces improves the fit to data to within a factor of 2 at nodes located on top of the glacial sediments. Our study shows that the anomalously large high‐frequency amplification recorded at and near NP.ARTY can be explained by a combination of topographic effects and near‐surface low‐velocity material with amplification effects on the high‐frequency ground motion by up to about 40% and an order of magnitude, respectively.

Alaska↗
USGS2025

Mine waste as a potential source of critical minerals and other commodities: Examples from the Four Corners states, USA

The growing demand for critical minerals and other mineral resources has raised concerns about possible supplies of these essential commodities. Mine waste is a potential source of these essential commodities. We compiled a geospatial database of publicly available data of the largest non-fuel mine waste features (>200,000 m 2 in areal extent) in the Four Corners states of the United States, where most of those features are from porphyry Cu deposits. The combined volume of those large porphyry Cu mine waste features is approximately 17 cubic kilometers, 60% of which is ore-related material such as tailings. Using publicly available data on density, grade, and previous recovery values, we estimate the contained endowments of Ag, As, Au, Bi, Cu, Mo, Re, S, Sb, Se, Te, and W. These estimates indicate endowments within ore-related mine wastes are collectively comparable to endowments of giant to supergiant deposits. If fully recovered, these commodities could meet current global demand from less than 1 year (Sb) to more than 200 years (Re), underscoring the enormous untapped resource potential of mine waste.

Arizona, Colorado, New Mexico, Utah↗
USGS2024

A comparative analysis of OpenET for evaluating evapotranspiration in California almond orchards

The almond industry in California faces water management challenges that are being exacerbated by droughts, climate change, and groundwater sustainability legislation. The Tree-crop Remote sensing of Evapotranspiration eXperiment (T-REX) aims to explore opportunities to improve precision irrigation management for woody perennial cropping systems. Almond orchards in the California Central Valley were equipped with eddy covariance flux measurements to evaluate satellite remote sensing-based evapotranspiration (RSET) models. OpenET provides high-resolution (30-m spatial and daily temporal) RSET data, synthesizing decades of research for practical water management. This study provides an evaluation of OpenET performance at six almond sites covering a large range in soils, age, and variety. It also compares OpenET ensemble evapotranspiration (ET) data with applied irrigation and precipitation records over an additional 148 almond orchards located in the Central Valley of California. Results show OpenET models, including the ensemble ET value, produced reasonable and actionable ET values, with overall coefficient of determination (R 2 ) and mean absolute error values of 0.73- and 0.95-mm d −1 at the daily time step, respectively. However, given the temporal sampling of Landsat (8-day revisit) and the interpolation methods used, the assessed ET models had difficulty in capturing short-term variability in almond ET; for example, the rapid decline in measured ET observed as a response to lack of irrigation preceding and during almond harvest. The study also drew attention to the spatial complexity in scenarios where irrigated orchards are surrounded by hot/dry areas, causing discrepancies between measured and modeled ET values. In comparison with irrigation records, OpenET ensemble ET was capable of quantifying water input (applied irrigation + precipitation) in almond orchards to within 13 % when evaluating monthly data. Initial results presented here reinforce the idea that RSET models, such as in OpenET, are powerful tools, yet their application requires nuanced understanding and careful consideration of local conditions.

California↗
USGS2025

Geochemistry and radiogenic isotopes constrain the mantle source region of the Mountain Pass Intrusive Suite, California

The Mountain Pass carbonatite stock is the largest rare earth element (REE) deposit and only active REE mine in the United States. The carbonatite intrusion and spatially associated alkaline silicate intrusions constitute the Mountain Pass Intrusive Suite, which is located within the Mojave Province in California. Both the carbonatite and the alkaline silicate rocks are enriched in large ion lithophile elements and light REEs, and less enriched to depleted in high field strength elements, indicating the mantle source region was metasomatically enriched in incompatible trace elements. The cause of this metasomatic mantle enrichment and the genetic relationship between the carbonatite and the alkaline silicate stocks are poorly understood. In this study, major and trace element geochemical data and isotopic (Rb-Sr, Sm-Nd, and Lu-Hf) data are presented to constrain genesis of the Mountain Pass Intrusive Suite, from mantle source region to the intrusion of the stocks. Our geochemical data are consistent with derivation of the alkaline silicate and carbonatite melts through partial melting from a shared mantle source region rather than through liquid immiscibility or fractional crystallization and separation of a carbothermal fluid. Although the Rb-Sr isotopic system in the Mountain Pass Intrusive Suite is disturbed at the whole-rock scale, the isotopic systems for whole-rock Sm-Nd (εNd i = ‐2.2 ± 0.8) and zircon Lu-Hf (εHf i = 0.1 ± 1.1) are robust and support mantle derivation of the magmas. Geochemical modeling using experimentally derived partition coefficients was used to identify possible causes of enrichment in incompatible elements through metasomatism in the mantle source region. Modeling of metasomatism by melts derived by partial melting of deeply subducted carbonated sediments approximates observed Mountain Pass Intrusive Suite trace element chemistry. Scattered εHf i in inherited zircon (2.8 ± 2.6) is consistent with derivation from an arc-related environment with substantial crustal contamination. Paleotectonic studies in the Mojave Province indicate that regional subduction preceded emplacement of the Mountain Pass Intrusive Suite by ∼300 Ma. Melting of the Mountain Pass source region may have been caused by post-collisional thermal relaxation and extension.

California↗
USGS1999

Estimation of magnitude and frequency of floods for streams in Puerto Rico: New empirical models

Flood-peak discharges and frequencies are presented for 57 gaged sites in Puerto Rico for recurrence intervals ranging from 2 to 500 years. The log-Pearson Type III distribution, the methodology recommended by the United States Interagency Committee on Water Data, was used to determine the magnitude and frequency of floods at the gaged sites having 10 to 43 years of record. A technique is presented for estimating flood-peak discharges at recurrence intervals ranging from 2 to 500 years for unregulated streams in Puerto Rico with contributing drainage areas ranging from 0.83 to 208 square miles. Loglinear multiple regression analyses, using climatic and basin characteristics and peak-discharge data from the 57 gaged sites, were used to construct regression equations to transfer the magnitude and frequency information from gaged to ungaged sites. The equations have contributing drainage area, depth-to-rock, and mean annual rainfall as the basin and climatic characteristics in estimating flood peak discharges. Examples are given to show a step-by-step procedure in calculating a 100-year flood at a gaged site, an ungaged site, a site near a gaged location, and a site between two gaged sites.

Puerto Rico↗