Abstract

Observations from the Oklahoma Mesonet and high resolution Weather Research and Forecasting model simulations were used to evaluate the effect that the dry line and large-scale atmospheric patterns had on drought evolution during 2011. Mesonet observations showed that a “dry” and “wet” pattern developed across Oklahoma due to anomalous atmospheric patterns. The location of the dry line varied due to this “dry” and “wet” pattern, with the average dry line location around 1.5° longitude further to the east than climatology. Model simulations were used to further quantify the impact of variable surface conditions on dry line evolution and convective initiation (CI) during April and May 2011. Specifically, soil moisture conditions were altered to depict “wet” and “dry” conditions across the domain by replacing the soil moisture values by each soil category’s porosity or wilting point value. Overall, the strength of the dry line boundary, its position, and subsequent CI were dependent on the modification of soil moisture. The simulations demonstrated that modifying soil moisture impacted the nature of the dry line and showed that soil moisture conditions during the first half of the warm season modified the dry line pattern and influenced the evolution and perpetuation of drought over Oklahoma.

1. Introduction

The Great Plains (GP) of the United States (US) is a region marked by periods of drought at multiple spatial and temporal scales ([1]; Christian et al., 2015). While drought impacts numerous regions of the United States [2], a lack of measurable precipitation in the GP, extended periods of anomalously high temperatures, misuse of natural resources, or a combination of the above yields vast impacts on the landscape and the ecosystem in the GP. For example, during the 1930s, intense decadal drought [3, 4] produced catastrophic socioeconomic impacts that led to substantial social migration out of the region due to the loss of crops and resources vital for the survival of livestock [5]. Even though current agricultural practices have improved and soil conservation practices have significantly reduced erosion in the GP, the region remains markedly susceptible to drought conditions.

Many contributing factors drive drought evolution in the GP including persistent synoptic patterns (e.g., Namais, 1983, and [69]), teleconnection patterns and sea surface temperature anomalies (e.g., [1016]), and surface-atmosphere coupling (e.g., [12, 17, 18]). Further, drought can persist at temporal scales spanning several years to decades, such as the “Dust Bowl” of the 1930s and intense drought during the 1950s, while also developing rapidly during the warm season over periods as short as a few weeks [1921].

As drought develops and evolves, surface characteristics change due to limited precipitation. For example, during the warm season, soils quickly dry due to increased evapotranspiration rates and vegetation may reach the wilting point without additional rainfall. Such vegetation changes impact the surface energy balance, as areas with less root-zone soil moisture partition the available net energy towards sensible heating while latent heat flux decreases. The exchange of heat and moisture between the surface and the atmosphere further impacts the structure and development of the planetary boundary layer (PBL [22, 23]) and such interactions between the surface conditions and the overlying atmosphere in the GP can modify atmospheric instability, convective precipitation, and regional precipitation coupling [2427].

At the same time, during the early warm season, precipitation and convective storm development in the Southern Great Plains (SGP) is influenced by the dry line [28, 29], a mesoscale boundary dividing dry, continental air to the west from moist, maritime air to the east (Hoch and Markowski, 2004, [28], and Schaefer, 1974a, 1974b). This is important for convective storm initiation during the early warm season and land-atmospheric coupling has been a focus of many dry line studies to deduce the effect of the land surface on dry line formation and development. Grasso [30] modeled the effect of differing soil moisture distributions on a specific dry line case and found that dry line propagation during the day is sensitive to changes in soil moisture. Further, when a constant soil moisture value was used over the domain, expansive dry line propagation did not occur. Ziegler et al. [31] modeled the sensitivity of dry line formation to differing soil moisture distributions using the Colorado State University Regional Atmospheric Modeling System (CSU-RAMS). The authors found that a west-to-east volumetric soil moisture gradient of 0.35 to 0.50 m3/m3 over 50 km was sufficient to initiate dry line formation, while a dry line was not produced when volumetric soil moisture was held constant at a value of either 0.35 or 0.50 m3/m3. The authors also found a connection between surface heat flux, boundary layer growth, and the strength of the surface moisture boundary. With less surface heat flux (their short grass and moderate soil moisture cases), Ziegler et al. [31] found that boundary layer growth was suppressed which further resulted in a weaker dry line than the control case. In addition, the study found that low-level convergence was less pronounced in these cases and that dry line formation was inhibited. Shaw et al. [32] showed that to accurately predict the dry line, the dry line morphology must be accurately tied to soil moisture and vegetation variability. The authors mention that a reduction of soil moisture over the entire domain caused an increase of sensible heat fluxes and a reduction of the west-to-east virtual temperature gradient. This allowed the moist air ahead of the surface boundary to mix more efficiently and thus the moisture gradient propagated further to the east with a slightly faster eastward movement. Shaw et al. [32] also theorized that moist soils ahead of the dry line could play a role in the maintenance of the strong moisture gradient at the boundary through faster transport of local moisture towards the dry line moisture gradient as opposed to transport from a source further away (e.g., the Gulf of Mexico). The authors conclude by noting that a fundamental relationship between the dry line and land surface conditions exists which yields an interesting question of cause and effect. For example, if this feature is influenced by the surface, how does the dry line’s influence on surface conditions further affect the dry line itself?

