SEARCH · USGS Science

Results for “Climate Dynamics”

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

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35Linked to original sources

USGS2022

Physiological and nutritional constraints on zooplankton productivity due to eutrophication and climate change predicted using a resource-based modeling approach

Emerging evidence suggests that zooplankton production is affected by physiological and nutritional constraints due to climate change and eutrophication, which in turn could have broad implications for food-web dynamics and fisheries production. In this study, we developed a resource-based zooplankton production dynamics model that causally links freshwater cladoceran and copepod daily production-to-biomass (P/B) ratios with water temperature, phytoplankton biomass and community composition, and zooplankton feeding selectivity. This model was used to evaluate constraints on zooplankton growth under four hypothetical scenarios: involving natural plankton community seasonal succession; lake fertilization to enhance fisheries production; eutrophication; and climatic warming. Our novel modeling approach predicts zooplankton production is strongly dependent on seasonal variation in resource availability and quality, which results in more complex zooplankton dynamics than predicted by simpler temperature dependent models. For mesotrophic and hypereutrophic lakes, our study suggests that the ultimate control over zooplankton P/B ratios shifts from physiological control during colder periods to strong resource control during warmer periods. Our resource-based model provided important insights into the nature of biophysical control of zooplankton under a changing climate that has crucial implications for food web energy transfer and fisheries production.

Canadian Journal of Fisheries and Aquatic Sciences↗

USGS2023

Comparison of ventifact orientations and recent wind direction indicators on the floor of Jezero crater, Mars

Wind-abraded rocks and aeolian bedforms have been observed at the Mars 2020 Perseverance landing site, providing evidence for recent and older wind directions. This study reports orientations of aeolian features measured in Perseverance images to infer formative wind directions. It compares these measurements with orbital observations, climate model predictions, and wind data acquired by the Mars Environmental Dynamics Analyzer. Three-dimensional orientations of flute textures on rocks, regolith wind tails extending from behind obstacles, and other aeolian features were measured using Digital Terrain Models (DTMs) derived from Mastcam-Z and navigation camera (Navcam) stereo images. Orientations of rock flutes measured in images acquired through Sol (martian day) 400 yielded a mean azimuth of 94° ± 7° (wind from the west). However, similar measurements of regolith wind tails indicate that recent sand-driving winds have been blowing from the east-southeast, nearly the opposite direction (mean azimuth = 285° ± 15°). Atmospheric modeling generally predicts net annual sand transport from the east-southeast at present, consistent with Perseverance regolith wind tail and orbital observations. The orientation of ventifact flutes thus suggests that they were formed under a different climate regime. Differences in orientations of recent and paleo-wind indicators have been noted at other Mars landing sites and may result from major orbital/axial changes that can cause significant changes in atmospheric circulation. Orientation differences between modern and older wind direction indicators at Jezero are useful clues to the climate history of the region.

JGR Planets↗

USGS2024

The noise is the signal: Spatio-temporal variability of production and productivity in high elevation meadows in the Sierra Nevada mountain range of North America

There are expectations that increasing temperatures will lead to significant changes in structure and function of montane meadows, including greater water stress on vegetation and lowered vegetation production and productivity. We evaluated spatio-temporal dynamics in production and productivity in meadows within the Sierra Nevada mountain range of North America by: (1) compiling Landsat satellite data for the Normalized Difference Vegetation Index (NDVI) across a 37-year period (1985–2021) for 8,095 meadows >2,500 m elevation; then, (2) used state-space models, changepoint analysis, geographically-weighted regression (GWR), and distance-decay analysis (DDA) to: (a) identify meadows with decreasing, increasing or no trends for NDVI; (b) detect meadows with abrupt changes (changepoints) in NDVI; and (c) evaluate variation along gradients of latitude, longitude, and elevation for eight indices of temporal dynamics in annual production (mean growing season NDVI; MGS) and productivity (rate of spring greenup; RSP). Meadows with no long-term change or evidence of increasing NDVI were 2.6x more frequent as those with decreasing NDVI (72% vs. 28%). Abrupt changes in NDVI were detected in 48% of the meadows; they occurred in every year of the study and with no indication that their frequency had changed over time. The intermixing of meadows with different temporal dynamics was a consistent pattern for monthly NDVI and, especially, the eight annual indices of MGS and RSP. The DDA showed temporal dynamics in pairs of meadow within a few 100 m of each other were often as different as those hundreds of kilometers apart. Our findings point strongly toward a great diversity of temporal dynamics in meadow production and productivity in the SNV. The heterogeneity in spatial patterns indicated that production and productivity of meadow vegetation is being driven by interplay among climatic, physiographic and biotic factors at basin and meadow scales. Thus, when evaluating spatio-temporal dynamics in condition for many high elevation meadow systems, what might often be considered “noise” may provide greater insight than a “signal” embedded within a large amount of variability.

