Climatic Trends of West Virginia: A Representative Appalachian Microcosm

During the late 19th and very early 20th centuries widespread deforestation occurred across the Appalachian region, USA. However, since the early 20th century, land cover rapidly changed from predominantly agricultural land use (72%; 1909) to forest. West Virginia (WV) is now the USA’s third most forested state by area (79%; 1989–present). It is well understood that land cover alterations feedback on climate with important implications for ecology, water resources, and watershed management. However, the spatiotemporal distribution of climatic changes during reforestation in WV remains unclear. To fill this knowledge gap, daily maximum temperature, minimum temperature, and precipitation data were acquired for eighteen observation sites with long periods of record (POR; ≥77 years). Results indicate an increasingly wet and temperate WV climate characterized by warming summertime minimum temperatures, cooling maximum temperatures year-round, and increased annual precipitation that accelerated during the second half (1959–2016) of the POR. Trends are elevation dependent and may be accelerating due to local to regional ecohydrological feedbacks including increasing forest age and density, changing forest species composition, and increasing globally averaged atmospheric moisture. Furthermore, results imply that excessive wetness may become the primary ecosystem stressor associated with climate change in the USA’s rugged and flood prone Appalachian region. The Appalachian region’s physiographic complexity and history of widespread land use changes makes climatic changes particularly dynamic. Therefore, mechanistic understanding of microto mesoscale climate changes is imperative to better inform decision makers and ensure preservation of the region’s rich natural resources.


Introduction
During the late 19th to early 20th century West Virginia (WV) was the frontier of American industrial capitalism where the rich natural resources of the Mountain State, especially timber and coal, were exploited for the benefit of aggressively populating eastern cities [1,2]. In 1880, it was estimated that two thirds of WV was covered by ancient growth hardwood forest [2], but by 1909 approximately 72% of WV was agriculture or pasture lands [3]. Widespread logging created a landscape characterized by vast acreages of dried slash that were easily ignited by passing locomotives [3] and burned with unusual intensity during dry periods [4]. By the 1920s, most lumbermen abandoned WV for lands further south and west [2] and extensive fire prevention efforts resulted in the redevelopment of vast hardwood forests [4]. When the first United States (US) Forest Service forest survey of WV was completed in 1949 an estimated 64% of WV lands were forested, which increased to 75% by the

Materials and Methods
Daily summaries of total precipitation, maximum temperature, and minimum temperature were acquired from the National Center for Environmental Information (NCEI), a branch of the National Oceanic and Atmospheric Administration (NOAA; [17]). Stations were selected with start dates as early as 1900, but no later than 1930, with a nearly continuous time series through the end of 2016. These temporal criteria ensured analyses included the rapid transition from primarily agricultural and pasture lands in the early 20th century to primarily forested land coverage by the middle 20th century and subsequent forest maturation [3,6]. All leap days (i.e., February 29th) were excluded from analyses and daily data were post-processed to ensure all missing dates were included resulting in a continuous time series with 365 days in each year. Annual averages of daily summaries of maximum temperature, minimum temperature, and total precipitation were estimated and counts of missing data points were tabulated annually for each of the 18 selected stations. The spatial distribution of the 18 selected stations is shown in Figure 1 and additional site-specific data are included in Table 2. It was assumed that observers and their equipment followed station siting guidelines established by the United States Weather Bureau in 1890 (Fiebrich 2009). The World Meteorological Organization's guidelines on missing data in the calculation of climatological normals recommends not calculating monthly averages when five or more daily values are missing [18], which scales to 60 or more daily values for annual averages. For the purposes of this work years with data gaps exceeding 15% of possible daily observations (>54 daily observations) were removed from analyses. Of 117 possible annual averages or totals, each observation location had between 77 and 110 annual averages for maximum and minimum temperature and between 84 and 110 annual precipitation totals.