During 2011, intense drought developed across the SGP [33], which yielded significant drought-related impacts across Texas and Oklahoma [34]. However, as drought evolved across Oklahoma an east-west gradient in drought severity developed in which areas to the east received significant precipitation while areas to the west received limited precipitation. As a result, a strong gradient in soil moisture conditions also developed across the region.

While the complex interactions between the dry line and the land surface have been extensively investigated, the connections between drought development, the dry line, and land-atmosphere interactions have not. Thus, this study examines the interrelationships between the dry line, convective initiation (CI), precipitation, and the evolution of drought across Oklahoma during 2011 using in situ observations and numerical simulations via the Weather Research and Forecasting Advanced Research (WRF-ARW) model. The overall goal was to determine the role of the dry line on surface conditions during the progression of drought and how these conditions influenced the dry line during the warm season.

2. Methodology

2.1. Oklahoma Mesonet

Oklahoma Mesonet (hereafter called Mesonet) data [35, 36] was used to analyze the soil and near-surface atmospheric conditions over the period spanning March through August of 2011. The Mesonet is an automated network of 121 stations (Figure 1) installed across the state of Oklahoma that collects observations of soil and atmospheric variables at 5–30-minute intervals. Each observing station includes a datalogger, solar panel, radio transceiver, lightning rod, and environmental sensors attached to, or surrounding, a 10 m tower [36]. To minimize the impacts of bare soil and fast growing vegetation, the stations are placed in areas with uniform, low-growing vegetation. The sensor layout is designed to obtain the highest quality data possible [37].

To gain a better understanding of the drought evolution, 2 m temperature and specific humidity and 10-m wind speed along with precipitation totals and soil moisture values collected by the Mesonet were analyzed. Anomalies were calculated from a station long-term average dataset spanning 1999–2014 if at least 80% of the observations from a given variable, period, and station were recorded in the dataset. This is especially important for the soil moisture anomalies given that the stations changed from 2000 to 2007, with 6 sites upgraded and 5 decommissioned [38]. The collected data were analyzed in terms of monthly averages. Daily averages and anomalies were investigated; however the monthly anomalies showed the same overall trends as the daily values and thus were used. A key point of caution in this analysis involves the use of monthly 5 cm soil moisture anomalies, as shallow soil moisture can change rapidly due to discrete recharge events [39].

2.1.1. Soil Moisture Data

Soil moisture is measured by the Campbell Scientific 229-L heat dissipation sensor at depths of 5, 25, and 60 cm [36, 38, 41]. Soil moisture measurements are not available at all stations at each depth due to site-specific limitations (e.g., shallow bedrock) and are measured to a depth of 5 cm at 108 sites, 25 cm at 106 sites, and 60 cm at 82 sites as of 2011 [38]. The variable used for this study is the Fractional Water Index (FWI), which is a normalized index for soil moisture calculated via [38]where is reference temperature difference, = 3.96°C, and = 1.38°C. The sensor is buried at the depth specified [41] and measures the response to a pulse of heat caused by a 50 mA current (for 21 seconds) passing through the resistor located next to the thermocouple itself located within a ceramic matrix [38]. The difference between the soil temperature after the pulse and the temperature before the pulse is used to derive , by using a linear regression to normalize the response of the individual sensor to an idealized sensor with a maximum of 3.96°C and a minimum of 1.38°C to derive sensor specific calibration coefficients, which are then used to calculate . This gives a maximum value of FWI at 1, no change in temperature between pulses ( = 1.38°C), and a minimum value of 0, a maximum change in temperature ( = 3.96°C). FWI is a useful quantitative measure because it is not impacted by differing soil types. It is also important to note that the majority of vegetation in Oklahoma will flourish when FWI is above 0.8, begin to wilt when FWI is around 0.5, and begin to perish when FWI values are less than 0.3 [40].

2.1.2. Specific Humidity Data

Specific humidity (kg kg−1) is calculated using where = 0.622, is pressure (hPa), and is vapor pressure (hPa) calculated from the Goff-Gratch [42] equation. Specific humidity is not measured directly, but rather computed from other measured quantities collected by the Mesonet (temperature, relative humidity, and air pressure).

2.2. Dry Line Classification

To determine the impact of the dry line on the environment, a dataset of dry line events was developed for March-August 2011. One of the defining features of Oklahoma is the east-to-west terrain height change. Use of dew point temperature and relative humidity would be impacted by this change of terrain height, so to remove this source of error the specific humidity field was used similar to that of Hoch and Markowski [43]. As such, a gradient of 3 g kg−1 (100 km)−1 in the specific humidity field determined using Mesonet data was required to identify a dry line.

