1. Introduction

The US is the fourth largest rice exporter in the world which is in part due to intensive production and the highest-yielding rice agroecosystems (FAOSTAT, 2015; Childs, 2019). Rice is an important economic crop in the US, with average national production of $371 per ha (McBride et al., 2018). Among the US major crops produced, rice has the highest average return on investment (USDA-NASS, 2016). The projected total rice export from 2019-2020 is 95.0 million hundredweight (cwt) (Childs, 2019). Major US rice producing states include Arkansas, California, Louisiana, Mississippi, Missouri and Texas with CA being the highest yielding and AR being the largest harvested areas in the country (USDA, 2018b). Total rice production in the US increased by 128% between 1960 and 2015 with the highest national mean yield reported in 2013 at 8617 kg∙ha⁻¹ for 1 million ha (USDA-NASS, 2016). These remarkable grain yield increases were driven by increased adoption of high yielding, fertilizer efficient, and disease-insect-resistant rice varieties and improved agronomic practices (i.e. irrigation, tillage). It has been projected that rice yield must increase at a rate of 1.2% - 2.4% from its current rate during the next decade to meet the future demand (Grassini et al., 2013). In an effort to sustain future US rice production, cropping systems will require significant yield gains per unit land area while conserving natural resources and reducing environmental degradation.

Initiatives for improving rice production system have been developed to meet the country’s increasing rice demand (Batres-Marquez et al., 2009; US Census 2011). Emerging strategies to increase grain yield potential include commercialization of hybrid rice, improved water management, smart fertilization, and stepped-up extension activities (McBride et al., 2018). While advancing rice production systems have been realized in the US for potential rice yield gains, crop production is highly impacted by weather fluctuations. One of the climatic factors affecting rice yields is the increasing air temperature. Specifically, the frequency of warm nights has increased for the past 60 years (Donat et al., 2013). It has been reported that each 1°C increase in night-time air temperature during the dry growing season led to 10% drop in rice yield (Peng et al., 2004). Other studies reported yield penalties due to high night-time air temperature (HNT) varied from 8% to 95% (Mohammed & Tarpley, 2009b; Mohammed & Tarpley, 2010; Zhang et al., 2013; Bahuguna et al., 2017). Differences in yield reduction were attributed to the rice varieties used and air temperature thresholds during the experiment. Overall, these reports provided evidence that HNT stress negatively impacts rice production.

Studies indicated that reduction in grain yield due to HNT is associated with lower spikelet fertility (Mohammed & Tarpley, 2010), lower percent pollen germination (Mohammed & Tarpley, 2009a), and increased night respiration rates (Mohammed & Tarpley, 2009b; Bahuguna et al., 2017). HNT also had direct effects on panicle development and panicle weight (Kanno & Makino, 2010). The starch content of rice grains was relatively lower among selected sensitive and high yielding rice varieties (Bahuguna et al., 2017). Grain quality was also affected under HNT through the development of chalky kernels, decreased grain weight and its dimensions (Cooper et al., 2008; Morita et al., 2005, Lanning et al., 2011; Bahuguna et al., 2017). The rationale of many studies behind understanding HNT impacts on grain yield has been mostly anchored on predicted global temperature trends based on climate model scenarios. According to IPCC’s current predictions, the global temperature may continue rising by 0.2°C per decade (IPCC, 2018).

Multiple weather databases are available for assessing temperature patterns across the global landscape. These weather databases (i.e. NOAA-NCDC, NASA-POWER) provide complete coverage of terrestrial surfaces with weather information derived typically from global computer models, remotely sensed surface satellite data, interpolated or actual weather station data. Despite various weather data sources and the need to assess climate variability effects on crop production, there have been no reports on the actual deviations of seasonal night-time air temperature during rice growth stages (i.e. R4 or flowering stage) that are critical for yield and grain quality (Counce et al., 2000). Studies on temperature patterns of specific regions during crop growth are rarely calculated and validated against recorded daily temperature data from well-maintained weather stations. In the study of heat stress, air temperatures are often set based on assumptions or estimations that do not have a ground basis for a specific region. Actual air temperature trends based on weather station data in a particular location or farm is important to agronomist, breeders and researchers to accurately design studies that properly address the level of heat stress in rice. Additionally, the use of field-based data accurately improves technological advancement in rice breeding and farming systems for current and future climate. To address local temperature changes that may influence rice yields, this study aimed to calculate the prevalent seasonal night-time air temperature in the four major US rice producing states during rice growth period. In addition, we assessed the accuracy of simulated air temperature from weather databases using actual air temperature data reported in weather stations. The overarching goal of this study was to provide evidence of significant trends in seasonal night air temperature changes in major US rice-growing regions.