Materials and Methods
Daily summaries of total precipitation, maximum temperature, and minimum temperature were acquired from the National Center for Environmental Information (NCEI), a branch of the National Oceanic and Atmospheric Administration (NOAA; [17]). Stations were selected with start dates as early as 1900, but no later than 1930, with a nearly continuous time series through the end of 2016. These temporal criteria ensured analyses included the rapid transition from primarily agricultural and pasture lands in the early 20th century to primarily forested land coverage by the middle 20th century and subsequent forest maturation [3,6]. All leap days (i.e., February 29th) were excluded from analyses and daily data were post-processed to ensure all missing dates were included resulting in a continuous time series with 365 days in each year. Annual averages of daily summaries of maximum temperature, minimum temperature, and total precipitation were estimated and counts of missing data points were tabulated annually for each of the 18 selected stations. The spatial distribution of the 18 selected stations is shown in Figure 1 and additional site-specific data are included in Table 2. It was assumed that observers and their equipment followed station siting guidelines established by the United States Weather Bureau in 1890 (Fiebrich 2009). The World Meteorological Organization's guidelines on missing data in the calculation of climatological normals recommends not calculating monthly averages when five or more daily values are missing [18], which scales to 60 or more daily values for annual averages. For the purposes of this work years with data gaps exceeding 15% of possible daily observations (>54 daily observations) were removed from analyses. Of 117 possible annual averages or totals, each observation location had between 77 and 110 annual averages for maximum and minimum temperature and between 84 and 110 annual precipitation totals. Investigating temporal trends of the central tendency of a data series is of interest in the context of identifying a changing climate [20,21]. Sen's slope estimator paired with the Mann-Kendall trend test is a robust, non-parametric method to estimate the linear trend and assess statistical significance [22] and was thus deemed appropriate for this work. Spatial characteristics of temporal trends were estimated and statistical significance (α = 0.05) assessed over the entire time series (1900-2016) and during the first  and second (1959-2016) halves of each time series at each observation location. Analyses for each half of the time series were performed because the data series was sufficiently long and the first half corresponded with reforestation whereas the second half Investigating temporal trends of the central tendency of a data series is of interest in the context of identifying a changing climate [20,21]. Sen's slope estimator paired with the Mann-Kendall trend test is a robust, non-parametric method to estimate the linear trend and assess statistical significance [22] and was thus deemed appropriate for this work. Spatial characteristics of temporal trends were estimated and statistical significance (α = 0.05) assessed over the entire time series  and during the first  and second (1959-2016) halves of each time series at each observation location. Analyses for each half of the time series were performed because the data series was sufficiently long and the first half corresponded with reforestation whereas the second half corresponded with forest maturation and globally averaged warming exceeding 0.65 • C [4][5][6]23]. Except for New Cumberland (Table 1, Figure 1), each observation location had a more complete time series of each variable during the second half (93.8%) rather than the first half (71.7%) of the POR suggesting early data gaps (1900-1930) could have influenced results. As a result, the data should not be considered stationary and readers should use caution when extrapolating results into the future. Changes in temporal characteristics (i.e., seasonality) were estimated for all 365 days using a continuous time series of daily data averaged across all eighteen locations for the entire POR and the second half of the POR (1959-2016). Seasonal analyses were not performed for the first half of the POR due to larger and more frequent data gaps. For plotting purposes, centered three-week moving averages of daily Sen's slope, upper and lower 95% confidence interval values were calculated to smooth day-to-day variability and better show seasonal changes of each variable. Centered moving averages (CMAs) represent a series of arithmetic means throughout a time series that is centered on the middle value (i.e., day 11 of the 21-day average). Three-week CMAs were selected because meteorological and astronomical definitions of seasonality are static in time and space and may be inadequate considering a swiftly changing climate [24]. Ultimately, the methods described in the current work are easily applied and provide valuable information about spatiotemporal changes in climate that are needed for effective decision making.

