Soil Moisture Conditions over Southeast Texas Prior to Hurricane Harvey

Soil Moisture Conditions over Southeast Texas Prior to Hurricane Harvey

As much-anticipated Hurricane Harvey approaches the southern and eastern coast of Texas today, it is worth examining the pre-existing soil moisture over the region to understand the capacity of the land surface to absorb the upcoming rainfall.  Granted, the amount of rainfall simulated by numerical guidance is off-the-charts high (e.g., today’s 0600 UTC initialized NAM model [Fig. 1] shows 84-hour maximum accumulated rainfall of over 60″ between Corpus Christie and Houston!!).  Thus, extreme flooding is anticipated, regardless of the amount that can be absorbed by the soils.


Figure 1.  The NCEP/NAM model 84-hour forecast of total accumulated precipitation (inches) over Southeastern Texas, from the simulation initialized at 0600 UTC 25 August 2017 [image courtesy of College of DuPage forecast page].

SPoRT manages a real-time simulation of the NASA Land Information System (hereafter, “SPoRT-LIS“), running over the Continental U.S. at ~3-km grid resolution.  The SPoRT-LIS product is a Noah land surface model climatological and real-time simulation over 4 model soil layers (0-10, 10-40, 40-100, and 100-200 cm).  The climatological simulation spans 1981-2013 and forms the basis for daily-updated total-column soil moisture percentiles (forthcoming in Fig. 3), in order to place current soil moisture values into historical context.  For real-time output, the Noah simulation is regularly updated four times per day as an extension of the long-term climatology simulation.  It includes NOAA/NESDIS daily global VIIRS Green Vegetation Fraction data, and the real-time SPoRT-LIS component also incorporates quantitative precipitation estimates (QPE) from the Multi-Radar Multi-Sensor (MRMS) gauge-corrected radar product.  The climatological SPoRT-LIS is based exclusively on atmospheric analysis input from the NOAA/NASA North American Land Data Assimilation System – version 2.

Relative Soil Moisture output from the SPoRT-LIS over the 0-100 cm layer is shown in Fig. 2 over Southeastern Texas and Louisiana at 1200 UTC this morning.  A marked gradient between very dry soils to the west and moist soils to the east occurs in the vicinity of the greater Houston metropolitan area.  The soils in the region bounded by Corpus Christi, San Antonio, Austin, and Houston (areas forecast to have the greatest rainfall from Hurricane Harvey) are extremely dry prior to Harvey’s landfall.  This dryness will help to some extent in absorbing the initial rainfall from Hurricane Harvey.  But with such excessive rainfall being forecast over a prolonged time period (3-5+ days), it won’t be long before the upper portions of the soil column saturates and widespread areal flooding occurs.  In addition, the high forecast rainfall rates could easily result in flash flooding (despite prevailing soil dryness), especially further inland where terrain plays a more important role in runoff and flash flooding.

The total column relative soil moisture percentile from 24 August shows that historically-speaking, the soil moisture is slightly drier than normal, particularly along the coastal plain between Corpus Christi and Houston (Fig. 3).  In this corridor, the soil moisture is generally between the 10th and 30th percentile compared to the 1981-2013 climatological distribution for 24 August.


Figure 2.  SPoRT-LIS relative soil moisture (RSM) distribution in the 0-1 meter layer across Southeastern Texas and Louisiana, valid 1200 UTC 25 August 2017.  RSM values of 0% represent wilting (vegetation cannot extract moisture from soil) and 100% represents saturation (subsequent rainfall becomes runoff).


Figure 3.  Total column (0-2 m) relative soil moisture percentile valid 24 Aug 2017, as compared to all 24 August soil moisture values from a 33-year climatological simulation of the SPoRT-LIS.

Finally, an hourly animation of the 1-day changes in 0-10 cm (top model layer) relative soil moisture show that the near-surface soils are quickly moistening between Corpus Christi and Houston, as the initial rainbands of Hurricane Harvey began impacting the coastal plain this morning.  As the soils continue to moisten rapidly from the top-down, subsequent rainfall will quickly lead to runoff and flooding.


Figure 4.  Hourly animation of 1-day change in top-layer (0-10 cm) relative soil moisture, for the time period spanning 0000-1400 UTC 25 August 2017.  Each hourly image is a simple difference in 0-10 cm relative soil moisture between the current and previous day at the same valid hour.  Line contours depict one-hour QPE from the MRMS product, as input to the real-time SPoRT-LIS.

Multiple Atmospheric Rivers Impact California in Early 2017

The state of California has been suffering from a multi-year drought that has severely depleted water resources and reservoir levels. Recent winters have failed to produce precipitation and mountain snows to replenish the losses during the dry summers. However, the situation has rapidly changed this winter, particularly in the past week when multiple atmospheric rivers have impacted the state.