California↗

USGS2022

Recent climate change has driven divergent hydrological shifts in high-latitude peatlands

High-latitude peatlands are changing rapidly in response to climate change, including permafrost thaw. Here, we reconstruct hydrological conditions since the seventeenth century using testate amoeba data from 103 high-latitude peat archives. We show that 54% of the peatlands have been drying and 32% have been wetting over this period, illustrating the complex ecohydrological dynamics of high latitude peatlands and their highly uncertain responses to a warming climate.

Nature Communications↗

USGS2009

Forecasting the combined effects of urbanization and climate change on stream ecosystems: from impacts to management options

Streams collect runoff, heat, and sediment from their watersheds, making them highly vulnerable to anthropogenic disturbances such as urbanization and climate change. Forecasting the effects of these disturbances using process-based models is critical to identifying the form and magnitude of likely impacts. Here, we integrate a new biotic model with four previously developed physical models (downscaled climate projections, stream hydrology, geomorphology, and water temperature) to predict how stream fish growth and reproduction will most probably respond to shifts in climate and urbanization over the next several decades. The biotic submodel couples dynamics in fish populations and habitat suitability to predict fish assemblage composition, based on readily available biotic information (preferences for habitat, temperature, and food, and characteristics of spawning) and day-to-day variability in stream conditions. We illustrate the model using Piedmont headwater streams in the Chesapeake Bay watershed of the USA, projecting ten scenarios: Baseline (low urbanization; no on-going construction; and present-day climate); one Urbanization scenario (higher impervious surface, lower forest cover, significant construction activity); four future climate change scenarios [Hadley CM3 and Parallel Climate Models under medium-high (A2) and medium-low (B2) emissions scenarios]; and the same four climate change scenarios plus Urbanization. Urbanization alone depressed growth or reproduction of 8 of 39 species, while climate change alone depressed 22 to 29 species. Almost every recreationally important species (i.e. trouts, basses, sunfishes) and six of the ten currently most common species were predicted to be significantly stressed. The combined effect of climate change and urbanization on adult growth was sometimes large compared to the effect of either stressor alone. Thus, the model predicts considerable change in fish assemblage composition, including loss of diversity. Synthesis and applications . The interaction of climate change and urban growth may entail significant reconfiguring of headwater streams, including a loss of ecosystem structure and services, which will be more costly than climate change alone. On local scales, stakeholders cannot control climate drivers but they can mitigate stream impacts via careful land use. Therefore, to conserve stream ecosystems, we recommend that proactive measures be taken to insure against species loss or severe population declines. Delays will inevitably exacerbate the impacts of both climate change and urbanization on headwater systems.

Journal of Applied Ecology↗

USGS1998

Drought-induced shift of a forest-woodland ecotone: Rapid landscape response to climate variation