2. Materials and Methods

2.1. Weather Data Sources

Air temperature data used in this study were obtained from meteorological stations that have historical records of regional and/or global weather data. Temperature data were taken from 1) National Aeronautics and Space Administration-Prediction of Worldwide Energy Resources database (NASA-POWER), 2) the National Oceanic and Atmospheric Administration-National Climatic Data Center (NOAA-NCDC), 3) California Irrigation Management Information Systems database (CIMIS) and 4) Integrated Agricultural Information and Management System Climatic database (iAIMS) (Table 1). These weather databases are publicly accessible, cover diverse spatial and temporal resolution, and are widely used by researchers for studying the effects of climate change in food systems (van Wart et al., 2013; Morell et al., 2016; Mourtzini et al., 2017). The NASA-POWER database is produced by the NASA Langley Research Center (LaRC) POWER Project funded through the NASA Earth Science/Applied Science Program. Air temperature in the NASA-POWER database is retrieved on a daily basis at 1 - 2 km² resolution and derived from satellite observations coupled with NASA’s Modern Era Retro-analysis for Research and Applications (MERRA-2), an integrated system of models (Goddard Earth Observing System Model) informed by data from modern satellite, ground stations and other sources. In contrast, the NOAA-NCDC database is a collection of daily historical observations from 40,000 weather stations around the world, in which data underwent several quality assurance procedures (NCDC, 2011). Most of these weather stations are located in airports or cities; hence, observations cannot accurately and reliably be used for agricultural applications. In contrast, CIMIS and iAIMS are weather databases that retrieved meteorological data from stations located mostly in agricultural areas. CIMIS is a program unit in the Water Use and Efficiency Branch, Division of Regional Assistance of the California

¹NOAA-NCDC, National Oceanic and Atmospheric Administration - National Climatic Data Center; NASA-POWER, National Aeronautics and Space Administration - Prediction of Worldwide Energy Resources database; CIMIS, California Irrigation Management Information Systems database; iAIMS, Integrated Agricultural Information and Management System Climatic database.