An atmospheric river is a concentrated channel of deep moisture that is transported from the tropical Pacific Oceanic regions to the West Coast of the United States.  These events are often associated with prodigious amounts of rainfall and mountain snows that lead to flooding, mudslides, and avalanches.  We have seen such events this past week impact California, especially the central and northern parts of the state.  CIRA’s total precipitable water product in Figures 1a and 1b depict two separate atmospheric rivers impinging on central California from 8 and 10 January 2017, respectively. The first wave transported a plume of tropical moisture from the south-southwest, which led to massive rainfall and high snow levels.  The second atmospheric river on the 10th was less directly connected to the tropics (coming in from the west-southwest), but nonetheless exhibited a well-focused transport of high moisture content.  Widespread flooding and mountain avalanches have resulted from these moisture plumes as the impacted California, as well as dramatic replenishment of reservoirs.


Figure 1.  CIRA total precipitable water product (inches) valid at (a) 2100 UTC 8 Jan 2017, and (b) 2100 UTC 10 Jan 2017.


SPoRT’s real-time instantiation of the Land Information System (aka “SPoRT-LIS”) has nicely depicted the substantial replenishment of the moisture content in the soils over California.  The SPoRT-LIS is an observations-driven, ~3-km resolution run of the Noah land surface model that consists of a 33-year climatology spanning 1981-2013, and real-time output at hourly intervals sent to select NOAA/NWS partnering forecast offices.  The one-year change in the SPoRT-LIS total column soil moisture at 1200 UTC 11 January (Fig. 2) shows large increases over most of California, particularly in the higher terrain (given by blue and purple shading).  [At the same time, annual degradation in soil moisture can be seen across the central and eastern U.S.]  Interestingly, a substantial portion of California’s annual soil moisture increases has occurred in just the past week (Fig. 3; SPoRT-LIS total column soil moisture change over the past week).  One can certainly see the important role that atmospheric rivers play in being “drought busters”!


Figure 2.  One-year change in the SPoRT-LIS total column relative soil moisture, valid 1200 UTC 11 January 2017.



Figure 3.  One-week change in the SPoRT-LIS total column relative soil moisture, valid 1200 UTC 11 January 2017.


A map of the SPoRT-LIS daily soil moisture percentiles from 11 January highlight the very wet anomaly over California relative to the 33-year soil moisture climatology (Fig. 4; similar to the pattern of annual soil moisture change from Fig. 2).  Blue shading denotes greater than or equal to the 98th percentile, thus indicating unusually wet soils on the tail end of the historical distribution.


Figure 4.  SPoRT-LIS total column relative soil moisture percentile from 11 January 2017.


Finally, SPoRT is acquiring and assimilating in real time the Soil Moisture Active Passive (SMAP) Level 2 (L2) retrievals produced by NASA/JPL into an experimental version of the SPoRT-LIS.  SPoRT is a SMAP Early Adopter and has a funded project to conduct soil moisture data assimilation experiments with LIS and evaluate impacts on land surface and numerical weather prediction models.  Figure 5 shows SMAP L2 retrievals of the evening overpasses from ~0000 UTC 11 January.  Panel (a) is the 36-km resolution radiometer product, while panel (b) shows the enhanced-resolution product, obtained from the SMAP radiometer by using Backus-Gilbert optimal interpolation techniques to provide data on a finer (9 km) grid.  The enhanced-resolution product provides much more detail of the wet soils in California, while retaining the same overall regional patterns as the original 36-km retrieval.  Given the loss of the active radar component of the SMAP mission, SPoRT plans to assimilate both the 36-km and 9-km products separately, and compare results on model accuracy.


Figure 5. SMAP Level 2 soil moisture retrievals for the evening overpasses from ~0000 UTC 11 January 2017; (a) 36-km resolution product; (b) enhanced 9-km resolution product.

U.S. Deep South flooding as depicted by NASA’s SMAP satellite and SPoRT-LIS

U.S. Deep South flooding as depicted by NASA’s SMAP satellite and SPoRT-LIS

A significant flooding event occurred over the U.S. Deep South from 8-10 March 2016 due to a slow-moving low pressure/front from southern Texas to the Mississippi River, combined with a deep tropical moisture plume.  Tremendous rainfall totals of 4-8″+ were depicted by the Multi-Radar Multi-Sensor (MRMS) 1-km gauge-corrected radar rainfall estimate product for consecutive 24-hour periods ending 1200 UTC 9 March and 10 March 2016 (Fig. 1).  The MRMS gridded product provides short-term input precipitation estimates to the real-time Land Information System (LIS) run at the NASA Short-term Prediction Research and Transition (SPoRT) Center.


Fig. 1.  Rainfall estimates from the Multi-Radar Multi-Sensor (MRMS) gauge-adjusted radar product for the 24-h period ending (a) 1200 UTC 9 Mar 2016, and (b) 1200 UTC 10 Mar 2016.