Climate during Period of Record
Monthly and daily averages of precipitation and maximum, minimum, and daily average temperatures are included in Table 1 and Figure 2, respectively. West Virginia's (WV) wettest month of the year was July when an average of 11.7 cm of precipitation was observed during the period of record (Table 1). July was also WV's warmest month of the year with a daily average temperature of 22.7 • C and average maximum and minimum temperatures of 29.4 • C and 16.0 • C, respectively. Alternatively, January was the coolest month with a daily average temperature of −0.3 • C and daily maximum and minimum temperatures of 5.2 • C and −5.8 • C, respectively. February, October, and November were similarly dry with each month averaging less than 7.5 cm of precipitation. At daily resolution, the warmest average maximum temperature was 29.9 • C on July 18th and the coolest average minimum temperature was −6.9 • C on January 28th ( Figure 2). Maximum daily average precipitation occurred on July 10th with 0.49 cm day −1 and minimum daily average precipitation occurred on November 21st with 0.14 cm day −1 (not shown). Between 1900 and 2016, WV's spatially averaged climate was cold with hot summers and no dry season (i.e., Köppen-Geiger Dfa; [25]).

Spatial Climatic Changes
The observation location, elevation, and POR for each selected site as well as trends in minimum temperature, maximum temperature, and total annual precipitation are summarized in Table 2 and shown spatially in Figure 3. Unless otherwise stated, estimated trends reported parenthetically represent trends averaged across observation sites with either positive or negative trends. Table 2. Eighteen observation stations in West Virginia, USA. Elevation (m), period of record (POR), and estimated linear slope for minimum temperature ( • C/decade), maximum temperature ( • C/decade), and total annual precipitation (cm/decade). Higher and lower elevation averages refer to above or below the median elevation (323 m), respectively. White SS = White Sulphur Springs, New Cumb. = New Cumberland, Elev = elevation, blue or brown shading = significant decrease (p < 0.05), red or green shading = significant increase (p < 0.05).

Spatial Climatic Changes
The observation location, elevation, and POR for each selected site as well as trends in minimum temperature, maximum temperature, and total annual precipitation are summarized in Table 2 and shown spatially in Figure 3. Unless otherwise stated, estimated trends reported parenthetically represent trends averaged across observation sites with either positive or negative trends.  Half of the eighteen observation sites in WV indicated a decreasing trend in minimum temperature (-0.05 °C/decade), while the other half indicated a greater increasing trend (0.08 °C/decade). Warming minimum temperatures were generally observed across the northern and southwestern portions of WV, but the two lowest elevation observation sites had the largest rate of warming (≥0.12 °C/decade) and trends were statistically significant (p < 0.00). Alternatively, thirteen observations sites spatially and topographically distributed across WV recorded decreasing trends (-0.14 °C/decade) in maximum temperatures. Three of the four largest negative trends (≤-0.20 °C/decade) were observed across north-central WV above the median station elevation (323 m) and all four of the largest negative trends were statistically significant (p < 0.00). Total annual precipitation increased (0.87 cm/decade) at twelve observation sites and decreased (-0.72 cm/decade) at the other six sites. The two sites with the largest increase (>1.5 cm/decade) and decrease (<-1 cm/decade) in annual precipitation were below and above the median elevation, respectively and all four trends were statistically significant (p < 0.03). In synthesis, annually averaged minimum temperatures increased more at lower elevations, maximum temperatures decreased more at higher elevations, and precipitation increased at lower elevations between 1900 and 2016 across WV.