In coming decades, global climate changes are expected to produce large shifts in vegetation distributions at unprecedented rates. These shifts are expected to be most rapid and extreme at ecotones, the boundaries between ecosystems, particularly those in semiarid landscapes. However, current models do not adequately provide for such rapid effects—particularly those caused by mortality—largely because of the lack of data from field studies. Here we report the most rapid landscape-scale shift of a woody ecotone ever documented: in northern New Mexico in the 1950s, the ecotone between semiarid ponderosa pine forest and piñon–juniper woodland shifted extensively (2 km or more) and rapidly (<5 years) through mortality of ponderosa pines in response to a severe drought. This shift has persisted for 40 years. Forest patches within the shift zone became much more fragmented, and soil erosion greatly accelerated. The rapidity and the complex dynamics of the persistent shift point to the need to represent more accurately these dynamics, especially the mortality factor, in assessments of the effects of climate change.

Proceedings of the National Academy of Sciences of↗

USGS2026

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing↗

USGS2023

Shifted sediment-transport regimes by climate change and amplified hydrological variability in cryosphere-fed rivers

Climate change affects cryosphere-fed rivers and alters seasonal sediment dynamics, affecting cyclical fluvial material supply and year-round water-food-energy provisions to downstream communities. Here, we demonstrate seasonal sediment-transport regime shifts from the 1960s to 2000s in four cryosphere-fed rivers characterized by glacial, nival, pluvial, and mixed regimes, respectively. Spring sees a shift toward pluvial-dominated sediment transport due to less snowmelt and more erosive rainfall. Summer is characterized by intensified glacier meltwater pulses and pluvial events that exceptionally increase sediment fluxes. Our study highlights that the increases in hydroclimatic extremes and cryosphere degradation lead to amplified variability in fluvial fluxes and higher summer sediment peaks, which can threaten downstream river infrastructure safety and ecosystems and worsen glacial/pluvial floods. We further offer a monthly-scale sediment-availability-transport model that can reproduce such regime shifts and thus help facilitate sustainable reservoir operation and river management in wider cryospheric regions under future climate and hydrological change.

Science Advances↗

USGS2016

Differential heating in the Indian Ocean differentially modulates precipitation in the Ganges and Brahmaputra basins

Indo-Pacific sea surface temperature dynamics play a prominent role in Asian summer monsoon variability. Two interactive climate modes of the Indo-Pacific—the El Niño/Southern Oscillation (ENSO) and the Indian Ocean dipole mode—modulate the amount of precipitation over India, in addition to precipitation over Africa, Indonesia, and Australia. However, this modulation is not spatially uniform. The precipitation in southern India is strongly forced by the Indian Ocean dipole mode and ENSO. In contrast, across northern India, encompassing the Ganges and Brahmaputra basins, the climate mode influence on precipitation is much less. Understanding the forcing of precipitation in these river basins is vital for food security and ecosystem services for over half a billion people. Using 28 years of remote sensing observations, we demonstrate that (i) the tropical west-east differential heating in the Indian Ocean influences the Ganges precipitation and (ii) the north-south differential heating in the Indian Ocean influences the Brahmaputra precipitation. The El Niño phase induces warming in the warm pool of the Indian Ocean and exerts more influence on Ganges precipitation than Brahmaputra precipitation. The analyses indicate that both the magnitude and position of the sea surface temperature anomalies in the Indian Ocean are important drivers for precipitation dynamics that can be effectively summarized using two new indices, one tuned for each basin. These new indices have the potential to aid forecasting of drought and flooding, to contextualize land cover and land use change, and to assess the regional impacts of climate change.

Remote Sensing↗

USGS1989

Velocities of antarctic outlet glaciers determined from sequential Landsat images

Approximately 91.0 percent of the volume of present-day glacier ice on Earth is in Antarctica; Greenland contains about another 8.3 percent of the volume. Thus, together, these two great ice sheets account for an estimated 99.3 percent of the total. Long-term changes in the volume of glacier ice on our planet are the result of global climate change. Because of the relationship of global ice volume to sea level (± 330 cubic kilometers of glacier ice equals ± 1 millimeter sea level), changes in the mass balance of the antarctic ice sheet are of particular importance. Whether the mass balance of the east and west antarctic ice sheets is positive or negative is not known. Estimates of mass input by total annual precipitation for the continent have been made from scattered meteorological observations (Swithinbank 1985). The magnitude of annual ablation of the ice sheet from calving of outlet glaciers and ice shelves is also not well known. Although the velocities of outlet glaciers can be determined from field measurements during the austral summer,the technique is costly, does not cover a complete annual cycle,and has been applied to just a few glaciers. To increase the number of outlet glaciers in Antarctica for which velocities have been determined and to provide additional data for under-standing the dynamics of the antarctic ice sheets and their response to global climate change, sequential Landsat image of several outlet glaciers were measured.