Department of Water Resources, and has 145 automated weather stations in California (https://cimis.water.ca.gov/). CIMIS was developed in 1982 and designed to assist crop growers and irrigators in managing their water resources more efficiently (CIMIS, 2020). The iAIMS climatic database is part of integrated soil, climate, and crop databases developed in 2007 by the Agroecosystems Research Group at the Texas A & M University Agricultural Research and Extension Center in Beaumont, Texas (Yang et al., 2010). The iAIMS climatic database is an archive of daily historical data obtained from five weather data sources namely: 1) Texas A & M University, Crop Weather Program Weather Station Network at Corpus Christi, TX (https://beaumont.tamu.edu/ClimaticData/), 2) NOAA’s National Weather Service Aviation Weather, Cooperative Observer Program (COOP) network, 3) Texas A&M University, Meteorological Aviation Report (METAR) Department of Atmospheric Sciences, 4) NCDC Global Surface Summary of Day Data that is based on data exchanged under the World Meteorological Organization (WMO) World Weather Watch Program according to WMO Resolution 40 (Cg-XII), and 5) NOAA NCDC Climate data online (https://www.ncdc.noaa.gov/cdo-web/) (Yang et al., 2010). The iAIMS is a web-based interface that summarizes climatic data, soil and cropland databases to provide dynamic access to consolidated and filtered data for studying cropping systems performance and management. It consists of climatic data from over 25,000 weather stations worldwide and serves as a support tool in data retrieval of information without duplication. This study used the climatic data component of iAIMS which consisted of temperature data from weather stations in TX.

In this study, hourly and daily air temperature data in rice producing states were mined from NOAA-NCDC, CIMIS and iAIMS climatic databases (Table 1). Air temperature data from NASA-POWER database were also used to compare simulated temperature data and observed weather station data collected from NOAA-NCDC, CIMIS, iAIMS databases. Air temperature data in these databases were archived on an hourly to daily basis and at high spatial resolution: 10 to 110 km (NOAA-NCDC) and 0.5˚ × 0.5˚ latitude and longitude global grid (NASA-POWER). In the analysis of night-time air temperature pattern, rice growing counties were identified in the rice producing states using the USDA-NASS crop website (http://www.nass.usda.gov/Surveys/Guide%20_to_NASS_%20Surveys/Crops_%20Stocks/index.php). More than 70 years of air temperature data were available for all rice states (1941-2018). The longest consecutive years available from NOAA-NCDC database for air temperature were CA. In order to maintain unbiased analysis of weather pattern, we considered equivalent time series with >20 years of data archived across all states. Air temperature parameters such as maximum, minimum, and mean air temperature data were not always available in all weather station databases (Table 1). In order to be accurate with long-term weather data analysis and because minimum air temperature data were not consistently monitored and recorded in all databases during the 76-year period, we used average hourly air temperature for all rice states and years in the trend analysis of air temperature.

2.2. Boundary Conditions

To calculate and summarize mean ambient night-time air temperature, the data retrieval was limited by the time of the day, month, and location. Air temperature data were collected from 20:00 to 06:00 within the rice growing period. The time of observation was limited to 10-hour period to capture the night field condition. According to Bahuguna et al. (2017), important physiological processes such as respiration which influence rice grain filling occurred at night. In this study, rice growing months were May to September for all locations and states. While flowering and reproductive stages are critical part in rice growth, under field growing condition, changes in mean air temperature do not occur only during specific period of growth but throughout the growing season. Grain yield of rice is not only affected by heat stress at those growth stages but mainly affected by crop and grain establishment that is influenced by heat stress at the beginning of the growing season and near maturity. We used the rice growing duration to capture the historical pattern of air temperature to which rice crop is exposed. Also, aiming for higher yield depends on increasing total biomass. According to Peng et al. (2004), total crop biomass is determined mainly by crop photosynthesis and respiration losses, and both are sensitive to temperature.

Four US rice states included in this study were: CA, AR, LA and TX accounting for about 90% of US rice area and production (USDA, 2018a). California is located in the US western region while TX, LA and AR are located in the South (Figure 1). Air temperature data from Mississippi (MS) and Missouri (MO) were not considered because of total rice harvested areas comprised only 5.7%

and 6.5% of the total US harvested ricearea, respectively. Also, available air temperature weather station data in these two states had <20 years record. According to USDA-NASS (2018c), AR covered approximately 80% of the total rice areas in US Delta Region. Overall, air temperature from 1940 to 2018 was retrieved from 23 locations across four major rice states (CA, AR, LA, TX), resulting in ca. (about) >1200 observations per site per year. For those sites that generated incomplete time-series due to erroneous entries or missing values (24% to 50%) we did not perform gap-filling or interpolation from air temperature from the two nearest weather stations to the weather station being estimated. To achieve real-time calculation of air temperature pattern, those sites-years with missing values were excluded in the trend analysis of night-time air temperature for each state.

2.3. Data Analysis

Average night-time air temperature was calculated from 1^st May to 30^th September from 20:00 to 06:00 for each location-state-year, separately for NOAA-NCDC, CIMIS and iAIMS weather station database and the NASA-POWER database. In each of the four US rice-growing states across all years and locations, linear regression was used to identify if a significant increase in night-time air temperature occurred in the past seven decades (SYSTAT, 2013). The air temperature data from NOAA-NCDC were regarded as the dependent variable while time in unit of the year was considered as the independent variable. The slope of regression was tested for significant difference at P < 0.05 level. We explored the accuracy of the interpolated air temperature from NASA-POWER database, a widely used weather database. Using the same location coordinates and years, simulated air temperature from NASA-POWER database was compared to actual air temperature data from active weather stations (NOAA-NCDC, CIMIS, iAIMS) in each rice-growing state. For a given state, accuracy and biases between interpolated and observed air temperature were assessed with the absolute mean error (ME), root mean square error (RMSE), and coefficient of variation (CV). These parameters were calculated as follows:

$\text{ME}=\frac{{\displaystyle {\sum }_{i=1}^n\left(Tem{p}_i^P-Tem{p}_i^O\right)}}{n}$ (1)

$\text{RMSE}=\sqrt{\frac{{\displaystyle {\sum }_{i=1}^n{\left(Tem{p}_i^P-Tem{p}_i^O\right)}^2}}{n}}$ (2)

$\text{CV}=\frac{\sigma }{\mu }$ (3)

where n is number of site-year, $Tem{p}_i^O$ is air temperature from active weather station for the ith site-year whereas $Tem{p}_i^P$ is the simulated/calculated night-time air temperature. The σ and μ are standard deviation and mean, respectively, of interpolated or observed air temperature. The mean error, ME, measures the mean magnitude and sign of the bias with air temperature from predicted or weather station data. Root mean square error, RMSE, measures the degree of agreement between weather data sources and quantifies average error on a relative basis compared to actual weather station data. CV measures the relative variability in a distribution.

3. Results and Discussion

3.1. Comparison between Estimated and Measured Air Temperature

The ability to accurately estimate temperature fluctuations under wide range of cropping environments depends on credible calculations and sources of weather data. In this study, hourly night-time air temperature data for some rice counties were fundamentally challenging to summarize due to daily time stamp were available for long-term weather data. In this study, we used temperature data from weather stations that reported hourly night air temperature to calculate changes in night air temperature. Prior to temperature trend analysis, we assessed the accuracy of night-time air temperature datasets from NASA-POWER database. The simulated datasets were compared to average daily measured night-time air temperature data from NOAA-NCDC, CIMIS and iAIMS databases. NASA-POWER air temperature datasets were different for most parts of US rice regions (CA, AR, LA) and deviations from measured values were consistent from 1981 to 2018 (Figure 2 and Figure 3). For CA and AR, NASA-POWER predicted values were approximately 2˚C - 4˚C lower compared to NOAA-NCDC’s observed values. In contrast, average hourly air temperature differences of roughly <1˚C were observed in LA and TX datasets and values between

predicted and observed tended to close the gap from 1980 to 2018 (Figure 2). The total number of observations in LA and TX datasets was 33% to 80% lower relative to observations used in CA and AR datasets (Figure 2). Deviations of predicted night air temperature from NASA-POWER database to night air temperature weather station data from various weather station sources were apparent (Figure 2 and Figure 3). To illustrate, predicted night air temperature in CA using NASA-POWER dataset underestimated night air temperature from NOAA-NCDC weather database (Figure 2). However, this relationship was reversed when NASA-POWER air temperature was compared to air temperature datasets from CIMIS database (Figure 3). The observed difference between NOAA-NCDC and CIMIS weather station data was probably due to the location of weather station managed by both organizations. NOAA-NCDC weather stations are generally located near airports or cities, while CIMIS weather stations are located near agricultural farms. These results showed relatively large uncertainties and a strong bias towards weather stations location and the total number of temperature data points for a given station. Overall, our results showed inconsistencies in predicting air temperature data and historical temperature trends when using the NASA-POWER data for night air temperature trend analysis in some US rice states.

3.2. Air Temperature from Weather Stations

The NOAA-NCDC database was used in the analysis of night-time air temperature patterns. All four states had weather stations with >20 years of recorded hourly air temperature and available for public use (Table 2). Following an exhaustive weather data collection, air temperature data were retrieved from 12, 5, 4, and 2 weather stations located near rice counties with 20 to 67 total recorded years for CA, AR, LA and TX, respectively (Table 2). California had the most number of observed air temperature data points due to weather stations availability with >20 years of data archived compared to the other three states (Table 2). Texas had the least number of historical air temperature data points on record with stations mostly located in the southeastern region of the state. Arkansas and LA had temperature data approximately 30% of CA’s temperature data owed to relatively fewer weather stations were available near rice counties.

Despite the extensive recorded years coverage of temperature data across rice states, very few well-maintained weather stations were located inside major rice producing counties. To illustrate these points, AR and LA are the largest southern rice states with total harvested rice area of 508,287 ha and 120,373 ha but there were only 4 to 5 weather stations with long-term data archived (Table 2). Several major producing rice counties within each state had no long-term air temperature data. This was true for Poinsett and Arkansas Counties in AR, Acadia and Jefferson Counties in LA, Wharton County in TX and Colusa and Sutter Counties in CA. Most weather stations found in all states were operating near

¹Total number of observations (air temperature data) was computed from the sum of total air temperature data in all locations and recorded years in each state.

regional and international airports. In fact, very few weather stations were located close to agricultural farms in the southern states. Due to sporadic measurements of air temperature in these regions, historical temperature trends may likely provide partial coverage of actual fluctuations of night-time air temperature in counties with fewer weather stations.

Within each state, historical weather station data varied among rice producing counties. In CA, >45 years of data were observed in all counties except in Butte with only 4 to 13 years of recorded data. Similarly, for AR, LA, and TX, the majority of temperature data archived were from 1 to 15 years. Gaps and many missing night-time air temperature data were common in rice counties particularly in AR and TX prior to 1980 (Figure 4). There were even isolated data in early 1940 with higher values compared to observations gathered after the gap. According to Schiermeier (2008), the post-war data anomaly characterized by a sudden measurement drop of night air temperatures recorded in early 1940’s was related mostly to the implementation of new methods to computing temperature as well as better instrumentation. Although more than 40 years of temperature data are publicly available, full temporal coverage of temperature in agricultural areas in the Southern states has not been achieved. Essentially, the underpinning causes for lack of daily temperature or missing data in major rice regions include; 1) diminishing number of active weather stations, 2) discontinued weather station data archive, 3) older records that are no longer accessible in real time, and 4) no physical weather stations existed on site. An increasing scarcity of weather stations near agricultural farms led many researchers to use weather data generated from satellite images or data interpolations from existing meteorological stations using derived empirical calculations. These predicted

