skip to main content

Title: Dry-to-Wet Soil Gradients Enhance Convection and Rainfall over Subtropical South America

Soil moisture–precipitation (SM–PPT) feedbacks at the mesoscale represent a major challenge for numerical weather prediction, especially for subtropical regions that exhibit large variability in surface SM. How does surface heterogeneity, specifically mesoscale gradients in SM and land surface temperature (LST), affect convective initiation (CI) over South America? Using satellite data, we track nascent, daytime convective clouds and quantify the underlying antecedent (morning) surface heterogeneity. We find that convection initiates preferentially on the dry side of strong SM/LST boundaries with spatial scales of tens of kilometers. The strongest alongwind gradients in LST anomalies at 30-km length scale underlying the CI location occur during weak background low-level wind (<2.5 m s−1), high convective available potential energy (>1500 J kg−1), and low convective inhibition (<250 J kg−1) over sparse vegetation. At 100-km scale, strong gradients occur at the CI location during convectively unfavorable conditions and strong background flow. The location of PPT is strongly sensitive to the strength of the background flow. The wind profile during weak background flow inhibits propagation of convection away from the dry regions leading to negative SM–PPT feedback whereas strong background flow is related to longer life cycle and rainfall hundreds of kilometers away from the CI location. Thus, the sign of the SM–PPT feedback is dependent on the background flow. This work presents the first observational evidence that CI over subtropical South America is associated with dry soil patches on the order of tens of kilometers. Convection-permitting numerical weather prediction models need to be examined for accurately capturing the effect of SM heterogeneity in initiating convection over such semiarid regions.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Hydrometeorology
Page Range / eLocation ID:
p. 1563-1581
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Site description. This data package consists of data obtained from sampling surface soil (the 0-7.6 cm depth profile) in black mangrove (Avicennia germinans) dominated forest and black needlerush (Juncus roemerianus) saltmarsh along the Gulf of Mexico coastline in peninsular west-central Florida, USA. This location has a subtropical climate with mean daily temperatures ranging from 15.4 °C in January to 27.8 °C in August, and annual precipitation of 1336 mm. Precipitation falls as rain primarily between June and September. Tides are semi-diurnal, with 0.57 m median amplitudes during the year preceding sampling (U.S. NOAA National Ocean Service, Clearwater Beach, Florida, station 8726724). Sea-level rise is 4.0 ± 0.6 mm per year (1973-2020 trend, mean ± 95 % confidence interval, NOAA NOS Clearwater Beach station). The A. germinans mangrove zone is either adjacent to water or fringed on the seaward side by a narrow band of red mangrove (Rhizophora mangle). A near-monoculture of J. roemerianus is often adjacent to and immediately landward of the A. germinans zone. The transition from the mangrove to the J. roemerianus zone is variable in our study area. An abrupt edge between closed-canopy mangrove and J. roemerianus monoculture may extend for up to several hundred meters in some locations, while other stretches of ecotone present a gradual transition where smaller, widely spaced trees are interspersed into the herbaceous marsh. Juncus roemerianus then extends landward to a high marsh patchwork of succulent halophytes (including Salicornia bigellovi, Sesuvium sp., and Batis maritima), scattered dwarf mangrove, and salt pans, followed in turn by upland vegetation that includes Pinus sp. and Serenoa repens. Field design and sample collection. We established three study sites spaced at approximately 5 km intervals along the western coastline of the central Florida peninsula. The sites consisted of the Salt Springs (28.3298°, -82.7274°), Energy Marine Center (28.2903°, -82.7278°), and Green Key (28.2530°, -82.7496°) sites on the Gulf of Mexico coastline in Pasco County, Florida, USA. At each site, we established three plot pairs, each consisting of one saltmarsh plot and one mangrove plot. Plots were 50 m^2 in size. Plots pairs within a site were separated by 230-1070 m, and the mangrove and saltmarsh plots composing a pair were 70-170 m apart. All plot pairs consisted of directly adjacent patches of mangrove forest and J. roemerianus saltmarsh, with the mangrove forests exhibiting a closed canopy and a tree architecture (height 4-6 m, crown width 1.5-3 m). Mangrove plots were located at approximately the midpoint between the seaward edge (water-mangrove interface) and landward edge (mangrove-marsh interface) of the mangrove zone. Saltmarsh plots were located 20-25 m away from any mangrove trees and into the J. roemerianus zone (i.e., landward from the mangrove-marsh interface). Plot pairs were coarsely similar in geomorphic setting, as all were located on the Gulf of Mexico coastline, rather than within major sheltering formations like Tampa Bay, and all plot pairs fit the tide-dominated domain of the Woodroffe classification (Woodroffe, 2002, "Coasts: Form, Process and Evolution", Cambridge University Press), given their conspicuous semi-diurnal tides. There was nevertheless some geomorphic variation, as some plot pairs were directly open to the Gulf of Mexico while others sat behind keys and spits or along small tidal creeks. Our use of a plot-pair approach is intended to control for this geomorphic variation. Plot center elevations (cm above mean sea level, NAVD 88) were estimated by overlaying the plot locations determined with a global positioning system (Garmin GPS 60, Olathe, KS, USA) on a LiDAR-derived bare-earth digital elevation model (Dewberry, Inc., 2019). The digital elevation model had a vertical accuracy of ± 10 cm (95 % CI) and a horizontal accuracy of ± 116 cm (95 % CI). Soil samples were collected via coring at low tide in June 2011. From each plot, we collected a composite soil sample consisting of three discrete 5.1 cm diameter soil cores taken at equidistant points to 7.6 cm depth. Cores were taken by tapping a sleeve into the soil until its top was flush with the soil surface, sliding a hand under the core, and lifting it up. Cores were then capped and transferred on ice to our laboratory at the University of South Florida (Tampa, Florida, USA), where they were combined in plastic zipper bags, and homogenized by hand into plot-level composite samples on the day they were collected. A damp soil subsample was immediately taken from each composite sample to initiate 1 y incubations for determination of active C and N (see below). The remainder of each composite sample was then placed in a drying oven (60 °C) for 1 week with frequent mixing of the soil to prevent aggregation and liberate water. Organic wetland soils are sometimes dried at 70 °C, however high drying temperatures can volatilize non-water liquids and oxidize and decompose organic matter, so 50 °C is also a common drying temperature for organic soils (Gardner 1986, "Methods of Soil Analysis: Part 1", Soil Science Society of America); we accordingly chose 60 °C as a compromise between sufficient water removal and avoidance of non-water mass loss. Bulk density was determined as soil dry mass per core volume (adding back the dry mass equivalent of the damp subsample removed prior to drying). Dried subsamples were obtained for determination of soil organic matter (SOM), mineral texture composition, and extractable and total carbon (C) and nitrogen (N) within the following week. Sample analyses. A dried subsample was apportioned from each composite sample to determine SOM as mass loss on ignition at 550 °C for 4 h. After organic matter was removed from soil via ignition, mineral particle size composition was determined using a combination of wet sieving and density separation in 49 mM (3 %) sodium hexametaphosphate ((NaPO_3)_6) following procedures in Kettler et al. (2001, Soil Science Society of America Journal 65, 849-852). The percentage of dry soil mass composed of silt and clay particles (hereafter, fines) was calculated as the mass lost from dispersed mineral soil after sieving (0.053 mm mesh sieve). Fines could have been slightly underestimated if any clay particles were burned off during the preceding ignition of soil. An additional subsample was taken from each composite sample to determine extractable N and organic C concentrations via 0.5 M potassium sulfate (K_2SO_4) extractions. We combined soil and extractant (ratio of 1 g dry soil:5 mL extractant) in plastic bottles, reciprocally shook the slurry for 1 h at 120 rpm, and then gravity filtered it through Fisher G6 (1.6 μm pore size) glass fiber filters, followed by colorimetric detection of nitrite (NO_2^-) + nitrate (NO_3^-) and ammonium (NH_4^+) in the filtrate (Hood Nowotny et al., 2010,Soil Science Society of America Journal 74, 1018-1027) using a microplate spectrophotometer (Biotek Epoch, Winooski, VT, USA). Filtrate was also analyzed for dissolved organic C (referred to hereafter as extractable organic C) and total dissolved N via combustion and oxidation followed by detection of the evolved CO_2 and N oxide gases on a Formacs HT TOC/TN analyzer (Skalar, Breda, The Netherlands). Extractable organic N was then computed as total dissolved N in filtrate minus extractable mineral N (itself the sum of extractable NH_4-N and NO_2-N + NO_3-N). We determined soil total C and N from dried, milled subsamples subjected to elemental analysis (ECS 4010, Costech, Inc., Valencia, CA, USA) at the University of South Florida Stable Isotope Laboratory. Median concentration of inorganic C in unvegetated surface soil at our sites is 0.5 % of soil mass (Anderson, 2019, Univ. of South Florida M.S. thesis via methods in Wang et al., 2011, Environmental Monitoring and Assessment 174, 241-257). Inorganic C concentrations are likely even lower in our samples from under vegetation, where organic matter would dilute the contribution of inorganic C to soil mass. Nevertheless, the presence of a small inorganic C pool in our soils may be counted in the total C values we report. Extractable organic C is necessarily of organic C origin given the method (sparging with HCl) used in detection. Active C and N represent the fractions of organic C and N that are mineralizable by soil microorganisms under aerobic conditions in long-term soil incubations. To quantify active C and N, 60 g of field-moist soil were apportioned from each composite sample, placed in a filtration apparatus, and incubated in the dark at 25 °C and field capacity moisture for 365 d (as in Lewis et al., 2014, Ecosphere 5, art59). Moisture levels were maintained by frequently weighing incubated soil and wetting them up to target mass. Daily CO_2 flux was quantified on 29 occasions at 0.5-3 week intervals during the incubation period (with shorter intervals earlier in the incubation), and these per day flux rates were integrated over the 365 d period to compute an estimate of active C. Observations of per day flux were made by sealing samples overnight in airtight chambers fitted with septa and quantifying headspace CO_2 accumulation by injecting headspace samples (obtained through the septa via needle and syringe) into an infrared gas analyzer (PP Systems EGM 4, Amesbury, MA, USA). To estimate active N, each incubated sample was leached with a C and N free, 35 psu solution containing micronutrients (Nadelhoffer, 1990, Soil Science Society of America Journal 54, 411-415) on 19 occasions at increasing 1-6 week intervals during the 365 d incubation, and then extracted in 0.5 M K_2SO_4 at the end of the incubation in order to remove any residual mineral N. Active N was then quantified as the total mass of mineral N leached and extracted. Mineral N in leached and extracted solutions was detected as NH_4-N and NO_2-N + NO_3-N via colorimetry as above. This incubation technique precludes new C and N inputs and persistently leaches mineral N, forcing microorganisms to meet demand by mineralizing existing pools, and thereby directly assays the potential activity of soil organic C and N pools present at the time of soil sampling. Because this analysis commences with disrupting soil physical structure, it is biased toward higher estimates of active fractions. Calculations. Non-mobile C and N fractions were computed as total C and N concentrations minus the extractable and active fractions of each element. This data package reports surface-soil constituents (moisture, fines, SOM, and C and N pools and fractions) in both gravimetric units (mass constituent / mass soil) and areal units (mass constituent / soil surface area integrated through 7.6 cm soil depth, the depth of sampling). Areal concentrations were computed as X × D × 7.6, where X is the gravimetric concentration of a soil constituent, D is soil bulk density (g dry soil / cm^3), and 7.6 is the sampling depth in cm. 
    more » « less
  2. null (Ed.)
    Abstract During the Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations-Cloud, Aerosol, and Complex Terrain Interactions (RELAMPAGO-CACTI) field experiments in 2018–19, an unprecedented number of balloon-borne soundings were collected in Argentina. Radiosondes were launched from both fixed and mobile platforms, yielding 2712 soundings during the period 15 October 2018–30 April 2019. Approximately 20% of these soundings were collected by highly mobile platforms, strategically positioned for each intensive observing period, and launching approximately once per hour. The combination of fixed and mobile soundings capture both the overall conditions characterizing the RELAMPAGO-CACTI campaign, as well as the detailed evolution of environments supporting the initiation and upscale growth of deep convective storms, including some that produced hazardous hail and heavy rainfall. Episodes of frequent convection were characterized by sufficient quantities of moisture and instability for deep convection, along with deep-layer vertical wind shear supportive of organized or rotating storms. A total of 11 soundings showed most unstable convective available potential energy (MUCAPE) exceeding 6000 J kg −1 , comparable to the extreme instability observed in other parts of the world with intense deep convection. Parameters used to diagnose severe-storm potential showed that conditions were often favorable for supercells and severe hail, but not for tornadoes, primarily because of insufficient low-level wind shear. High-frequency soundings also revealed the structure and evolution of the boundary layer leading up to convection initiation, convectively generated cold pools, the South American low-level jet (SALLJ), and elevated nocturnal convection. This sounding dataset will enable improved understanding and prediction of convective storms and their surroundings in subtropical South America, as well as comparisons with other heavily studied regions such as the central United States that have not previously been possible. 
    more » « less
  3. This study analyzes an ensemble of numerical simulations of a heavy rainfall event east of Taiwan on 9 June 2020. Heavy rainfall was produced by quasi-stationary back-building mesoscale convective systems (MCS) associated with a mei-yu front. Global model forecast skill was poor in location and intensity of rainfall. The mesoscale ensemble showed liminal conditions between heavy rainfall or little to no rainfall. The two most accurate and two least accurate ensemble members are selected for analysis via validation against radar-estimated rainfall observations. All members feature moist soundings with low levels of free convection (LFC) and sufficient instability for deep convection. We find that stronger gradients in 100-m θe and θv in the most accurate members associated with a near-surface frontal boundary focus the lifting mechanism for deep, moist convection and enhanced rainfall. As the simulations progress, stronger southerly winds in the least accurate members advect drier mid-level air into the region of interest and shift the near-surface boundary further north and west. Analysis of the verification ensemble mean analysis reveals a strong near-surface frontal boundary similarly positioned as in the most accurate members and dry air aloft more similar to that in the least accurate members, suggesting that the positioning of the frontal boundary is more critical to accurately reproducing rainfall patterns and intensity in this case. The analyses suggest that subtle details in the simulation of frontal boundaries and mesoscale flow structures can lead to bifurcations in producing extreme or almost no rainfall. Implications for improved probabilistic forecasts of heavy rainfall events will be discussed. 
    more » « less
  4. Abstract

    Frequent deep convective thunderstorms and mesoscale convective systems make the Córdoba region, near the Sierras de Córdoba mountain range, one of the most active areas on Earth for hail activity. Analysis of hail observations from trained observers and social media reports cross-referenced with operational radar observations identified the convective characteristics of hail-producing convective systems in central Argentina over a 6-month period divided into early (October–December 2018) and late seasons (January–March 2019). Reflectivity and dual-polarization characteristics from the Córdoba operational radar [Radar Meteorológico Argentina (RMA1)] were used to identify the convective modes of convective cells at time of positive hail indicators. Analysis of ERA5 upper-air and surface data examined convective environments of hail events and identified representative dynamic and thermodynamic environments. A majority of early season hail-producing cells were classified as discrete convection, while discrete and multicell occurrence evened out in the late season. Most hail-producing cells initiated directly adjacent to the Sierras in the late season, while cell initiation and hail production is further spread out in the early season. Dividing convective events into dynamic/thermodynamic regimes based on values of 1000 J kg−1of CAPE and vertical wind shear of 20 m s−1results in most early season events reflecting shear-dominant characteristics (low CAPE, high shear) and most late-season events exhibiting CAPE-dominant characteristics (high CAPE, low shear). Strength and placement of low-level temperature and moisture anomalies/advection and upper-level jets largely defined the differences in the dominant regimes.

    Significance Statement

    This study used regional radar data alongside hail reports from trained observers and social media to better understand the types and timing of storms identified as producing hail, given the lower resolution of satellite studies. Dividing the hail season (October–December; January–March) showed that within hail season, early season storms tended to be singular storms that formed across the region in environments with strong vertical winds and weak instability. Late-season storms were a mix of singular storms and multicellular storm systems focused on the mountains in weak vertical winds and strong instability. These results show differences from satellite studies and identify key representative hail-producing radar features and environmental regimes for this region, which could guide hail risk analysis within the severe-weather season.

    more » « less
  5. null (Ed.)
    Abstract Subtropical South America (SSA) east of the Andes Mountains is a global hotspot for mesoscale convective systems (MCSs). Wide convective cores (WCCs) are typically embedded within mature MCSs, contribute over 40% of SSA’s warm-season rainfall, and are often associated with severe weather. Prior analysis of Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR) data identified WCCs in SSA and associated synoptic conditions during austral summer. As WCCs also occur during the austral spring, this study uses the 16-year TRMM PR dataset and ERA5 reanalysis to compare anomalies in environmental conditions between austral spring (SON) and summer (DJF) for the largest and smallest WCCs in SSA. During both seasons, large WCCs are associated with an anomalous mid-level trough that slowly crosses the Andes Mountains and a northerly South American low-level jet (SALLJ) over SSA, though the SON trough and SALLJ anomalies are stronger and located farther northeastward than in DJF. A synoptic pattern evolution resembling large WCC environments is illustrated through a multi-day case during the RELAMPAGO field campaign (10-13 November 2018). Unique high-temporal resolution soundings showed strong mid-level vertical wind shear associated with this event, induced by the juxtaposition of the northerly SALLJ and southerly near-surface flow. It is hypothesized that the Andes help create a quasi-stationary trough/ridge pattern such that favorable synoptic conditions for deep convection persist for multiple days. For the smallest WCCs, anomalously weaker synoptic-scale forcing was present compared to the largest events, especially for DJF, pointing to future work exploring MCS formation under weaker synoptic conditions. 
    more » « less