Antarctic Journal of the United States↗

USGS2012

Hotspot of accelerated sea-level rise on the Atlantic coast of North America

Climate warming does not force sea-level rise (SLR) at the same rate everywhere. Rather, there are spatial variations of SLR superimposed on a global average rise. These variations are forced by dynamic processes, arising from circulation and variations in temperature and/or salinity, and by static equilibrium processes, arising from mass redistributions changing gravity and the Earth's rotation and shape. These sea-level variations form unique spatial patterns, yet there are very few observations verifying predicted patterns or fingerprints. Here, we present evidence of recently accelerated SLR in a unique 1,000-km-long hotspot on the highly populated North American Atlantic coast north of Cape Hatteras and show that it is consistent with a modelled fingerprint of dynamic SLR. Between 1950–1979 and 1980–2009, SLR rate increases in this northeast hotspot were ~ 3–4 times higher than the global average. Modelled dynamic plus steric SLR by 2100 at New York City ranges with Intergovernmental Panel on Climate Change scenario from 36 to 51 cm (ref. 3); lower emission scenarios project 24–36 cm (ref. 7). Extrapolations from data herein range from 20 to 29 cm. SLR superimposed on storm surge, wave run-up and set-up will increase the vulnerability of coastal cities to flooding, and beaches and wetlands to deterioration.

Atlantic Coast, Cape Cod, Cape Hatteras, North Ame↗

USGS2018

Microspatial ecotone dynamics at a shifting range limit: plant–soil variation across salt marsh–mangrove interfaces

Ecotone dynamics and shifting range limits can be used to advance our understanding of the ecological implications of future range expansions in response to climate change. In the northern Gulf of Mexico, the salt marsh–mangrove ecotone is an area where range limits and ecotone dynamics can be studied in tandem as recent decreases in winter temperature extremes have allowed for mangrove expansion at the expense of salt marsh. In this study, we assessed aboveground and belowground plant–soil dynamics across the salt marsh–mangrove ecotone quantifying micro-spatial patterns in horizontal extent. Specifically, we studied vegetation and rooting dynamics of large and small trees, the impact of salt marshes (e.g. species and structure) on mangroves, and the influence of vegetation on soil properties along transects from underneath the mangrove canopy into the surrounding salt marsh. Vegetation and rooting dynamics differed in horizontal reach, and there was a positive relationship between mangrove tree height and rooting extent. We found that the horizontal expansion of mangrove roots into salt marsh extended up to eight meters beyond the aboveground boundary. Variation in vegetation structure and local hydrology appear to control mangrove seedling dynamics. Finally, soil carbon density and organic matter did not differ within locations across the salt marsh-mangrove interface. By studying aboveground and belowground variation across the ecotone, we can better predict the ecological effects of continued range expansion in response to climate change.

Oecologia↗

USGS2007

Forest dynamics in Oregon landscapes: Evaluation and application of an individual-based model