This analysis determined the largest distance between any two points across the identified boundary that contained a 3 g kg−1 gradient and then matched it against the prescribed dry line qualification criteria to determine if the identified boundary was a dry line. Analysis was completed at 0000 UTC to compare results with the Hoch and Markowski [43] study. During the analysis several active dry line cases were identified that had been overtaken by surface fronts prior to 0000 UTC. To include these dry line cases they were identified at 2100 UTC and are noted as such in the results. Dry line longitudes were determined following the Hoch and Markowski [43] methodology. As such, longitudes were found by taking the longitude of the easternmost point on the furthest east specific humidity contour that was analyzed as part of the dry line moisture gradient.

To determine if CI occurred along any identified dry line boundary, remotely sensed observations were used. The National Center for Atmospheric Research (NCAR) provides an archive of weather data that is publicly available at the Warm Season Precipitation Episodes (WSPE: http://www2.mmm.ucar.edu/episodes/) website. The archive includes visible satellite images and composite radar data for each day from 1996 to present. Dry line location was determined via Mesonet data and satellite images were used to identify convective clouds along the boundary. The longitude of CI was determined by the furthest west convective cloud that developed on or near the dry line via manual inspection of satellite data, with convective cloud formation determined using archived visible band Geostationary Operational Environmental Satellite (GOES) data from the NCAR WSPE website. Three classes of dry lines were determined from this analysis:(1)A dry line that did not produce any convective clouds or precipitation(2)A dry line that produced convective clouds but no precipitation(3)A dry line that produced both convective clouds and precipitation

Production of precipitation was determined using Mesonet precipitation data in conjunction with the dry line position; however it is possible that dry line precipitation was not identified due to the limitations of the Mesonet observations (i.e., that a station was not in the area that received precipitation). The longitude of precipitation for each dry line event that produced precipitation was found by analyzing Mesonet data and determining the westernmost station that received rainfall (greater than or equal to 1 mm) that could be traced back to storms that initiated along the dry line boundary. This was done so as not to artificially inflate the statistics that determine the furthest east station which received dry line produced precipitation as convection could occur some distance away from the actual dry line boundary. However, typical dry lines are quasi-linear boundaries [43] so with the method used to define the dry line longitude, rainfall longitudes may appear further to the west than the “mean” dry line that initiated the rainfall.

2.3. WRF-ARW Modeling Configuration

The WRF-ARW was used as the main investigation tool for the modeling portion of this study. A two-way nested grid configuration (Figure 2) was used with the outer domain covering 810 × 710 km with 10 km horizontal grid spacing, and the inner domain covering 592 × 452 km with 2 km horizontal grid spacing. Each domain had 114 vertical levels, with a model top of 50 hPa, using a mass-based terrain following vertical coordinate system. With regard to boundary layer vertical grid spacing a total of 13 vertical levels are below 2 km with a total of 18 vertical levels below 3 km, the typical boundary layer height found in most of the simulations. The fine vertical resolution in the PBL promotes a realistic simulation of processes important to dry line formation and evolution. Boundary conditions over both domains were set using the North American Mesoscale (NAM) Forecast System 12 km, 6 hourly update output data. Both domains utilized the Rapid Radiative Transfer Model (RRTM) longwave [44] and Dudhia [45] shortwave radiation schemes. The Unified NCEP/NCAR/AFWA Noah Land Surface Model was utilized to simulate surface conditions through all simulations with 4 soil layers and accounts for frozen soil and multilayer snow. The Yonsei University PBL Scheme (YSU [46]) was chosen given the results of Gibbs et al. [47] which showed that the YSU scheme more accurately reproduced the observed surface heat fluxes over portions of Oklahoma compared to a local closure scheme. The Kain-Fritsch [48] cumulus scheme was used on the outer domain whereas no cumulus scheme was used on the inner domain. The Morrison 2-moment microphysics scheme [49] was used over both domains. Comparison of simulations completed with different microphysics schemes in the inner domain showed that the Morrison double moment scheme outperformed the WDM-6 scheme [50] in the test case with regard to the location of precipitation (not shown).

2.3.1. Initialization of Surface Conditions and Sensitivity Analysis Design

To develop the initialization of soil conditions, the High Resolution Land Data Assimilation System (HRLDAS [51]) was used. This system is based on the Noah Land Surface Model (LSM [52, 53]) and the Noah LSM’s predecessor the Oregon State University (OSU) LSM [5457] which uses input data to build a more accurate background of soil moisture and temperature. To achieve this, a 5-year spin-up period (2006–2011) was completed using HRLDAS and North American Regional Reanalysis (NARR) data for initialization and forcing. This duration was chosen due to analyses by Monroe et al. 2014 who noted large differences between short spin-up (<2 years) simulations and simulations using a 5-year surface spin-up via HRLDAS. Cosgrove et al. [58] also noted that in the central Great Plains that spin-up for 49 months was required to reach a percent cutoff threshold (how long it takes for yearly changes in monthly averaged model output to decrease to a certain threshold) of 0.01%, which represents fine scale model equilibrium with root-zone soil moisture. Thus, a longer spin-up duration of 5 years was required to reach a more representative surface state.

To quantify the impact of surface conditions on the dry line several simulations were completed in which the initial soil moisture conditions were altered. Instead of using homogenous (wet or dry) soil moisture conditions [31, 32] or a fractional (−%, +%) change [30, 59] the sensitivity runs were initialized with soil moisture conditions based on the grid point soil type and maximum/minimum soil “wetness.” As such, the magnitude of soil moisture was replaced by the wilting point (Dry) and soil porosity (Moist) of each soil type at each grid point throughout the soil depth as determined by the WRF Postprocessing System based on United States Geological Survey (USGS) soil categories. A description of these values is available in Table 1. Physically, this type of modification is a more realistic approach than using a homogenous modification of the magnitude of soil moisture given the available water in the soil is largely determined by soil texture.

3. Drought and Dry Line Analysis

3.1. Drought Analysis

During the months of April and May 2011 drought evolved in a distinct west-to-east pattern over the state of Oklahoma (Figure 3). Drought was established in southern Oklahoma at the beginning of April (Figure 3(a)) and conditions during April and May modified the evolution of drought towards a more east-west gradient of drought conditions (Figure 3(b)). To determine the surface environmental conditions that led to this drought evolution, Mesonet observations were analyzed. Precipitation was anomalously high in eastern Oklahoma during April and May, with western Oklahoma receiving much less rainfall (Figure 4). The precipitation anomalies depict the “dry and wet” pattern that occurred during the two months. This result was also depicted by Su and Dickinson [60] who modeled the SGP region in order to understand the spatial patterns of soil moisture-precipitation feedback during the month of April 2011. They found that atmospheric conditions were such that deep convection was inhibited over western Oklahoma and enhanced by moisture convergence in eastern Oklahoma. However, the main driver for the anomalous rainfall in April was due to synoptic features rather than the dry line. In contrast, the anomalous rainfall that occurred in May (4/9 rainfall events driven by dry lines) was produced more so by dry line events when compared to April (2/11 rainfall events driven by dry lines). This pattern of precipitation over the state provided the conditions that led to easing of drought conditions across eastern Oklahoma and drought intensification over western Oklahoma. Climatologically, the precipitation that normally occurs in this period is abundant over eastern Oklahoma while being crucial in western Oklahoma to stave off drought development as the warm season progresses. Illston et al. [40] noted that during phase II (transitional drying phase) the large fluctuations of FWI that are occurring, especially in western Oklahoma, are due to the variable rainfall that occurs in the region and recharges the soil moisture (Figure 5). As such, precipitation from March to June is important to the evolution of early warm season drought in Oklahoma. This is in contrast to the other phases (I, plateau phase; III, enhanced drying phase; IV, recharge phase), which show a more consistent pattern. Thus, with little precipitation over western Oklahoma during these months, the drought continued to strengthen, while the copious amounts of precipitation in eastern Oklahoma eased drought conditions there. An area of interest during May is the region of anomalously high precipitation in west central Oklahoma (Figure 4(b)). These positive anomalies were caused by isolated dry line convection and other surface boundary events. An easing of drought occurred due to these rainfall events but overall did not have a significant impact on the overall evolution of drought conditions.

The surface temperature anomaly patterns were not similar in the months of April and May (Figure 6). In April, temperatures were largely above normal across the entire state. These conditions aided the development of drought across western Oklahoma as high temperatures along with a lack of precipitation dried out the surface. In eastern Oklahoma, the impact of these higher than average temperatures was lessened due to the copious amounts of rainfall. In May, however, a different pattern emerged. Temperature anomalies over central and eastern Oklahoma were below normal, by 0.2 to 1.3°C, while western Oklahoma was above normal by 0.2 to 0.9°C. This pattern would further enhance drought in extreme western Oklahoma while aiding in the easing of drought in the eastern portions of the state.

Analysis of soil moisture anomalies during April and May (Figure 7) depicts the effect that the anomalous precipitation and temperature patterns had on the surface. During April, the 5 cm FWI anomalies show near average (−0.25 to 0.10) soil moisture in eastern Oklahoma, with below average (0.30 to 0.50) values in western Oklahoma. Even with the anomalous precipitation in eastern Oklahoma soil moisture values were not substantially above normal. This is likely due to the anomalous temperature pattern seen over Oklahoma that would increase evapotranspiration (ET) and thus decrease the effect the precipitation would have on soil moisture, along with increased surface runoff due to the excessive amounts of rainfall. In May, a different pattern emerged, with 5 cm FWI anomalies above normal over central and eastern Oklahoma while remaining below normal in western Oklahoma. While large areas of positive anomalies in central and eastern Oklahoma were due to the above average precipitation, areas of positive 5 cm FWI anomalies in central Oklahoma cannot be attributed to above average rainfall. This likely shows that precipitation was percolating deeper into the soil more effectively than would occur during a typical month. This could be due to different soil characteristics in this area or by decreased evaporation in this region of lower rainfall anomalies. With lower surface temperatures in this portion of the state at the time it can be seen that evapotranspiration would be reduced and thus soil moisture values would be increased. While soil moisture observations are collected at 5, 25, and 60 cm, only the 5 cm FWI anomalies were included due to analysis of the other depths showing similar trends during these two months.

3.2. Dry Line Analysis

From March to August a total of 36 dry lines were analyzed with four occurring in March, six in April, seventeen in May, eight in June, zero in July, and one in August. The approximate location of each of these dry lines is shown in Figure 8. This figure indicates two areas of preferred dry line locations, over the central to west central portion of the state and in the northwest portion of the state into the panhandle. Coffer et al. [61] mentioned in their analysis of model dry line position errors that 2011 was an active year with frequent passages of shortwave troughs across the central Plains. Hane [62] states that dry lines that form in active synoptic regimes propagate further to the east during the day than those that develop during quiescent periods because the increased westerly momentum associated with the synoptic disturbance aids the vertical turbulent mixing in shifting the boundary towards the east. While it is difficult to remove the influence of the synoptic pattern on dry line location, a distinction between dry line synoptic classifications can be made in order to gain better understanding of the influence the drought had on dry line location.

Results of this analysis are shown in Table 2. The results revealed that the average dry line position was 2.5° (2.4°) longitude further to the east during April (May) when compared to climatology. Even with this eastward shift in the dry line, the overall seasonal progression of the average dry line location towards the west was similar to climatology. Linear trend analyses for each type of dry line show that dry lines with or without precipitation transitioned towards the west through the 3-month period, while dry lines with CI and no precipitation did not exhibit any discernable shift in average longitude through the season. With respect to differing dry line synoptic classifications, active dry lines were on average located further to the west than quiescent dry lines at 0000 UTC (Table 2). During April, active dry lines were on average located at 97.7°W while quiescent dry lines were located at 96.9°W. This trend continued into May with the active dry lines average longitude at 99.0°W and the quiescent dry lines at 98.3°W. This is counter to Hane [62] which mentions that active dry lines propagate further to the east than quiescent dry lines during the daily cycle. These results were not expected, as propagation of synoptically active dry lines is due to surface conditions and the synoptic flow pattern, with propagation of quiescent dry lines due to surface conditions alone. However, with synoptically active dry lines being influenced by the frontal boundary trailing beyond the moisture gradient these boundaries would have a finite distance they could propagate before being overtaken by the surface front, while quiescent dry lines would not have that limitation. Thus, with surface conditions conducive to more rapid eastward propagation when compared to normal, it would not be unexpected to see quiescent dry lines at 0000 UTC being further to the east on average compared to synoptically active dry lines.

To determine if land-atmosphere interactions played a role in dry line CI, average longitudes of dry lines with and without coincident precipitation were compared. Through modification of the prestorm environment, surface conditions have a large influence on whether CI can occur. Ford et al. [63] found that storms over Oklahoma preferentially initiate over drier soils compared with wet soils but further mentioned that storm initiation is not determined solely by soil moisture as storms were analyzed to initiate above wetter soils as well. Recalling the “dry and wet” pattern displayed in the observations (Figures 4 and 7), dry lines located further to the east would be located closer to areas in which soils were wetter than average. With dry lines being closer to these wetter areas their results would suggest that CI would be less likely to occur. However, a wetter than average soil also allows for higher convective available potential energy (CAPE) due to increased low-level moisture along with a decrease in convective inhibition (CIN), the level of free convection (LFC), and surface temperature [63]. Thus, while the potential for deep convection is increased due to higher CAPE, lower CIN, and a lower LFC, the potential for deep convection is also limited because surface convective forcing is diminished when surface temperatures decrease over wetter soils due to enhanced evapotranspiration.

The average dry line longitudes for precipitating dry lines indicate that this “dry and wet” pattern influenced dry line CI. Average dry lines that produced precipitation during April and May were approximately 1.5° (1.5° and 1.4°, resp.) longitude further to the east than dry lines that produced CI but did not produce precipitation at any Mesonet Stations. However, analysis on the longitude of CI and precipitation due to the dry line boundaries tends to disagree. Average longitudes for dry line produced CI (April 98.6°W; May 99.3°W) were further to the west than the average dry line longitude (April 97.3°W; May 98.6°W). This shows that CI that occurred due to the dry line boundaries originated behind the dry line boundaries and thus was unlikely to be directly caused by the surface conditions gradient owing to the “dry and wet” pattern. Dry line produced precipitation, however, agrees with the hypothesis that the pattern impacted the boundaries. Average rainfall longitudes were further to the east or at the monthly average dry line longitude. For April, dry line produced precipitation on average was at 97.5°W while during May this was at 98.1°W. This appears to show that rainfall was on average behind the dry line in April, but this was an effect of the few number of precipitating dry lines in April and the manner in which rainfall longitude for each event was determined.

When taking into account synoptic influence, an important result is seen. During the month of May active dry lines (99.0°W) were further to the west than quiescent dry lines (98.3°W) on average (Table 2). This indicates that even with synoptic forcing dry lines that were influenced by the surface alone developed and propagated further to the east. Precipitation from these dry lines was influenced as well, as rainfall that occurred from quiescent dry lines was 1° longitude further to the east than rainfall initiated by active dry lines. Conversely, during April, the “dry and wet” pattern originated due to anomalous rainfall caused by a large number of synoptic waves along with the few precipitating dry lines. Thus, the influence of the pattern would not have been significant during this month.

4. WRF-ARW Modeling Case Studies and Sensitivity Simulations

From observations, an apparent shift of the dry line boundary towards eastern Oklahoma was seen in conjunction with a soil moisture anomaly gradient over the central portion of the state. To quantify the effect of the soil moisture gradient on specific dry line cases during this season, simulations were completed using control and modified soil moisture fields. Using the WRF-ARW modeling system, four case studies were investigated, two from April (08 and 14 April) and two from May (08 and 21 May). To quantify the effect that the surface (i.e., soil moisture) had on the dry line boundary under varying synoptic conditions, two of these cases were synoptically forced (14 April and 21 May) while the other two occurred during quiescent periods (08 April and 08 May). Simulations using HRLDAS soil moisture values (Control) well represented each case with minor differences when compared to Mesonet observations. Dry line boundaries were shifted slightly to the east (in Control) compared to observations; however this is in line with a study by Coffer et al. [61] that found an eastward bias of the dry line in the operational NSSL-WRF system. Surface fields were also well represented by the Control simulations, with each simulation incurring small perturbations to the surface temperature and moisture fields. Simulated reflectivity and precipitation totals showed that the WRF-ARW was able to accurately depict the nature of each dry line event, from CI to the later portions of storm evolution.

Overall, modifying soil moisture in the WRF simulations had a definitive impact on dry line structure and subsequent CI. The Dry simulations were drier and warmer than the Control with a dry line that was shifted to the east (Figures 9, 10, 11, and 12). The PBL produced in a Dry simulation was typically deeper than the Control and much drier while CI would occur earlier in time than in the Control simulations. For the Moist simulations, the surface moisture was greater with cooler surface temperatures. The dry line boundary was typically shifted to the west (Figures 9, 10, 11, and 12), and like the Dry cases more diffuse than Control. Further, the PBL depth was shallower, with more moisture within the layer which further produced lower LCL heights. Of note, the depth of the PBL with relation to LCL heights appeared to play a significant role in determining whether CI occurred (April 8, May 8, and May 21 Moist simulations), along with the effect the cooler near-surface temperatures would have in stabilizing the environment.

Because the focus of this modeling study was to quantify the impact of the existing soil moisture regime (“dry and wet” regime) on the dry line boundary, the location of the dry line boundary was the primary focus within the case studies. The modification of soil moisture to the May dry line cases (Figures 9 and 10) produced more significant changes to the dry line boundary than the same changes made to the April dry line cases (Figures 11 and 12). The removal of soil moisture on the May 21 dry line event removed the dry line boundary and changed the nature of the event altogether. This can be seen in Figure 9 as the Dry simulation proxy dry line boundary (10 g/kg specific humidity contour; value chosen because this contour well represented the dry line boundary in the Control simulation) was shifted nearly outside of the domain. Adding soil moisture over the domain in the May 21 dry line case caused a large shift in the dry line boundary towards the west (Figure 9). The addition of soil moisture to the May 8 dry line case shifted the dry line from central Oklahoma to the Oklahoma/Texas panhandle border (Figure 10). Further, the addition of soil moisture in this case caused a large increase in atmospheric moisture in the PBL over western Oklahoma that was not seen in Control as soil moisture conditions were particularly dry over western Oklahoma at that time. The removal of soil moisture on the May 8 dry line was also significant. The structure of the surface moisture gradient was similar to that of Control; however drying caused the shifting of this gradient towards lower values of specific humidity and hence the difference seen in Figure 10. The changes to the PBL were typical for a Dry simulation, with an increase of PBL height and overall being drier than the Moist simulation. In the April 08 and April 14 simulations (Figures 11 and 12), these shifts in the dry line boundary and the changes in structure noted in the May simulations did not occur. The modifications to soil moisture did have an impact on dry line location and structure in the April 08 simulation (Figure 11), an expected result given that it was a quiescent dry line. However, the impacts were not as significant as with the May 08 dry line simulations, the other simulated quiescent dry line. The April 14 simulated dry line boundary did incur shifts to the east (DRY) and west (MOIST); however they did not change the surface structure of the dry line boundary (Figure 12). These results show that land-atmosphere interactions were more impactful on May dry lines (Figures 9 and 10) than on April dry lines (Figures 11 and 12) within our WRF simulations. While synoptic influences on dry line location and evolution were similar in April and May, the impact the surface imparted on dry line cases during May was increased compared to the influences the surface had on the April dry lines.

5. Discussion and Conclusions

This study was completed to investigate the role of the dry line in the evolution of the drought over Oklahoma during 2011. An analysis of in situ observations from the Oklahoma Mesonet provided background knowledge into the drought intensity and evolution along with dry line characteristics during the period. Using the observations, a criteria based dry line identification method was designed similar to that of the Hoch and Markowski [43] dry line identification method to develop a dataset of dry line statistics. Lastly, a modeling study using the WRF system was used to quantify the impact of soil conditions on the dry line in select cases with the goal of determining the effect of varying soil conditions on the dry line during 2011 and modifications to the dry line pattern and subsequent precipitation.

Analysis of Mesonet observations during the period provided several key results about the evolution of the 2011 drought. During the spring, conditions were ideal for drought to progress from southern Oklahoma into the northern portions of the state yet during April precipitation (Figure 4(a)) in eastern Oklahoma eased drought conditions. Su and Dickinson [60] found that synoptic conditions were favorable for CI to occur in eastern Oklahoma while being inhibited in western Oklahoma during the month of April. Further, during May, conditions were such that drought did not progress in eastern Oklahoma due to below normal temperatures (Figure 6(b)) and above normal precipitation (Figure 4(b)). This anomalous precipitation during May in eastern and portions of west central Oklahoma was partially due to dry line events (8 out of 17 rainfall events). However, in western Oklahoma extreme drought continued to intensify during April and May (Figure 4). As such, the intensification of drought in eastern Oklahoma was delayed from April to June while across western Oklahoma drought continued to intensify and subsequently produced a “dry and wet” pattern (dry in the west, wet in the east) during this period.

Given the “dry and wet” pattern, the role of the surface conditions due, in part to drought, on dry line development was investigated. Using the Hoch and Markowski [43] dry line climatology for comparison, results demonstrated that average dry line longitudes during the period were anomalously located towards the east (Table 2); average monthly dry line longitudes typically occur near 100°W, while during 2011 the investigation showed that average monthly dry line longitudes were near 98°W. Analysis of CI showed that dry lines in April and May which produced precipitating storms were located further to the east than dry lines that did not initiate any storms with precipitation. Thus, the dry line was influenced by the “dry and wet” pattern with a mean position located near the gradient between drought and nondrought conditions. The mean location of the numerous (6 during April, 17 during May) dry line events would modify surface conditions owing to the atmospheric pattern that defines a dry line event. West of these dry lines a drier air mass would be conducive for higher temperatures, enhanced ET, and greater overall atmospheric demand thus enhancing the existing drought conditions in these areas. Conversely, to the east of these boundaries enhanced moisture advection from the south and reduced vapor pressure deficit inhibited ET and limited drought development. Additionally, convective precipitation that develops along the dry line boundary propagates eastward and further limits drought development in this portion of the “dry and wet” pattern. As such, the results of this study show that there was a likely enhancement of the “dry and wet” pattern over central Oklahoma owing to the numerous dry line events that occurred within the month of May and the subsequent precipitation/ET pattern that these events caused.

Separating dry lines into synoptic classifications revealed an important result: synoptically active dry lines were on average further to the west than quiescent dry lines. This is especially true for the month of May in which quiescent dry lines that produced precipitation (as well as the associated precipitation and CI longitudes) were all further to the east when compared to active dry lines. As such, the more eastward placement of quiescent precipitating dry lines is important due to the fact that quiescent dry lines are more susceptible to surface conditions given that surface fluxes are the main driver for the turbulent mixing that propagate dry lines eastward [62]. Thus, these results show that surface conditions were impacting dry line position. It appears that when synoptic patterns influenced dry line boundaries during May that propagation towards the east was hindered compared to when surface conditions were the main driver of dry line evolution.

With dry line positions being much further to the east than the climatological normal, the meteorological patterns that are associated with these features would also be shifted towards the east. This would drive precipitation patterns associated with the dry line further to the east, along with the influences of the dry air mass that forms behind the dry line. Thus, the features impact surface conditions in such a way that would strengthen the “dry and wet” pattern over the state, as the dry air mass would continue to dry out the western extent of the state, while the moist air mass to the east would continue to provide the eastern portions of the state with lower temperatures, less surface evaporation, and a higher chance of deep convection and precipitation. As such, the pattern of drought that developed during late April and May was likely influenced by atmospheric conditions that developed due to the large number of dry line boundaries which were further to the east than the climatological normal.

A numerical modeling study into the effects of soil moisture conditions on specific dry line cases during the 2011 period was completed and four cases were chosen: two quiescent dry lines (April 8 and May 8) and two active dry lines of varying synoptic forcing strength (strongly forced April 14 and weakly forced May 21). The Control simulations demonstrated that the WRF produced reasonable atmospheric moisture gradients and precipitation patterns compared to observations. While differences were noted, these were typical to the WRF modeling system shown by previous studies [61]. Sensitivity simulations modifying soil moisture were completed to quantify the impacts of soil moisture on the dry line location during the study period. Results showed that during April these modifications led to typical results, with the addition or removal of soil moisture modifying the dry line events in such a way that did not influence the overall event that was simulated in the Control simulations. However, during May this was not the case. The modification of soil moisture in the WRF simulations led to changes compared to their Control simulations (Figures 9, 10, 11, and 12). In the Dry (Moist) simulation the dry line was shifted to the east (west) and had a different structure than that observed in Control. Thus, land-atmosphere interactions impacting dry line evolution were more significant during the month of May (Figures 9 and 10) than during April (Figures 11 and 12) within our WRF simulations. These extreme cases and the lack of such results from the April dry line simulations show there is a clear difference between the interactions between the dry line and surface conditions during these two months. Thus, the results show that if a different soil moisture pattern would have developed during April, then the May dry line would have had a different influence on the overall drought development. Shifting the dry line back towards the climatological normal (~100°W) would have shifted the moist air mass back towards the west and allowed for a higher chance of rainfall in the western extent of Oklahoma and aided in drought recovery, while shifting the dry line to the east would have reduced the rainfall over the entire state and would have aided in the development of the drought further to the east. However, the extent that drying over the entire state would impact overall dry line position is a difficult question in and of itself, as dry lines cannot propagate indefinitely to the east given extremely dry soil conditions over the entire region.

Overall the results show that surface conditions which developed during April 2011 affected dry line evolution during May. Because the dry line impacted precipitation development on the moist side of the wet-dry pattern, it further influenced the evolution of the drought during the months of April and May by abating the evolution of drought towards eastern Oklahoma until later in the warm season. These results provide evidence to a drought feedback process over the SGP. Surface processes begin to dry out soils through enhanced ET which leads to soil moisture anomalies which then influence the propagation of future dry line boundaries that can further strengthen or modify existing soil moisture anomalies through CI and precipitation. Thus, the evolution of surface conditions can enhance surface boundaries that cause the enhancement or modification of established surface anomalies (drought conditions). Future work should investigate this possible feedback loop further, using different drought cases during different years. Understanding of these processes could aid in the forecast of rapid drought evolution within the SGP, a result similarly noted in Otkin et al. [19].

Finally, it is important to note that numerous simulations were completed to determine the final configuration used in the numerical analyses. However, varying initial conditions and model configurations may lead to varying final results and possibly different overall conclusions. Future analyses are needed to determine the overall sensitivity of the dry line to varying drought evolution scenarios. While only a few modeling cases were shown, the cases were representative of the broader overall conditions during the period. Thus, the results increased the overall understanding of the dry line to rapidly varying surface conditions while further pointing to the need for further investigation across multiple drought years and large-scale environments driving the variability of surface conditions. Further, because in situ observations from the Mesonet were used to identify the dry line boundaries, cases that occurred outside of Oklahoma were not included. Even so, the numerous cases identified, and above average number of dry line events compared to the Hoch and Markowski [43] climatology, show that such a limitation was minor as the climatological numbers found in Hoch and Markowski [43] include dry lines that occurred over the entire SGP region.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported, in part, by the NOAA Climate Program Office’s Sectoral Applications Research Program (SARP) Grant NA13OAR4310122 and the Agriculture and Food Research Initiative Competitive Grant 2012-02355 from the USDA National Institute of Food and Agriculture (Grant nos. 2013-69002-23146 and 2016-68002-24967). Some of the computing for this project was performed at the OU Supercomputing Center for Education & Research (OSCER) at the University of Oklahoma (OU). Oklahoma’s taxpayers fund the Oklahoma Mesonet through the Oklahoma State Regents for Higher Education and the Oklahoma Department of Public Safety. The U.S. Drought Monitor is jointly produced by the National Drought Mitigation Center at the University of Nebraska-Lincoln, the United States Department of Agriculture, and the National Oceanic and Atmospheric Administration, Map Courtesy of NDMC-UNL.