weather parameters are commonly used to assess climate change and agricultural-related research such as decision-support tool systems to assist in crop and farm management (i.e. Foley et al., 2005; Lobell et al., 2011; Morell et al., 2016). However, several studies have compared interpolated weather data to actual weather station data and determined weather data simulation accuracy perform when used as input for crop model estimates (Van wart et al., 2013; Mourtzinis et al., 2017). Findings from these studies indicated that simulated weather station data could not replace measured weather data as a basis for field-scale agricultural applications i.e. yield, crop growth stage, water balance. Given the uncertainties of simulated weather data for crop research applications, calculations of night-time air temperature patterns in this study were based only on measured weather station data.

3.3. Night-Time Air Temperature Pattern in the Four Rice Producing States

On average, measured night-time air temperature data across all rice growing areas and years were above 23˚C (Figure 4, Table 3) and ranged from 22.6˚C to 29.5˚C during the rice growing period. The magnitude of mean seasonal night-time air temperature was highest in CA and lowest in TX (Figures 3(a)-(d)). Measured night-time air temperature fluctuated between ±0.2˚C and ±4˚C in all states with the most frequent in CA and intense fluctuation occurred in AR (Figures 4(e)-(h)). The minimal intensity of temperature fluctuations occurred in LA (Figure 3(g)). In all states, the rate of increase of night-time air temperature ranged from 0.01˚C to 0.02˚C per year since 1941. During the last 67-year period, night-time air temperature increased by 1.12^°C, 0.53˚C, and 0.30˚C in CA, AR, TX, respectively and no increase was observed in LA (0.0˚C) (Table 3). On average, calculated seasonal night-time air temperature using the NOAA-NCDC weather data were in close agreement with actual seasonal night-time