Entire Time Series (1900-2016)
Half of the eighteen observation sites in WV indicated a decreasing trend in minimum temperature (-0.05 • C/decade), while the other half indicated a greater increasing trend (0.08 • C/decade). Warming minimum temperatures were generally observed across the northern and southwestern portions of WV, but the two lowest elevation observation sites had the largest rate of warming (≥0.12 • C/decade) and trends were statistically significant (p < 0.00). Alternatively, thirteen observations sites spatially and topographically distributed across WV recorded decreasing trends (-0.14 • C/decade) in maximum temperatures. Three of the four largest negative trends (≤-0.20 • C/decade) were observed across north-central WV above the median station elevation (323 m) and all four of the largest negative trends were statistically significant (p < 0.00). Total annual precipitation increased (0.87 cm/decade) at twelve observation sites and decreased (-0.72 cm/decade) at the other six sites. The two sites with the largest increase (>1.5 cm/decade) and decrease (<-1 cm/decade) in annual precipitation were below and above the median elevation, respectively and all four trends were statistically significant (p < 0.03). In synthesis, annually averaged minimum temperatures increased more at lower elevations, maximum temperatures decreased more at higher elevations, and precipitation increased at lower elevations between 1900 and 2016 across WV.

First Half (1900-1958)
Twelve observation sites were characterized by increasing minimum temperatures (0.17 • C/decade) and the remaining six sites had decreasing minimum temperatures (−0.12 • C/decade). Just two sites below the median elevation had decreasing minimum temperatures and the two largest decreases (≤0.14 • C/decade), both statistically significant (p < 0.04), were found at higher elevations (604 and 722 m). Six of the twelve sites with positive minimum temperature trends were statistically significant (p < 0.02) and four of the six sites, including the largest positive trend (0.40 • C/decade), were below the median elevation. Ten of the observation sites indicated warming maximum temperature trends (0.20 • C/decade), the remaining eight sites indicated cooling (−0.15 • C/decade). Half of sites with observed warming and cooling maximum temperatures were above or below the median elevation. Seven observation sites indicated increasing trends in total annual precipitation (1.85 cm/decade) and eleven indicated decreasing trends (−2.70 cm/decade). Three of the eleven sites with decreasing trends were significant (p < 0.03) and all were above the median elevation. In summary, annually averaged minimum temperatures increased particularly at lower elevations, maximum temperatures increased, and precipitation decreased in WV between 1900 and 1958.

Second Half (1959-2016)
Fourteen climate observation sites indicated increasing minimum temperatures (0.22 • C/decade) and the remaining four sites indicated decreasing minimum temperatures (−0.10 • C/decade). Eight observation sites indicated increasing trends in annually averaged maximum temperatures (0.15 • C/decade) and ten observation sites indicated decreasing trends (−0.17 • C/decade). Half of all sites with observed warming or cooling minimum and maximum temperatures were above or below the median elevation, but sites with warming minimum (maximum) temperatures warmed 0.04 • C/decade (0.05 • C/decade) slower at higher elevation stations. Additionally, sites with cooling minimum (maximum) temperatures cooled −0.10 • C/decade (−0.08 • C/decade) faster at lower (higher) elevations. All eighteen observation sites indicated increasing trends (2.50 cm/decade) in total annual precipitation, but lower elevation sites had trends 1.09 cm/decade larger than higher elevation sites. Additionally, statistically significant (p < 0.05) increases in total annual precipitation were confined to observation locations below 400 meters in elevation across the western half of WV. As a result, minimum temperatures increased particularly at lower elevations, maximum temperatures decreased particularly at higher elevations, and precipitation increased state-wide with significant (p < 0.05) increases at lower elevations across western WV between 1959 and 2016.