The FORCLIM model of forest dynamics was tested against field survey data for its ability to simulate basal area and composition of old forests across broad climatic gradients in western Oregon, USA. The model was also tested for its ability to capture successional trends in ecoregions of the west Cascade Range. It was then applied to simulate present and future (1990-2050) forest landscape dynamics of a watershed in the west Cascades. Various regimes of climate change and harvesting in the watershed were considered in the landscape application. The model was able to capture much of the variation in forest basal area and composition in western Oregon even though temperature and precipitation were the only inputs that were varied among simulated sites. The measured decline in total basal area from tall coastal forests eastward to interior steppe was matched by simulations. Changes in simulated forest dominants also approximated those in the actual data. Simulated abundances of a few minor species did not match actual abundances, however. Subsequent projections of climate change and harvest effects in a west Cascades landscape indicated no change in forest dominance as of 2050. Yet, climate-driven shifts in the distributions of some species were projected. The simulation of both stand-replacing and partial-stand disturbances across western Oregon improved agreement between simulated and actual data. Simulations with fire as an agent of partial disturbance suggested that frequent fires of low severity can alter forest composition and structure as much or more than severe fires at historic frequencies. ?? 2007 by the Ecological Society of America.

Ecological Applications↗

USGS2014

Simulating future residential property losses from wildfire in Flathead County, Montana

Wildfire damages to private residences in the United States and elsewhere have increased as a result of expansion of the wildland-urban interface (WUI) and other factors. Understanding this unwelcome trend requires analytical frameworks that simulate how various interacting social, economic, and biophysical factors influence those damages. A methodological framework is developed for simulating expected residential property losses from wildfire [E(RLW)], which is a probabilistic monetary measure of wildfire risk to residential properties in the WUI. E(RLW) is simulated for Flathead County, Montana for five, 10-year subperiods covering the period 2010-2059, under various assumptions about future climate change, economic growth, land use policy, and forest management. Results show statistically significant increases in the spatial extent of WUI properties, the number of residential structures at risk from wildfire, and E(RLW) over the 50-year evaluation period for both the county and smaller subareas (i.e., neighborhoods and parcels). The E(RLW) simulation framework presented here advances the field of wildfire risk assessment by providing a finer-scale tool that incorporates a set of dynamic, interacting processes. The framework can be applied using other scenarios for climate change, economic growth, land use policy, and forest management, and in other areas.

Montana↗

USGS2018

Assessing historical and projected carbon balance of Alaska: A synthesis of results and policy/management implications

We summarize the results of a recent interagency assessment of land carbon dynamics in Alaska, in which carbon dynamics were estimated for all major terrestrial and aquatic ecosystems for the historical period (1950–2009) and a projection period (2010–2099). Between 1950 and 2009, upland and wetland (i.e., terrestrial) ecosystems of the state gained 0.4 Tg C/yr (0.1% of net primary production, NPP), resulting in a cumulative greenhouse gas radiative forcing of 1.68 × 10 −3 W/m 2 . The change in carbon storage is spatially variable with the region of the Northwest Boreal Landscape Conservation Cooperative (LCC) losing carbon because of fire disturbance. The combined carbon transport via various pathways through inland aquatic ecosystems of Alaska was estimated to be 41.3 Tg C/yr (17% of terrestrial NPP). During the projection period (2010–2099), carbon storage of terrestrial ecosystems of Alaska was projected to increase (22.5–70.0 Tg C/yr), primarily because of NPP increases of 10–30% associated with responses to rising atmospheric CO 2 , increased nitrogen cycling, and longer growing seasons. Although carbon emissions to the atmosphere from wildfire and wetland CH 4 were projected to increase for all of the climate projections, the increases in NPP more than compensated for those losses at the statewide level. Carbon dynamics of terrestrial ecosystems continue to warm the climate for four of the six future projections and cool the climate for only one of the projections. The attribution analyses we conducted indicated that the response of NPP in terrestrial ecosystems to rising atmospheric CO 2 (~5% per 100 ppmv CO 2 ) saturates as CO 2 increases (between approximately +150 and +450 ppmv among projections). This response, along with the expectation that permafrost thaw would be much greater and release large quantities of permafrost carbon after 2100, suggests that projected carbon gains in terrestrial ecosystems of Alaska may not be sustained. From a national perspective, inclusion of all of Alaska in greenhouse gas inventory reports would ensure better accounting of the overall greenhouse gas balance of the nation and provide a foundation for considering mitigation activities in areas that are accessible enough to support substantive deployment.

Alaska↗

USGS2015