¹ME: Mean error was calculated using Equation (1). ²RMSE: Root mean square error was calculated using Equation (2). ³CV: coefficient of variations was calculated from Equation (3). ⁴Air temperature observations from weather stations. ⁵Air temperature was calculated using linear regression model.

air temperature across four states with modest estimation of night air temperature data (ME = 0.6˚C to 1.3˚C) (Table 3). Similarly, %RMSE for calculated and observed night-time air temperature based on weather database were 8.68%, 4.21%, 2.25% and 1.83% for CA, AR, LA, and TX, respectively.

The shift in seasonal night-time air temperature was significant in CA and AR (P-level ≤ 0.01) while LA and TX had no significant difference (P-level ≤ 0.10). The overall increase in night-time air temperature in CA was in the range of 2 to 22 times greater than the increase in night-time air temperature in the Southern states (Figure 4). The temperature rise observed in this study was unlikely due to local air pollution, as warming trends are consistent with temperature increases found elsewhere (Pathak et al., 2003; Peng et al., 2004). In fact, temperature increase in CA was consistent with the increase in minimum temperature of 1.13˚C reported by Peng et al. (2004) during the 25-year period of temperature trend analysis in International Rice Research Institute (IRRI) Farm, Philippines. Lobell and Ortiz-Monasterio (2007) analyzed a 15 to 24-year dataset to determined variation of regional spring wheat yields and temperatures for three irrigated sites in western North America and found that in San Luis-Mexicali Valley of Mexico and in the Imperial Valley of California, historical wheat yields were associated with minimum and maximum temperature but projected responses to maximum temperature were three times larger than minimum temperature. Wheat yields in these regions were reduced by 8% to 10% for 1˚C rise of minimum or maximum temperature based on CERES model simulations.