Entire Time Series (1900-2016)
Centered three-week moving averages (CMAs) of state-averaged, daily precipitation and temperature trends are indicative of seasonal trends in WV (Figure 4). Daily precipitation trends were generally small (92% between ±0.0005 cm year −1 ), except for increases during May, decreases during June, and increases during July when 12 of the 17 significant (p < 0.05) daily trends occurred (Figure 4a). All but four of the 23 significant daily trends (p < 0.05) in average daily temperature were negative and the significant positive trends occurred in June, August, and December (Figure 4b). Daily average temperature trends were most negative during January (−0.021 • C year −1 ), October (−0.008 • C year −1 ), and March (−0.008 • C year −1 ) whereas April (0.006 • C year −1 ) and December (0.003 • C year −1 ) were most positive. Nineteen of the 26 significant (p < 0.05) daily minimum temperature trends were positive, especially during June, July, and August when sixteen significantly positive trends occurred (p < 0.05; Figure 4c). All seven significantly decreasing (p < 0.05) daily minimum temperature trends occurred between January 22nd and April 8th. All significant daily maximum temperature trends (n = 73) were negative, but CMAs were briefly positive during April, late November, and early December (Figure 4d). Between May and October, WV's growing season, 79% of the significant (p < 0.05) negative maximum temperature trends occurred. In summary, daily precipitation trends were largest during May, June, and July; average temperatures decreased most during January; minimum temperatures increased most during June, July, and August; and maximum temperatures decreased year-round between 1900 and 2016.   . Centered three-week moving averages of (a) daily precipitation, (b) average temperature, (c) minimum temperature, and (d) maximum temperature trends averaged across all eighteen stations between 1900 and 2016. Grey shading represents the 95% confidence interval. Brown and green (red and blue) shading represents statistically significant (p < 0.05) daily precipitation (temperature) trends.

Second Half (1959-2016)
Between 1959 and 2016, 312 of the 365 (85%) daily precipitation trends were between ±0.001 cm year −1 (not shown) and twelve days (four negative; eight positive) were characterized by significant trends (p < 0.05; Figure 5a). Average temperatures cooled between late January and early March when five of the nine (56%) significantly negative (p < 0.05) daily trends occurred (Figure 5b). Average temperatures warmed at an average rate of 0.015 • C year −1 during June, July, and August (JJA) when eleven of the sixteen (69%) significantly positive (p < 0.05) daily trends occurred. Similarly, minimum temperatures warmed an average of 0.022 • C year −1 during JJA when thirteen of the seventeen (76%) significantly positive (p < 0.05) daily trends occurred (Figure 5c). Alternatively, maximum temperatures cooled an average of −0.018 • C year −1 between January 1st and July 1st when twenty of the twenty three (87%) significantly negative trends (p < 0.05) occurred. Between May and October, WV's growing season, average trends in maximum and minimum temperatures were −0.016 • C year −1 and 0.014 • C year −1 respectively. In summary, daily precipitation trends were largest between March and August; average temperatures increased most during JJA; minimum temperatures increased most during JJA; and maximum temperatures decreased between January 1st and July 1st between 1959 and 2016.
Water 2019, 11, x FOR PEER REVIEW 10 of 15 Figure 5. Centered three-week moving averages of (a) daily precipitation, (b) average temperature, (c) minimum temperature, and (d) maximum temperature trends averaged across all eighteen stations between 1959 and 2016. Grey shading represents the 95% confidence interval. Brown and green (red and blue) shading represents statistically significant (p < 0.05) daily precipitation (temperature) trends.

Discussion
Extensive LULC changes occurred across West Virginia (WV) and the broader Appalachian region between 1900 and 2016 that greatly contrast global LULC changes [26]. WV's dominant LULC transitioned from agricultural land (72%) to native forest cover (64%) between 1909 and 1949 [5]. Forest cover continued to increase in areal coverage (79%) and density (i.e., timber volume) through 2012 [6]. LULC changes influence biogeophysical regulation of climate associated with changes in albedo, turbulent fluxes of energy, mass, and momentum at the land-atmosphere interface, and the hydrologic cycle [26][27][28][29]. Simulations of coupled land-atmosphere interactions indicated that temperate forests are characterized by lower albedo and evapotranspiration values relative to Figure 5. Centered three-week moving averages of (a) daily precipitation, (b) average temperature, (c) minimum temperature, and (d) maximum temperature trends averaged across all eighteen stations between 1959 and 2016. Grey shading represents the 95% confidence interval. Brown and green (red and blue) shading represents statistically significant (p < 0.05) daily precipitation (temperature) trends.

