Abstract
Forest mortality caused by convective storms (windthrow) is a major disturbance in the Amazon. However, the linkage between windthrows at the surface and convective storms in the atmosphere remains unclear. In addition, the current Earth system models (ESMs) lack mechanistic links between convective wind events and tree mortality. Here we find an empirical relationship that maps convective available potential energy, which is well simulated by ESMs, to the spatial pattern of large windthrow events. This relationship builds connections between strong convective storms and forest dynamics in the Amazon. Based on the relationship, our model projects a 51 ± 20% increase in the area favorable to extreme storms, and a 43 ± 17% increase in windthrow density within the Amazon by the end of this century under the high-emission scenario (SSP 585). These results indicate significant changes in tropical forest composition and carbon cycle dynamics under climate change.
Similar content being viewed by others
Introduction
Convective storms account for ~50% of the tree mortality across Amazonia at scales ranging from individual trees1 to large forest gaps exceeding thousands of hectares2,3,4,5, resulting in uprooted and broken trees (windthrow6). Despite its importance, a detailed understanding of the mechanisms of convective storms producing windthrow events remains unclear. Given the sensitivity of convective storms to climate change, a better understanding of the interactions between convective storms and windthrow events is needed to explain changes in forest composition and ecosystem processes under the elevated disturbance regimes expected with continued climate system warming.
Convective storms are deep moist convective cells, formally known as Mesoscale Convective Systems (MCSs). In the Amazon basin, 64% of MCSs last 3–5.5 hours, while 0.7% of MCSs have a lifespan equal to or over 24 hours7. Some of the most extreme convective systems, such as squall lines, also occur over the Amazon region8. The occurrence of MCS is most common between afternoon and evening in the local time9,10. The development of these convective cells includes an updraft of moist air, which forms the cumulonimbus cloud, and a downdraft, which is associated with heavy precipitation, strong surface winds, hail, and cloud-to-ground lightning9,11. Warm and moist near-surface air in the Amazon region provides a favorable environment for the development of MCSs7. MCSs are driven by moist convective instability, which can be well represented by convective available potential energy (CAPE). CAPE determines the speed limit of upward convection. CAPE is also fundamental to our understanding of future convective storms12. CAPE increases robustly with warming across the tropics in ESMs, suggesting an increasing likelihood of favorable environments for most severe storms11,13,14,15. In addition to CAPE, vertical windshear physically displaces downdraft air from the storm, which prolongs and intensifies storms, and also further influences the development of further convections11,16.
In the Amazon, MCSs have significant impacts on forest ecosystem processes, community structure and composition, and the regional carbon balance5,17,18,19. MCSs produce 50–90% of the annual precipitation in the tropics, shaping key aspects of tropical forest function and structure11. Severe winds and heavy rainfall associated with MCSs also cause widespread windthrows in the Amazon rainforest10,20. At the continental scale, Nelson et al.5 mapped 330 large windthrows in the Brazilian Amazon using satellite images and explored the size and spatial pattern of windthrow events. Field studies of windthrows in the Amazon show there is a positive correlation between the satellite-derived disturbance intensity metrics and field-measured tree mortality within windthrows2,17. The study of a windthrow event in the Manaus region documented a high regional mortality from a single squall line17. In addition to wind and rainfall, lightning associated with convective storms can also be a major agent of canopy disturbance and large tree mortality21,22.
Previous estimates of windthrow events in the Amazon relied on measuring extreme rainfall events and cloud top temperature. The highest frequencies of windthrows were found to coincide with the frequency of heaviest rainfall days2,5,23 and strong convective activities measured from cloud top temperature10,18. Storms that caused windthrows were the most intense storms, capable of ascending furthest into the upper troposphere or lower stratosphere; therefore, these storms come with the lowest cloud top temperature and heaviest rainfall8,10,24,25. However, extreme rainfall events are poorly projected over the tropical regions in ESMs26, and realistic storm cloud top projections are not represented in the models. Moreover, it is still unclear how the cloud physics of MCSs will change under a warming climate. Therefore, a reliable climate variable with robust model projection in the future that correlates with windthrows is lacking, and the current suite of ESMs27 is unable to project future forest windthrows dynamics.
This study seeks to investigate the spatial pattern of windthrow events that are larger than 25,000 m2 (hereafter ‘large windthrows’) and their correlation with CAPE over 30 years while demonstrating that CAPE is an appropriate proxy for estimating the spatial pattern of large windthrow events in the Amazon. CAPE has been well studied in ESMs13,14,16,28. Compared to the relatively uniform geographic variability of weak windshear11 over the Amazon region (Supplementary Fig.1), CAPE shows a greater spatial variation. Furthermore, the frequency of severe thunderstorm environments increases with warming across the tropics, mainly due to changes in the distribution of CAPE13. Therefore, CAPE is selected for this study as the only climate variable for the spatial pattern analysis of windthrows. A model is also developed using CAPE to project future windthrow events, and it projects a large increase of windthrow events in the Amazon region by the end of this century.
In this study, we manually map 1012 large windthrow events encompassing 30 years from 1990–2019 (Fig. 1a) using images from Landsat 8 OLI/TIRS collection (USGS) over the entire Amazon region (see Methods for the details). Windthrow density is generated using all of these windthrow events. This study uses only large windthrow events that can be identified using Landsat satellite images. These large windthrows are produced by extreme convective storms10, such as squall lines17. Although this approach omits the vast majority of small windthrow events (≤25,000 m2)19, the general spatial pattern of large windthrows occurrence is used as a tractable proxy for all windthrows. Hourly CAPE is obtained from ERA5 reanalysis data29. The means of CAPE in the local time 13:00–19:00 (UTC 17:00–23:00) over 1990–2019 are calculated (hereafter mean afternoon CAPE).
Results
The spatial pattern of windthrows and CAPE
As seen in Fig. 1, the spatial pattern of mean afternoon CAPE matches with the density of large windthrows. High densities of large windthrow events are in the northwest and central Amazon regions (Fig. 1b) that experience the highest CAPE (Fig. 1d). Few large windthrow events are found on the northeast side of the Amazon region, where CAPE is low, indicating less severe storms30. The region with the highest windthrow density over 12 events per 10,000 km2 is in the state of Amazonas in Brazil, in between Rio Negro and the Amazon River. Another dense region with windthrow density over 7 events per 10,000 km2 is near the east border of Columbia. Both regions occur in locations where the mean afternoon CAPE exceeds 1000 J kg−1 (Fig. 1c). In addition to mean afternoon CAPE, mean daily CAPE over 30 years have a similar spatial pattern (Supplementary Fig. 2), but it exhibits high values over rivers due to the nighttime CAPE intensification (Supplementary Fig. 3).
A relationship maps CAPE to windthrow density
We found an empirical relationship using mean afternoon CAPE to reproduce the spatial variability of windthrows and project the frequency of future windthrows in the Amazon (Fig. 2a). The relationship uses the locations of windthrow events shown in Fig. 1a and the distribution of mean afternoon CAPE shown in Fig. 1c. The values of CAPE at the locations of the 1012 windthrow events are sorted and split into six equal-sized groups (of 168 or 169 each). This determines the boundaries of the CAPE bins used in Fig. 2a, which plots the current windthrow density, i.e., the number of windthrow events (168 or 169) divided by the area of the Amazon with that CAPE (see Methods for the details). Fig 2a shows that a CAPE threshold at 1023 J kg−1 can be used to define if the environment is favorable for high density of windthrows. A similar threshold is defined by a mechanistically-based approach using the data shown in Fig. 1b, d (Supplementary Fig. 4) and is supported by convection theory, numerical models, and observations in previous studies31,32. 62% of the Amazon area has a relatively low windthrow density where CAPE is lower than 1023 J kg−1. 38% of the Amazon region has a mean CAPE over 1023 J kg−1 (Fig. 2c), which provides a storm-favorable condition with high windthrow density. In general, windthrows occur nearly 8 times more frequently in the favorable environments (CAPE > 1023 J kg−1, hereafter storm-favorable environments) than in the non-favorable environments (CAPE ≤ 1023 J kg−1). The mean uncertainty (SD) of these values, estimated via bootstrapping, is 0.2. The relationship shows that CAPE is a strong predictor of the spatial pattern of large windthrow events in the Amazon region and provides a step function that can model windthrow density using CAPE values.
Future increase in storm-favorable environments and windthrows
We developed a model based on the CAPE and windthrow density relationship to assess how climate change will affect CAPE and then affect the density of future large windthrows in the Amazon region. Outputs from 10 ESMs in CMIP6 were analyzed and incorporated into our model. Figure 2b shows the distribution of area by the current and future CAPE values over the Amazon region. The comparison between the spatial pattern of current and future CAPE can be found in Supplementary Fig. 5a, b. The current CAPE distribution was generated from mean afternoon ERA5 CAPE. For future CAPE in 2070–2099, we calculated the mean of CAPE values from 10 CMIP6 ESMs and adjusted the values to make them comparable with ERA 5 CAPE values (see Methods for the details). Table 1 lists the fractional changes in storm-favorable CAPE for the ESMs in our ensemble (see Supplementary Table 1 for details of the ESMs) under the highest emission scenario SSP585. All 10 ESMs agree on an increasing trend of CAPE over the Amazon region, with a mean increase of 26 ± 9% (mean ± SD) from the past 30 years (1990–2015) to the last 30 years of the 21st century (Fig. 2d). There is a high level of agreement between models on the spatial pattern and magnitude of the CAPE increase (Fig. 2d). The magnitude of future CAPE indicates strong atmospheric instability that could produce more frequent storms. While the fraction of area that is storm-favorable in the current mean CAPE map is 38%, ESMs project that 58 ± 8% of the Amazon region becomes storm-favorable by the end of the 21st century. The area with CAPE over 1023 J kg−1 increases 51 ± 20% by the end of the 21st century (Table 1), indicating a much larger area with storm-favorable conditions within the Amazon. An area of ~1,390,000 km2 in the northwestern, southern, and central Amazon, including forested areas in central Columbia, the northern part of Bolivia, the south-eastern part of Peru, and central Brazil, where windthrows are infrequent now, are projected to experience an eightfold or greater increase in windthrow events by 2100 (Fig. 2c). The northwestern and western Amazon, where large windthrow density is relatively high now (2.4 windthrows/10,000 km2), are estimated to experience a 33–50% increase in the number of windthrow events in the next 50–80 years (Supplementary Fig. 5c).
Using the model, all ESMs in our ensemble predict a significant increase in future windthrow events. Overall, windthrow density in the Amazon region is estimated to increase 43 ± 17% over the 21st century (Table 1).
Discussion
Tree mortality plays important roles in determining forest carbon balance across Amazonia, and variable disturbance regimes increase uncertainty in the tropical forest carbon sink capacity over time. Extreme convective storms are important drivers of tree mortality in the Amazon region and affect the biomass patterns and function composition of Amazonian forests33. In this study, we provide a framework for representing the coupling between forest mortality on the land surface and windstorms in the atmosphere. Previous studies showed that tree mortality in the central Amazon was higher in wet months as a result of extreme storms, even in drought years34,35. Elevated forest mortality related to extreme storms in the central Amazon was found in the La Niña wet year of 1990, La Niña year of 1999 followed by prolonged warmer temperatures and drought in the previous years, and in the drought year of 200534. The relationship between extreme storms and tree mortality also implies that increasing frequency of convective storms36 may contribute to the observed increase in tree mortality and a weakening of the carbon sink across Amazonia37,38,39.
This analysis highlights the potential for predicting the rate of future windstorm-driven tree mortality, a driver of tree mortality that currently is not included in ESMs, and emphasizes the need to improve land-atmosphere coupling in ESMs. The results in this paper generated using ESM outputs and the model indicate that over half of the Amazon region is projected to experience average conditions favorable for intense convective storms with a large increase in windthrow events by the end of this century under the highest emission scenario SSP585. The projected increase in the number of large windthrows potentially facilitate the establishment of forest species with fast turnover rate2,33,34. These species are vulnerable to disturbances, which leads to greater forest mortality and in turn create more forest gaps. This positive feedback loop speeds up forest dynamics with faster carbon-water and nutrient cycling and less carbon storage2,33,34. The projected significant increase in tree mortality associated with windthrow density in the Amazon by the end of the 21st century has the potential to re-shape tropical forest structure and composition, with sizable impacts on regional carbon balance.
This first step in projecting the impact of future extreme wind events on Amazon forests has several caveats, which we discuss here to facilitate the interpretation of our results. Potential biases can exist in disturbance identification and sampling effort. 2.8% of the Amazon region was not included in this study due to the lack of clear remote sensing images (clouds <20%). Windthrows without a characteristic shape were eliminated from this study. We only considered the spatial pattern of large windthrows in this study while omitting the vast majority of windthrows that are smaller than 25,000 m2. Since the focus of this study is on the spatial pattern of large windthrow events, the underestimation of windthrow events will not have a significant impact on the results. In addition, the large windthrow dataset in this study provides a comprehensive view over the entire Amazon which would be difficult to achieve using small-scale windthrow datasets.
There are also biases arising from regional variation in vegetation composition or vegetation dynamics. Compared with the central Amazon, forests in the northwestern Amazon are more vulnerable to windthrow mortality and recover faster due to a number of factors, including tree stem density, soil type, root architecture, tree diversity, microclimate condition, topographic exposure2,40. A previous study showed that windthrows were spatially more common in the northwestern Amazon than in the central Amazon2. With high windthrow tree mortality, forest in the northwestern amazon took ~20 years to recover to 90% of pre-disturbance biomass, and forests in the central Amazon took ~40 years to recover, and the biomass recovery depends on the windthrow severity and time since disturbance33,40. Therefore, fewer large windthrows may be identified in this study in the northwestern Amazon due to relatively faster recovery. More accurate windthrow identification could be done on a yearly basis. Such analysis, however, is beyond the scope of the current study.
The function of CAPE and large windthrow density developed in this study is based on the assumption that increases in CAPE will result in increases in the frequency of extreme convective storms, which create more large windthrows. This assumption is supported by convection-permitting models that project an increase in extreme storm frequency with a warming climate41. However, the assumption may not be applied to weak or moderate convective storms, which are projected to decrease in the future41. Studies on small windthrows (<25,000 m2) and the corresponding convective storms can be carried out to explore their relationships. Future studies could involve convection-permitting models42 that are able to simulate MSCs to reduce the uncertainty in storm process behavior when exploring changes in storm frequency and intensity in the Amazon region.
The study refrains from discussing the relationship between storm intensification and forest structure changes. How the increasing CAPE in the Amazon region affects precipitation rate43 and wind gust potential44, which are the major causes of tree mortality10,20, has not been well explored. We emphasize the need for additional field observation datasets of severe convective winds and extreme rainfall rate associated with convective storms in the Amazon region. Such observations would improve the understanding of how forest dynamics link with convective storm intensity.
Methods
Identification of windthrow events
Landsat images from January 1st 2018 to December 31st 2019 were filtered on 20% or less of cloud coverage, and only the least cloudy image at each location was selected to make an image composite covering the entire Amazon region. In total, 395 least cloudy Landsat 8 images within the Amazon boundary during 2018–2019 were selected and displayed in false color (red: shortwave infrared band, green: near-infrared band, blue: red band) on Google Earth Engine for windthrow events identification (Supplementary Fig. 6). Hollow regions on Supplementary Fig. 6 (2.8% of the total area of the Amazon region) indicated that no clear images with <20% cloud were found in these areas. Spectral mixture analysis was applied on each image to detect potential windthrows (similar methods have been used in previous windthrow studies5,17). Each pixel was unmixed to fractions of image-derived endmembers, including green vegetation (GV), non-photosynthetic vegetation (NPV), and shade. GV and NPV fraction images were normalized without shade and then used to identify windthrows. Windthrows were identified manually as large fan-shape areas with high NPV fraction. Each potential windthrow was then visually checked using false color images and evaluated based on authors’ 15-year experience working with windthrow and remote sensing2,3,17,19. “New” windthrows that occurred within 1 year were spectrally more visible based on their clear fan-shape5,10 (diverging from a central area with small pixels scattered at the tail) and their relatively distinguishable reddish colors (Supplementary Fig. 7a, due to high reflectance in shortwave infrared band from woody biomass), while “old” windthrows (Supplementary Fig. 7b) occurred >1 year before the identification were displayed in bright green colors (due to reflectance in near-infrared band from the pioneer species). “Old” windthrows account for ~80% of total identified windthrows, and they were verified using historical Landsat images that can go as far as 1984 (when Landsat 5 was launched). “Old” windthrows were validated once they were found with clear shape and more distinguish color on the historical Landsat images (Supplementary Fig. 7c). 10–15% of “old” windthrows without fan-shape were eliminated from this study because it was hard to identify if they were windthrows or other types of forest disturbance. The minimum size of windthrows identified in this study was 25,000 m2. This process generated the location and rough size of 1012 visible (both “old” and “new”) windthrow scars with fan-shaped patch, scattered small disturbance pixels tails, and an area of over 25,000 m2 (Supplementary Fig. 8). Based on a gap-size probability distribution function that simulates the entire disturbance gradient from all sizes of windthrows19, the proportion of total tree mortality represented by large windthrows (>25,000 m2) identified in this study is 0.5–1.1%.
Among 1012 visual identified windthrows, the occurrence year of 125 windthrows were identified using Landsat 5,7,8, MODIS, and TRMM dataset (Supplementary Table 2), and 38 windthrows from these 125 windthrows had clear remote sensing evidence to validate their occurring date (Supplementary Table 3). It is difficult to get the accurate year and date of occurrence of all identified windthrows. Previous studies showed that windthrows in the northwestern Amazon took ~20 years to recover to 90% of “pre-disturbance” biomass from all damage classes while forests in the central Amazon took ~40 years to recover40. The biomass recovery depends on the windthrow severity and time since disturbance33. Based on the recovery time (20–40 years) and the time of windthrow identification (2018–2019), we estimated that these 1012 windthrows most likely occurred within 30 years (between 1990 and 2019), and the estimated occurrence period was validated using the range of the occurrence year (1986–2019) of 125 windthrow cases.
Windthrow density data
The windthrow density shown in Fig. 1b was generated using 1012 windthrow points in QGIS45. We created a 2.5° by 2.5° grid map, and the windthrow density was calculated by counting the number of windthrows in each grid. These values were then converted to a density with units of counts of windthrows per 10,000 km2. We chose 2.5 degrees to aggregate the data to make sure that over 50% of grids have at least 1 windthrow event while still preserving the spatial distribution of mean afternoon CAPE over the Amazon. The contour lines displayed in Fig. 1c were generated using the “Contour” function on the windthrow density map in QGIS.
Meteorological data
To derive the correlation between windthrow density and meteorological variables, we used ERA 5 global reanalysis hourly CAPE on single levels from 1979 to present at 0.25° × 0.25° resolution provided by the European Center for Medium-Range Weather Forecasts. ERA 5 CAPE was computed by considering parcels of air departing at different pressure levels below the 35 kPa level, with maximum–unstable algorithm under a pseudo-adiabatic assumption46. Afternoon mean CAPE map was calculated as the average of hourly CAPE data from 17:00–23:00 UTC (13:00–19:00 local time in Amazon) over all the months between 1990 and 2019. We chose to average CAPE over 30 years because these windthrow events occurred in these 30 years and calculating the average can help capture the overall spatial pattern of CAPE and minimize the influence of interannual climate variability on windthrow events.
To project future windthrow density in the Amazon for the end of the 21st century, we analyzed meteorological output from 10 ESMs that participated in CMIP6 (https://www.wcrp-climate.org/wgcm-cmip/wgcm-cmip6). The models used in this research were listed in Supplementary Table 1. We extracted daily surface temperature (tas), specific humidity (huss), surface pressure (ps), temperature (ta) from these models to calculate daily nondilute, near-surface-based, adiabatic CAPE. CMIP6 CAPE was calculated by considering the buoyancy of a near-surface parcel lifted adiabatically to a series of discrete pressure levels (100 kPa to 10 kPa in increments of 10 kPa). CMIP6 CAPE is calculated as follows:
Where \({{{{{\rm{d}}}}}}p\) = 10 kPa, H is the Heaviside unit step function, and \({b}_{i}=\frac{1}{{\rho }_{i}}-\frac{1}{{\rho }_{e,i}}\), with \({\rho }_{i}\) being the parcel density at pressure level i and \({\rho }_{e,i}\) being the environmental density at pressure level i.
The future projections in our analysis were based on SSP585, a high-emission scenario with high radiative forcing by the end of the century. We calculated mean daily CAPE over 1990–2015 as current CMIP6 CAPE and mean daily CAPE over 2070–2099 as future CMIP6 CAPE. Since different approaches were used to calculate ERA 5 CAPE and CMIP6 CAPE47, the absolute CAPE values of the two datasets are not comparable. Therefore, for each ESM model, we scaled future CMIP6 CAPE by multiplying, grid-wise, the delta CAPE generated from an individual model in CMIP6 with the ERA 5 current mean afternoon CAPE (Fig. 1c) as follows:
The delta CAPE indicated the projected increase in CAPE from 1990–2015 to the end of the 21st century. In this way, a scaled CMIP6 future CAPE map was generated for each model, and an ensemble-mean scaled CMIP6 CAPE map over 10 ESM models can be found in Supplementary Fig. 5b. The scaled CMIP6 future CAPE values were within plausible range compared to the ERA 5 current mean afternoon CAPE values, and both current and future CAPE maps were used to produce the increase in area with high CAPE values (>1023 J kg−1) in Table 1. However, it is worth noting that the scaling with relative changes in delta CAPE (%) is more sensitive to CMIP historical baseline conditions than absolute changes of CAPE (J kg−1), which will likely introduce a larger scaled spread (min/max CAPE changes).
The increase in area with storm-favorable environments was calculated as follows:
Where areacurrent is the area of CAPE > 1023 J kg−1 for current ERA 5 CAPE, and areafuture is the area of CAPE > 1023 J kg−1 for the scaled CMIP6 future CAPE.
A model of windthrow density
We developed a model based on the relationship between satellite-derived windthrow density and mean afternoon CAPE from the ERA 5 reanalysis over 1990–2019. The non-parametric model provides a look-up table of windthrow density as a function of CAPE within the range of observations. Counts of observed windthrow events and Amazon’s area were separately binned by CAPE using the same bins, producing two histograms of CAPE. The ratio of the former to the latter gives the density of windthrow events (windthrow events per area) as a function of CAPE. To avoid noise at the tails of the histograms, the six CAPE bins were chosen such that each bin would have about the same number of windthrow events (either 168 or 169). The total number of windthrow events is given by the sum over bins of the product of windthrow density and area. The minimum and maximum of current ERA 5 mean afternoon CAPE was 42 and 1549. The minimum CAPE value of the first bin was extended to 0 and the maximum CAPE value of the last bin was extended to infinity under the assumption that the windthrow density is similar for neighboring values. Based on the windthrow density and CAPE relationship used in the model, it is the increase in the area with high CAPE that then leads to an increase in the number of windthrow events.
It is worth noting that the future windthrow density produced by models may be underestimated because the windthrow observations within regions with high CAPE were incomplete due to high cloud coverage. Moreover, the non-parametric model makes the conservative assumption that the windthrow density does not increase at higher, as-yet unobserved values of mean afternoon CAPE.
Future projections of windthrow density
We combined the non-parametric relationship (Fig. 2a) with the future CAPE map generated from the ten CMIP6 ESMs (adjusted by ERA 5 mean CAPE values) to estimate the changes in windthrow density at the end of the 21st century. We estimated uncertainties for windthrow density projections by combining information about model-to-model differences. The analysis yielded a set of 10 estimates. The overall windthrow density increase and uncertainty were estimated using the mean increase and one standard deviation from the ensemble of the 10 models.
Data availability
The windthrow location datasets are available from the corresponding authors upon request. Current CAPE (1990–2019) datasets were retrieved from ERA5 hourly data on single levels from 1959 to present48 (available at https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview). Future CAPE datasets were calculated using daily surface temperature (tas), specific humidity (huss), surface pressure (ps), temperature (ta) variables from 10 models in CMIP6 (available at https://esgf-node.llnl.gov/search/cmip6/). Source datasets, including windthrow density datasets, processed mean ERA 5 CAPE datasets, and processed mean scaled CMIP6 future CAPE, and processed datasets to produce Fig. 2a and future windthrows are available at https://doi.org/10.15486/ngt/188360449.
Code availability
The codes used in generating this paper’s results are available at https://doi.org/10.15486/ngt/1883604.
References
Esquivel-Muelbert, A. et al. Tree mode of death and mortality risk factors across Amazon forests. Nat. Commun. 11, 5515 (2020).
Negrón-Juárez, R. I. et al. Vulnerability of Amazon forests to storm-driven tree mortality. Environ. Res. Lett. 13, 054021 (2018).
Negrón-Juárez, R. I. et al. Detection of subpixel treefall gaps with landsat imagery in Central Amazon forests. Remote Sens. Environ. 115, 3322–3328 (2011).
Chambers, J. Q., Negrón-Juárez, R. I., Hurtt, G. C., Marra, D. M. & Higuchi, N. Lack of intermediate-scale disturbance data prevents robust extrapolation of plot-level tree mortality rates for old-growth tropical forests. Ecol. Lett. 12, 22–25 (2009).
Nelson, B. W. et al. Forest disturbance by large blowdowns in the Brazilian Amazon. Ecology 75, 853–858 (1994).
Mitchell, S. J. Wind as a natural disturbance agent in forests: a synthesis. Forestry 86, 147–157 (2013).
Rehbein, A., Ambrizzi, T. & Mechoso, C. R. Mesoscale convective systems over the Amazon basin. Part I: climatological aspects. Int. J. Climatol. 38, 215–229 (2018).
Garstang, M., Harold, L., Massie, J., Halverson, J., Greco, S. & Scala, J. Amazon coastal squall lines, part i: structure and kinematics. Mon. Weather Rev. 122, 608–622 (1994).
Nunes, A. M. P., Silva Dias, M. A. F., Anselmo, E. M. & Morales, C. A. Severe convection features in the Amazon basin: a trmm-based 15-year evaluation. Front. Earth Sci. 4, 1–14 (2016).
Garstang, M., White, S., Shugart, H. H. & Halverson, J. Convective cloud downdrafts as the cause of large blowdowns in the Amazon rainforest. Meteorol. Atmos. Phys. 67, 199–212 (1998).
Schumacher, R. S. & Rasmussen, K. L. The formation, character and changing nature of mesoscale convective systems. Nat. Rev. Earth Environ. 1, 300–314 (2020).
Seeley, J. T. & Romps, D. M. Why does tropical convective available potential energy (CAPE) increase with warming? Geophys. Res. Lett. 42, 10429–10437 (2015).
Singh, M. S., Kuang, Z., Maloney, E. D., Hannah, W. M. & Wolding, B. O. Increasing potential for intense tropical and subtropical thunderstorms under global warming. Proc. Natl Acad. Sci. USA 114, 11657–11662 (2017).
Diffenbaugh, N. S., Scherer, M. & Trapp, R. J. Robust increases in severe thunderstorm environments in response to greenhouse forcing. Proc. Natl Acad. Sci. USA 110, 16361–16366 (2013).
Romps, D. M. Clausius-clapeyron scaling of CAPE from analytical solutions to RCE. J. Atmos. Sci. 73, 3719–3737 (2016).
Seeley, J. T. & Romps, D. M. The effect of global warming on severe thunderstorms in the United States. J. Clim. 28, 2443–2458 (2015).
Negrón-Juárez, R. I. et al. Widespread Amazon forest tree mortality from a single cross-basin squall line event. Geophys. Res. Lett. 37, 1–5 (2010).
Espírito-Santo, F. D. B. et al. Size and frequency of natural forest disturbances and the Amazon forest carbon balance. Nat. Commun. 5, 1–6 (2014).
Chambers, J. Q. et al. The steady-state mosaic of disturbance and succession across an old-growth central Amazon forest landscape. Proc. Natl Acad. Sci. USA 110, 3949–3954 (2013).
Lugo, A. E. Visible and invisible effects of hurricanes on forest ecosystems: an international review. Austral Ecol. 33, 368–398 (2008).
Yanoviak, S. P. et al. Lightning is a major cause of large tree mortality in a lowland neotropical forest. N. Phytol. 225, 1936–1944 (2020).
Gora, E. M., Burchfield, J. C., Muller-Landau, H. C., Bitzer, P. M. & Yanoviak, S. P. Pantropical geography of lightning-caused disturbance and its implications for tropical forests. Glob. Chang. Biol. 26, 5017–5026 (2020).
Espírito-Santo, F. D. B. et al. Storm intensity and old-growth forest disturbances in the Amazon region. Geophys. Res. Lett. 37, 1–6 (2010).
Scofield, R. The NESDIS operational convective precipitation-estimation technique. Mon. Weather Rev. 115, 1773–1793 (1987).
Romatschke, U. & Houze, R. A. Extreme summer convection in South America. J. Clim. 23, 3761–3791 (2010).
Hagos, S. M. et al. The relationship between precipitation and precipitable water in CMIP6 simulations and implications for tropical climatology and change. J. Clim. 34, 1587–1600 (2021).
Eyring, V. et al. Overview of the coupled model intercomparison project phase 6 (CMIP6) experimental design and organization. Geosci. Model. Dev. 9, 1937–1958 (2016).
Chen, J., Dai, A., Zhang, Y. & Rasmussen, K. L. Changes in convective available potential energy and convective inhibition under global warming. J. Clim. 33, 2025–2050 (2020).
Hersbach, H. et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 146, 1999–2049 (2020).
Taszarek, M., Allen, J. T., Marchio, M. & Brooks, H. E. Global climatology and trends in convective environments from ERA5 and rawinsonde data. npj Clim. Atmos. Sci. 4, 1–11 (2021).
Rennó, N. & Ingersoll, A. Natural convection as a heat engine: a theory for CAPE. J. Atmos. Sci. 53, 572–585 (1996).
Klemp, M. L. & Weisman, J. B. The dependence of numerically simulated convective storms on vertical wind shear and buoyancy. Mon. Weather Rev. 110, 504–520 (1982).
Magnabosco Marra, D. et al. Windthrows control biomass patterns and functional composition of Amazon forests. Glob. Chang. Biol. 24, 5867–5881 (2018).
Aleixo, I. et al. Amazonian rainforest tree mortality driven by climate and functional traits. Nat. Clim. Chang. 9, 384–388 (2019).
Fontes, C. G., Chambers, J. Q. & Higuchi, N. Revealing the causes and temporal distribution of tree mortality in Central Amazonia. Ecol. Manag. 424, 177–183 (2018).
Harel, M. & Price, C. Thunderstorm trends over Africa. J. Clim. 33, 2741–2755 (2020).
Brienen, R. J. W. et al. Long-term decline of the Amazon carbon sink. Nature 519, 344–348 (2015).
Gatti, L. V. et al. Amazonia as a carbon source linked to deforestation and climate change. Nature 595, 388–393 (2021).
Baccini, A. et al. Tropical forests are a net carbon source based on aboveground measurements of gain and loss. Science 358, 230–234 (2017).
David Urquiza Muñoz, J. et al. Recovery of forest structure following large-scale windthrows in the northwestern amazon. Forests 12, 1–15 (2021).
Rasmussen, K. L., Prein, A. F., Rasmussen, R. M., Ikeda, K. & Liu, C. Changes in the convective population and thermodynamic environments in convection-permitting regional climate simulations over the United States. Clim. Dyn. 55, 383–408 (2020).
Prein, A. F. et al. Simulating North American mesoscale convective systems with a convection-permitting climate model. Clim. Dyn. 55, 95–110 (2020).
Fitzpatrick, R. G. J. et al. What drives the intensification of mesoscale convective systems over the West African Sahel under climate change? J. Clim. 33, 3151–3172 (2020).
Brown, A. & Dowdy, A. Severe convective wind environments and future projected changes in Australia. J. Geophys. Res. Atmos. 126, 1–17 (2021).
QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. http://qgis.osgeo.org. (2020).
Groenemeijer, P., Púčik, T., Tsonevsky, I. & Bechtold, P. An overview of convective available potential energy and convective inhibition provided by NWP models for operational forecasting. ECMWF Tech. Memo. 1–19 (2019).
Chen, Y. et al. Future increases in Arctic lightning and fire risk for permafrost carbon. Nat. Clim. Chang. 11, 404–410 (2021).
Hersbach, H. et al. ERA5 hourly data on single levels from 1959 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). https://doi.org/10.24381/cds.adbb2d47 (2018).
Feng, Y., Negron-Juarez, R., Romps, D. & Chambers, J. Amazon windthrow disturbances are likely to increase with storm frequency under global warming: data and codes. 1.0. NGEE tropics data Collection. (dataset). https://doi.org/10.15486/ngt/1883604 (2022).
Acknowledgements
This research was supported as part of the Next Generation Ecosystem Experiments-Tropics, funded by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under contract number DE-AC02-05CH11231. We acknowledge the World Climate Research Program, which, through its Working Group on Coupled Modeling, coordinated and promoted CMIP6. We thank the climate modeling groups for producing and making available their model output, the Earth System Grid Federation (ESGF) for archiving the data and providing access, and the multiple funding agencies that support CMIP6 and ESGF.
Author information
Authors and Affiliations
Contributions
R.I.N. and J.Q.C. conceived the study. Y.F. and D.M.R. designed the research and developed the methods. R.I.N. collected the datasets of 1012 large windthrows. Y.F. and D.M.R. collected and calculated CAPE datasets from ERA 5 and CMIP6. Y.F. performed the analyses and wrote the manuscript. R.I.N., J.Q.C., D.M.R., and Y.F. have reviewed and helped to revise the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Feng, Y., Negrón-Juárez, R.I., Romps, D.M. et al. Amazon windthrow disturbances are likely to increase with storm frequency under global warming. Nat Commun 14, 101 (2023). https://doi.org/10.1038/s41467-022-35570-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-35570-1