The majority of the wettest cities in the United States are located in the Southeast. The major city with the most precipitation is New Orleans, Louisiana, which receives an average of 1592 millimeters (62.7 inches) of precipitation every year, based on an average between 1981 and 2010.

This statistic shows the ten major U.S. cities with the most rainy days per year between 1981 and 2010. Rochester, New York, had an average of about 167 days per year with precipitation. The sunniest city in the U.S. was Phoenix, Arizona, with an average of 85 percent of sunshine per day.

In 2024, Louisiana recorded ***** inches of precipitation. This was the highest precipitation within the 48 contiguous U.S. states that year. On the other hand, Nevada was the driest state, with only **** inches of precipitation recorded. Precipitation across the United States Not only did Louisiana record the largest precipitation volume in 2024, but it also registered the highest precipitation anomaly that year, around 14.36 inches above the 1901-2000 annual average. In fact, over the last decade, rainfall across the United States was generally higher than the average recorded for the 20th century. Meanwhile, the driest states were located in the country's southwestern region, an area which – according to experts – will become even drier and warmer in the future. How does global warming affect precipitation patterns? Rising temperatures on Earth lead to increased evaporation which – ultimately – results in more precipitation. Since 1900, the volume of precipitation in the United States has increased at an average rate of **** inches per decade. Nevertheless, the effects of climate change on precipitation can vary depending on the location. For instance, climate change can alter wind patterns and ocean currents, causing certain areas to experience reduced precipitation. Furthermore, even if precipitation increases, it does not necessarily increase the water availability for human consumption, which might eventually lead to drought conditions.

Brevard (North Carolina) recorded the wettest precipitation anomaly in 2019 with around ** inches of precipitation more than average. Wet years were seen in cities that were hit by heavy rain events such as tropical storms.

Historical changes of annual temperature and precipitation indices at selected 210 U.S. cities

This dataset provide:

Annual average temperature, total precipitation, and temperature and precipitation extremes calculations for 210 U.S. cities.

Historical rates of changes in annual temperature, precipitation, and the selected temperature and precipitation extreme indices in the 210 U.S. cities.

Estimated thresholds (reference levels) for the calculations of annual extreme indices including warm and cold days, warm and cold nights, and precipitation amount from very wet days in the 210 cities.

Annual average of daily mean temperature, Tmax, and Tmin are included for annual average temperature calculations. Calculations were based on the compiled daily temperature and precipitation records at individual cities.

Temperature and precipitation extreme indices include: warmest daily Tmax and Tmin, coldest daily Tmax and Tmin , warm days and nights, cold days and nights, maximum 1-day precipitation, maximum consecutive 5-day precipitation, precipitation amounts from very wet days.

Number of missing daily Tmax, Tmin, and precipitation values are included for each city.

Rates of change were calculated using linear regression, with some climate indices applied with the Box-Cox transformation prior to the linear regression.

The historical observations from ACIS belong to Global Historical Climatological Network - daily (GHCN-D) datasets. The included stations were based on NRCC’s “ThreadEx” project, which combined daily temperature and precipitation extremes at 255 NOAA Local Climatological Locations, representing all large and medium size cities in U.S. (See Owen et al. (2006) Accessing NOAA Daily Temperature and Precipitation Extremes Based on Combined/Threaded Station Records).

Resources:

See included README file for more information.

Additional technical details and analyses can be found in: Lai, Y., & Dzombak, D. A. (2019). Use of historical data to assess regional climate change. Journal of climate, 32(14), 4299-4320. https://doi.org/10.1175/JCLI-D-18-0630.1

Other datasets from the same project can be accessed at: https://kilthub.cmu.edu/projects/Use_of_historical_data_to_assess_regional_climate_change/61538

ACIS database for historical observations: http://scacis.rcc-acis.org/

GHCN-D datasets can also be accessed at: https://www.ncei.noaa.gov/data/global-historical-climatology-network-daily/

Station information for each city can be accessed at: http://threadex.rcc-acis.org/

• 2024 August updated -

  Annual calculations for 2022 and 2023 were added.

  Linear regression results and thresholds for extremes were updated because of the addition of 2022 and 2023 data.

  Note that future updates may be infrequent.

• 2022 January updated -

  Annual calculations for 2021 were added.

  Linear regression results and thresholds for extremes were updated because of the addition of 2021 data.