Unlike in the study of Peng et al. (2004), the magnitude of air temperature increase in this study was based on average night air temperature. Long-term trend analysis for minimum air temperature was not performed in this study due to minimum air temperature data were inconsistently recorded in many weather stations prior to 1980. Generally, the similar increase in air temperature trends observed in this study and other parts of the globe implicitly assumed that temperature changes occurred in minimum and mean night air temperature in rice growing regions. These results may reflect a strong greenhouse gas-induced warming during the last half of the 20^th century. The shift in seasonal night-time air temperature was less evident in the Southern states than in CA (Figures 4(a)-(d)). In this case, a possible reason for the small temperature increase is the small number of observations in Southern rice counties relative to observations in CA. The total number of observations was 3 to 6 times lower in southern states than in CA. Another possible reason for the minimal temperature rise is likely due to the weak warming trend of surface air temperature in the Southern region caused by changes in sea surface temperature. Wang et al. (2009) reported that climate trend analysis over the Southern US from 1950 to 2000 showed a slight to no increase in air temperature pattern during late spring and summer seasons and trends were attributed to the decadal climate variability pattern in the Pacific and Atlantic oceans.

Results of this paper confirm the observed warming of global surface air temperature during the last part of the 20^th century (IPCC, 2014; Ray et al. 2015; Hansen et al. 2010; Easterling et al., 1997). According to IPCC (2014), the expected changes in surface air temperature by the end of the 21^th century are predicted to be in the range of 0.3˚C to 4.8˚C. Studies on crop growth models showed that predictions from climate models of significant warming of air temperature might lead to crop yield reductions (Rosenzweig & Parry, 1994; Wheeler et al., 2000; Lobell et al., 2011). Recent research assessing the potential impacts of high night air temperature on rice yield had focused on constant air temperature treatments that were set to produce significant outcomes (i.e. Tamaki et al., 1989; Cooper et al., 2008; Coast et al., 2019). These studies had, for the most part, treated each trait of rice yield potential without relation to actual changes in air temperature under local rice environments. Although these adaptive experiments provided better understanding of rice growth mechanism and the components of grain yield under a set of constant air temperature thresholds, results may have limited implications for potential future sustainable intensive crop production for a specific region. In order to achieve reliable field crop responses to heat stress, heat/temperature stress thresholds should be based on local temperature difference relative to actual ambient temperature in the region. Accurate representation of field growing temperature reliably identifies cause-effect relationships of rice response tolerance, productivity and heat stress. The information discussed in our paper is often overlooked in climate articles because most authors focused on the data analysis and the validity of mathematical models and algorithms used. With the vulnerability of agriculture to deterioration due to climate change, we need results that have direct use to crop growers, breeders, and researchers which will aid them in developing resilient cultivars and farming systems that will sustain food for current and future human population.

4. Conclusion

Results in this study showed seasonal night-time air temperature change occurred in major US rice producing states during the last 67-years. Average seasonal night-time air temperature has ranged from 0.0˚C to 1.12˚C with greater fluctuations in CA and AR. While many open-accessed weather databases completely provide climatological data, accuracy and reliability of modeled data to predict local weather patterns relied on capturing variability across specific locations and time. Without the presence of weather stations in major rice counties in the US particularly in the southern states, research on temperature effects or heat stress on crop response remains uncertain, and it may lead to a strong bias towards non-agricultural areas. This study justifies the need to establish more weather stations in agricultural areas with strong emphasis in the US Southern States where frequent and intense weather fluctuations are more likely to happen. Currently, crop growers in some parts of the US have started maintaining their own local weather stations. Local or federal support to implement more privately owned weather stations could potentially improve accuracy of weather monitoring for the region. This paper highlights the observed warming of surface air temperature during the last part of the 20^th century. Although the overall results of our study occur in relatively small portion of the US agricultural areas, they strongly support the realization that warming of seasonal night air temperature will likely continue in the future. Our study provides direct evidence that air temperature changes varied in specific crop region and representing a constant and estimated growth temperature value to crop-specific region may not accurately and reliably reflect actual crop growing environments.

Acknowledgements

This work is funded by the Wheat and Rice Center for Heat Resilience (WRCHR; http://wrchr.org/) National Science Foundation Award Number #1736192. Additional support was provided by the Arkansas Biosciences Institute, Arkansas State University and the Office of National Programs through the HQ 0500 Air Quality funds United States Department of Agriculture, Agricultural Research Services. We also thank Dr. Ahmed Hashem of Arkansas State University for his assistance in mapping rice counties in the US.