Discussion
Extensive LULC changes occurred across West Virginia (WV) and the broader Appalachian region between 1900 and 2016 that greatly contrast global LULC changes [26]. WV's dominant LULC transitioned from agricultural land (72%) to native forest cover (64%) between 1909 and 1949 [5]. Forest cover continued to increase in areal coverage (79%) and density (i.e., timber volume) through 2012 [6]. LULC changes influence biogeophysical regulation of climate associated with changes in albedo, turbulent fluxes of energy, mass, and momentum at the land-atmosphere interface, and the hydrologic cycle [26][27][28][29]. Simulations of coupled land-atmosphere interactions indicated that temperate forests are characterized by lower albedo and evapotranspiration values relative to irrigated croplands resulting in significantly warmer temperatures [29]. However, near surface eddy-covariance observations were compared between forest cover and an adjacent, non-irrigated grass field in nearby North Carolina, USA [30]. Results of Reference [30] showed the cooling effect of greater bulk aerodynamic conductance in forests was about three times larger than warming attributable to albedo effects. Therefore, the rapid transition from non-irrigated agricultural to forested land uses in WV may have resulted in a net cooling effect.
Net cooling in WV was characterized by decreasing state-average maximum temperatures (−1.0 • C) and increasing minimum temperatures (+0.4 • C) between 1900 and 2016 [9], but spatial and seasonal changes were not addressed. At the regional scale, changes in maximum and minimum temperatures contrasted the reported warming of 0.64 • C and 0.94 • C across the broader Northeast US region that included WV, respectively [31]. However, temperature trends in WV were more consistent with the Southeast US region during the period of study where maximum and minimum temperatures increased 0.16 • C and 0.76 • C, respectively [31]. In WV, minimum temperatures increased more at lower elevation locations (Table 1) where old-growth hardwood forests were harvested, conceivably with relatively less effort than WV's rugged mountainous terrain [2]. Between 1900 and 2016, increasing daily minimum temperature trends were consistently positive during JJA when 84% of significant (p < 0.05) increasing daily trends occurred. Similarly, 76% of the seventeen significantly positive (p < 0.05) daily average trends between 1959 and 2016 occurred during JJA. Warming summertime minimum temperatures in WV's humid subtropical climate suggests warmer dew point temperatures that may alter dewfall dynamics [32] and contribute to WV's increasingly humid climate [9]. Maximum temperatures cooled significantly (p < 0.05) on 72 days between 1900 and 2016 ( Figure 4) and annual trends were larger at higher elevations between 1959 and 2016 ( Table 2) where forest cover may be less disturbed and more mature [3]. Cooling high-elevation maximum temperatures may be explained by topographically induced circulations focusing afternoon cloud development along high-elevation ridges [33] that could be more frequent in an increasingly humid climate [9]. Further investigation into the mechanism(s) responsible for differential temperature trends at low and high elevation locations is needed, but warmer minimum and cooler maximum temperatures are consistent with attenuation of the diurnal temperature range due to forest cover [34] or cloud cover [35]. In addition to LULC changes in WV, increased forest coverage and maturity, and agricultural intensification (i.e., fertilization and irrigation; [36]) upwind of WV (i.e., Midwestern US) likely contributed to WV's increasingly humid WV climate [37]. Additional contributing mechanisms include (but are not limited to) increases in globally averaged atmospheric moisture associated with warming air and sea surface temperatures [38] or regional precipitation recycling [39]. Irrespective of the forcing mechanism(s), increased cloud cover aligns with an increasingly wet, temperate, and humid WV climate [35] and suggests regional climate forcing associated with coupled land-atmosphere interactions exceeded radiative forcing in WV (i.e., greenhouse gases; [40]).
Extrapolating annual trends over the entire POR then averaging across all eighteen locations indicated total annual precipitation increased by 3.8%, which is consistent with the national average reported in the Fourth National Climate Assessment [31]. However, precipitation increased at all locations between 1959 and 2016 when linear trends averaged across all sites indicated a 14.2 cm (13.2%) increase in precipitation relative to the epoch average. This precipitation trend was more than triple the rate of increase during the entire POR suggesting precipitation trends may be accelerating as a result of global scale increases in water vapor content [41] and super Clausius-Clapeyron scaling (SCCS) of precipitation extremes [42]. Significant (p < 0.05) increasing precipitation trends between 1959 and 2016 were estimated at lower elevation sites across western WV where forest harvesting, regrowth, and maturation was widespread [2] implying greater feedbacks associated with LULC change. The methods documented in the current work could be used to address the need for data-driven assessments of feedbacks between climate and LULC change in many locations, globally [43]. Understanding of LULC feedbacks on precipitation are particularly important in WV where some of the most extreme rainfall accumulations in the world were observed at time scales less than 6 h [44]. Thus, WV may be particularly vulnerable to changes in the frequency and magnitude (i.e., SCCS) of precipitation extremes and the impact of LULC changes (e.g., reforestation) on precipitation extremes and SCCS should be investigated.
The combination of an increasingly wet and temperate climate should help to secure WV's water resource quantity, but water quality problems exist in WV including (but not limited to) LULC change [45], acid mine drainage [46], and pathogenic water contamination [47]. Additionally, WV's increasingly temperate climate combined with decreasing temperature variance [9] may prolong episodes of excessive soil wetness and humid conditions that increases the vulnerability of terrestrial ecosystems to pathogens and fungal-like oomycetes such as Phytophtora and Pythium [48]. Similarly, increasing precipitation variance [9] indicates wet years are becoming more extreme, which may exacerbate vulnerabilities associated with pests and pathogens [49][50][51][52][53]. As a result, further investigation into changing climate variability may help to better explain observed changes in ecosystem biodiversity, carbon dynamics, and hydrology at a range of spatiotemporal scales [54]. However, increasing average precipitation (i.e., freshwater availability) signals potential for increasing agricultural productivity to ameliorate the food desert crisis in WV through conservation practices sponsored by the United States Department of Agriculture's (USDA) Environmental Quality Incentives Program (EQIP; [55]). In particular EQIP's High Tunnel System Initiative allows farmers to more efficiently deliver water and nutrients to plants, better protect plants from pests and pathogens, and minimize LULC change by increasing productivity per unit land area. Ultimately, climatic changes in WV may be distinct from the Northeast US climate region and Appalachia's complex topography may therefore necessitate refined climatic regions based on climate trends observed over a long POR and physiographic properties [56], including (but not limited to) elevation, slope and aspect, historic LULC changes, and forest age and composition.