• 2021 January updated -

  Annual calculations for 2020 were added.

  Linear regression results and thresholds for extremes were updated because of the addition of 2020 data.

• 2020 January updated -

  Annual calculations for 2019 were added.

  Linear regression results and thresholds for extremes were updated because of the addition of 2019 data.

  Thresholds for all 210 cities were combined into one single file – Thresholds.csv.

• 2019 June updated -

  Baltimore was updated with the 2018 data (previously version shows NA for 2018) and new ID to reflect the GCHN ID of Baltimore-Washington International AP. city_info file was updated accordingly.

  README file was updated to reflect the use of "wet days" index in this study. The 95% thresholds for calculation of wet days utilized all daily precipitation data from the reference period and can be different from the same index from some other studies, where only days with at least 1 mm of precipitation were utilized to calculate the thresholds. Thus the thresholds in this study can be lower than the ones that would've be calculated from the 95% percentiles from wet days (i.e., with at least 1 mm of precipitation).

Extreme rainfall events can lead to catastrophic consequences such as widespread property damage or loss of lives. In engineering design, we take into account different characteristics of extreme events, such as average recurrence interval, duration and volume of rainfall delivered. We have so far relied on historical rainfall observations to estimate these characteristics, however, there is evidence that rainfall patterns in the future will look different from the past. Global Climate Models can be used to understand to some extent, the magnitude of these changes. In this repository, we provide updated rainfall volumes for 24-hour duration storms using 5 different downscaled climate model datasets and for selected cities in the United States. For more information on the update process, please refer to the associated publication of this dataset.

U.S. 15 Minute Precipitation Data is digital data set DSI-3260, archived at the National Climatic Data Center (NCDC). This is precipitation data. The primary source of data for this file is approximately 2,000 mostly U.S. weather stations operated or managed by the U.S. National Weather Service. Stations are primary, secondary, or cooperative observer sites that have the capability to measure precipitation at 15 minute intervals. This dataset contains 15-minute precipitation data (reported 4 times per hour, if precip occurs) for U.S. stations along with selected non-U.S. stations in U.S. territories and associated nations. It includes major city locations and many small town locations. Daily total precipitation is also included as part of the data record. NCDC has in archive data from most states as far back as 1970 or 1971, and continuing to the present day. The major parameter is precipitation amounts at 15 minute intervals, when precipitation actually occurs.

The Baltimore radar rainfall dataset was developed from a multi-sensor analysis combining radar rainfall estimates from the Sterling, VA WSR88D radar (KLWX) with measurements from a collection of ground based rain gages. The archived data have a 15-minute time resolution and a grid resolution of 0.01 degree latitude/longitude (approximately 1 km x 1 km); 15-minute rainfall accumulations for each grid are in mm. The dataset spans 22 years, 2000-2021, and covers an area of approximately 4,900 km^2 (70 by 70 grids, each with approximate area of 1 km^2) surrounding the Baltimore, MD metropolitan area (Figure 1). The rainfall data cover the six months from April to September of each year. This is the period with most intense sub-daily rainfall and the period for which radar measurements are most accurate. Figure 1 illustrates the climatological analyses of mean annual frequency of days with at least 1 hour rainfall exceeding 25 mm. The striking spatial variability of convective rainfall is illustrated in Figure 2 by the April-September climatology of annual lightning strikes.

As with many long-term environmental data sets, sensor technology has changed during the time period of the archive. The Sterling, VA WSR88D radar underwent a hardware upgrade from single polarization to dual polarization in 2012. Prior to the upgrade, rainfall was estimated using a conventional radar-reflectivity algorithm (HydroNEXRAD) which converts reflectivity measurements in polar coordinates from the lowest sweep to rainfall estimates on a 0.01 degree latitude-longitude grid at the surface (see Seo et al. 2010 and Smith et al. 2012 for details on the algorithm). The polarimetric upgrade introduced new measurements into the radar-rainfall algorithm. In addition to reflectivity, the operational rainfall product, Digital Precipitation Rate (DPR), directly uses differential reflectivity and specific differential phase shift measurements to estimate rainfall (https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C00708; see also Giangrande and Ryzhkov 2008). Details of the algorithm structure and parameterization for the DPR radar-rainfall estimates have been modified during the 10-year period of the data set.