The soil moisture response to the heavy rainfall was captured nicely by NASA’s Soil Moisture Active-Passive (SMAP) satellite Level 2 retrieval product of 0-5 cm volumetric soil moisture for the early morning overpass across the Central U.S from 9 March (Fig. 2a).  Very high retrieved volumetric soil moisture of 0.45 or higher is seen from eastern Texas across much of northern Louisiana and southern/central Arkansas, aligned well with the areas that received the heaviest rainfall.  The corresponding SPoRT-LIS modeled soil moisture analysis for the 0-10 cm layer of the Noah land surface model (Fig. 2b) shows a reasonable agreement with the SMAP retrieved soil moisture but with slightly less dynamic range than the SMAP data — a reasonable result given that the SMAP retrieval is valid over a shallower, near-surface layer than the LIS top model layer.


Fig. 2.  (a) Soil Moisture Active-Passive (SMAP) 0-5 cm volumetric soil moisture retrieval valid 1223 UTC 9 Mar 2016, and (b) corresponding SPoRT-Land Information System (LIS) 0-10 cm volumetric soil moisture analysis valid 1200 UTC 9 Mar 2016.

The SPoRT-LIS total column relative soil moisture (RSM) has been found to qualitatively correlate with areas of river/areal flooding when exceeding ~60% across northern Alabama.  The depiction of total column RSM from the early morning of 10 March shows a large area exceeding 65% (blue  shading in Fig. 3a) across eastern Texas, northern Louisiana, and southern/central Arkansas.  The weekly change in total column RSM (Fig. 3b) highlights the regions that experienced the largest increases in soil moisture in response to the heavy rainfall.  As the soil approaches saturation/field capacity, most new rainfall goes directly to runoff, thereby exacerbating the flooding situation.


Fig. 3.  SPoRT-LIS analysis valid 1200 UTC 10 Mar 2016 of (a) total column relative soil moisture (RSM), and (b) one-week change in total column RSM.

Another product available from the real-time SPoRT-LIS is the total column RSM percentile, which shows how anomalous the current soil moisture conditions are relative to a historical, 33-year climatological database of modeled soil moisture. The percentile product from the morning of 10 March (Fig. 4) shows that areas of eastern Texas and northwestern Louisiana (and a few other spots) exceed the 98th percentile for the current day — indicating that these present-day soil moisture values are about the most moist it has been in the last 33+ years for 10 March.

Notice that there are a few anomalous “bulls-eyes” of dry percentiles in southern/eastern Arkansas.  These areas are caused by corrupt input precipitation data driving the SPoRT-LIS land surface model simulation.  The issue is currently being corrected by the operational organizations managing the rain gauge input.  However, it should be noted that SPoRT is working to implement near real-time data assimilation of the SMAP Level 2 soil moisture retrievals into its LIS simulation.  Routine assimilation of SMAP satellite soil moisture will help correct anomalies caused by poor input precipitation, thereby resulting in more robust soil moisture analyses for situational awareness and disaster-response applications.


Fig. 4. SPoRT-LIS total column RSM percentile valid 1200 UTC 10 Mar 2016.

Finally, the areas of active flooding at U.S. Geological stream gauges and the NOAA/NWS flood watch/warning map are given in Fig. 5 for the afternoon of 10 March.  The axis of flash flood warnings (dark red color) aligns quite well with the total column RSM percentiles exceeding the 98th percentile in Fig. 4, whereas the broader footprint of all flood warnings/watches corresponds closely with the total column RSM above the 65% value (blue shading in Fig. 3a).


Fig. 5.  Plot of U.S. Geological Survey stream gauges with active flooding (top), and NOAA/NWS active [flash] flood watches and warnings (bottom-right) for the afternoon of 10 Mar 2016.

SFR Product Verifies Snow Coverage over Four Corners

The NESDIS Snowfall Rate (SFR) product assessment is in full swing at NWS Albuquerque and forecasters are already capturing some good cases over data sparse regions. The first week of January 2016 was very active across New Mexico as back to back winter storm systems crossed the area. The second system in the series crossed over the Four Corners region on 4 January 2016, producing light to moderate snowfall rates for several hours. The forecaster on shift noted the observation at Farmington, NM (KFMN) indicated light snow with a visibility of 5 statute miles. A quick glance at the SFR procedure used in Figure 1a shows the extent of any precipitation echoes well to the east of KFMN at 0000 UTC 5 January 2016. The nearest radar (KABX, not shown) is located roughly 150 miles southeast of KFMN near Albuquerque, NM. The arrival of a SFR product at 0010 UTC 5 January 2016 showed the extent of the precipitation was much greater with the merged POES image overlaid on the radar data (Figure 1b). Sampled liquid equivalent values in the light green areas to the east of KFMN were near 0.03″/hour.