Conclusions
Long term changes in West Virginia's LULC associated with widespread deforestation in the early 20th century, rapid forest regrowth through the mid-20th century, and subsequent forest maturation resulted in climatic changes distinct from the broader Northeast US climate region. Overall climatic trends averaged across the state suggest an increasingly temperate climate with increasing summertime minimum temperatures, particularly at lower elevations, and decreasing maximum temperatures year-round, particularly at higher elevations. Total annual precipitation is increasing and trends appear to be accelerating, particularly across the lower elevations of western WV where forest regrowth and maturation was particularly widespread. The mechanisms forcing observed changes in WV climate require much more investigation. This article therefore serves as an alert for greatly needed climate investigations in the Appalachian region to provide science-based guidance for future land-use/ land management activities in the region. This is critical since a wetter and more temperate climate may improve some ecosystems services (i.e., water security) while others may become increasingly stressed by increasing vulnerability to invasive pests and pathogens. Additionally, observed climatic changes support increasing agricultural productivity to address the food desert crisis through sustainable agricultural practices that minimize LULC changes (i.e., high tunnel systems). Therefore, climate observations over a long period of record should be used to identify spatiotemporal patterns of dynamic climatic changes at local to regional scales to better understand climatic changes, inform decision makers, and ensure preservation of WV's rich natural resources.