A storm-based (daily) multiplicative mean field bias has been applied to both datasets. The mean field bias is computed as the ratio of daily rain gage rainfall at a point to daily radar rainfall for the bin that contains the gage. The rain gage dataset is compiled from rain gages in the Baltimore metropolitan region and surrounding areas and includes gages acquired from both Baltimore City and Baltimore County, and the Global Historical Climatology Network daily (GHCNd). Mean field bias improves rainfall estimates and diminishes the impacts of changing measurement procedures.

The dataset has been archived in 2 formats: netCDF gridded rainfall, 1 file for each 15-minute time period, and csv or excel format point rainfall (1 point at the center of each grid) in a timeseries format with 1 file per calendar month covering the entire 70x70 domain. The csv files are in folders organized by calendar year. The first five columns in each file represent year, month, day, hour, and minute and can be combined to generate a unique date-time value for each time step. Each additional column is a complete time series for the month and represents data from one of the 1-km2 grid cells in the original data set.

The latitude and longitude coordinates for each pixel in the grid are provided. The latitude and longitude represent the centroid of the cell, which is square when represented in latitude and longitude coordinates and rectangular when represented in other distance-based coordinate systems such as State Plane or Universal Transverse Mercator. There are 4900 pixels in the domain. In order to visualize the data using GIS or other software, the user needs to associate each column in the annual rainfall file with the latitude and longitude values for that grid cell number.

These data may be subject to modest revision or reformatting in future versions. The current version is version 2.0 and is being offered to users who wish to explore the data. We will revise this document as needed.

These daily weather records were compiled from a subset of stations in the Global Historical Climatological Network (GHCN)-Daily dataset. A weather record is considered broken if the value exceeds the maximum (or minimum) value recorded for an eligible station. A weather record is considered tied if the value is the same as the maximum (or minimum) value recorded for an eligible station. Daily weather parameters include Highest Min/Max Temperature, Lowest Min/Max Temperature, Highest Precipitation, Highest Snowfall and Highest Snow Depth. All stations meet defined eligibility criteria. For this application, a station is defined as the complete daily weather records at a particular location, having a unique identifier in the GHCN-Daily dataset. For a station to be considered for any weather parameter, it must have a minimum of 30 years of data with more than 182 days complete in each year. This is effectively a 30-year record of service requirement, but allows for inclusion of some stations which routinely shut down during certain seasons. Small station moves, such as a move from one property to an adjacent property, may occur within a station history. However, larger moves, such as a station moving from downtown to the city airport, generally result in the commissioning of a new station identifier. This tool treats each of these histories as a different station. In this way, it does not thread the separate histories into one record for a city. Records Timescales are characterized in three ways. In order of increasing noteworthiness, they are Daily Records, Monthly Records and All Time Records. For a given station, Daily Records refers to the specific calendar day: (e.g., the value recorded on March 7th compared to every other March 7th). Monthly Records exceed all values observed within the specified month (e.g., the value recorded on March 7th compared to all values recorded in every March). All-Time Records exceed the record of all observations, for any date, in a station's period of record. The Date Range and Location features are used to define the time and location ranges which are of interest to the user. For example, selecting a date range of March 1, 2012 through March 15, 2012 will return a list of records broken or tied on those 15 days. The Location Category and Country menus allow the user to define the geographic extent of the records of interest. For example, selecting Oklahoma will narrow the returned list of records to those that occurred in the state of Oklahoma, USA. The number of records broken for several recent periods is summarized in the table and updated daily. Due to late-arriving data, the number of recent records is likely underrepresented in all categories, but the ratio of records (warm to cold, for example) should be a fairly strong estimate of a final outcome. There are many more precipitation stations than temperature stations, so the raw number of precipitation records will likely exceed the number of temperature records in most climatic situations.

In 2011, Buffalo, New York was the major city in the United States with the most partial to heavy cloud cover, with 311 days of clouds in that year. Seattle, Pittsburg, Rochester, and Cleveland rounded out the top five cities.

Buffalo’s climate

Buffalo, New York, is located on the eastern end of Lake Erie and is the origin point of the Niagara River, and its location on Lake Erie helps to regulate the city’s climate. However, between 1981 and 2010, it had an average of 167 days with more than 0.01 inches of rainfall per year, and also had an average wind speed of 11.8 miles per hour.

The second largest city in New York

Buffalo is the second largest city in New York state, with a metro area population of 1.13 million in 2017. The city is not far from Niagara Falls, which was listed as one of the most expensive summer destinations in the state of New York.

