Philippa Kaur

The shell of the endangered hawksbill sea turtle (pictured) is prized for trinkets and jewellery.Credit: Reinhard Dirscherl/SPL

Poaching is less of a threat to the survival of sea turtles than it once was, a new analysis suggests1. Illegal sea-turtle catch has dropped sharply since 2000, with most of the current exploitation occurring in areas where turtle populations are relatively healthy.This study is the first worldwide estimate of the number of adult sea turtles moved on the black market. According to the analysis, more than one million sea turtles were illegally harvested between 1990 and 2020. But the researchers also found that the illegal catch from 2010 to 2020 was nearly 30% lower than that in the previous decade.“The silver lining is that, despite the seemingly large illegal take, exploitation is not having a negative impact on sea-turtle populations on a global scale. This is really good news,” says co-author Jesse Senko, a marine conservation scientist at Arizona State University in Tempe. The research was published 7 September in Global Change Biology.Turtles for trinketsFor millennia, humans have used both adult sea turtles and their eggs as a food source and for cultural practices. In the past 200 years, however, many sea turtle populations declined steeply as hunting rose to meet a growing demand for turtle-based goods. In Europe, North America and Asia, sea-turtle shells were used to make combs, jewelry and furniture inlays. Turtles were also hunted for meat and for use in traditional medicine.The rise in turtle hunting meant that, by 2014, an estimated 42,000 sea turtles were legally harvested every year, and an unknown number of sea turtles were sold on the black market. Today, six of the seven sea-turtle species found around the globe are endangered owing to a deadly combination of habitat destruction, poaching and accidental entanglement in fishing gear.To pin down how many sea turtles were illegally harvested, Senko and his colleagues surveyed sea-turtle specialists and sifted through 150 documents, including reports from non-governmental organizations, papers in peer-reviewed journals and news articles.

Source: Ref. 1

By combining this information, the researchers made a conservative estimate that around 1.1 million sea turtles were illegally caught between 1990 and 2020. Nearly 90% of these turtles were funneled into China and Japan, largely from a handful of middle- and low-income countries (see ‘Long-distance turtle transport’). Of the species that could be identified, the most frequently exploited were the endangered green turtles (Chelonia mydas), hunted for meat, and the critically endangered hawksbill turtles (Eretmochelys imbricata), prized for their beautiful shells.However, the data also showed that the number of illegally caught turtles decreased from around 61,000 each year between the start of 2000 and the end of 2009 to around 44,000 in the past decade (see ‘More sea turtles swim free’). And, although there were exceptions, most sea turtles were taken from relatively robust populations that were both large and genetically diverse.

Source: Ref. 1

Although sea turtles seem to be doing well globally, this doesn’t mean that threats to regional populations can be ignored, says Emily Miller, an ecologist at the Monterey Bay Aquarium Research Institute in California. The study pins down where — and for whom — sea turtles are being exploited, which could help conservationists to target communities for advocacy, she says.Overall, the numbers signal that conservation efforts could be working, says Senko. “Contrary to popular belief, most sea-turtle populations worldwide are doing quite well,” he says. “The number of turtles being exploited is a shocker, but the ocean is big, and there are a lot of turtles out there.” More