Figure 1a. Liquid equivalent values of the merged SFR product valid 0000 UTC 5 January 2016. KFMN is denoted by the white circle. Note the extent of the radar coverage is well east of KFMN.

Figure 1b. Liquid equivalent values of the merged SFR product valid 0010 UTC 5 January 2016. KFMN is denoted by the white circle. Note the extent of the snowfall coverage is much greater with the addition of the POES image.

The Terminal Aerodrome Forecast (TAF) issued for KFMN shortly before the receipt of this image indicated temporary fluctuations in the visibility to 1 statute mile with light snow and an overcast ceiling near 1,200 ft between 0000 UTC and 0400 UTC (Instrument Flight Rules, IFR). It is not clear whether any operational changes occurred based on the receipt of the merged SFR product or whether the product increased confidence on the IFR forecast. However, it is entirely possible given the improvement in product latency compared to the 2015 assessment that the imagery could be used in this way.

The webcam available at San Juan College just a short distance from the KFMN observation showed significant decreases in the visibility between 330pm and shortly after sunset (Figure 2a and 2b). The two images below show the decrease in surface visibility as well as notable accumulations on grassy surfaces in front of the college. An observer 3 miles southeast of Farmington did report a total accumulation of 1″ from this event. The merged SFR product did in fact show higher rates immediately to the east of KFMN. The last image in the series shows the impact on travel conditions noted by the NM Department of Transportation web page (Figure 3). The areal coverage of the difficult travel impacts (yellow highlights) was greater than that depicted by what can be seen based on poor radar coverage.

Figure 2a. Webcam at San Juan College around 330pm. Note the light snowfall beginning to develop over the distant mesas behind the college.

Figure 2a. Webcam at San Juan College around 330pm. Note the light snowfall beginning to develop over the distant mesas behind the college.

Figure 2b. Webcam at San Juan College shortly after sunset. Note the dramatic decrease in visibility and light snow accumulations on grassy surfaces in front of the college.

Figure 2b. Webcam at San Juan College shortly after sunset. Note the dramatic decrease in visibility and light snow accumulations on grassy surfaces in front of the college.

Figure 3. Screen capture of NM DOT web page showing areal coverage of difficult travel conditions (yellow highlights) and some text summaries detailing the impacts.

Heavy morning snow on 03Mar2014 with some freezing rain far southeast



During the weekend of Mar 02-03 2014, several weather features moved northeast across the area. The precipitation started out as rain across West Virginia with some freezing rain, sleet and snow across portions of southeast Ohio. Colder air began to filter into the region and as it did, the precipitation changed from rain to freezing rain to sleet and finally to snow. By 603 AM, the precipitation had turned to snow across all of West Virginia, but for portions of the extreme southeast counties.

I have attached two images from around 6 AM on Mar 3rd. The first image showed the radar data from KRLX at 603 am while the second image was the 607 am SFR product and 6 AM surface observations. When comparing these images, the “best” SFR signal for heavy snow was located along a line where the precipitation transitioned from freezing rain to snow. The heaviest signal in the SFR images was actually located over Mingo County where a total of 8 inches of snow was reported from the storm.

Freezing rain was falling across portions of extreme southeast West Virginia. Bluefield WV (KBLF) is located southeast of the “SFR” heaviest snow signal in an area where the SFR product is not showing anything. The SFR product did a great job across that area as the 6 AM KBLF observation indicated freezing rain was falling at that time.

SFR performance at temperatures below 22 DegF


On Jan 26 2014, an upper level shortwave caused an area of light snow across Ohio, western Pennsylvania and the northern counties of West Virginia. Surface temperatures were quite cold with readings generally in the teens. Even at these cold temperatures, the SFR product did indicate snowfall across the far northern counties of our forecast area.

The maximum snowfall rates indicated on the 1605 UTC product was about 0.3 to 0.4 inches per hour. Based on reports, these numbers appear to be representative of what actually was occurring.

While this is just one case, the SFR product appears to work reasonably well at temperatures below 22 DegF.

Snow Advisory Event on Jan 25 2014


On Jan 25 2014, a mid-level shortwave moved across the region generating light to moderate snow. I have included screen captures of the 1118 UTC regional radar mosaic and surface observations…along with a 1120 UTC Snowfall Rate Product and surface observations.

It looks like the SFR product did not detect all of the snow that was falling around 11 UTC. But the misses can generally be described as either (1) the surface temperatures being too cold or (2) the probabilistic model, that is part of the calulations, indicating probabilities that were too low to determine if there was snow.

Once you know all of the details on how the product is calculated, I think this product did a good job at detailing where the snowfall was occurring.

The highest snowfall rates indicated by this image was around 0.3 to 0.5 inches which seems to be representative of what was occurring.