On the continental scale, climate is an important determinant of the distributions of plant taxa and ecoregions. To quantify and depict the relations between specific climate variables and these distributions, we placed modern climate and plant taxa distribution data on an approximately 25-kilometer (km) equal-area grid with 27,984 points that cover Canada and the continental United States (Thompson and others, 2015). The gridded climatic data include annual and monthly temperature and precipitation, as well as bioclimatic variables (growing degree days, mean temperatures of the coldest and warmest months, and a moisture index) based on 1961-1990 30-year mean values from the University of East Anglia (UK) Climatic Research Unit (CRU) CL 2.0 dataset (New and others, 2002), and absolute minimum and maximum temperatures for 1951-1980 interpolated from climate-station data (WeatherDisc Associates, 1989). As described below, these data were used to produce portions of the "Atlas of relations between climatic parameters and distributions of important trees and shrubs in North America" (hereafter referred to as "the Atlas"; Thompson and others, 1999a, 1999b, 2000, 2006, 2007, 2012a, 2015). Evolution of the Atlas Over the 16 Years Between Volumes A & B and G: The Atlas evolved through time as technology improved and our knowledge expanded. The climate data employed in the first five Atlas volumes were replaced by more standard and better documented data in the last two volumes (Volumes F and G; Thompson and others, 2012a, 2015). Similarly, the plant distribution data used in Volumes A through D (Thompson and others, 1999a, 1999b, 2000, 2006) were improved for the latter volumes. However, the digitized ecoregion boundaries used in Volume E (Thompson and others, 2007) remain unchanged. Also, as we and others used the data in Atlas Volumes A through E, we came to realize that the plant distribution and climate data for areas south of the US-Mexico border were not of sufficient quality or resolution for our needs and these data are not included in this data release. The data in this data release are provided in comma-separated values (.csv) files. We also provide netCDF (.nc) files containing the climate and bioclimatic data, grouped taxa and species presence-absence data, and ecoregion assignment data for each grid point (but not the country, state, province, and county assignment data for each grid point, which are available in the .csv files). The netCDF files contain updated Albers conical equal-area projection details and more precise grid-point locations. When the original approximately 25-km equal-area grid was created (ca. 1990), it was designed to be registered with existing data sets, and only 3 decimal places were recorded for the grid-point latitude and longitude values (these original 3-decimal place latitude and longitude values are in the .csv files). In addition, the Albers conical equal-area projection used for the grid was modified to match projection irregularities of the U.S. Forest Service atlases (e.g., Little, 1971, 1976, 1977) from which plant taxa distribution data were digitized. For the netCDF files, we have updated the Albers conical equal-area projection parameters and recalculated the grid-point latitudes and longitudes to 6 decimal places. The additional precision in the location data produces maximum differences between the 6-decimal place and the original 3-decimal place values of up to 0.00266 degrees longitude (approximately 143.8 m along the projection x-axis of the grid) and up to 0.00123 degrees latitude (approximately 84.2 m along the projection y-axis of the grid). The maximum straight-line distance between a three-decimal-point and six-decimal-point grid-point location is 144.2 m. Note that we have not regridded the elevation, climate, grouped taxa and species presence-absence data, or ecoregion data to the locations defined by the new 6-decimal place latitude and longitude data. For example, the climate data described in the Atlas publications were interpolated to the grid-point locations defined by the original 3-decimal place latitude and longitude values. Interpolating the data to the 6-decimal place latitude and longitude values would in many cases not result in changes to the reported values and for other grid points the changes would be small and insignificant. Similarly, if the digitized Little (1971, 1976, 1977) taxa distribution maps were regridded using the 6-decimal place latitude and longitude values, the changes to the gridded distributions would be minor, with a small number of grid points along the edge of a taxa's digitized distribution potentially changing value from taxa "present" to taxa "absent" (or vice versa). These changes should be considered within the spatial margin of error for the taxa distributions, which are based on hand-drawn maps with the distributions evidently generalized, or represented by a small, filled circle, and these distributions were subsequently hand digitized. Users wanting to use data that exactly match the data in the Atlas volumes should use the 3-decimal place latitude and longitude data provided in the .csv files in this data release to represent the center point of each grid cell. Users for whom an offset of up to 144.2 m from the original grid-point location is acceptable (e.g., users investigating continental-scale questions) or who want to easily visualize the data may want to use the data associated with the 6-decimal place latitude and longitude values in the netCDF files. The variable names in the netCDF files generally match those in the data release .csv files, except where the .csv file variable name contains a forward slash, colon, period, or comma (i.e., "/", ":", ".", or ","). In the netCDF file variable short names, the forward slashes are replaced with an underscore symbol (i.e., "_") and the colons, periods, and commas are deleted. In the netCDF file variable long names, the punctuation in the name matches that in the .csv file variable names. The "country", "state, province, or territory", and "county" data in the .csv files are not included in the netCDF files. Data included in this release: - Geographic scope. The gridded data cover an area that we labelled as "CANUSA", which includes Canada and the USA (excluding Hawaii, Puerto Rico, and other oceanic islands). Note that the maps displayed in the Atlas volumes are cropped at their northern edge and do not display the full northern extent of the data included in this data release. - Elevation. The elevation data were regridded from the ETOPO5 data set (National Geophysical Data Center, 1993). There were 35 coastal grid points in our CANUSA study area grid for which the regridded elevations were below sea level and these grid points were assigned missing elevation values (i.e., elevation = 9999). The grid points with missing elevation values occur in five coastal areas: (1) near San Diego (California, USA; 1 grid point), (2) Vancouver Island (British Columbia, Canada) and the Olympic Peninsula (Washington, USA; 2 grid points), (3) the Haida Gwaii (formerly Queen Charlotte Islands, British Columbia, Canada) and southeast Alaska (USA, 9 grid points), (4) the Canadian Arctic Archipelago (22 grid points), and (5) Newfoundland (Canada; 1 grid point). - Climate. The gridded climatic data provided here are based on the 1961-1990 30-year mean values from the University of East Anglia (UK) Climatic Research Unit (CRU) CL 2.0 dataset (New and others, 2002), and include annual and monthly temperature and precipitation. The CRU CL 2.0 data were interpolated onto the approximately 25-km grid using geographically-weighted regression, incorporating local lapse-rate estimation and correction. Additional bioclimatic variables (growing degree days on a 5 degrees Celsius base, mean temperatures of the coldest and warmest months, and a moisture index calculated as actual evapotranspiration divided by potential evapotranspiration) were calculated using the interpolated CRU CL 2.0 data. Also included are absolute minimum and maximum temperatures for 1951-1980 interpolated in a similar fashion from climate-station data (WeatherDisc Associates, 1989). These climate and bioclimate data were used in Atlas volumes F and G (see Thompson and others, 2015, for a description of the methods used to create the gridded climate data). Note that for grid points with missing elevation values (i.e., elevation values equal to 9999), climate data were created using an elevation value of -120 meters. Users may want to exclude these climate data from their analyses (see the Usage Notes section in the data release readme file). - Plant distributions. The gridded plant distribution data align with Atlas volume G (Thompson and others, 2015). Plant distribution data on the grid include 690 species, as well as 67 groups of related species and genera, and are based on U.S. Forest Service atlases (e.g., Little, 1971, 1976, 1977), regional atlases (e.g., Benson and Darrow, 1981), and new maps based on information available from herbaria and other online and published sources (for a list of sources, see Tables 3 and 4 in Thompson and others, 2015). See the "Notes" column in Table 1 (https://pubs.usgs.gov/pp/p1650-g/table1.html) and Table 2 (https://pubs.usgs.gov/pp/p1650-g/table2.html) in Thompson and others (2015) for important details regarding the species and grouped taxa distributions. - Ecoregions. The ecoregion gridded data are the same as in Atlas volumes D and E (Thompson and others, 2006, 2007), and include three different systems, Bailey's ecoregions (Bailey, 1997, 1998), WWF's ecoregions (Ricketts and others, 1999), and Kuchler's potential natural vegetation regions (Kuchler, 1985), that are each based on distinctive approaches to categorizing ecoregions. For the Bailey and WWF ecoregions for North America and the Kuchler potential natural vegetation regions for the contiguous United States (i.e.,

Display and/or analyses requiring spatially distributed monthly or annual precipitation for the climatological period 1961-90. Parameter-elevation Regressions on Independent Slopes Model (PRISM) derived raster data is the underlying data set from which the polygons and vectors were created. PRISM is an analytical model that uses point data and a digital elevation model (DEM) to generate gridded estimates of annual, monthly and event-based climatic parameters. Detailed descriptions of the PRISM raster data can be found with the accompanying raster data sets. For further information, the online PRISM homepage can be found at