Study areaSoutheastern Chongqing, China (107° 14′–109° 19′ E, 28° 9′–30° 32′ N), has an area of about 19,800 km2 (Fig. 1). The study area has a subtropical monsoon climate. And the area has four distinct seasons, with an annual average temperature of 16.2 °C and abundant rainfall, with an average annual rainfall of 1209 mm. This region is located in the central part of the Wuling mountains, which is characterized by medium and low mountainous landforms, with an average altitude of greater than 1000 m. The water system (the Wujiang River system) in the study area is well developed, with a large drainage area and rich groundwater resources. The soil is dominated by yellow soil and limestone soil, and the sensitivity to soil erosion is high. The district exhibits the typical ecological fragility of karst areas, with barren soil, fragmented surfaces, a single community, and a low ecological carrying capacity. The area includes six counties: Qianjiang district, Shizhu Tujia Autonomous county, Xiushan Tujia and Miao Autonomous county, Youyang Tujia and Miao Autonomous county, Wulong district, and Pengshui Miao and Tujia Autonomous county. The coverage rates of the carbonatite layers in these counties are 42.11, 67.77, 25.70, 34.80, 59.70 and 88.46%, respectively38, and the average coverage of the carbonatite layers is 53.09%, making this a representative area of karst rocky desertification.Data and image pre-processingIn the study, the remote sensing data were obtained from the United States Geological Survey (USGS, https://earthexplorer.usgs.gov/), including landsat-5 thematic mapper (TM) images acquired in 2001, 2006 and 2011 and Landsat-8 operational land imager (OLI) images obtained in 2016 and 2021 (Table 1). The spatial resolution is 30 m. In order to ensure the comparability of spectral characteristics, the data collection was conducted from May to September when the vegetation grew better. In order to meet the usage requirements, the cloud cover of each image used is below 10%. For the images with poor quality, the adjacent years were selected for replacement. The difference in ecological quality between adjacent years in the same region was not particularly large. In order to represent the actual situation of the ecological environment quality in the target year as much as possible, we tried to minimize the replaced part in each target year. A total of 20 images were collected in this study. The images downloaded were all L1T products, which had undergone systematic radiometric correction and geometric correction, so precise geometric correction was no longer performed. Before the subsequent processing, all 20 images were preprocessed by radiometric calibration, atmospheric correction, image mosaicking and cropping. Then these images were calculated to obtain NDVI, WET, NDBSI, LST and RI. And based on the preprocessed Landsat images, support vector machine classification was performed to obtain the land use (LU) status.Table 1 Information of images used in this study.Full size tableThe topographical data included the elevation (EV) and slope (SP) data. Among them, the elevation data was provided by the official website of the United States Geological Survey (USGS, https://earthexplorer.usgs.gov/). And the slope data was calculated from the elevation data. The meteorological data, including the monthly average temperature (MT), monthly mean precipitation (PR), monthly even relative humidity (RH), and monthly total sunshine hours (SH) from May to September of the target year, were got from the China Meteorological Data Network (http://data.cma.cn/). In addition, socioeconomic data, including the population density (PD) and gross domestic product (GDP), were obtained from the statistical yearbooks of each district and county in the study area. The nighttime light (NTL) data were obtained from the National Oceanic and Atmospheric Administration (NOAA, https://www.noaa.gov/). The above data and LU were used as the influencing factors of ecological quality to analyze the reasons for the change of local ecological environment quality. The statistical data and monitoring data of each evaluation index used to construct the EI come from the statistical yearbooks, water resources bulletin and soil and water conservation bulletin of each district and county.MethodologyStudy frameworkA framework was developed for evaluating the ecological quality in southeastern Chongqing from 2001 to 2021 in the study. And the framework included three parts: data preparation, construction of the MRSEI, and the analysis of the ecological status in the region. Figure 2 presents the detailed information about the framework. The operations of band calculation, normalization and PCA were all carried out using the ENVI 5.3 software (https://www.harrisgeospatial.com).Figure 2The study framework.Full size imageIndicators used in MRSEIThe greenness, humidity, heat, dryness, and degree of rocky desertification were used to construct the MRSEI. The NDVI39 was chosen to characterize the greenness. The humidity component acquired from the tasseled cap transformation (WET)40 was selected to represent the humidity. The LST41 was used to represent the heat, the normalized difference build-up soil index (NDBSI)42 was used to characterize the dryness. The RI was applied to characterize the degree of rocky desertification.The NDVI is an important indicator for monitoring the physical and chemical properties of vegetation, and it can be employed to calculate the vegetation coverage, leaf area index, and so on19. In addition, it eliminates some radiation errors and has a stronger response to surface vegetation. It has been widely used in vegetation remote sensing monitoring. The equation for calculating the NDVI is as follows39:$$ {text{NDVI}} = {{(uprho }}_{{{text{NIR}}}} – {uprho }_{{{text{Red}}}} {)}/{{(uprho }}_{{{text{NIR}}}} {{ + uprho }}_{{{text{Red}}}} ), $$
(1)
where ({uprho }_{{{text{NIR}}}}) is the reflectance of the near-infrared band and ({uprho }_{{{text{Red}}}}) refers to the reflectance of the red band corresponding to each image.The WET can effectively reflect the humidity conditions of the surface vegetation, water, and soil, and can reveal the changes in the ecological environment, such as soil degradation. Therefore, it is commonly used in ecological environment monitoring43. The WET can be expressed as40,43:$$ {text{WET}}_{{{text{TM}}}} { = 0}{{.3102uprho }}_{{{text{Red}}}} { + 0}{{.2021uprho }}_{{{text{Green}}}} { + 0}{{.0315uprho }}_{{{text{Blue}}}} { + 0}{{.1594uprho }}_{{{text{NIR}}}} – {0}{{.6806uprho }}_{{{text{SWIR1}}}} – {0}{{.6109uprho }}_{{{text{SWIR2}}}} , $$
(2)
$$ {text{WET}}_{{{text{OLI}}}} { = 0}{{.3283uprho }}_{{{text{Red}}}} { + 0}{{.1972uprho }}_{{{text{Green}}}} { + 0}{{.1511uprho }}_{{{text{Blue}}}} { + 0}{{.3407uprho }}_{{{text{NIR}}}} – {0}{{.7117uprho }}_{{{text{SWIR1}}}} – {0}{{.4559uprho }}_{{{text{SWIR2}}}} , $$
(3)
where ({uprho }_{{text{i}}} ,) is the reflectance of band i.The NDBSI is expressed as the average of two indicators, the bare soil index (SI)44 and the index-based built-up index (IBI)45. It can be applied to characterize the dryness. The calculation formulas are44,45:$$ {text{IBI }} = {text{ }}left[ {2uprho _{{{text{SWIR1}}}} /left( {uprho _{{{text{SWIR1}}}} + {text{ }}uprho _{{{text{NIR}}}} } right) – uprho _{{{text{NIR}}}} /(uprho _{{{text{NIR}}}} + {text{ }}uprho _{{{text{Red}}}} } right) – uprho _{{{text{Green}}}} /(uprho _{{{text{Green}}}} + {text{ }}uprho _{{{text{SWIR1}}}} )]/[2uprho _{{{text{SWIR1}}}} /left( {uprho _{{{text{SWIR1}}}} + {text{ }}uprho _{{{text{NIR}}}} } right) + {text{ }}uprho _{{{text{NIR}}}} /(uprho _{{{text{NIR}}}} + {text{ }}uprho _{{{text{Red}}}} ) + {text{ }}uprho _{{{text{Green}}}} /(uprho _{{{text{Green}}}} + {text{ }}uprho _{{{text{SWIR1}}}} )], $$
(4)
$$ {text{SI = }}left[ {{uprho }_{{{text{SWIR1}}}} {{ + uprho }}_{{{text{red}}}} – left( {{uprho }_{{{text{Blue}}}} {{ + uprho }}_{{{text{NIR}}}} } right)} right]/left[ {{uprho }_{{{text{SWIR1}}}} {{ + uprho }}_{{{text{red}}}} { + }left( {{uprho }_{{{text{Blue}}}} {{ + uprho }}_{{{text{NIR}}}} } right)} right], $$
(5)
$$ {text{NDBSI = (IBI + SI)/2,}} $$
(6)
where ({uprho }_{{text{i}}} ,) is the reflectance of band i.The LST is closely related to natural processes and human phenomena such as crop yield, vegetation growth and distribution, surface water cycle, etc. It can well reflect the state of the surface ecological environment. The atmospheric correction method is used to invert the LST here46,47, it can be expressed as:$$ {text{L = gain}} times {text{DN + bias,}} $$
(7)
$$ {text{T = K}}_{{2}} /{text{ln}}left( {frac{{{text{K}}_{{1}} }}{{text{L}}}{ + 1}} right){,} $$
(8)
$$ {text{LST = T}}/left[ {{1 + }left( {frac{{{lambda T}}}{{upalpha }}} right){{lnvarepsilon }}} right]{,} $$
(9)
where L is the radiation value in the thermal infrared band, DN is the gray value, gain and bias is the gain value and offset value of the L-band, which was got from the image header file. And T is the temperature value at the sensor; K1 and K2 are calibration parameters respectively (for TM, K1 = 607.76 W/(m2 sr μm), K2 = 1260.56 K; for TIRS, K1 = 774.89 W/(m2 sr μm), K2 = 1321.08 K); λ is the central wavelength of thermal infrared band; α = 1.438 × 10−2 m K. ε is the surface emissivity and the value is estimated by the vegetation index mixture model48,49. It is calculated as follows:$$ {text{VFC = }}frac{{{text{NDVI}} – {text{NDVI}}_{{{text{Soil}}}} }}{{{text{NDVI}}_{{{text{Veg}}}} – {text{NDVI}}_{{{text{Soil}}}} }}, $$
(10)
$$ {text{d}}_{{upvarepsilon }} { = }left( {{1} – {upvarepsilon }_{{text{s}}} } right){{ times (1}} – {text{VFC) }}times text{F} times upvarepsilon _{{text{v}}} , $$
(11)
$$ {{upvarepsilon = upvarepsilon }}_{{text{v}}} times {text{ VFC}} + varepsilon _{{text{s}}} {{ times }}left( {{1} – {text{FVC}}} right){text{ + d}}_{{upvarepsilon }} , $$
(12)
where VFC is the vegetation fractional cover, ({text{NDVI}}_{{{text{Veg}}}}) is the NDVI of the pixel covered by full vegetation and the pixels with NDVI  > 0.72 are regarded as pure vegetation pixels; ({text{NDVI}}_{{{text{Soil}}}}) is the NDVI of the bare pixel and the pixels with NDVI  More

Berendsen RL, Pieterse CMJ, Bakker PAHM. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012;17:478–86.CAS 
PubMed 
Article 

Google Scholar 
Haas D, Défago G. Biological control of soil-borne pathogens by fluorescent pseudomonads. Nat Rev Microbiol. 2005;3:307–19.CAS 
PubMed 
Article 

Google Scholar 
Whipps JM. Microbial interactions and biocontrol in the rhizosphere. J Exp Bot. 2001;52:487–511.CAS 
PubMed 
Article 

Google Scholar 
Mendes R, Kruijt M, De Bruijn I, Dekkers E, Van Der Voort M, Schneider J, et al. Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science 2011;332:1097–100.CAS 
PubMed 
Article 

Google Scholar 
Jousset A, Becker J, Chatterjee S, Karlovsky P, Scheu S, Eisenhauer N. Biodiversity and species identity shape the antifungal activity of bacterial communities. Ecology 2014;95:1184–90.PubMed 
Article 

Google Scholar 
Becker J, Eisenhauer N, Scheu S, Jousset A. Increasing antagonistic interactions cause bacterial communities to collapse at high diversity. Ecol Lett. 2012;15:468–74.PubMed 
Article 

Google Scholar 
Hu J, Wei Z, Friman VP, Gu SH, Wang XF, Eisenhauer N, et al. Probiotic diversity enhances rhizosphere microbiome function and plant disease suppression. mBio. 2016;7:e01790–16.CAS 
PubMed 
PubMed Central 

Google Scholar 
Mehrabi Z, McMillan VE, Clark IM, Canning G, Hammond-Kosack KE, Preston G, et al. Pseudomonas spp. diversity is negatively associated with suppression of the wheat take-all pathogen. Sci Rep. 2016;6:1–10.Article 
CAS 

Google Scholar 
Ma Z, Geudens N, Kieu NP, Sinnaeve D, Ongena M, Martins JC, et al. Biosynthesis, chemical structure, and structure-activity relationship of orfamide lipopeptides produced by Pseudomonas protegens and related species. Front Microbiol. 2016;7:1–16.
Google Scholar 
Yan Q, Philmus B, Chang JH, Loper JE. Novel mechanism of metabolic co-regulation coordinates the biosynthesis of secondary metabolites in Pseudomonas protegens. Elife 2017;6:e22835.PubMed 
PubMed Central 
Article 

Google Scholar 
Ramette A, Moënne-Loccoz Y, Défago G. Prevalence of fluorescent pseudomonads producing antifungal phloroglucinols and/or hydrogen cyanide in soils naturally suppressive or conducive to tobacco black root rot. FEMS Microbiol Ecol. 2003;44:35–43.CAS 
PubMed 
Article 

Google Scholar 
Raaijmakers JM, Weller DM. Natural Plant Protection by 2,4-Diacetylphloroglucinol-Producing Pseudomonas spp. in Take-All Decline Soils. Mol Plant-Microbe Interact. 1998;11:144–52.CAS 
Article 

Google Scholar 
Murata K, Suenaga M, Kai K. Genome Mining Discovery of Protegenins A–D, Bacterial Polyynes Involved in the Antioomycete and Biocontrol Activities of Pseudomonas protegens. ACS Chem Biol. 2021. https://pubs.acs.org/doi/10.1021/acschembio.1c00276. Online ahead of print.Achkar J, Xian M, Zhao H, Frost JW. Biosynthesis of Phloroglucinol. J Am Chem Soc. 2005;127:5332–3.CAS 
PubMed 
Article 

Google Scholar 
Bangera MG, Thomashow LS. Identification and Characterization of a Gene Cluster for Synthesis of the Polyketide Antibiotic 2,4-Diacetylphloroglucinol from Pseudomonas fluorescens Q2-87. J Bacteriol. 1999;181:3155–63.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bottiglieri M, Keel C. Characterization of PhlG, a hydrolase that specifically degrades the antifungal compound 2,4-diacetylphloroglucinol in the biocontrol agent Pseudomonas fluorescens CHA0. Appl Environ Microbiol. 2006;72:418–27.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Yan X, Yang R, Zhao R-X, Han J-T, Jia W-J, Li D-Y, et al. Transcriptional Regulator PhlH Modulates 2,4-Diacetylphloroglucinol Biosynthesis in Response to the Biosynthetic Intermediate and End Product. Appl Environ Microbiol. 2017;83:e01419–17.CAS 
PubMed 
PubMed Central 

Google Scholar 
Dorrestein PC, Yeh E, Garneau-Tsodikova S, Kelleher NL, Walsh CT. Dichlorination of a pyrrolyl-S-carrier protein by FADH2- dependent halogenase PltA during pyoluteorin biosynthesis. Proc Natl Acad Sci USA. 2005;102:13843–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Thomas MG, Burkart MD, Walsh CT. Conversion of L-proline to pyrrolyl-2-carboxyl-S-PCP during undecylprodigiosin and pyoluteorin biosynthesis. Chem Biol. 2002;9:171–84.CAS 
PubMed 
Article 

Google Scholar 
Schnider-Keel U, Seematter A, Maurhofer M, Blumer C, Duffy B, Gigot-Bonnefoy C, et al. Autoinduction of 2,4-diacetylphloroglucinol biosynthesis in the biocontrol agent Pseudomonas fluorescens CHA0 and repression by the bacterial metabolites salicylate and pyoluteorin. J Bacteriol. 2000;182:1215–25.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Brodhagen M, Henkels MD, Loper JE. Positive autoregulation and signaling properties of pyoluteorin, an antibiotic produced by the biological control organism Pseudomonas fluorescens Pf-5. Appl Environ Microbiol. 2004;70:1758–66.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Maurhofer M, Baehler E, Notz R, Martinez V, Keel C. Cross Talk between 2,4-Diacetylphloroglucinol-Producing Biocontrol Pseudomonads on Wheat Roots. Appl Environ Microbiol. 2004;70:1990–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Clifford JC, Buchanan A, Vining O, Kidarsa TA, Chang JH, McPhail KL, et al. Phloroglucinol functions as an intracellular and intercellular chemical messenger influencing gene expression in Pseudomonas protegens. Environ Microbiol. 2016;18:3296–308.CAS 
PubMed 
Article 

Google Scholar 
Kidarsa TA, Goebel NC, Zabriskie TM, Loper JE. Phloroglucinol mediates cross-talk between the pyoluteorin and 2,4-diacetylphloroglucinol biosynthetic pathways in Pseudomonas fluorescens Pf-5. Mol Microbiol. 2011;81:395–414.CAS 
PubMed 
Article 

Google Scholar 
Hassan KA, Johnson A, Shaffer BT, Ren Q, Kidarsa TA, Elbourne LDH, et al. Inactivation of the GacA response regulator in Pseudomonas fluorescens Pf-5 has far-reaching transcriptomic consequences. Environ Microbiol. 2010;12:899–915.CAS 
PubMed 
Article 

Google Scholar 
Dubuis C, Haas D. Cross-species GacA-controlled induction of antibiosis in pseudomonads. Appl Environ Microbiol. 2007;73:650–4.CAS 
PubMed 
Article 

Google Scholar 
Hansen ML, He Z, Wibowo M, Jelsbak L. A Whole-Cell Biosensor for Detection of 2,4- Diacetylphloroglucinol (DAPG)-Producing Bacteria from Grassland Soil. Appl Environ Microbiol. 2021;87:e01400–e01420.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hesse C, Schulz F, Bull CT, Shaffer BT, Yan Q, Shapiro N, et al. Genome‐based evolutionary history of Pseudomonas spp. Environ Microbiol. 2018;20:2142–59.CAS 
PubMed 
Article 

Google Scholar 
Lozano-Andrade CN, Strube ML, Kovács ÁT. Complete genome sequences of four soil-derived isolates for studying synthetic bacterial community assembly. Microbiol Resour Announc. 2021;10:e00848–21.CAS 
PubMed Central 
Article 

Google Scholar 
Le Roux M, Kirkpatrick RL, Montauti EI, Tran BQ, Brook Peterson S, Harding BN, et al. Kin cell lysis is a danger signal that activates antibacterial pathways of Pseudomonas aeruginosa. Elife. 2015;2015:1–65.
Google Scholar 
Tyc O, van den Berg M, Gerards S, van Veen JA, Raaijmakers JM, de Boer W, et al. Impact of interspecific interactions on antimicrobial activity among soil bacteria. Front Microbiol. 2014;5:1–10.
Google Scholar 
Qi SS, Bogdanov A, Cnockaert M, Acar T, Ranty-Roby S, Coenye T, et al. Induction of antibiotic specialized metabolism by co-culturing in a collection of phyllosphere bacteria. Environ Microbiol. 2021;23:2132–51.CAS 
PubMed 
Article 

Google Scholar 
Cornforth DM, Foster KR. Competition sensing: The social side of bacterial stress responses. Nat Rev Microbiol. 2013;11:285–93.CAS 
PubMed 
Article 

Google Scholar 
LeRoux M, Peterson SB, Mougous JD. Bacterial danger sensing. J Mol Biol. 2015;427:3744–53.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Westhoff S, van Wezel GP, Rozen DE. Distance-dependent danger responses in bacteria. Curr Opin Microbiol. 2017;36:95–101.PubMed 
Article 

Google Scholar 
Davies J, Spiegelman GB, Yim G. The world of subinhibitory antibiotic concentrations. Curr Opin Microbiol. 2006;9:445–53.CAS 
PubMed 
Article 

Google Scholar 
Garbeva P, Silby MW, Raaijmakers JM, Levy SB, Boer WDE. Transcriptional and antagonistic responses of Pseudomonas fluorescens Pf0-1 to phylogenetically different bacterial competitors. ISME J. 2011;5:973–85.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Abrudan MI, Smakman F, Grimbergen AJ, Westhoff S, Miller EL, Van Wezel GP, et al. Socially mediated induction and suppression of antibiosis during bacterial coexistence. Proc Natl Acad Sci USA. 2015;112:11054–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kehe J, Ortiz A, Kulesa A, Gore J, Blainey PC, Friedman J. Positive interactions are common among culturable bacteria. Sci Adv. 2021;7:1–10.Article 
CAS 

Google Scholar 
Yang KM, Kim JS, Kim HS, Kim YY, Oh JK, Jung HW, et al. Lactobacillus reuteri AN417 cell-free culture supernatant as a novel antibacterial agent targeting oral pathogenic bacteria. Sci Rep. 2021;11:1–16.Article 
CAS 

Google Scholar 
Dubern JF, Lugtenberg BJJ, Bloemberg GV. The ppuI-rsaL-ppuR quorum-sensing system regulates biofilm formation of Pseudomonas putida PCL1445 by controlling biosynthesis of the cyclic lipopeptides putisolvins I and II. J Bacteriol. 2006;188:2898–906.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wellington S, Peter Greenberg E. Quorum sensing signal selectivity and the potential for interspecies cross talk. mBio. 2019;10:e00146–19.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Duffy BK, Défago G. Zinc Improves Biocontrol of Fusarium Crown and Root Rot of Tomato by Pseudomonas fluorescens and Represses the Production of Pathogen Metabolites Inhibitory to Bacterial Antibiotic Biosynthesis. Phytopathology. 1997;87:1250–7.CAS 
PubMed 
Article 

Google Scholar 
Li W, Estrada-de los Santos P, Matthijs S, Xie G-L, Busson R, Cornelis P, et al. Promysalin, a Salicylate-Containing Pseudomonas putida Antibiotic, Promotes Surface Colonization and Selectively Targets Other Pseudomonas. Chem Biol. 2011;18:1320–30.CAS 
PubMed 
Article 

Google Scholar 
Parnell JJ, Berka R, Young HA, Sturino JM, Kang Y, Barnhart DM, et al. From the lab to the farm: An industrial perspective of plant beneficial microorganisms. Front Plant Sci. 2016;7:1–12.Article 

Google Scholar 
Berendsen RL, van Verk MC, Stringlis IA, Zamioudis C, Tommassen J, Pieterse CMJ, et al. Unearthing the genomes of plant-beneficial Pseudomonas model strains WCS358, WCS374 and WCS417. BMC Genomics. 2015;16:1–23.CAS 
Article 

Google Scholar 
Niu B, Paulson JN, Zheng X, Kolter R. Simplified and representative bacterial community of maize roots. Proc Natl Acad Sci USA. 2017;114:E2450–E2459.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zhuang L, Li Y, Wang Z, Yu Y, Zhang N, Yang C, et al. Synthetic community with six Pseudomonas strains screened from garlic rhizosphere microbiome promotes plant growth. Micro Biotechnol. 2021;14:488–502.CAS 
Article 

Google Scholar 
Zobel S, Benedetti I, Eisenbach L, De Lorenzo V, Wierckx N, Blank LM. Tn7-Based Device for Calibrated Heterologous Gene Expression in Pseudomonas putida. ACS Synth Biol. 2015;4:1341–51.CAS 
PubMed 
Article 

Google Scholar 
Van Gestel J, Weissing FJ, Kuipers OP, Kovács ÁT. Density of founder cells affects spatial pattern formation and cooperation in Bacillus subtilis biofilms. ISME J. 2014;8:2069–79.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hmelo LR, Borlee BR, Almblad H, Love ME, Randall TE, Tseng BS, et al. Precision-engineering the Pseudomonas aeruginosa genome with two-step allelic exchange. Nat Protoc. 2015;10:1820–41.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Yang L, Hengzhuang W, Wu H, Damkiær S, Jochumsen N, Song Z. et al. Polysaccharides serve as scaffold of biofilms formed by mucoid Pseudomonas aeruginosa. FEMS Immunol Med Microbiol. 2012;65:366–76.CAS 
PubMed 
Article 

Google Scholar  More

Andela, N. et al. A human-driven decline in global burned area. Science 356, 1356 (2017).CAS 
Article 

Google Scholar 
Sloan, S., Locatelli, B., Wooster, M. J. & Gaveau, D. L. A. Fire activity in Borneo driven by industrial land conversion and drought during El Niño periods, 1982–2010. Glob. Environ. Change 47, 95–109 (2017).Article 

Google Scholar 
Kelley, D. I. et al. How contemporary bioclimatic and human controls change global fire regimes. Nat. Clim. Change 9, 690–96 (2019).Article 

Google Scholar 
Jolly, W. M. et al. Climate-induced variations in global wildfire danger from 1979 to 2013. Nat. Commun. 6, 7537 (2015).Ward, D. S., Shevliakova, E., Malyshev, S. & Rabin, S. Trends and variability of global fire emissions due to historical anthropogenic activities. Glob. Biogeochem. Cycles 32, 122–42 (2018).CAS 
Article 

Google Scholar 
Earl, N. & Simmonds, I. Spatial and temporal variability and trends in 2001–2016 global fire activity. J. Geophys. Res. Atmos. 123, 2524–36 (2018).Article 

Google Scholar 
Giglio, L., Randerson, J. T. & van der Werf, G. R. Analysis of daily, monthly, and annual burned area using the fourth-generation Global Fire Emissions Database (GFED4). J. Geophys. Res. Biogeosci. 118, 317–28 (2013).Article 

Google Scholar 
Doerr, S. H. & Santín, C. Global trends in wildfire and its impacts: perceptions versus realities in a changing world. Philos. Trans. R. Soc. B Biol. Sci. 371, 20150345 (2016).Article 

Google Scholar 
van Lierop, P., Lindquist, E., Sathyapala, S. & Franceschini, G. Global forest area disturbance from fire, insect pests, diseases and severe weather events. Forest Ecol. Manag. 352, 78–88 (2015).Article 

Google Scholar 
Zheng, B. et al. Increasing forest fire emissions despite the decline in global burned area. Sci. Adv. 7, eabh2646 (2021).Article 

Google Scholar 
Andela, N. & van der Werf, G. R. Recent trends in African fires driven by cropland expansion and El Niño to La Niña transition. Nat. Clim. Change 4, 791–95 (2014).Article 

Google Scholar 
Van der Werf, G. R. et al. Global fire emissions and the contribution of deforestation, savanna, forest, agricultural, and peat fires (1997–2009). Atmos. Chem. Phys. 10, 11707–35 (2010).Article 
CAS 

Google Scholar 
Balch, J. K. et al. Negative fire feedback in a transitional forest of southeastern Amazonia. Glob. Change Biol. 14, 2276–87 (2008).Article 

Google Scholar 
Cochrane, M. A. & Laurance, W. F. Synergisms among fire, land use, and climate change in the Amazon. Ambio 37, 522–27 (2008).Article 

Google Scholar 
Gaveau, D. L. A. et al. Major atmospheric emissions from peat fires in Southeast Asia during non-drought years: evidence from the 2013 Sumatran fires. Sci. Rep. 4, 6112 (2014).Vadrevu, K. P. et al. Trends in vegetation fires in South and Southeast Asian countries. Sci. Rep. 9, 7422 (2019).Article 
CAS 

Google Scholar 
Sloan, S., Tacconi, L. & Cattau, M. E. Fire prevention in managed landscapes: recent successes and challenges in Indonesia. Mitig. Adapt. Strateg. Glob. Change 26, Article 32 (2021).Article 

Google Scholar 
Gaveau, D. L. A., Descales, A., Salim, M. A., Shields, D. & Sloan, S. Refined burned-area mapping protocol using Sentinel-2 data increases estimate of 2019 Indonesian burning. Earth Syst. Sci. Data, https://doi.org/10.5194/essd-2021-113, (2021).Randerson, J. T., Chen, Y., van der Werf, G. R., Rogers, B. M. & Morton, D. C. Global burned area and biomass burning emissions from small fires. J. Geophys. Res. Biogeosci. 117, G04012 (2012).Article 
CAS 

Google Scholar 
Field, R. D., van der Werf, G. R. & Shen, S. S. P. Human amplification of drought-induced biomass burning in Indonesia since 1960. Nat. Geosci. 2, 185–88 (2009).CAS 
Article 

Google Scholar 
Huijnen, V. et al. Fire carbon emissions over maritime Southeast Asia in 2015 largest since 1997. Sci. Rep. 6, 26886 (2016).CAS 
Article 

Google Scholar 
Tacconi, L. Preventing fires and haze in Southeast Asia. Nat. Clim. Change 6, 640–43 (2016).Article 

Google Scholar 
Koplitz, S. N. et al. Public health impacts of the severe haze in Equatorial Asia in September–October 2015: demonstration of a new framework for informing fire management strategies to reduce downwind smoke exposure. Environ. Res. Lett. 11, 094023 (2016).Article 

Google Scholar 
Kiely, L. et al. Air quality and health impacts of vegetation and peat fires in Equatorial Asia during 2004–2015. Environ. Res. Lett.15, 094054 (2020).Article 

Google Scholar 
Crippa, P. et al. Population exposure to hazardous air quality due to the 2015 fires in Equatorial Asia. Sci. Rep. 6, 37074 (2016).CAS 
Article 

Google Scholar 
Glauber, A. J. & Gunawan, I. The Cost of Fire: An Economic Analysis of Indonesia’s 2015 Fire Crisis. (The World Bank, Washington, D.C., (2016).Tan, Z. D., Carrasco, L. R. & Taylor, D. Spatial correlates of forest and land fires in Indonesia. Int. J. Wildland Fire 29, 1088–99 (2020).Article 

Google Scholar 
Marlier, M. E. et al. Fire emissions and regional air quality impacts from fires in oil palm, timber, and logging concessions in Indonesia. Environ. Res. Lett. 10, 085005 (2015).Article 
CAS 

Google Scholar 
Vetrita, Y. & Cochrane, M. A. Fire frequency and related land-use and land-cover changes in Indonesia’s peatlands. Remote Sens. 12, 5 (2020).Nikonovas, T., Spessa, A., Doerr, S. H., Clay, G. D. & Mezbahuddin, S. Near-complete loss of fire-resistant primary tropical forest cover in Sumatra and Kalimantan. Commun. Earth Environ. 1, 65 (2020).Article 

Google Scholar 
Field, R. Biomass burning in Indonesia: Signs of Progress in 2019?, http://www.columbia.edu/~rf2426/index_files/20200128.Field.GSFC.NoOz.pdf, January, NASA Goddard Space Flight Center, (2019).Watts, J. et al. Incentivising compliance: evaluating the effectiveness of targeted village incentives for reducing forest and peat fires. Forest Policy Econ. 108, 101956 (2019).Wijedasa, L. et al. Carbon emissions from peat forests will continue to increase despite emission-reduction schemes. Glob. Change Biol. 24, 4598–613 (2018).Article 

Google Scholar 
Sloan, S., Meyfroidt, P., Rudel, T. K. & Bongers, F. & Chazdon Robin, L. The forest transformation: Planted tree cover and regional dynamics of tree gains and losses. Glob. Environ. Change 59, 101988 (2019).Article 

Google Scholar 
Albar, I., Jaya, I. N. S., Saharjo, B. H., Kuncahyo, B. & Vadrevu, K. P. Spatio-temporal analysis of land and forest fires in Indonesia using MODIS active fire dataset, in Land-Atmospheric Research Applications in South and Southeast Asia (eds K P Vadrevu et al.), p. 105-27 (Springer International Publishing, 2018).Miettinen, J., Shi, C. & Liew, S. C. Fire distribution in Peninsular Malaysia, Sumatra and Borneo in 2015 with special emphasis on peatland fires. Environ. Manage. 60, 747–57 (2017).Article 

Google Scholar 
Fanin, T. & van der Werf, G. R. Precipitation–fire linkages in Indonesia (1997–2015). Biogeosciences 14, 3995–4008 (2017).Article 

Google Scholar 
Wiggins, E. B. et al. Smoke radiocarbon measurements from Indonesian fires provide evidence for burning of millennia-aged peat. Proc. Natl. Acad. Sci. USA 115, 12419 (2018).CAS 
Article 

Google Scholar 
Page, S. E. et al. The amount of carbon released from peat and forest fires in Indonesia during 1997. Nature 420, 61–65 (2002).CAS 
Article 

Google Scholar 
Lohberger, S., Stängel, M., Atwood, E. C. & Siegert, F. Spatial evaluation of Indonesia’s 2015 fire-affected area and estimated carbon emissions using Sentinel-1. Glob. Change Biol. 24, 644–54 (2018).Article 

Google Scholar 
van der Werf, G. R. et al. Global fire emissions estimates during 1997–2016. Earth Syst. Sci. Data 9, 697–720 (2017).Article 

Google Scholar 
Field, R. D. et al. Indonesian fire activity and smoke pollution in 2015 show persistent nonlinear sensitivity to El Niño-induced drought. Proc. Natl Acad. Sci. USA 113, 9204–09 (2016).CAS 
Article 

Google Scholar 
Austin, K. G. et al. Shifting patterns of oil palm driven deforestation in Indonesia and implications for zero-deforestation commitments. Land Use Policy 69, 41–48 (2017).Article 

Google Scholar 
Pan, X., Chin, M., Ichoku, C. & Field, R. Connecting Indonesian fires and drought with the type of El Niño and phase of the Indian Ocean Dipole during 1979–2016. J. Geophys. Res. Atmos. 123, (2018).van der Werf, G. R. et al. Climate regulation of fire emissions and deforestation in Equatorial Asia. Proc. Natl Acad. Sci. USA 105, 20350–55 (2008).Article 

Google Scholar 
Wooster, M. J., Roberts, G., Perry, G. L. W. & Kaufman, Y. J. Retrieval of biomass combustion rates and totals from fire radiative power observations: FRP derivation and calibration relationships between biomass consumption and fire radiative energy release. J. Geophys. Rese. Atmos. 110, (2005).Spessa, A. et al. Seasonal forecasting of fires over Kalimantan, Indonesia. Nat. Hazards Earth Syst. Sci. 15, 429–42 (2015).Article 

Google Scholar 
Siegert, F., Ruecker, G., Hinrichs, A. & Hoffmann, A. A. Increased damage from fires in logged forests during droughts caused by El Niño. Nature 414, 437–40 (2001).CAS 
Article 

Google Scholar 
Fernandes, K. et al. Heightened fire probability in Indonesia in non-drought conditions: the effect of increasing temperatures. Environ. Res. Lett. 12, 054002 (2017).Article 

Google Scholar 
Herawati, H. & Santoso, H. Tropical forest susceptibility to and risk of fire under changing climate: a review of fire nature, policy and institutions in Indonesia. Forest Policy Econ. 13, 227–33 (2011).Article 

Google Scholar 
Nepstad, D. et al. Slowing Amazon deforestation through public policy and interventions in beef and soy supply chains. Science 344, 1118–23 (2014).CAS 
Article 

Google Scholar 
Dennis, R. A Review of Fire Projects In Indonesia, 1982-1998. (CIFOR, Bogor, Indonesia, 1999).de Groot, W. J., Field, R. D., Brady, M. A., Roswintiarti, O. & Mohamad, M. Development of the Indonesian and Malaysian fire danger rating systems. Mitig. Adapt. Strateg. Glob. Change 12, 165 (2006).Article 

Google Scholar 
Clough, Y. et al. Land-use choices follow profitability at the expense of ecological functions in Indonesian smallholder landscapes. Nat. Commun. 7, 13137 (2016).CAS 
Article 

Google Scholar 
Bissonnette, J.-F. & De Koninck, R. The return of the plantation? Historical and contemporary trends in the relation between plantations and smallholdings in Southeast Asia. J. Peasant Stud. 44, 918–38 (2017).Article 

Google Scholar 
Gaveau, D. L. A. et al. Slowing deforestation in Indonesia follows declining oil palm expansion and lower oil prices. PLOS ONE 17, e0266178 (2022).Svatoňová, T., Herák, D. & Kabutey, A. Financial profitability and sensitivity analysis of palm oil plantation in Indonesia. Acta Univ. Agric. Silvic. Mendelianae Brunensis 63, 1365–73 (2015).Article 

Google Scholar 
Gaveau, D. L. A. et al. Rapid conversions and avoided deforestation: examining four decades of industrial plantation expansion in Borneo. Scientific Reports 6, (2016).Simamora, A. P. Govt says no to converting peatland into plantations, The Jakarta Post. August (2010).Satriastanti, F. E. Jokowi bans new oil palm and mining concessions, Mongabay.com April (2016).Sloan, S., Edwards, D. P. & Laurance, W. F. Does Indonesia’s REDD+ moratorium on new concessions spare imminently-threatened forests? Conserv. Lett. 5, 222–31 (2012).Article 

Google Scholar 
Busch, J. et al. Reductions in emissions from deforestation from Indonesia’s moratorium on new oil palm, timber, and logging concessions. Proc. Natl Acad Sci USA 112, 1328–33 (2015).CAS 
Article 

Google Scholar 
Forsyth, T. Public concerns about transboundary haze: a comparison of Indonesia, Singapore, and Malaysia. Glob. Environ. Change 25, 76–86 (2014).Article 

Google Scholar 
Carbon Conservation. Fire Free Village Program – Review 2017. (Carbon Conservation, Singapore, (2017).Gaveau, D. L. A. et al. Overlapping land claims limit the use of satellites to monitor no-deforestation committments and no-burning compliance. Conserv. Lett. 10, 257–64 (2017).Article 

Google Scholar 
EarthData. MODIS Collection 6 Active-Fire Detections standard scientific data (MCD14ML), NASA EarthData, https://earthdata.nasa.gov/firms (2019).Giglio, L., Schroeder, W. & Justice, C. O. The Collection 6 MODIS active fire detection algorithm and fire products. Remote Sens. Environ. 178, 31–41 (2016).Article 

Google Scholar 
Sloan, S., Cattau, M.E. Discrete Fire Events, their Severity, and their Ignitions, as Derived from MODIS MCD 14ML Active-Fire Detection Data for Indonesia, 2002-2019. Sean Sloan and Megan E. Cattau, Datadryad.org. (2022).Cattau, M. E. et al. Sources of anthropogenic fire ignitions on the peat-swamp landscape in Kalimantan, Indonesia. Glob. Environ. Change 39, 205–19 (2016).Article 

Google Scholar 
Wooster, M. J., Perry, G. L. W. & Zoumas, A. Fire, drought and El Niño relationships on Borneo during the pre-MODIS era (1980–2000). Biogeosciences 9, 317–40 (2012).Article 

Google Scholar 
Tansey, K., Beston, J., Hoscilo, A., Page, S. E. & Paredes Hernández, C. U. Relationship between MODIS fire hot spot count and burned area in a degraded tropical peat swamp forest in Central Kalimantan, Indonesia. J. Geophys. Res. 113, (2008).Oom, D., Silva, P. C., Bistinas, I. & Pereira, J. M. C. Highlighting biome-specific sensitivity of fire size distributions to time-gap parameter using a new algorithm for fire event individuation. Remote Sens. 8, 663 (2016).Schroeder, W. et al. Validation of GOES and MODIS active fire detection products using ASTER and ETM plus data. Remote Sens. Environ. 112, 2711–26 (2008).Article 

Google Scholar 
Hantson, S., Padilla, M., Corti, D. & Chuvieco, E. Strengths and weaknesses of MODIS hotspots to characterize global fire occurrence. Remote Sens. Environ. 131, 152–59 (2013).Article 

Google Scholar 
Tanpipat, V., Honda, K. & Nuchaiya, P. MODIS hotspot validation over Thailand. Remote Sens. 1, 1043–54 (2009).Article 

Google Scholar 
Liew, S. C., Shen, C., Low, J., Lim, A. & Kwoh, L. K. The 24th Asian Conference on Remote Sensing and 2003 International Symposium on Remote Sensing (ACRS2003). p. 671-73 (Asian Association on Remote Sensing), November 3–7.Fornacca, D., Ren, G. & Xiao, W. Performance of three MODIS fire products (MCD45A1, MCD64A1, MCD14ML), and ESA Fire_CCI in a mountainous area of northwest Yunnan, China, characterized by frequent small fires. Remote Sens. 9, 1131 (2017).Article 

Google Scholar 
Schroeder, W., Oliva, P., Giglio, L. & Csiszar, I. A. The New VIIRS 375m active fire detection data product: algorithm description and initial assessment. Remote Sens. Environ. 143, 85–96 (2014).Article 

Google Scholar 
Giglio, L., Boschetti, L., Roy, D. P., Humber, M. L. & Justice, C. O. The Collection 6 MODIS burned area mapping algorithm and product. Remote Sens. Environ. 217, 72–85 (2018).Article 

Google Scholar 
Roy, D. P., Boschetti, L., Justice, C. O. & Ju, J. The Collection 5 MODIS burned area product — Global evaluation by comparison with the MODIS active fire product. Remote Sens. Environ. 112, 3690–707 (2008).Article 

Google Scholar 
Miettinen, J., Langner, A. & Siegert, F. Burnt area estimation for the year 2005 in Borneo using multi-resolution satellite imagery. Int. J. Wildland Fire 16, 45–53 (2007).Luo, R., Hui, D., Miao, N., Liang, C. & Wells, N. Global relationship of fire occurrence and fire intensity: a test of intermediate fire occurrence-intensity hypothesis. J. Geophys. Res. Biogeosci. 122, 1123–36 (2017).Article 

Google Scholar 
Andela, N. et al. The Global Fire Atlas of individual fire size, duration, speed, and direction. Earth Syst. Sci. Data 11, 529–52 (2019).Article 

Google Scholar 
Andela, N., Morton, D. C., Giglio, L. & Randerson, J. T. Global Fire Atlas with Characteristics of Individual Fires, 2003-2016, ORNL Distributed Active Archive Center, https://doi.org/10.3334/ORNLDAAC/1642, https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1642 (2019).Field, R. D. & Shen, S. S. P. Predictability of carbon emissions from biomass burning in Indonesia from 1997 to 2006. J. Geophys. Res. Biogeosci. 113, G04024 (2008).Article 
CAS 

Google Scholar 
Fuller, D. O. & Murphy, K. The ENSO-fire dynamic in insular Southeast Asia. Clim. Change 74, 435–55 (2006).Article 

Google Scholar 
Field, R. D. et al. Development of a global fire weather database. Nat. Hazards Earth Syst. Sci. 15, 1407–23 (2015).Article 

Google Scholar 
Huffman, G. J. GPM IMERG Final Precipitation gridded data, L3 1 month 0.1 degree x 0.1 degree, version 06B. NASA Precipitation Processing System, Goddard Earth Sciences Data and Information Services Center (GES DISC). https://storm-pps.gsfc.nasa.gov/storm/; https://pmm.nasa.gov/data-access/downloads/gpm (2019).Huffman, G. J. et al. The TRMM Multisatellite Precipitation Analysis (TMPA): quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydrometeorol. 8, 38–55 (2007).Article 

Google Scholar 
Funk, C. et al. The climate hazards infrared precipitation with stations—a new environmental record for monitoring extremes. Sci. Data 2, 150066 (2015).Article 

Google Scholar 
Hsu, J., Huang, W.-R., Liu, P.-Y. & Li, X. Validation of CHIRPS precipitation estimates over taiwan at multiple timescales. Remote Sens. 13, 254 (2021).Rozante, J. R., Vila, D. A., Barboza Chiquetto, J., Fernandes, A. D. A. & Souza Alvim, D. Evaluation of TRMM/GPM blended daily products over Brazil. Remote Sens. 10, 882 (2018).Prakash, S., Mitra, A. K., Pai, D. S. & AghaKouchak, A. From TRMM to GPM: how well can heavy rainfall be detected from space? Adv. Water Resour. 88, 1–7 (2016).Article 

Google Scholar 
Ma, Q. et al. Performance evaluation and correction of precipitation data using the 20-year IMERG and TMPA precipitation products in diverse subregions of China. Atmos. Res. 249, 105304 (2021).Article 

Google Scholar 
Nwachukwu, P. N., Satge, F., Yacoubi, S. E., Pinel, S. & Bonnet, M.-P. From TRMM to GPM: how reliable are satellite-based precipitation data across Nigeria? Remote Sens. 12, 3964 (2020).Popovych, V. F. & Dunaieva, I. A. Assessment of the GPM IMERG and CHIRPS precipitation estimations for the steppe part of the Crimea. Meteorol. Hydrol. Water Manage 9, (2021).Navarro, A. et al. Assessment of IMERG precipitation estimates over Europe. Remote Sens. 11, 2470 (2019).Dezfuli, A. K. et al. Validation of IMERG precipitation in Africa. J. Hydrometeorol. 18, 2817–25 (2017).Article 

Google Scholar 
Efron, B. & Tibshirani, R. J. An Introduction to the Bootstrap. (Chapman and Hall, Boca Raton, FL, USA, 1993).Pérez-Hoyos, A., Rembold, F., Kerdiles, H. & Gallego, J. Comparison of global land cover datasets for cropland monitoring. Remote Sens. 9, 1118 (2017).ESA. Annual land-cover product, 1992 to 2019/present, based on MERIS 300-m and ancillary SPOT, AVHRR, Sentinel-3 and PROB-V satellite data. European Space Agency (ESA) European Centre for Medium-Range Weather Forecasts (ECMFW) Copernicus Climate Change Service (C3S) Climate Change Initiative (CCI), https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=overview; http://maps.elie.ucl.ac.be/CCI/viewer/download.php; http://www.esa-landcover-cci.org/ (2020).Defourny, P. Product User Guide and Specification: ICDR Land Cover 2016 to 2019 (Version 2.1.1 of ESA Coperninus Climate Change Intitiative Annual 300-m Land-Cover Classifications). (Universitie Catholique du Lovain, Louvain, Belgium, (2020).Vetrita, Y. & Cochrane, M. A. Annual Burned Area from Landsat, Mawas, Central Kalimantan, Indonesia, 1997-2015, ORNL Distributed Active Archive Center, https://doi.org/10.3334/ORNLDAAC/1708, https://daac.ornl.gov/CMS/guides/Annual_Burned_Area_Maps.html; https://daac.ornl.gov/cgi-bin/dataset_lister.pl?p=33 (2019). More

Transforming our World: The 2030 Agenda for Sustainable Development (UN, 2015).May, R. M. Will a large complex system be stable? Nature 238, 413–414 (1972).CAS 
PubMed 
Article 

Google Scholar 
Pimm, S. L. The complexity and stability of ecosystems. Nature 307, 321–326 (1984).Article 

Google Scholar 
McCann, K. S. The diversity–stability debate. Nature 405, 228–233 (2000).CAS 
PubMed 
Article 

Google Scholar 
Ives, A. R. & Carpenter, S. R. Stability and diversity of ecosystems. Science 317, 58–62 (2007).CAS 
PubMed 
Article 

Google Scholar 
Donohue, I. et al. Navigating the complexity of ecological stability. Ecol. Lett. 19, 1172–1185 (2016).PubMed 
Article 

Google Scholar 
Yachi, S. & Loreau, M. Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis. Proc. Natl Acad. Sci. USA 96, 1463–1468 (1999).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Isbell, F. et al. Biodiversity increases the resistance of ecosystem productivity to climate extremes. Nature 526, 574–577 (2015).CAS 
PubMed 
Article 

Google Scholar 
Tilman, D., Reich, P. B. & Knops, J. M. Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature 441, 629–632 (2006).CAS 
PubMed 
Article 

Google Scholar 
Loreau, M. et al. Biodiversity as insurance: from concept to measurement and application. Biol. Rev. 96, 2333–2354 (2021).PubMed 
Article 

Google Scholar 
Thibaut, L. M. & Connolly, S. R. Understanding diversity–stability relationships: towards a unified model of portfolio effects. Ecol. Lett. 16, 140–150 (2013).PubMed 
Article 

Google Scholar 
Xu, Q. et al. Consistently positive effect of species diversity on ecosystem, but not population, temporal stability. Ecol. Lett. 24, 2256–2266 (2021).PubMed 
Article 

Google Scholar 
Hector, A. et al. General stabilizing effects of plant diversity on grassland productivity through population asynchrony and overyielding. Ecology 91, 2213–2220 (2010).CAS 
PubMed 
Article 

Google Scholar 
Hautier, Y. et al. Anthropogenic environmental changes affect ecosystem stability via biodiversity. Science 348, 336–340 (2015).CAS 
PubMed 
Article 

Google Scholar 
Cardinale, B. J. et al. Biodiversity loss and its impact on humanity. Nature 486, 59–67 (2012).CAS 
PubMed 
Article 

Google Scholar 
Tilman, D., Isbell, F. & Cowles, J. M. Biodiversity and ecosystem functioning. Annu. Rev. Ecol. Evol. Syst. 45, 471–493 (2014).Article 

Google Scholar 
Isbell, F. et al. Linking the influence and dependence of people on biodiversity across scales. Nature 546, 65–72 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gonzalez, A. et al. Scaling-up biodiversity–ecosystem functioning research. Ecol. Lett. 23, 757–776 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Hooper, D. U. et al. Effects of biodiversity on ecosystem functioning: a consensus of current knowledge. Ecol. Monogr. 75, 3–35 (2005).Article 

Google Scholar 
Wang, S. & Loreau, M. Ecosystem stability in space: alpha, beta and gamma variability. Ecol. Lett. 17, 891–901 (2014).PubMed 
Article 

Google Scholar 
Wang, S. & Loreau, M. Biodiversity and ecosystem stability across scales in metacommunities. Ecol. Lett. 19, 510–518 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Wang, S. et al. Biotic homogenization destabilizes ecosystem functioning by decreasing spatial asynchrony. Ecology 102, e03332 (2021).Zhang, Y., He, N., Loreau, M., Pan, Q. & Han, X. Scale dependence of the diversity–stability relationship in a temperate grassland. J. Ecol. 106, 1227–1285 (2018).PubMed 
Article 

Google Scholar 
Wang, S., Lamy, T., Hallett, L. M. & Loreau, M. Stability and synchrony across ecological hierarchies in heterogeneous metacommunities: linking theory to data. Ecography 42, 1200–1211 (2019).Article 

Google Scholar 
Hautier, Y. et al. General destabilizing effects of eutrophication on grassland productivity at multiple spatial scales. Nat. Commun. 11, 5375 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Liang, M., Liang, C., Hautier, Y., Wilcox, K. R. & Wang, S. Grazing-induced biodiversity loss impairs grassland ecosystem stability at multiple scales. Ecol. Lett. 24, 2054–2064 (2021).PubMed 
Article 

Google Scholar 
Qiao, X. et al. Spatial asynchrony matters more than alpha stability in stabilizing ecosystem productivity in a large temperate forest region. Glob. Ecol. Biogeogr. 31, 1133–1146 (2022).Article 

Google Scholar 
Catano, C. P., Fristoe, T. S., LaManna, J. A. & Myers, J. A. Local species diversity, beta-diversity and climate influence the regional stability of bird biomass across North America. Proc. R. Soc. B 287, 20192520 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Patrick, C. J. et al. Multi‐scale biodiversity drives temporal variability in macrosystems. Front. Ecol. Environ. 19, 47–56 (2021).Article 

Google Scholar 
Wilcox, K. R. et al. Asynchrony among local communities stabilises ecosystem function of metacommunities. Ecol. Lett. 20, 1534–1545 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Zhang, Y. et al. Nitrogen addition does not reduce the role of spatial asynchrony in stabilising grassland communities. Ecol. Lett. 22, 563–571 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Garcia, R. A., Cabeza, M., Rahbek, C. & Araujo, M. B. Multiple dimensions of climate change and their implications for biodiversity. Science 344, 1247579 (2014).PubMed 
Article 
CAS 

Google Scholar 
García-Palacios, P., Gross, N., Gaitán, J. & Maestre, F. T. Climate mediates the biodiversity–ecosystem stability relationship globally. Proc. Natl Acad. Sci. USA 115, 8400–8405 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Hillebrand, H. On the generality of the latitudinal diversity gradient. Am. Nat. 163, 192–211 (2004).PubMed 
Article 

Google Scholar 
Qian, H. & Ricklefs, R. E. A latitudinal gradient in large-scale beta diversity for vascular plants in North America. Ecol. Lett. 10, 737–744 (2007).PubMed 
Article 

Google Scholar 
Kraft, N. J. B. et al. Disentangling the drivers of β diversity along latitudinal and elevational gradients. Science 333, 1755–1758 (2011).CAS 
PubMed 
Article 

Google Scholar 
Ma, Z. et al. Climate warming reduces the temporal stability of plant community biomass production. Nat. Commun. 8, 15378 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Song, J. et al. A meta-analysis of 1,119 manipulative experiments on terrestrial carbon-cycling responses to global change. Nat. Ecol. Evol. 3, 1309–1320 (2019).PubMed 
Article 

Google Scholar 
Valencia, E. et al. Synchrony matters more than species richness in plant community stability at a global scale. Proc. Natl Acad. Sci. USA 117, 24345–24351 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gilbert, B. et al. Climate and local environment structure asynchrony and the stability of primary production in grasslands. Glob. Ecol. Biogeogr. 29, 1177–1188 (2020).Article 

Google Scholar 
Hallett, L. M. et al. Biotic mechanisms of community stability shift along a precipitation gradient. Ecology 95, 1693–1700 (2014).PubMed 
Article 

Google Scholar 
Hong, P. et al. Biodiversity promotes ecosystem functioning despite environmental change. Ecol. Lett. 25, 555–569 (2022).PubMed 
Article 

Google Scholar 
. Plant presence and percent cover, RELEASE-2021. NEON (National Ecological Observatory Network) https://doi.org/10.48443/abge-r811 (2021).Barnett, D. T. et al. The plant diversity sampling design for The National Ecological Observatory. Netw. Ecosphere 10, e02603 (2019).
Google Scholar 
Lasky, J. R., Uriarte, M. & Muscarella, R. Synchrony, compensatory dynamics, and the functional trait basis of phenological diversity in a tropical dry forest tree community: effects of rainfall seasonality. Environ. Res. Lett. 11, 115003 (2016).Article 

Google Scholar 
Inchausti, P. & Halley, J. Investigating long-term ecological variability using the global population dynamics database. Science 293, 655–657 (2001).CAS 
PubMed 
Article 

Google Scholar 
Luo, M. et al. The effects of dispersal on spatial synchrony in metapopulations differ by timescale. Oikos 130, 1762–1772 (2021).Article 

Google Scholar 
Pimm, S. L. & Redfearn, A. The variability of population densities. Nature 334, 613–614 (1988).Article 

Google Scholar 
Craven, D. et al. Multiple facets of biodiversity drive the diversity–stability relationship. Nat. Ecol. Evol. 2, 1579–1587 (2018).PubMed 
Article 

Google Scholar 
Peet, R. K., Wentworth, T. R. & White, P. S. A flexible, multipurpose method for recording vegetation composition and structure. Castanea 63, 262–274 (1998).
Google Scholar 
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).Article 

Google Scholar 
Loreau, M. & de Mazancourt, C. Species synchrony and its drivers: neutral and nonneutral community dynamics in fluctuating environments. Am. Nat. 172, E48–E66 (2008).PubMed 
Article 

Google Scholar 
Lefcheck, J. S. piecewiseSEM: piecewise structural equation modelling inr for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579 (2016).Article 

Google Scholar 
Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & R Core Team. nlme: Linear and nonlinear mixed effects models. R package v.3.1–152 https://CRAN.R-project.org/package=nlme (2021).R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019). More

IPCC, 2014, Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change.Bindoff, N. L. et al. Changing Ocean, Marine Ecosystems, and Dependent Communities. IPCC Spec. Rep. Ocean Cryosph. a Chang. Clim. [H.-O. Pörtner, D.C. Roberts, V. Masson-Delmotte, P. Zhai, M. Tignor, E. Poloczanska, K. Mintenbeck, A. Alegría, M. Nicolai, A. Okem, J. Petzold, B. Rama, N.M. Weyer (eds.)]. Press 447–588 (2019).Cheng, L. et al. Upper Ocean Temperatures Hit Record High in 2020. Adv. Atmos. Sci. 38, 523–530 (2021).Article 

Google Scholar 
Smale, D. A. et al. Marine heatwaves threaten global biodiversity and the provision of ecosystem services. Nat. Clim. Chang. 9, 306–312 (2019).Article 

Google Scholar 
Hobday, A. J. et al. A hierarchical approach to defining marine heatwaves. Prog. Oceanogr. 141, 227–238 (2016).Article 

Google Scholar 
Garrabou, J. et al. Mass mortality in Northwestern Mediterranean rocky benthic communities: Effects of the 2003 heat wave. Glob. Chang. Biol. 15, 1090–1103 (2009).Article 

Google Scholar 
Frölicher, T. L. & Laufkötter, C. Emerging risks from marine heat waves. Nat. Commun. 9, 2015–2018 (2018).Article 

Google Scholar 
Oliver, E. C. J. et al. Longer and more frequent marine heatwaves over the past century. Nat. Commun. https://doi.org/10.1038/s41467-018-03732-9 (2018).Article 

Google Scholar 
Garcia-Herrera, R., Díaz, J., Trigo, R. M., Luterbacher, J. & Fischer, E. M. A review of the european summer heat wave of 2003. Crit. Rev. Environ. Sci. Technol. 40, 267–306 (2010).Article 

Google Scholar 
Marbà, N., Jordà, G., Agustí, S., Girard, C. & Duarte, C. M. Footprints of climate change on Mediterranean Sea biota. Front. Mar. Sci. 2, 56 (2015).Holbrook, N. J. et al. Keeping pace with marine heatwaves. Nat. Rev. Earth Environ. https://doi.org/10.1038/s43017-020-0068-4 (2020). in press.Article 

Google Scholar 
Oliver, E. C. J., Wernberg, T., Benthuysen, J., Chen, K. & Eds. Advances in Understanding Marine Heatwaves and Their Impacts. Lausanne: Frontiers Media SA. vol. 7 (2020).Smale, D. A. & Wernberg, T. Satellite-derived SST data as a proxy for water temperature in nearshore benthic ecology. Mar. Ecol. Prog. Ser. 387, 27–37 (2009).Article 

Google Scholar 
Schlegel, R. W., Oliver, E. C. J., Wernberg, T. & Smit, A. J. Nearshore and offshore co-occurrence of marine heatwaves and cold-spells. Prog. Oceanogr. 151, 189–205 (2017).Article 

Google Scholar 
Rutgersson, A., Jaagus, J., Schenk, F. & Stendel, M. Observed changes and variability of atmospheric parameters in the Baltic Sea region during the last 200 years. Clim Res. 61, 177–190 (2014).Liblik, T. & Lips, U. Stratification has strengthened in the baltic sea – an analysis of 35 years of observational data. Front. Earth Sci. 7, 1–15 (2019).Article 

Google Scholar 
Reusch, T. B. H. et al. The Baltic Sea as a time machine for the future coastal ocean. Sci. Adv. 4, eaar8195 (2018).Hu, S. et al. Observed strong subsurface marine heatwaves in the tropical western Pacific Ocean. Environ. Res. Lett. 16, 104024 (2021).Scannell, H. A., Johnson, G. C., Thompson, L., Lyman, J. M. & Riser, S. C. Subsurface Evolution and Persistence of Marine Heatwaves in the Northeast Pacific. Geophys. Res. Lett. 47, 1–10 (2020).Article 

Google Scholar 
Schaeffer, A. & Roughan, M. Subsurface intensification of marine heatwaves off southeastern Australia: The role of stratification and local winds. Geophys. Res. Lett. 44, 5025–5033 (2017).Article 

Google Scholar 
WMO, Guide to Climatological Practices. (2018).Hobday, A. J. et al. Categorizing and naming marine heatwaves. Oceanography 31, 162–173 (2018).Article 

Google Scholar 
Zanna, L., Khatiwala, S., Gregory, J. M., Ison, J. & Heimbach, P. Global reconstruction of historical ocean heat storage and transport. Proc. Natl. Acad. Sci. U. S. A. 116, 1126–1131 (2019).CAS 
Article 

Google Scholar 
Reynolds, R. W. et al. Daily high-resolution-blended analyses for sea surface temperature. J. Clim. 1, 5473–5496 (2007).Veneranta, L., Vanhatalo, J. & Urho, L. Detailed temperature mapping–Warming characterizes archipelago zones. Estuar. Coast. Shelf Sci. 182, 123–135 (2016).Article 

Google Scholar 
Merkouriadi, I. & Leppäranta, M. Long-term analysis of hydrography and sea-ice data in Tvärminne, Gulf of Finland, Baltic Sea. Clim. Change 124, 849–859 (2014).CAS 
Article 

Google Scholar 
Woolway, R. I. et al. Lake heatwaves under climate change. Nature 589, 402–407 (2021).CAS 
Article 

Google Scholar 
Frölicher, T. L., Fischer, E. M. & Gruber, N. Marine heatwaves under global warming. Nature 560, 360–364 (2018).Article 

Google Scholar 
Rey, J., Rohat, G., Perroud, M., Goyette, S. & Kasparian, J. Shifting velocity of temperature extremes under climate change. Environ. Res. Lett. 15, 034027 (2020).Oliver, E. C. J. et al. Marine Heatwaves. Ann. Rev. Mar. Sci. 13, 313–342 (2021).Article 

Google Scholar 
Bennett, J. M. et al. The evolution of critical thermal limits of life on Earth. Nat. Commun. 1–9 (2021) https://doi.org/10.1038/s41467-021-21263-8.Holbrook, N. J. et al. A global assessment of marine heatwaves and their drivers. Nat. Commun. 10, 1–13 (2019).CAS 
Article 

Google Scholar 
Kniebusch, M., Meier, H. E. M., Neumann, T. & Börgel, F. Temperature variability of the baltic sea since 1850 and attribution to atmospheric forcing variables. J. Geophys. Res. Ocean. 124, 4168–4187 (2019).Article 

Google Scholar 
Merkouriadi, I. & Leppäranta, M. Influence of sea ice on the seasonal variability of hydrography and heat content in Tvärminne, Gulf of Finland. Ann. Glaciol. 56, 274–284 (2015).Article 

Google Scholar 
Haapala, J. Upwelling and its influence on nutrient concentration in the coastal area of the Hanko Peninsula, entrance of the Gulf of Finland. Estuarine, Coastal and Shelf Science 38, 507–521 (1994).CAS 
Article 

Google Scholar 
Sorte, C. J. B., Fuller, A. & Bracken, M. E. S. Impacts of a simulated heat wave on composition of a marine community. Oikos 119, 1909–1918 (2010).Article 

Google Scholar 
Pansch, C. et al. Heat waves and their significance for a temperate benthic community: A near-natural experimental approach. Glob. Chang. Biol. 24, 4357–4367 (2018).Article 

Google Scholar 
Morón Lugo, S. C. et al. Warming and temperature variability determine the performance of two invertebrate predators. Sci. Rep. 10, 1–14 (2020).Article 

Google Scholar 
Humborg, C. et al. High emissions of carbon dioxide and methane from the coastal Baltic Sea at the end of a summer heat wave. Front. Mar. Sci. 6, 1–14 (2019).Article 

Google Scholar 
Laakso, L. et al. 100 Years of atmospheric and marine observations at the Finnish Utö Island in the Baltic Sea. Ocean Sci. 14, 617–632 (2018).Article 

Google Scholar 
Høyer, J. L. & Karagali, I. Sea surface temperature climate data record for the North Sea and Baltic Sea. J. Clim. 29, 2529–2541 (2016).Article 

Google Scholar 
Schlegel, R. W. & Smit, A. J. heatwaveR: A central algorithm for the detection of heatwaves and cold-spells. J. Open Source Softw. 3, 821 (2018).Article 

Google Scholar 
Schlegel, R. W., Oliver, E. C. J., Hobday, A. J. & Smit, A. J. Detecting Marine Heatwaves With Sub-Optimal Data. Front. Mar. Sci. 6, 1–14 (2019).Article 

Google Scholar  More

MaterialsInformation on wood volume and the physical environment of the plots were obtained from the US Forest Service Forest Inventory and Analysis (USFS-FIA)22. The FIA database categorizes each plot into one of 33 forest groups, but 23 groups do not have sufficient data in the control period (before 1990) to enable robust matching and so were dropped from this study. As a result, several western forest groups (e.g., Douglas-fir) were not included in our study. The following ten forest groups [(1) Loblolly/Shortleaf Pine, (2) Slash/Shortleaf Pine, (3) White/Red/Jack Pine, (4) Spruce/Fir, (5) Elm/Ash/Cottonwood, (6) Maple/Beech/Birch, (7) Oak/Hickory, (8) Oak/Gum/Cypress, (9) Aspen/Birch, and (10) Oak/Pine] all had more than 5000 observations and large numbers of observations both from before 1990 and from 2000 on. Data for the 48 conterminous states from evaluation years between 1968 and 2018 were included in the study. We limited our analysis to plots with trees from 1 to 100 years of age, resulting in trees that had been planted somewhere between 1869 and 2018—a period during which atmospheric CO2 increased from roughly 287 to more than 406 ppm32,33,34. The geographic distribution of the ten forest groups presented in Fig. 2 shows in orange all counties in which the USFS recorded in at least one year between 1968 and 2018 the presence of a plot of the respective forest group that met the age requirements for inclusion in this study. Precipitation and temperature data were obtained from the PRISM Climate Group41.MethodsResults in Tables 1 and 2 are based on estimated exponential tree-volume functions of the generalized form shown in Eq. 1. The left-hand side is the natural log of the volume per hectare in the central stem of trees on each plot in cubic meters. Volume is assumed to be a function of age, the logged cumulative lifetime concentration of CO2, and other variables, including plot-specific variables that vary across plots but not time (Xi), weather variables that vary across plots and time (Wit), and time-specific fixed effects that vary across time but not plots (Et).$${{{{mathrm{Ln}}}}},{left(frac{{{{{{rm{Volume}}}}}}}{{{{{{rm{Hectare}}}}}}}right)}_{it}= ,alpha+{beta }_{0}frac{1}{{{{{{{rm{Age}}}}}}}_{{{{{{rm{it}}}}}}}}+{beta }_{1},{{{{mathrm{Ln}}}}}({{{{{rm{CumCO}}}}}}2{{{{{{rm{Life}}}}}}}_{{{{{{rm{t}}}}}}})\ +{beta }_{2}{{{{{{rm{X}}}}}}}_{{{{{{rm{i}}}}}}}+{beta }_{3}{{{{{{rm{W}}}}}}}_{{{{{{rm{it}}}}}}}+{beta }_{4}{{{{{{rm{E}}}}}}}_{{{{{{rm{t}}}}}}}+{varepsilon }_{it}$$
(1)
The nonparametric smearing estimate method was used to transform logged-volume results into a volume in cubic meters per hectare42. The climate variables, obtained from the PRISM Climate Group41 and described in Supplementary Table 1, enter as cubic polynomials of the lifetime seasonal temperature and precipitation averages that a plot of a given age at a given time experienced.The variable for atmospheric carbon was constructed as the logarithmic transformation of the sum of yearly atmospheric CO2 exposures over the lifetime of the stand. Other site-specific covariates were obtained from the FIA data (Supplementary Table 2), such as the availability of water, the quality of the soil, the photoperiod of the plot, whether disturbances had impacted the land, and whether the land was publicly or privately owned43,44.The time-specific fixed effects (Et) in the model control for episodic factors like nitrogen deposition and invasive species, which are correlated with time but cannot be observed over space for the whole time period. These time-dummy variables account for underlying, unobservable systematic differences between the 21st-century period when atmospheric CO2 was higher and the pre-period when levels were much lower. Controlling for these factors aids the identification of the impact of elevated CO2, which varies annually.A potential concern is that wood volume changes over time could be related to an increased number of trees per hectare rather than increased wood volume of the trees. To assess whether controls for the stocking condition were needed, we examined data on the number of trees per acre of each forest type. First, we looked at a group of southern states (Supplementary Table 3) and found double-digit percentage changes in tree stocking between 1974 and 2017 for seven of the nine forest groups. However, the changes were mixed, with four having increased tree density and five decreasing tree density. The FIA data do not record the Aspen/Birch forest group as present in these southern states in these evaluations.Examination of a group of northern states involved a comparison of the average stocking conditions around 1985 with those in 2017. The changes in tree density for these forest types (Supplementary Table 4) were also split with four showing increased stocking and five having less dense stocking. The change for Loblolly/Shortleaf pine was relatively large, with stocking density increasing by 27.2%. Slash/Longleaf was not recorded as present in these states in these evaluations.Next, we analyzed changes, over the period from around 1985 to 2017, in all states east of the 100th meridian, as those states comprised the bulk of the data in our study (Supplementary Table 5). Results for seven of the ten forest groups showed a less dense composition. Loblolly/Shortleaf pine again was shown to have become more densely stocked, with an increase of 13.2%.The last check included all of the 48 conterminous states and compared changes in stocking conditions from years around 1985 to 2017 (Supplementary Table 6). Seven of the ten forest groups showed decreased stocking density over time. Not surprisingly (because most Loblolly/Shortleaf is located in the Eastern US), the change in Loblolly/Shortleaf pine density is the same for this check as was shown in the results in Supplementary Table 5. Based on the results from all these comparisons and given that stocking density has changed over time, we controlled for it both in the matching and in the multivariate-regression analysis.Genetic matching (GM), the primary approach used for this analysis, combines propensity score matching and Mahalanobis matching techniques45. The choice of GM was made after initially considering other approaches, such as nearest-neighbor propensity score matching with replacement and a non-matching, pooled regression approach. These three options were tested on the samples for Loblolly/Shortleaf pine and Oak/Hickory, and the regression results are presented in Supplementary Data 3-4.The results across these different approaches were quite similar, suggesting that the results are not strongly driven by methodological choice. We focused on matching rather than a pooled regression approach to help reduce bias and provide estimates closer to those that would be obtained in a randomized controlled trial. When choosing the specific matching approach, we considered that standard matching methods are equal percent bias reducing (EPBR) only in the unlikely case that the covariate distributions are all roughly normal46 and that EPBR may not be desirable, as in the case where one of two covariates has a nonlinear relationship with the dependent variable16. We also noted that GM is a matching algorithm that at each step minimizes the largest bias distance of the covariates24 and that GM has been shown to be a more efficient estimator than other methods like the inverse probability of treatment weighting and one-to-one greedy nearest-neighbor matching24,47,48,49. Additionally, when the distributions of covariates are non-ellipsoidal, this nonparametric method has been shown to minimize bias that may not be captured by simple minimization of mean differences50. Lastly, as sample size increases, this approach will converge to a solution that reduces imbalance more than techniques like full or greedy matching48,51,52. Given the support that this choice has in the literature, we decided to employ GM to create all the matched data used in this study using R software53.Artificial regeneration of forest stands, noted as planting throughout the text is used as the main proxy for the impact of forest management. The other indicator of management activity is what can be described as interventions, which are a range of human on-site activities that the USFS details22. We define unmanaged land as stands with natural regeneration and where no interventions occurred on the plot.To create Table 1, we first excluded all plots on which there had been either planting activities or some type of human intervention. Then, we created treatment and control groups by forming two time periods separated by an intervening period of ten years to ensure a more than a marginal difference between the groups in terms of lifetime exposure to atmospheric CO2. The control period used forest plot data sampled between 1968 and 1990, and the treatment period used forest plots sampled between 2000 and 2018. Note that even though the earlier period contains more years, there are fewer overall observations.Matches were then made to balance the treatment and control groups based on the following observable covariates: (1) Seasonal Temperature, (2) Seasonal Precipitation, (3) Stocking Condition, (4) Aspect, (5) Age, (6) Physiographic Class, and (7) Site Class. The propensity score was defined as a logit function of the above covariates to generate estimates of the probability of treatment. Calipers with widths less than or equal to 0.2 standard deviations of the propensity score were also employed to remove at least 98% of bias49.Balance statistics for the primary covariates are presented in Supplementary Data 1–2 and show a strong balance for all covariates across all forest groups. Thus for each forest group, our sample of plots includes control plots (pre-1990) and treatment plots (post-2000) that are comparable (balanced) in climate and other biophysical attributes.After trimming our sample using this matching process and obtaining strongly balanced matches, we turned to regression analysis, where we employed Stata software54. To confirm that we had the most appropriate model structure, tests of the climate and atmospheric carbon variables were undertaken using various polynomial forms, and the main variable of interest, atmospheric carbon, was tested both using a linear lifetime cumulative CO2 variable and a logarithmic transformation of that variable. Results (Supplementary Data 5–10) show that the climate variables were not improved with complexity beyond cubic form. Moreover, selection tools, like the Akaike and Bayesian information criterion, favored the cubic choice, and so we utilized the cubic formulation throughout this study. Results for the CO2 variable were similar in both sign and significance for the linear and logged form. We use the logged form as it allows easier interpretation of the effect, suppresses heteroscedasticity, and removes the assumption that each unit increase in CO2 exposure will have a linear (constant) effect on volume.The estimated effect of CO2 exposure for each forest group (Supplementary Data 12–21) was estimated using alternate specifications of the independent variables included in Eq. 1. For each forest type, the Model (1) specification (Eq. 2) is the basis for the results presented in Table 1. The β0 coefficient details the impact on the volume of the main variable of interest, atmospheric carbon.$${{{{mathrm{Ln}}}}}left(frac{volume}{hectare}right)= alpha+{beta }_{0},{{{{mathrm{Ln}}}}}({{{{{{rm{Lifetime}}}}}}{{{{{rm{CO}}}}}}}_{2})+{beta }_{1}frac{1}{{{{{{rm{Age}}}}}}}+{beta }_{2}{{{{{rm{Site}}}}}},{{{{{rm{Class}}}}}}\ +{beta }_{3}{{{{{rm{Seasonal}}}}}},{{{{{rm{Temperature}}}}}}+{beta }_{4}{{{{{rm{Seasonal}}}}}},{{{{{{rm{Temp}}}}}}}^{2}+{beta }_{5}{{{{{rm{Seasonal}}}}}},{{{{{{rm{Temp}}}}}}}^{3}\ +{beta }_{6}{{{{{rm{Seasonal}}}}}},{{{{{rm{Precipitation}}}}}}+{beta }_{7}{{{{{rm{Seasonal}}}}}},{{{{{{rm{Precip}}}}}}}^{2}+{beta }_{8}{{{{{rm{Seasonal}}}}}},{{{{{{rm{Precip}}}}}}}^{3}\ +{beta }_{9}{{{{{rm{Stocking}}}}}}+{beta }_{10}{{{{{rm{Disturbances}}}}}}+{beta }_{11}{{{{{rm{Physiographic}}}}}},{{{{{rm{Class}}}}}}+{beta }_{12}{{{{{rm{Aspect}}}}}}\ +{beta }_{13}{{{{{rm{Slope}}}}}}+{beta }_{14}{{{{{rm{Elevation}}}}}}+{beta }_{15}{{{{{rm{Latitude}}}}}}+{beta }_{16}{{{{{rm{Longitude}}}}}}+{beta }_{17}{{{{{rm{Ownership}}}}}}\ +{beta }_{18}{{{{{rm{Time}}}}}},{{{{{rm{Dummies}}}}}}+{beta }_{19}{{{{{rm{Seasonal}}}}}},{{{{{rm{Vapor}}}}}},{{{{{rm{Pressure}}}}}},{{{{{rm{Deficit}}}}}}\ +{beta }_{20}{{{{{rm{Length}}}}}},{{{{{rm{of}}}}}},{{{{{rm{Growing}}}}}},{{{{{rm{Season}}}}}}+{{{{{rm{varepsilon }}}}}}$$
(2)
After estimating Eq. 2 for each forest type individually (Supplementary Data 12–21), all plots were pooled across forest groups, with additional forest-group dummy variables, to estimate a general tree-volume function (Supplementary Data 22).Our main Model (1) results are provided in Supplementary Data 12–22, along with three additional models that assess the robustness of the elevated CO2 effect to different specifications. The simplest specification, Model (4), included only stand age, CO2 exposure, and a time-dummy variable. Model (3) took the Model (4) base and added in an array of site-specific variables, including those for the climate. Model (2) was similar to Model (1) in that it included the impact of vapor pressure deficit and the length of the growing season on the variables included in Model (3), but it differed from Model (1) in that it tested an alternate approach to capturing the impact of underlying, unobservable systematic differences like nitrogen deposition.Using the estimated coefficients from the preferred Model (column 1) specification (Eq. 2), the estimated change in growing-stock volume between two CO2 exposure scenarios was calculated at ages 25, 50, and 75. The first scenario examined CO2 exposure up to 1970 (that is, when calculating growing-stock volume for a 25-year-old stand, the CO2 exposure would have the summation of the yearly values for the years from 1946 to 1970 [310 to 326 ppm CO2]). The second scenario examined CO2 exposure up to 2015 (that is, when calculating growing-stock volume for a 25-year-old stand, the CO2 exposure was the summation of the yearly values for the years from 1991 to 2015 [347 to 401 ppm CO2])32,33,34. In both scenarios, climate variables were maintained at their 1970 exposure levels, covering the same historical years (e.g., for a 25-year-old stand, 1946 to 1970 were the years of interest), while using seasonal, not annual values and calculating average values, not lifetime summations.Forest dynamics in the Western US differ from those in the East (e.g., generally drier conditions; greater incidence of large wildfires) and as most of the observations for this study are of forest groups located in the 33 states that the USFS labels as comprising the Eastern US, robustness tests were conducted to assess whether results would differ were only eastern observations utilized. Three forest groups [(1) Loblolly/Shortleaf pine, (2) Oak/Gum/Cypress, and (3) Slash/Longleaf pine] have no observations in the Western US. A fourth, White/Red/Jack Pine, has a slight presence in a few Western states, but no western observations were selected in the original matching process (Supplementary Data 2). For the other six forest groups, all observations from Western US states were dropped. As can be seen from Fig. 2, this had the biggest impact on Aspen/Birch and Elm/Ash/Cottonwood. With this data removed, the GM matching algorithm was again used. Balance statistics are presented in Supplementary Data 23 and again show a strong balance for all covariates across all forest groups. With matches made, the average treatment effect on the treated was estimated using the Model (1) specification used to create Table 1. Regression results are presented in Supplementary Data 24,25, and a revised version of Table 1 for just the observations from the Eastern US is presented as Supplementary Table 7.As an additional robustness check on the results in Table 1, we tested an alternative functional form of the volume function. This alternative volume function is shown in Eq. 3. It has a similar shape as the function used for the main results in the paper, however, this equation cannot be linearized with logs in a similar way. Thus, it was estimated with nonlinear least squares, using the matched samples of naturally regenerated forests for individual forest groups, as well as the aggregated sample.$$frac{{{{{{mathrm{Volume}}}}}}}{{{{{{mathrm{Hectare}}}}}}}=a/(b+exp (-c,ast ,{{{{{rm{Age}}}}}}))$$
(3)
We began by estimating two separate growth functions, one for the pre-1990 (low CO2) period and one for the post-2000 (high CO2) period using Eq. 3. That is, observations from the pre-1990 (low CO2) control period and from the post-2000 (high CO2) treatment period were handled in separate regressions. For this initial analysis with the nonlinear volume function, we did not control for CO2 concentration or other factors that could influence volume across sites (e.g., weather, soils, slope, aspect), and thus, results likely show the cumulative impact of these various factors. Using the regression results (Supplementary Data 26), we calculated the predicted volume for the pre-1990 and post-2000 periods and compared the predicted volumes (Supplementary Table 8).Next, we tested this yield function on the combined sample (containing both control and treatment observations) and all forest groups. Here the model was expanded to better identify the impact of elevated CO2 by including all covariates. Instead of using a dummy variable for each forest group, though, a single dummy variable was used to differentiate hardwoods from softwoods. Once again, the equation was logarithmically transformed for ease of comparison with the results presented in Table 1. All covariates were originally input, but those which were not significant were removed. That process yielded the functional form shown in Eq. 4. Results for the regression are presented in Supplementary Data 27. The predicted change in volume due to CO2 fertilization from 1970 to 2015 is shown in Supplementary Table 9.$$frac{{{{{{mathrm{Volume}}}}}}}{{{{{{mathrm{Hectare}}}}}}}= big(a0+a1,ast ,{{{{{rm{Time}}}}}},{{{{{rm{Dummy}}}}}}+a2,ast ,{{{{mathrm{Ln}}}}}({{{{{rm{LifetimeCO}}}}}}2)+a{3},ast ,{{{{mathrm{Ln}}}}}({{{{{rm{Seasonal}}}}}},{{{{{rm{Temperature}}}}}})\ +a{4},ast ,{{{{mathrm{Ln}}}}}({{{{{rm{Seasonal}}}}}},{{{{{rm{Precipitation}}}}}})+a{5},ast ,{{{{{rm{Site}}}}}},{{{{{rm{Class}}}}}}\ +a6,ast ,{{{{{rm{Physiographic}}}}}},{{{{{rm{Dummy}}}}}}+a{7},ast ,{{{{{rm{Aspect}}}}}},{{{{{rm{Dummy}}}}}}+a{8},ast ,{{{{{rm{Stocking}}}}}},{{{{{rm{Code}}}}}}\ +a9,ast ,{{{{{rm{Disturbances}}}}}}+a{10},ast ,{{{{{rm{Hardwood}}}}}}/{{{{{rm{Softwood}}}}}},{{{{{rm{Dummy}}}}}}left.right) /left(right.b{0}+b{1},ast ,{{{{{rm{Time}}}}}},{{{{{rm{Dummy}}}}}}\ +b{2},ast ,{{{{mathrm{Ln}}}}}({{{{{rm{Lifetime}}}}}},C{O}_{2})+b3,ast ,{{{{mathrm{Ln}}}}}({{{{{rm{Seasonal}}}}}},{{{{{rm{Temperature}}}}}})\ +b{4},ast ,{{{{mathrm{Ln}}}}}({{{{{rm{Seasonal}}}}}},{{{{{rm{Precipitation}}}}}})+b5,ast ,{{{{{rm{Site}}}}}},{{{{{rm{Class}}}}}}\ +b6,ast ,{{{{{rm{Physiographic}}}}}},{{{{{rm{Dummy}}}}}}+b{7},ast ,{{{{{rm{Aspect}}}}}},{{{{{rm{Dummy}}}}}}+b8,ast ,{{{{{rm{Stocking}}}}}},{{{{{rm{Code}}}}}}\ +b9,ast ,{{{{{rm{Disturbances}}}}}}+b{10},ast ,{{{{{rm{Hardwood}}}}}}/{{{{{rm{Softwood}}}}}},{{{{{rm{Dummy}}}}}}\ +exp left(right.-left(right.c{0}+c{1},ast ,{{{{{rm{Time}}}}}},{{{{{rm{Dummy}}}}}}+c{2},ast ,{{{{{rm{Lifetime}}}}}},{{{{{{rm{CO}}}}}}}_{2}\ +c{3},ast ,{{{{mathrm{Ln}}}}}({{{{{rm{Seasonal}}}}}},{{{{{rm{Temperature}}}}}})+c{4},ast ,{{{{mathrm{Ln}}}}}({{{{{rm{Seasonal}}}}}},{{{{{rm{Precipitation}}}}}})+c{5},ast ,{{{{{rm{Site}}}}}},{{{{{rm{Class}}}}}}\ +c{6},ast ,{{{{{rm{Physiographic}}}}}},{{{{{rm{Dummy}}}}}}+c{7},ast ,{{{{{rm{Aspect}}}}}},{{{{{rm{Dummy}}}}}}+c{8},ast ,{{{{{rm{Stocking}}}}}},{{{{{rm{Code}}}}}}\ +c{9},ast ,{{{{{rm{Disturbances}}}}}}+c{10},ast ,{{{{{rm{Hardwood}}}}}}/{{{{{rm{Softwood}}}}}},{{{{{rm{Dummy}}}}}}left.right),ast ,{{{{{rm{Age}}}}}}left.right)left.right)$$
(4)
As the results using the nonlinear volume functions were similar in sign and magnitude to the multivariate-regression results and as the practice of matching and then running a multivariate-regression represents a doubly robust econometric approach that has been shown to yield results that are robust to misspecification in either the matching or the regression model47,55,56,57, the main text results are based on estimations utilizing multivariate-regression analysis post-matching.To develop Table 2, which compares naturally regenerated stands with planted stands, we used the same general approach as was used to create Table 1. The analysis and comparison of planted and naturally regenerated stands was conducted only for stands with enough observations of both to make a comparison: White/Red/Jack, Slash/Longleaf, and Loblolly/Shortleaf pine. We followed the same matching and regression procedures as above, but conducted the matching separately for naturally regenerated and planted stands. We also limited the data to stands less than or equal to 50 years of age, as there are few planted stands of older ages due to the economics of rotational forestry35,36,37,38,39,40. Balance statistics for the matched samples are presented in Supplementary Data 28–30. Again, the matching process resulted in a good balance in observable plot characteristics, which implies that we achieved comparable treatment and control plots.Using the matched data, we estimated the same regression as in Eq. 2. Estimation results, which use the Model (2) specification from Supplementary Data 19–21 that was used with the data for these three forest groups from ages 1–100, are presented in Supplementary Data 30–32. A comparison of the parameter estimates on the natural log of lifetime CO2 exposure between the results for ages 1–50 (from Supplementary Tables 31–33) and those for ages 1–100 (from Supplementary Data 19–21) is presented in Supplementary Table 10.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Beer, C. et al. Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate. Science 329, 834 (2010).Article 

Google Scholar 
Anderegg, W. R. L. et al. Pervasive drought legacies in forest ecosystems and their implications for carbon cycle models. Science 349, 528–532 (2015).Article 

Google Scholar 
Ciais, P. et al. Europe-wide reduction in primary productivity caused by the heat and drought in 2003. Nature 437, 529–533 (2005).Article 

Google Scholar 
Peñuelas, J. et al. Shifting from a fertilization-dominated to a warming-dominated period. Nat. Ecol. Evol. 1, 1438–1445 (2017).Article 

Google Scholar 
Morin, X. et al. Temporal stability in forest productivity increases with tree diversity due to asynchrony in species dynamics. Ecol. Lett. 17, 1526–1535 (2014).Article 

Google Scholar 
Isbell, F. et al. Biodiversity increases the resistance of ecosystem productivity to climate extremes. Nature 526, 574 (2015).Article 

Google Scholar 
De Boeck, H. J. et al. Patterns and drivers of biodiversity–stability relationships under climate extremes. J. Ecol. 106, 890–902 (2018).Article 

Google Scholar 
Grossiord, C. Having the right neighbors: how tree species diversity modulates drought impacts on forests. N. Phytol. 228, 42–49 (2020).Article 

Google Scholar 
O’Brien, M. J. et al. Resistance of tropical seedlings to drought is mediated by neighbourhood diversity. Nat. Ecol. Evol. 1, 1643–1648 (2017).Article 

Google Scholar 
Gazol, A. & Camarero, J. J. Functional diversity enhances silver fir growth resilience to an extreme drought. J. Ecol. 104, 1063–1075 (2016).Article 

Google Scholar 
Pretzsch, H., Schütze, G. & Uhl, E. Resistance of European tree species to drought stress in mixed versus pure forests: evidence of stress release by inter-specific facilitation. Plant Biol. 15, 483–495 (2013).Article 

Google Scholar 
Grossiord, C. et al. Tree diversity does not always improve resistance of forest ecosystems to drought. P. Natl Acad. Sci. USA 111, 14812–14815 (2014).Article 

Google Scholar 
Grossiord, C. et al. Does drought influence the relationship between biodiversity and ecosystem functioning in boreal forests. Ecosystems 17, 394–404 (2014).Article 

Google Scholar 
Loreau, M., Mouquet, N. & Gonzalez, A. Biodiversity as spatial insurance in heterogeneous landscapes. P. Natl Acad. Sci. USA 100, 12765 (2003).Article 

Google Scholar 
Lloret, F. et al. Woody plant richness and NDVI response to drought events in Catalonian (northeastern Spain) forests. Ecology 88, 2270–2279 (2007).Article 

Google Scholar 
He, Q. & Bertness, M. D. Extreme stresses, niches, and positive species interactions along stress gradients. Ecology 95, 1437–1443 (2014).Article 

Google Scholar 
Hafner, B. D. et al. Hydraulic redistribution under moderate drought among English oak, European beech and Norway spruce determined by deuterium isotope labeling in a split-root experiment. Tree Physiol. 37, 950–960 (2017).Article 

Google Scholar 
Forrester, D. I. & Bauhus, J. A review of processes behind diversity–productivity relationships in forests. Curr. For. Rep. 2, 45–61 (2016).Article 

Google Scholar 
Vitali, V., Forrester, D. I. & Bauhus, J. Know your neighbours: drought response of Norway spruce, silver fir and Douglas fir in mixed forests depends on species identity and diversity of tree neighbourhoods. Ecosystems 21, 1215–1229 (2018).Article 

Google Scholar 
The State of the World’s Forests 2020: Forests, Biodiversity and People (FAO and UNEP, 2020).Zhang, J., Fu, B., Stafford-Smith, M., Wang, S. & Zhao, W. Improve forest restoration initiatives to meet Sustainable Development Goal 15. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-020-01332-9 (2020).Schulze, K., Malek, Ž. & Verburg, P. H. Towards better mapping of forest management patterns: a global allocation approach. For. Ecol. Manage. 432, 776–785 (2019).Article 

Google Scholar 
Harris, N. L. et al. Global maps of twenty-first century forest carbon fluxes. Nat. Clim. Change https://doi.org/10.1038/s41558-020-00976-6 (2021).Maxwell, S. L. et al. Area-based conservation in the twenty-first century. Nature 586, 217–227 (2020).Article 

Google Scholar 
Williams, A. P. et al. Large contribution from anthropogenic warming to an emerging North American megadrought. Science 368, 314–318 (2020).Article 

Google Scholar 
Blackman, C. et al. Leaf hydraulic vulnerability is related to conduit dimensions and drought resistance across a diverse range of woody angiosperms. N. Phytol. 188, 1113–1123 (2010).Article 

Google Scholar 
Liang, J. et al. Positive biodiversity–productivity relationship predominant in global forests. Science 354, aaf8957 (2016).Article 

Google Scholar 
Wieczynski, D. J. et al. Climate shapes and shifts functional biodiversity in forests worldwide. P. Natl Acad. Sci. USA 116, 587–592 (2019).Article 

Google Scholar 
Tomppo, E. et al. National Forest Inventories: Pathways for Common Reporting (Springer, 2010).Chirici, G. et al. National Forest Inventories: Contributions to Forest Biodiversity Assessments (Springer, 2011).Magnussen, S., Smith, B. & Uribe, S. National Forest inventories in North America for monitoring forest tree species diversity. Plant Biosyst. 141, 113–122 (2007).Article 

Google Scholar 
Lesiv, M. et al. Global forest management data for 2015 at a 100 m resolution. Sci. Data 9, 199 (2022).Article 

Google Scholar 
Vicente-Serrano, S. M., Beguería, S. & López-Moreno, J. I. A. Multiscalar drought index sensitive to global warming: the Standardized Precipitation Evapotranspiration Index. J. Clim. 23, 1696–1718 (2010).Article 

Google Scholar 
Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850 (2013).Article 

Google Scholar 
Forest Resources Assessment 2015 (FAO, 2015).Lyapustin, A. et al. Scientific impact of MODIS C5 calibration degradation and C6+ improvements. Atmos. Meas. Tech. 7, 4353–4365 (2014).Article 

Google Scholar 
Didan, K. & Brreto, A. NASA MEaSUREs Vegetation Index and Phenology (VIP) Phenology EVI2 Yearly Global 0.05Deg CMG, NASA EOSDIS Land Processes DAAC, https://doi.org/10.5067/MEaSUREs/VIP/VIPPHEN_EVI2.004 (2016).Olson, D. M. et al. Terrestrial Ecoregions of the World: A New Map of Life on Earth: a new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity. Bioscience 51, 933–938 (2001).Article 

Google Scholar 
Kline, T. J. B. Sample issues, methodological implications, and best practices. Can. J. Behav. Sci. 49, 71–77 (2017).Article 

Google Scholar 
Gourlet-Fleury, S. et al. Tropical forest recovery from logging: a 24 year silvicultural experiment from Central Africa. Phil. Trans. R. Soc. B 368, 20120302 (2013).Article 

Google Scholar 
Obiang, N. L. E. et al. Spatial pattern of central African rainforests can be predicted from average tree size. Oikos 119, 1643–1653 (2010).Article 

Google Scholar 
Plotkin, J. B. et al. Predicting species diversity in tropical forests. P. Natl Acad. Sci. USA 97, 10850–10854 (2000).Article 

Google Scholar 
Graham, M. H. Confronting multicollinearity in ecological multiple regression. Ecology 84, 2809–2815 (2003).Article 

Google Scholar 
Tukey, J. W. Exploratory Data Analysis (Addison-Wesley,1977). More