Rangeland monitoring reveals long-term plant responses to precipitation and grazing at the landscape scale

Managers of rangeland ecosystems require methods to track the condition of natural resources over large areas and long periods of time as they confront climate change and land use intensification. We demonstrate how rangeland monitoring results can be synthesized using ecological site concepts to understand how climate, site factors, and management actions affect long-term vegetation dynamics at the landscape-scale. Forty-six years of rangeland monitoring conducted by the Bureau of Land Management (BLM) on the Colorado Plateau reveals variable responses of plant species cover to cool-season precipitation, land type (ecological site groups), and grazing intensity. Dominant C 3 perennial grasses (Achnatherum hymenoides , Hesperostipa comata) , which are essential to support wildlife and livestock on the Colorado Plateau, had responses to cool-season precipitation that were at least twice as large as the dominant C 4 perennial grass (Pleuraphis jamesii) and woody vegetation. However, these C 3 perennial grass responses to precipitation were reduced by nearly one-third on grassland ecological sites with fine- rather than coarse-textured soils, and there were no detectable C 3 perennial grass responses to precipitation on ecological sites dominated by a dense-growing shrub, Coleogyne ramosissima . Heavy grazing intensity further reduced the responses of C 3 perennial grasses to cool-season precipitation on ecological sites with coarse-textured soils and surprisingly reduced the responses of shrubs as well. By using ecological site groups to assess rangeland condition, we were able to improve our understanding of the long-term relationships between vegetation change and climate, land use, and site characteristics, which has important implications for developing landscape-scale monitoring strategies.

Arizona, Colorado, New Mexico, Utah↗

USGS2019

Modeling spatially and temporally complex range dynamics when detection is imperfect

Species distributions are determined by the interaction of multiple biotic and abiotic factors, which produces complex spatial and temporal patterns of occurrence. As habitats and climate change due to anthropogenic activities, there is a need to develop species distribution models that can quantify these complex range dynamics. In this paper, we develop a dynamic occupancy model that uses a spatial generalized additive model to estimate non-linear spatial variation in occupancy not accounted for by environmental covariates. The model is flexible and can accommodate data from a range of sampling designs that provide information about both occupancy and detection probability. Output from the model can be used to create distribution maps and to estimate indices of temporal range dynamics. We demonstrate the utility of this approach by modeling long-term range dynamics of 10 eastern North American birds using data from the North American Breeding Bird Survey. We anticipate this framework will be particularly useful for modeling species’ distributions over large spatial scales and for quantifying range dynamics over long temporal scales.

Scientific Reports↗

USGS2021

Climatic aridity shapes post-fire interactions between Ceanothus spp. and Douglas-fir (Pseudotsuga menziesii) across the Klamath Mountains

Climate change is leading to increased drought intensity and fire frequency, creating early-successional landscapes with novel disturbance–recovery dynamics. In the Klamath Mountains of northwestern California and southwestern Oregon, early-successional interactions between nitrogen (N)-fixing shrubs ( Ceanothus spp.) and long-lived conifers (Douglas-fir) are especially important determinants of forest development. We sampled post-fire vegetation and soil biogeochemistry in 57 plots along gradients of time since fire (7–28 years) and climatic water deficit (aridity). We found that Ceanothus biomass increased, and Douglas-fir biomass decreased with increasing aridity. High aridity and Ceanothus biomass interacted with lower soil C:N more than either factor alone. Ceanothus biomass was initially high after fire and declined with time, suggesting a large initial pulse of N-fixation that could enhance N availability for establishing Douglas-fir. We conclude that future increases in aridity and wildfire frequency will likely limit post-fire Douglas-fir establishment, though Ceanothus may ameliorate some of these impacts through benefits to microclimate and soils. Results from this study contribute to our understanding of the effects of climate change and wildfires on interspecific interactions and forest dynamics. Management seeking to accelerate forest recovery after high-severity fire should emphasize early-successional conifer establishment while maintaining N-fixing shrubs to enhance soil fertility.

California, Oregon↗