sustainability

Article Accumulation of Urban Insect Pests in : 50 Years’ Observations on ( camphora)

Zhiyuan Xiang 1,2, Meifang Zhao 1,2,* and U. S. Ogbodo 2

1 Faculty of Life Science and Technology, Central South University of Forestry and Technology, 410004, China; [email protected] 2 National Engineering Laboratory for Applied Forest Ecological Technology in Southern China, Changsha 410004, China; [email protected] * Correspondence: [email protected]; Tel.: +86-0731-8562-3868

 Received: 2 December 2019; Accepted: 19 February 2020; Published: 20 February 2020 

Abstract: Since China experienced a rapid and unprecedented process of urbanization and climate change from 1978 onwards, pest outbreaks were frequently reported on urban forests, which reflects a significant imbalance between natural regulation and human control. Based on information extracted from all journal articles and reports about insect pests on camphor tree (Cinnamomum camphora) in urban China, we characterized historical patterns and trends in pest outbreaks over large areas. Our results suggested that (1) most distribution areas of C. camphora in urban China had pest records (14 provinces) over the last 50 years, especially at the south-eastern coastal areas; (2) pests on camphor tree in urban China showed an accelerated growth since the 1990s; and (3) pests on camphor tree in urban China were characterized by native and leaf-feeding . Urbanization seems to positively correlate with urban pest outbreaks. Changes of urban pest outbreaks could largely be described by synchronic changes of socio-economic indicators, of which CO2 emissions as metric tons per capita is the most significant predictor, followed by GDP and human population. Thus, managers and city planners should allocate resources to socio-economic-related pest outbreaks for a sustainable ecosystem.

Keywords: Cinnamomum camphora; insect pest; outbreak; urban tree; climate change

1. Introduction The rapid change in pest population densities in forests, known as forest pest outbreaks, and their dynamics represent critical systematic evidence for ecosystem stability [1]. For example, pest outbreaks occur when are less phenotypically plastic in their ability to cope with conditional changes compared to insect pests [2]. In addition, pest outbreaks can reflect increased urbanization, uncertainties in socio-economic requirements, climatic change impacts, and insect pest sensitivities. Numerous studies have focused on pest outbreaks at community or local scales under climate impacts [3–5], which can contribute to forest dynamics modeling, theoretical validating, and calibrating [6]. Results can also be valuable reference information for informed prediction of potential future changes under global climate change, development of adequate management strategies [6], and characterization of life in changing environments, such as biodiversity [7] and adaptive evolution mechanisms [8]. Urban forests are essential green infrastructure, which are recognized for generating a range of benefits and providing valuable ecosystem services (e.g., controlling pollution). They are highly dynamic, inherently vulnerable, and commonly subject to high levels of stress and disturbances [9–11]. As a result, managing urban forests for sustainable urban development is becoming increasingly complex

Sustainability 2020, 12, 1582; doi:10.3390/su12041582 www.mdpi.com/journal/sustainability Sustainability 2020, 12, 1582 2 of 15 under accelerated urbanization, owing to large uncertainties in socio-economic requirements [12], climatic change impacts (higher temperature, higher CO2 concentrations, increased nitrogen deposition rates, and longer growing seasons), and insect pests’ sensitivities (e.g., natural enemies, competitors and vectors of taxa, feeding types, invasive intensities, etc.). In this context of rapid and complex change, baseline characterization of historical patterns and trends in pest outbreaks over large areas builds a necessary foundation for urban planning and urban forestry, public health risk identification, and integrated pest management, as well as for global change vulnerability monitoring, modelling, and mapping [6,13]. China has experienced a rapid and unprecedented process of urbanization since 1978 (see Appendix 1, Figure S1, Table S1 in Supplementary Material for maps and detailed descriptions of the tempo-spatial distribution of urbanization in China), which entails increasing urban tree cultivation [14]. Cinnamomum camphora, also known as camphor tree or camphor laurel, is a medium to large broadleaved fast-growing tree species with a dense crown. This species is a typical indigenous street and park tree in China and also a famous non-indigenous tree in the rest of the world that is widely cultivated for timber production, carbon sequestration, and provision of critical ecosystem services and economic benefits in urban areas (Figure1b, see Figure 1b, see Appendix 3, Table S3 in Supplementary Material for a detailed view of its global distribution). It is universally native to East Asia and Pacific West Bank, especially from coastal areas to the heartland in China (20 provinces and municipalities) and adjacent islands, including and Hainan (Figure1c, see Figure 1c, see Appendix 3, Table S3 in Supplementary Material for a detailed view of its distribution in China). The area for cultivating camphor trees in Shanghai, the global financial center and transport hub, for example, accounts for 23.39% of all street/park trees according to the National Forestry Database from the China Forestry Science Data Center (http://www.cfsdc.org/). In Changsha, the core city of the Changsha–Zhuzhou–Xiangtan urban agglomeration, camphor trees comprise more than half of the total street/park tree populations. However, there have been frequent reports of pest outbreaks of C. camphora in China (see Appendix 2, Figure S2, Table S2 in Supplementary Material for maps and detailed descriptions of C. camphora-hosted insect pest damages). For example, 91.1% of C. camphora had experienced termite (e.g., Odontotermes formosanus) invasion in Wuhan, one of the central cities in China [15]. They were also seriously damaged by a newly discovered weevil pest in Shanghai [16]. Because environmental factors and insect outbreaks are distributed unevenly or randomly across a heterogeneous urban forest landscape. Detection of insect pest outbreaks is mostly made based on annual aerial surveys [17]. Since aerial overview survey coverage is limited in China (Figure1a), many forests that are infected by insect pests, especially those in urban landscapes, unfortunately go uncounted. There are also strong (though poorly studied) correlations among urban forest pest outbreak characteristics and socio-economic indicators. Which factors are responsible for the observed patterns of long-term insect pest outbreaks, and what possible drivers could be used to explain the magnitude of impact of forest insect pests, are also often difficult to determine. The aim of this study is to explore ecosystem stability of urban forests based on the characterization of historic patterns and trends in pest outbreaks. The objectives were to (1) determine the diversity of insect pests on C. camphora and (2) correlate and explain the urban forest pest outbreaks with indirect factors. Sustainability 2020, 12, 1582 3 of 15

Sustainability 2020, 12, x FOR PEER REVIEW 3 of 16

Figure 1. GeographicalFigure 1. Geographical distribution distribution of of ((a) studies studies and andinsect insecttypes according types to according Senf et al. [17], to (b Senf) C. et al. [17], (b) C. camphoracamphorathroughout throughout the the world, (c) ( C.c) camphoraC. camphora in China,in and China, (d) reported and ( doutbreaks) reported of urban outbreaks insect of urban pests of C. camphora in China over past half-century (see Appendix 2, Figure S2, Table S2 in insect pests of C. camphora in China over past half-century (see Appendix 2, Figure S2, Table S2 Supplementary Material for photos and detailed descriptions of C. camphora-hosted insect pest in Supplementarydamages). Material Panel (a) represents for photos patterns and of the detailed current state descriptions of forest insect of disturbanceC. camphora reviewed-hosted from insect pest damages). Panelthe literature (a) represents on remot patternse sensing. of Different the current insect types: state barkof forest beetle insect (olive triangles),disturbance broadleaved reviewed from the literature ondefoliator remote sensing. (green triangles), Different coniferous insect defoliator types: (orange beetle triangles) (olive and triangles), other fluid broadleaved feeders (red defoliator triangles). Panel (b) shows C. camphora is a widely distributed urban tree species in the world (see (green triangles),Appendix coniferous 2, Table S2 defoliator in Supplementary (orange Material), triangles) i.e., the most and popular other fluidnative feederstree species (red in tropical triangles). Panel (b) shows C.and camphora subtropicalis aChina widely (green distributed circles), and has urban been introduced tree species by other in the countries world (cyan (see circles). Appendix Photo 2, Table S2 in Supplementaryin the lower Material), right corner i.e., is the C. camphora most popular in the urban native street ( treeChangsha, species Hunan in tropicalprovince, China and subtropical). Panel China (d) indicates native C. camphora has experienced rapid growth of pest outbreaks (see Appendix 4, (green circles),Table and S4 in has Supplementary been introduced Material for by included other reports) countries, pink (cyan= prevalent circles). area, low Photo frequency in the; red lower right corner is C. camphora= severely damagedin the urban area, high street frequency. (Changsha, Hunan province, China). Panel (d) indicates native C. camphora has experienced rapid growth of pest outbreaks (see Appendix 4, Table S4 in Supplementary 2. Materials and Methods Material for included reports), pink = prevalent area, low frequency; red = severely damaged area, high frequency.2.1. The C. camphora-Hosted Insect Pest Dataset

2. Materials2.1.1. and Data Methods Collection We assembled a database of urban trees and insect pests coupled population time series in China 2.1. The C. camphora-Hostedfrom sources including Insect Google Pest Dataset scholar, Scopus, Web of Science, China Knowledge Resource Integrated Database (http://www.cnki.net), and bibliographic databases. Adhering to the PRISMA 2.1.1. Data Collectionguidelines [18], we conducted a systematic literature search using these data sources (initial database query 8 June 2017), for all published journal articles, newspaper reports, conference proceeding We assembledpapers, and a dissertations database that of were urban on insect trees pest and disturbance insect of pests C. camphora coupled trees in population urban areas with time series in China from sourcesno limitation including on publication Google dates. scholar,We identified Scopus, specific Webkeyword of Science,strings to search China relevant Knowledge studies Resource or reports, including “Chinese/China”, “urban/city”, “pest/insect”, “outbreak”, “camphor Integrated Database (http://www.cnki.net), and bibliographic databases. Adhering to the PRISMA guidelines [18], we conducted a systematic literature search using these data sources (initial database query 8 June 2017), for all published journal articles, newspaper reports, conference proceeding papers, and dissertations that were on insect pest disturbance of C. camphora trees in urban areas with no limitation on publication dates. We identified specific keyword strings to search relevant studies or reports, including “Chinese/China”, “urban/city”, “pest/insect”, “outbreak”, “camphor tree/camphor laurel/C. camphora”, “leaf/stem/root”, “habit/life history”, “record”, and “damage/injury/death/wilt/defoliate/tilt/broken/uproot/overthrown”. We also scanned the reference section of the studies or reports found in this research for other relevant articles. All titles and abstracts were screened and assessed for relevance using the following entrance criteria: (1) reported at least one insect/pest species that damage at least one organ/tissue of leaves, stem/bark/trunk, seeds/flowers, Sustainability 2020, 12, 1582 4 of 15 and/or roots; (2) acquired data from field investigations or reports; (3) had similar methods and specific information on site, tree, and insect/pest properties; (4) is a complete article from which all data and results were made available for analysis; and (5) would be possible to calculate the synthetic statistics and address meta-analytic questions. In addition, articles with the following characteristics were excluded from our analysis: (1) repeated publications or duplicated reports; (2) studies without clear pest records (e.g., physiological studies, morphological studies, or field experiment); (3) studies with poor data quality (data could not be analyzed) and/or quantity (incomplete data), such as studies with population records of insect/pests but without the corresponding records about the host trees; (4) studies with low information availability, from which the diversity of insect pests on C. camphora and tempo-spatial information for pest outbreaks could not be calculated for statistical analysis (attempts have been made to contact the corresponding author(s) to obtain the missing data whenever possible); (5) reports on other tree species of Cinnamomum, such as Cinnamomum porrectum (Roxb.) Kosterm., Cinnamomum subavenium Miq., and Cinnamomum micranthum (Hayata) were excluded; (6) studies in which any specific compounds (or management methods) had been applied to facilitate or inhibit pest outbreaks; and (7) articles that were just summaries of previous works. The literature search resulted in a total of 72 unique publications that met the strict criteria for our data compilation (see Appendix 4, Table S4 in Supplementary Material for 48 journal articles and 24 newspaper reports).

2.1.2. Data Filtering and Compilation Data from all qualifying studies were compiled into datasets. The following information was extracted from each article: (1) first author name, observation sites, geographic coordinates, sampling site description, recording time, and areas of sampled fields; (2) plant taxa and pest species; (3) tree growth characteristics including stand age, stand density, monoculture/mixture, tree height/size/architecture, etc.; (4) information description about the possible history of tree cultivate; (5) description about damaging plant organ/parts and symptoms; and (6) extent and duration of insect pest disturbances and method of detection. Various attributes determine the distribution and population density of insect pests [19]. Thus to facilitate the diversity and comparison analysis, we classified the complete extracted data sets into several major categories using the variable information of the respective publications. For extracted information, we synthesized pest information based on three categories: (1) tempo-spatial properties; (2) the taxonomic information and habitat characteristics; and (3) symptoms, causes, and types of feeders (see Appendix 5, Tables S5–S7 in Supplementary Material for detailed metadata descriptions and verification criteria). (1) Tempo-spatial properties Temporal properties included period, year, and month. Spatial properties included province, city/county, and location. One pest occurrence in our list referred to one pest event with a specific location and date [19]. (2) The taxonomic information and habitat characteristics Insect pests were classified into species, genus, family and order. We rectified the Chinese name according to the Chinese Animal Scientific Database (http://www.zoology.csdb.cn/) and the National digital-museum of animal specimens (http://museum.ioz.ac.cn/). We also recorded Latin names and taxonomic information in these national databases. We excluded pest records whose Chinese names were unspecific or cannot be found. We then checked the Latin name by related global databases (e.g., http://www.biolib.cz/en/main/). We applied Pivot Tables in Excel to sort the records and find mismatch taxonomic information based on the similarities in Latin names (see Appendix 6, Table S8 in Supplementary Material for the rectified results from the original records). Insect pests were classified into native species or non-native species according to locations where the species was first identified. Sustainability 2020, 12, 1582 5 of 15

(3) Types of feeders, symptoms, and causes. We group pest species into three types of feeders according to Cummins (1973)’s trophic functional groups [20]: leaf feeders, stem-feeders, and other-feeders, which differed significantly in structure and function [20,21] (see Appendix 5, Table S7 in Supplementary Material for our cluster). Leaf feeders, also known as foliage feeders or defoliators, are defoliating insects or leaf-eating insects feeding on the foliage of trees, including guilds such as leaf rollers, leaf miners, leaf strip feeder, leaf pit feeders, and leaf sap feeders, skeletonizers, and whole leaf feeders. Stem feeders are stem-eating insects that feed on the stem of trees, including bark beetles, borers, sap-sucking insects, and meristem feeders. Other feeders are insects that feed on the , seed, or root of trees and those feeding on mixed plant organs. For symptoms, we corrected the original description of infestation based on the handbook of the major forest pests in southeast Europe from the Food and Agriculture Organization (http: //www.fao.org/3/a-i4084e.pdf) and analytical approaches for diagnosing plant problems from the US Forest Service and CAL Fire (http://caforestpestcouncil.org/wp-content/uploads/2008/06/Insect-and- Disease-Training-Manual.pdf). Then we categorized symptoms as follows: leaves (12 symptoms), stems (8 symptoms), seeds (5 symptoms), seedlings and cuttings (4 symptoms), flowers (4 symptoms), fruit (6 symptoms), and root system (5 symptoms) (see Appendix 5, Tables S6 and S7 in Supplementary Material for recognizing infestation and damaging symptoms and signs). For causes of pest damage, we verified records according to insect pest analysis from the United States Department of Agriculture (https://www.fs.fed.us/rm/pubs/rmrs_gtr204/rmrs_gtr204_045_048. pdf) and Food and Agriculture Organization (http://www.fao.org/docrep/006/Y4837E/y4837e06.htm), and we then classified causes into four categories (http://www.eagri.org/eagri50/ENTO232/lec11.pdf). In addition, we recorded mean annual temperature, mean annual precipitation, and annual relative humidity of pest-infested areas according to the China Meteorological Administration (http: //data.cma.cn/). We recorded the range of temperature and precipitation, or humidity if the climate information for pest occurrence was vague. We also recorded tree age and categorized the host trees into young trees, mature trees, and old trees according to the United States Department of Agriculture (https://www.fs.fed.us/psw/publications/bingham/bingham3.pdf).

2.2. Pest cumulative Characteristics Natural and anthropogenic factors have resulted in altered environmental conditions that influence changes in abundance and diversity of insect pests [22], which then possibly affect the intensity and frequency of pest outbreaks on host [23]. The C. camphora-hosted insect pest outbreak dataset that was compiled are multi-proxy datasets of pest complexes on camphor trees. In order to explain the dataset records to understand the long-term patterns of community change as well as the effects of the diversity of insect pests on the stability of forest functioning, we used cumulative amount and frequency of insect pest detection over time as measures of intensity and frequency of pest outbreaks on camphor trees as used in Senf et al. [17] and Aukema et al. [19]. The cumulative amount and frequency of insect pest detection are usually analyzed by computing total number of detected establishments of insect pests for certain periods or a specific year on the basis of absolute numbers of detections involved in an outbreak across the study sites at any given time. We used linear and non-linear regressions to explore the cumulative detections of insect pest outbreaks and frequency of insect pest detection over time (cumulative frequency or yearly detection) in different categories (through 20 determinative variables), including total occurrence, taxonomic ranks (species, genera, family, and order), habitat characteristics (native and non-native), and feeding preferences (leaf, stem, and others): (1) We evaluated the distribution patterns of pest outbreak in ArcGIS v 9.3 (ESRI, Redlands, CA, USA). The charts were plotted in software OriginPro v 8.0 (OriginLab Co., Northampton, MA, USA); and (2) in order to assess the long-term patterns of pest complex on camphor trees, we regressed pest outbreak frequency (as response variables) against all calendar years (nonlinear terms for cumulative frequency and linear for frequency). We tested for correlations and covariances between the calendar year and Sustainability 2020, 12, 1582 6 of 15 variance in intensity and frequency of pest outbreaks on camphor trees using Pearson’s correlation coefficients (r). We calculated p-values (t-tests) to determine the significance (at a level of 0.05) of these relationships.

2.3. Socio-Economic Indicators The relationship between pest outbreak and calendar year can be complex and synthetic, including indirect drivers [9], such as higher urbanization intensities, unreasonable urban planning, and management, which results in a harsh growing environment and makes it difficult for trees to survive and reach maturity. Therefore, to examine the possible relationship between socioeconomic characteristics and emerging urban pest outbreaks, we regressed a synchronization term between random economic (or population-driven) processes against calendar years. We used degree of urbanization, fossil fuels’ share of energy consumption, composition of gross domestic product (GDP), changes of nitrous oxide (N2O) emissions, and carbon dioxide (CO2) emissions as proxies for socio-economic indicators. We extracted 22 indicators with clearly increasing or decreasing trends according to economic or population-driven topics in the World Development Indicators (https://www.eui.eu/Research/Library/ResearchGuides/Economics/Statistics/DataPortal/WDI), which were (1) three indicators for degree of urbanization (percentage of urban population of total population, percentage of population living in metropolitan areas that had a population of more than one million people, and the percentage of urban population living in the largest metropolitan area—Chongqing); (2) three indicators for fossil fuels’ share of energy consumption (percentage of fossil fuels energy consumption out of total energy consumption, fossil fuels’ share of total net electricity production, and coal share of total electricity production); (3) eight indicators for GDP (contribution of forest resources, services, industry, agriculture, trade, imports, exports, and external balance to overall GDP); (4) one indicator for changes of nitrous oxide (N2O) emissions (the percentage change in N2O from 1990 levels); and (5) seven indicators for CO2 emissions (CO2 emissions as metric tons per capita and from solid fuel, manufacturing industries and construction, residential buildings and commercial and public services, electricity and heat production, liquid fuel, and gaseous fuel consumption). Pearson’s correlation coefficients between pest outbreaks and the composite indices of socio-economic factors were calculated. Multiple linear regression analysis was carried out using all the composite socio-economic indices as independent variables, and pest outbreaks during 1964–2017 as dependent variables. A step-wise procedure was used in the decision to include significant independent variables. We here focused on overall accuracy as a measure of the coefficient of determination (R2) to estimate model performances (IBM SPSS Statistics v 21.0, IBM Corp., Armonk, NY, USA). We chose regressive slopes and descriptive variable domains to test for the presence of possible indirect causes.

3. Results

3.1. Insect Pest Spectrum on C. camphora Most distribution areas of C. camphora in urban China had pest records (14 provinces) over the last 50 years (Figure1d), resulting in a total of 772 pest occurrences (7 orders, 58 families, 225 genus, 233 species) (see Appendix 6, Table S8 in Supplementary Material for all recorded insect pest lists). Aphids (314 occurrences, 86 species) dominate these ecosystems, followed by moths (255 occurrences, 74 species) and beetles (116 occurrences, 58 species). Insect pests on C. camphora in urban China were characterized by native leaf feeders. Of the 772 pest records, 73% (n = 563) were native species (see Appendix 8, Figure S4, Table S10 in Supplementary Material for all recorded insect pests categorized by native vs. non-native). Leaf-feeding insects account for 73% of all species, followed by stem-feeding insects (105 species, 450 occurrences) and others-feeding insects (19 species, 81 occurrences) (see Appendix 9, Figure S5, Table S11 in Supplementary Material for feeding compositions). Sustainability 2020, 12, 1582 7 of 15 Sustainability 2020, 12, x FOR PEER REVIEW 7 of 16

Insect pests on C. camphora in urban China were characterized by native leaf feeders. Of the 772 3.2. Spatial Aggregationpest records, 73% (n = 563) were native species (see Appendix 8, Figure S4, Table S10 in Supplementary Material for all recorded insect pests categorized by native vs. non-native). Leaf- Pests werefeeding found insects in account Jiangxi for Province 73% of all species, for the followed first time, by stem which-feeding is located insects (105 in species, the southern 450 band of the middle andoccurrences) lower and others- River.feeding insects Severely (19 species, damaged 81 occurrences) areas were (see Appendix most numerous 9, Figure S5 in, Table coastal regions, which wereS11 characterized in Supplementary by M aaterial higher for feeding intensity compositions). and frequency of pest outbreaks (see Appendix 7,

Figure S3, Table3.2. Spatial S9 in Aggregation Supplementary Material for spatial aggregations).

3.3. Temporal TrendsPests were found in Jiangxi Province for the first time, which is located in the southern band of the middle and lower Yangtze River. Severely damaged areas were most numerous in coastal regions, Figure2which demonstrated were characterized that bothby a higher intensity intensity and and frequency frequency of of pest pest outbreaks outbreaks (see Appendix on C. camphora7, were Figure S3, Table S9 in Supplementary Material for spatial aggregations). monotonically increased since the 1960s. The cumulative frequency was described well in a power function, while3.3. Temporal frequency Trends of insect pest detection over time was fitted in a linear model (Figure2 and Appendix 10, TableFigure S12 2 demo inn Supplementarystrated that both intensity Material), and frequency with higherof pest outbreaks cumulative on C. camphora frequency were observed in native leaf-feedingmonotonically species. increased Specifically, since the 1960s frequency. The cumulati ofve pest frequency outbreaks was described on C. camphorawell in a powerin urban China function, while frequency of insect pest detection over time was fitted in a linear model (Figure 2 and was observedAppendix at a yearly 10, Table increase S12 in Supplementary rate of approximately Material), with 12higher occurrence, cumulative 8frequency species, observed 7 genera, in 4 families, and 1 ordernative (Figure leaf2-feedinga, that species. is a, ySpecifically,= ax + b,frequency a=0.766, of pest 0.694, outbreaks 0.438, on 0.075,C. camphora separately), in urban China while 9 native species, 6 non-nativewas observed species at a yearly (Figure increase2 rateb, thatof approximately is a, y = ax12 occurrence,+ b, a = 1.231,8 species, 0.851, 7 genera, 0.562, 4 families, separately), or 9 and 1 order (Figure 2a, that is a, y = ax + b, a=0.766, 0.694, 0.438, 0.075, separately), while 9 native leaf feeders,species, 2 stem 6 non feeders,-native species and 1 ( otherFigure 2b, feeder that is (Figurea, y = ax 2+ c,b, a that = 1.231, isa, 0.851, y = 0.562ax, +separately)b, a = 0.872,, or 9 0.19, 0.079 separately) perleaf feeders, decade. 2 stem feeders, and 1 other feeder (Figure 2c, that is a, y = ax + b, a = 0.872, 0.19, 0.079 separately) per decade.

Figure 2. Temporal patterns of cumulative frequency and yearly pest outbreaks (1964–2017) as classified Figure 2. Temporal patterns of cumulative frequency and yearly pest outbreaks (1964–2017) as by four traditionalclassified taxonomic by four traditional ranks taxonomic (species, ranks genera, (species, family, genera, andfamily, order) and order) (a), ( pooleda), pooled as as whole occurrence or binned byoccurrence two types or binnedof habitat by two characteristics types of habitat (nativecharacteristics and non-native)(native and non (b-native)), or three (b), or groups three of feeding preferencesgroups (leaf, stem,of feeding and pre others)ferences ((leaf,c). Statistical stem, and others) parameters (c). Statistical (r, Slope, parameters and (r, P) Slope, indicate and P) how general indicate how general temporal models fit the data (see Appendix 10, Table S12 in Supplementary temporal modelsMaterial fit for the non data-linear(see simulation Appendix and linear 10, regression Table S12 results). in Supplementary Cumulative power Materialmodels (bold for non-linear simulation andsolid linearlines in left regression panels): occurrence results). (black), Cumulative species (olive), power genera models (dark green), (bold family solid (green), lines order in left panels): occurrence (black),(light green), species native (olive), (blue), non genera-native (dark(cyan), green),leaf feeders family (red), stem (green), feeders order (orange), (light other green), feedersnative (blue), non-native (cyan), leaf feeders (red), stem feeders (orange), other feeders (light magenta). Linear growth models (fitted red lines with a shaded area of 95% confident intervals in right panels): gray circles represent the frequency of detection for each recorded year. Unrepresentative data from early reports or province-wide surveys were excluded (open triangles). Pearson correlation was used to evaluate the correlation coefficient between variables.

3.4. Causes and Controller Attribution We chose temporal trends (Figure3), synchronic covariant coe fficients (Figure4), and multi-regressive model simulations (Table1) to test for the presence of possible indirect causes for pest outbreaks. Sustainability 2020, 12, 1582 8 of 15 Sustainability 2020, 12, x FOR PEER REVIEW 9 of 16

Figure 3. Historical trends of socio-economic indicators. (A) Degree of urbanization: (a) urban Figure 3. Historical trends of socio-economic indicators. (A) Degree of urbanization: (a) urban population (%population of total (% population), of total population (b),) ( populationb) population in in urban urban agglomerations agglomerations of more than of more one million than one million (% of total population),(% of total population), and (c and) population (c) population in in thethe largest largest city city (% of (% urban ofurban population) population).. (B) Energy (B) Energy consumptionconsumption of fossil of fuels: fossil fuels (a): fossil(a) fossil fuel fuel energy energy consumption consumption (% of total (% energy of total consumption), energy ( consumption),b) electricity production from fossil fuels (% of total electricity production), and (c) electricity production b c ( ) electricityfrom production coal sources from (% of fossil total electricity fuels (% production). of total electricity (C) Gross production),domestic product and components ( ) electricity (% of production from coal sourcesGDP) as (% derived of total from electricity (a) forest rents, production). (b) services ( valueC) Gross-added, domestic (c) industry product value-added, components (d) (% of GDP) as derivedagriculture from value (a) forest-added, rents,(e) trade, (b ()f) servicesimports of value-added,goods and services, (c )(g industry) exports of value-added,goods and services, (d ) agriculture and (h) external balance on goods and services. (D) Nitrous oxide emissions (% change from 1990). value-added, (e) trade, (f) imports of goods and services, (g) exports of goods and services, and (E) CO2 emissions as (a) metric tons per capita, from (b) solid fuel consumption (% of total CO2 (h) external balanceemissions), on (c) goodsmanufacturing and services. industries (andD) construction Nitrous oxide (% of total emissions fuel combustion), (% change (d) res fromidential 1990). (E) CO2 emissions asbuildings (a) metric and commercial tons per and capita, public from services (b )(% solid of total fuel fuel consumption combustion), (e) (%electricity of total and COheat2 emissions), (c) manufacturingproduction, industries total (% of and total construction fuel combustion), (% (f of) liquid total fuel fuel consumption combustion), (% of total (d) CO residential2 emissions), buildings and and (g) gaseous fuel consumption (% of total CO2 emissions). commercial and public services (% of total fuel combustion), (e) electricity and heat production, total

(% of total fuel combustion), (f) liquid fuel consumption (% of total CO2 emissions), and (g) gaseous

fuel consumption (% ofSustainability total 20 CO20, 12,2 x FORemissions). PEER REVIEW 10 of 16

Figure 4. Heatmap illustratingFigure 4. Heatmap Pearson’s illustrating Pearson’s correlation correlation coefficients coeffi betweencients socio-economic between indicators socio-economic indicators and pest complex on camphor trees. Blue and red levels correspond to strengths of the correlations, and pest complex on camphorstrong positive trees. (r = 0.6~1.0), Blue moderate and positive red (r = 0.3~0.6), levels weak correspondpositive (r = 0.0~0.3), weak to negative strengths of the correlations, (r = −0.3~0.0), moderate negative (r = −0.6 ~−0.3), and strong negative (r = −1.0 ~−0.6). Statistical strong positive (r = 0.6~1.0),significance moderate is indicated by positive ● p < 0.001, ●(r p < 0.01= ,0.3~0.6), 〇 p < 0.05. weak positive (r = 0.0~0.3), weak negative (r = 0.3~0.0), moderate negative (r = 0.6 ~ 0.3), and strong negative (r = 1.0 ~ 0.6). Statistical − The linear relationship between− the− random economic or population-driven indicators and− − significance is indicatedcalendar by yearsp < was0.001, highly significantp < 0.01, (all p < 0.001p expect< 0.05. p for CO2 emissions from solid fuel consumption• and CO2 emissions• from liquid fuel consumption, which was 0.01 and 0.03 respectively). Most indicators were more significantly positively# related to the calendar year than insect pest detection over the same time (r for 13 indicators > 0.6). Regression slope of pest outbreaks and socio- economic indicators over time indicated that there was comparatively higher convergence between pest outbreaks and some socio-economic indicators (r for 14 indicators > 0.0 and p < 0.05). The correlation analysis indicated that pest outbreak dynamics response to the urban population (% of total population), population in urban agglomeration of more than one million (% of total population), fossil fuel energy consumption (% of total energy consumption), GDP compositions (including forest rents, services, agriculture, trade, and imports and exports of goods and services), changes of N2O emissions, and CO2 emissions (Figure 4). Overall, cumulative amount of pest outbreaks showed a greater response to socio-economic indicators than yearly pest outbreaks. The outbreak of leaf feeder was significantly correlated with socio-economic indicators (p < 0.05), although other types of feeders were nearly unaffected. Multi-regression analysis (Table 1) indicates that most of historical trends of the pest outbreak metrics could be explained by socio-economic indicators (p for 17/20 models < 0.05). Changes of cumulative amount of pest outbreaks could largely be described by synchronic changes of socio- economic indicators (R2 for 10/10 models > 0.93), of which CO2 emissions as metric tons per capita is Sustainability 2020, 12, 1582 9 of 15

Table 1. The multiple linear models for pest outbreaks as regressed by socio-economic predictors.

Parameter * Yi R2 F Value p bij ai X2 X6 X8 X10 X13 X14 X16 Y 65.21 10.21 - - - - - 4.27 1.49 39.11 1.87 0.969 239.22 <0.001 1 ± ± ± Y 43.18 9.55 - - - - - 3.82 1.39 36.85 1.75 0.970 241.42 <0.001 2 ± ± ± Y 2.75 6.33 ------28.1 1.32 0.968 450.6 <0.001 3 ± ± Y 45.62 6.87 1.65 0.49 - 1.8 0.32 - - - 8.96 0.5 0.995 1067.36 <0.001 4 − ± − ± ± ± Y 19.02 12.80 - - - - 1.19 0.48 - - 0.260 6.27 <0.05 5 − ± ± Y 16.75 11.62 - - - - 1.08 0.43 - - 0.260 6.26 <0.05 6 − ± ± Y 10.68 7.07 - - - - 0.75 0.26 - - 0.321 8.08 <0.05 7 − ± ± Y 6.37 1.04 - - - 0.26 0.08 - - - 0.412 11.51 <0.01 8 ± − ± Y 53.09 16.25 - - - - - 6.02 2.37 63.27 2.98 0.970 244.62 <0.001 9 ± ± ± Y 110.04 55.68 - - 4.28 1.6 - - 5.19 1.54 39.26 3.18 0.977 209.21 <0.001 10 − ± ± ± ± Y 20.21 4.60 ------17.48 0.96 0.956 329.62 <0.001 11 ± ± Y 25.24 20.37 - - - - 1.74 0.76 - - 0.220 5.24 <0.05 12 − ± ± Y13 ------>0.05 Y 234.00 102.71 - 3.08 1.3 - - - - - 0.337 5.57 <0.05 14 − ± ± Y 13.82 10.42 - - - - - 3.41 1.52 39.68 1.91 0.969 232.28 <0.001 15 ± ± ± Y 31.09 24.95 - - 2.02 0.72 - - 2.5 0.69 18.23 1.42 0.978 226.31 <0.001 16 − ± ± ± ± Y 2.04 0.73 ------2.23 0.15 0.933 211.39 <0.001 17 ± ± Y 17.23 11.77 - - - - 1.15 0.44 - - 0.280 6.834 <0.05 18 − ± ± Y19 ------>0.05 Y20 ------>0.05 Model: (Y ) = (a ) + (Σ(b X )), i = 1~20 and j = 1~22. (Y ) donate the 20 outcome variables of the regression equation, i.e., the observed cumulative and yearly pest outbreaks as indicated i i ij× j i in Figure2:Y 1 ~ Y8 = yearly cumulative amount (Y1~Y4) and yearly(Y5~Y8) frequency of pest outbreaks by taxonomic species, genera, family, and order, respectively, Figure2a; Y 9~Y14 = yearly cumulative amount (Y9~Y11) and yearly (Y12~Y14) frequency of pest outbreaks pooled as whole occurrence or binned by two types of habitat characteristics (native and non-native), respectively, Figure2b; Y 15~Y20 = yearly cumulative amount (Y15~Y17) and yearly (Y18~Y20) frequency of pest outbreaks binned by three groups of feeding preferences (leaf, stem, and others), respectively, Figure2c. While the (X j) refer the included (Step-wise method, p < 0.05) seven predictors in the model selected from the 22 socio-economic indicators shown in Figure3:X 2 = population in urban agglomeration of more than one million (% of total population), Figure3Ab; X 6 = electricity production from coal sources (% of total electricity production), Figure3Bc; X 8,X10,X13,X14 = percentage of total GDP as derived from services value added (X8, Figure3Cb), agriculture value-added (X 10, Figure3Cd), exports of goods and services (X13, Figure3Cg), and external balance on goods and services (X 14, Figure3Ch), respectively; X 16 = CO2 emissions as metric tons per capita (Figure3Da). Sustainability 2020, 12, 1582 10 of 15

The linear relationship between the random economic or population-driven indicators and calendar years was highly significant (all p < 0.001 expect p for CO2 emissions from solid fuel consumption and CO2 emissions from liquid fuel consumption, which was 0.01 and 0.03 respectively). Most indicators were more significantly positively related to the calendar year than insect pest detection over the same time (r for 13 indicators > 0.6). Regression slope of pest outbreaks and socio-economic indicators over time indicated that there was comparatively higher convergence between pest outbreaks and some socio-economic indicators (r for 14 indicators > 0.0 and p < 0.05). The correlation analysis indicated that pest outbreak dynamics response to the urban population (% of total population), population in urban agglomeration of more than one million (% of total population), fossil fuel energy consumption (% of total energy consumption), GDP compositions (including forest rents, services, agriculture, trade, and imports and exports of goods and services), changes of N2O emissions, and CO2 emissions (Figure4). Overall, cumulative amount of pest outbreaks showed a greater response to socio-economic indicators than yearly pest outbreaks. The outbreak of leaf feeder was significantly correlated with socio-economic indicators (p < 0.05), although other types of feeders were nearly unaffected. Multi-regression analysis (Table1) indicates that most of historical trends of the pest outbreak metrics could be explained by socio-economic indicators (p for 17/20 models < 0.05). Changes of cumulative amount of pest outbreaks could largely be described by synchronic changes of 2 socio-economic indicators (R for 10/10 models > 0.93), of which CO2 emissions as metric tons per capita is the most significant predictor, followed by the contribution of external balance on goods and services in total GDP. Synchronic changes of socio-economic indicators only explain 26%~41.2% of the changes for yearly pest outbreaks (R2 for 7/10 models = 0.26~0.412), and 70% of the models were only significantly simulated from the contribution of exports of goods and services in total GDP (p for 5/7 models < 0.05). The cumulative amount of outbreaks for pest species = 65.21 + 39.11 CO × 2 emissions as metric tons per capita + 4.27 contribution of external balance on goods and services in × total GDP (R2 = 0.969, p < 0.001).

4. Discussion

4.1. Improved Reference Data for Assessing Insect Pest Outbreak

4.1.1. Field Observational Database Except for large-scale monitoring like remote sensing, a high-quality database based on field observation in a large tempo-spatial scale is especially important to characterize the historical patterns and trends for pest accumulations. Existing databases provide limited information for pest species. For example, Senf et al. [17] reviewed the worldwide literature on the remote sensing of insect outbreaks (Figure1a), resulting in 27 insect species and attributes that include feeding preference, habitat characteristics, host species, and spatial properties; Aukema et al. [19] compiled a list for nonindigenous forest insects with taxonomic information and feeding preference information established in the United States between 1860 to 2006, which resulted in more than 450 insects; Liebhold et al. [24] assembled 62 species of damaging invasive forest insects established in the continental USA, with feeding preference, date of detection, and distribution information recorded. Besides, there is no database based on field studies concentrated in urban areas. We present the first dataset of insect pests on native C. camphora in urban China over the past half-century based on all relevant pest attributes, which included taxonomic ranks, types of feeders, habitat characteristics, symptoms, signs, damage, and causes (see Appendix 11, Table S14 in Supplementary Material for the dataset). Compiled pest attributes and pest records in a large tempo-spatial scale could help us to collect and integrate data temporally and spatially, to analyze the causes and progress of the pest accumulation and to predict future infestation patterns. Sustainability 2020, 12, 1582 11 of 15

4.1.2. Habitat Characteristics Both the native species and non-native species showed a decreasing trend from southeast coastal areas to northwest inland areas, which reflected the similarity in their need for environmental conditions. Studies had shown that water potential remained positively associated with pest outbreaks of both native species and non-native species in urban areas [25]. In general, native species have maintained higher fitness than non-native species on C. camphora in China over the past half-century, with a faster accumulation process than non-native species. Related studies also indicated there are more obstacles to the assemblage of non-native species, such as their association with invasion pathways and variation in pest life history [26]. In addition, a lot of studies provided evidence that warming is significantly related to native species in urban areas [27,28], while the environmental effects on non-native species need to be further studied. Studies had further revealed the complexity and specificity in the establishment processes for non-native species, such as the more significant impacts of connectivity of host plant patches, the number of international tourist arrivals, exotic vegetation cover, and regional GDP on than on native species [29,30]. However, phenotypic plasticity and resource availability can allow non-native pest species to succeed across heterogeneous environments, such as urban forest ecosystems. Since 2010, urban forests in China may experience more frequent fluctuations in allocation, growth, and reproduction.

4.1.3. Types of Feeders Regional pest outbreaks were characterized by co-occurrence patterns of insect pests in different types of feeders, especially the pest aggregation of leaf feeders, stem feeders, and other feeders, or an aggregation of leaf feeders and stem feeders. The pest accumulation and detection in all of the types of feeders also showed a decreasing trend from eastern coastal areas to western interior regions. Related studies proved there exists similarity among insects in different feeding-types in response to conditional changes. For instance, the positive impacts of temperature and tree water potential were observed both in leaf feeders and stem feeders in cities [31]. In addition, there was an association between urban land-use and species richness and abundance both in leaf feeders and stem feeders [32]. Leaf feeders played a more significant role in the overall insect pest outbreaks of C. camphora in China, followed by stem feeders, then other feeders. Studies had indicated that warming and CO2 had strong positive effects on foliar consumption for insect pests in urban areas [27,33]. Insects also responded differently to experimentally induced water deficit in plants classified by sap-feeding insects and chewing insects [34]. However, the response of stem feeders and other feeders to conditional changes in urban areas remains to be studied.

4.2. Accelerated Urban Pest Outbreaks in China Since the 1990s Urbanization over the last century has been accompanied by accelerated tree growth [35,36], increased pest fitness and abundance [37], and enhanced pest suitability and plasticity [38]. Related studies put emphasis on improper plantations when discussing the causes for pest damage (see Table1 for detailed causes). However, the random economic or population-driven process can explain the pest damage in some cases: Researchers have realized the acceleration of insect pest outbreaks. For example, an accumulation of non-native pest insects (8 species per 10 years) was detected on in [9]. An accumulation of nonindigenous forest insects (25 species per 10 years) was also observed around the United States [19]. In our research, pest outbreaks of C. camphora in urban China showed an increasing trend with a well-fitted power function over time. Especially, pest outbreaks of C. camphora in urban China accelerated since the 1990s. Research also indicated that global urbanization has increased rapidly with a nearly tripled number of mega-cities and a doubled population living in medium-sized cities since the 1990s [39]. Global emissions of all major greenhouse gases have increased considerably since 1990, such as net emissions of CO2 that increased by 42% according to the United Sustainability 2020, 12, 1582 12 of 15

States Environmental Protection Agency (https://www.epa.gov/climate-indicators/climate-change- indicators-global-greenhouse-gas-emissions). In addition, increasing drought under global warming since the 1980s has been observed [40]. However, there was a significant reduction in the pest abundance since 1990s in some studies [41,42] about croplands, forests, and grassland, which may also indicate that pest outbreaks may not always be explained by the direct impacts of environmental factors. Insect pests are widely distributed in growing regions of C. camphora in urban China. In particular, the coastal areas had been seriously disturbed with a higher cumulative number and frequency of pest detection, which agreed with the research that showed a highly concentrated pest distribution in the north-eastern region of USA [24]. It was also proved that the increase of pest records was most likely due to economic development and urban expansion since the 1970s [43]. The study further showed that industrialization and urbanization in China were spatially concentrated in eastern coastal cities [44]. On the other hand, relatively high temperatures and suitable humidity of eastern areas in China [45] would motivate pest growth and reproduction. Especially, extreme climatic events like droughts were observed at a significantly increasing trend since the 1960s [46,47]. Particularly, Jiangxi had the highest frequency and diversity of insect pests on C. camphora when compared to other provinces. The changeable climate in Jiangxi is characterized by frequent droughts [48] and increasing extreme precipitation events since the 1960s [49]. Especially in 1964, pests on C. camphora were first detected in Jiangxi, and extreme weather phenomena represented by drought and heavy rain occurred sequentially before the pest outbreak [50]. Similarly, there was pest detection in urban areas in Montes Claros, Brazil, with the tropical savanna climate [51]. In addition, the increased pest damage was also observed on heavily watered trees in urban areas [21]. It may highlight the sensibility, suitability, and plasticity of insect pests when with conditional changes [52].

5. Conclusions It’s essential to estimate the ecosystem stability for urban trees to manage urban sustainability. In this study, we mapped the tempo-spatial changes of vulnerability for C. camphora in urban China between 1964 and 2017 based on pest records from field observation. As a result, a total of 772 pest occurrences (7 orders, 58 families, 225 genera, and 233 species) were recorded, with native leaf feeders dominating the ecosystems in terms of species richness, abundance, and accumulation rate. Pest outbreaks were highly focused in the southeastern coastal areas, with the largest variation and highest frequency over the last 50 years. We also revealed the most likely indirect causes (random economic or population-driven process) for temporal patterns of intensity and frequency of pest outbreaks on camphor trees. Our results can provide detailed information for pest identification in the future and reference information for urban planning and urban forestry, integrative pest management, and data monitoring, as well as for forest dynamics modeling, theoretical validating, and calibrating.

Supplementary Materials: The following are available online at http://www.mdpi.com/2071-1050/12/4/1582/s1, Figure S1: Tempo-spatial distribution of urbanization in China, Figure S2: Photos showing for recorded pest outbreaks in native C. camphora among groups of leaf feeding (a–k), stem feeding (l–p) and others (q–r), Figure S3: Regional patterns and provincial distributions of total (a) occurrence, (b) species, (c) genera, (d) family, and (e) order for C. camphora hosted pest outbreaks over past half Century, Figure S4: Tempo-spatial distribution of pest outbreaks on the recording of overall (a) and provincial (b-o) insect pest disturbances by types of insect pest habitat characteristics, Figure S5: Tempo-spatial distribution of pest outbreaks on the recording of overall (a) and provincial (b–o) insect pest disturbances by types of feeding types, Table S1: Panel captions for Figure S1, Table S2: Panel captions for Figure S2, Table S3: Worldwide distribution of C. camphora shown in Figure1b, Table S4: Selected studies for data extraction, Table S5: Meta data for our pest database, Table S6: Recognizing common signs and symptoms, Table S7: Description of types of feeders clusters, Table S8: List of insect pests found on native urban tree species C. camphora in China in the period 1964–2017, Table S9: Summary of reported outbreaks (Frequency) of insect pests of C. camphora in native geographic range in Chinese urban area, Table S10: Summary of reported native vs. non-native insect pests of C. camphora in native geographic range in China among different province, Table S11: Summary of reported insect pests of C. camphora in native geographic range in China among different province by three feeding types, Table S12: Non-linear simulation and linear regression parameters for Figure2, Table S13: Linear regression parameters for Figure3, Table S14: the dataset. Sustainability 2020, 12, 1582 13 of 15

Author Contributions: M.Z. performed analyses, developed theory, and wrote the paper; Z.X. compiled data and wrote the paper; U.S.O. reviewed the paper. All authors have read and agreed to the published version of the manuscript. Funding: This study was financially supported by the Natural Science Foundation of China (Nos. 31971456 and 31600355) and the Program for Science and Technology Innovative Research Team in Higher Educational Institutions of Hunan Province. Acknowledgments: The authors would like to thank Xi Peng and Sheng Zhang for their suggestions and comments. We thank the editor and two anonymous reviewers for their careful reading of our manuscript and constructive comments. Conflicts of Interest: The authors declare no conflict of interest.

References

1. Bradshaw, C.J.A.; Leroy, B.; Bellard, C.; Roiz, D.; Albert, C.; Fournier, A.; Barbet-Massin, M.; Salles, J.-M.; Simard, F.; Courchamp, F. Massive yet grossly underestimated global costs of invasive insects. Nat. Commun. 2016, 7, 12986. [CrossRef][PubMed] 2. Kant, M.R.; Jonckheere, W.; Knegt, B.; Lemos, F.; Liu, J.; Schimmel, B.C.J.; Villarroel, C.A.; Ataide, L.M.S.; Dermauw, W.; Glas, J.J.; et al. Mechanisms and ecological consequences of plant defence induction and suppression in herbivore communities. Ann. Bot. 2015, 115, 1015–1051. [CrossRef][PubMed] 3. Boyd, I.L.; Freer-Smith, P.H.; Gilligan, C.A.; Godfray, H.C.J. The consequence of tree pests and diseases for ecosystem services. Science 2013, 342, 1235773. [CrossRef][PubMed] 4. Kurz, W.A.; Dymond, C.C.; Stinson, G.; Rampley, G.J.; Neilson, E.T.; Carroll, A.L.; Ebata, T.; Safranyik, L. Mountain beetle and forest carbon feedback to climate change. Nature 2008, 452, 987–990. [CrossRef] [PubMed] 5. Kozlov, M.V.; Stekolshchikov, A.V.; Söderman, G.; Labina, E.S.; Zverev, V.; Zvereva, E.L. Sap-feeding insects on forest trees along latitudinal gradients in northern Europe: A climate-driven patterns. Glob. Chang. Biol. 2015, 21, 106–116. [CrossRef][PubMed] 6. Kautz, M.; Meddens, A.J.H.; Hall, R.J.; Arneth, A. Biotic disturbances in Northern Hemisphere forests—A synthesis of recent data, uncertainties and implications for forest monitoring and modelling. Glob. Ecol. Biogeogr. 2017, 26, 533–552. [CrossRef] 7. Beudert, B.; Bässler, C.; Thorn, S.; Noss, R.; Schröder, B.; Dieffenbach-Fries, H.; Foullois, N.; Müller, J. Bark beetles increase biodiversity while maintaining drinking water quality: Bark beetles, biodiversity and drinking water. Conserv. Lett. 2015, 8, 272–281. [CrossRef] 8. Johnson, M.T.J.; Munshi-South, J. Evolution of life in urban environments. Science 2017, 358, eaam8327. [CrossRef] 9. Paine, T.D.; Millar, J.G.; Daane, K.M. Accumulation of pest insects on Eucalyptus in California: Random process or smoking gun. J. Econ. Entomol. 2010, 103, 1943–1949. [CrossRef] 10. Meineke, E.K.; Dunn, R.R.; Frank, S.D. Early pest development and loss of biological control are associated with urban warming. Biol. Lett. 2014, 10, 20140586. [CrossRef] 11. Kowarik, I.; Hiller, A.; Planchuelo, G.; Seitz, B.; von der Lippe, M.; Buchholz, S. Emerging urban forests: Opportunities for promoting the wild side of the urban green infrastructure. Sustainability 2019, 11, 6318. [CrossRef] 12. Van Grunsven, R.H.A.; Becker, J.; Peter, S.; Heller, S.; Hölker, F. Long-term comparison of attraction of flying insects to streetlights after the transition from traditional light sources to light-emitting diodes in urban and peri-urban settings. Sustainability 2019, 11, 6198. [CrossRef] 13. Seidl, R.; Rammer, W. Climate change amplifies the interactions between wind and bark beetle disturbances in forest landscapes. Landsc. Ecol. 2017, 32, 1485–1498. [CrossRef][PubMed] 14. Zhang, K.H.; Song, S. Rural–urban migration and urbanization in China: Evidence from time-series and cross-section analyses. China Econ. Rev. 2003, 14, 386–400. [CrossRef] 15. Appel, A.G.; Hu, X.P.; Zhou, J.; Qin, Z.; Zhu, H.; Chang, X.; Wang, Z.; Liu, X.; Liu, M. Observations of the biology and ecology of the black-winged termite, Odontotermes formosanus Shiraki (Termitidae: Isoptera), in Camphor, Cinnamomum camphora (L.) (). Psyche J. Entomol. 2012, 2012, 123102. 16. Huang, H.; Wu, S.; Huang, J.; Wu, S.; Gao, L.; Yang, L. Diagnosis and damage of weevil pest Pagiophloeus tsushimanus on camphor tree (in Chinese). J. Zhejiang A&F Univ. 2014, 31, 764–767. Sustainability 2020, 12, 1582 14 of 15

17. Senf, C.; Seidl, R.; Hostert, P. Remote sensing of forest insect disturbances: Current state and future directions. Int. J. Appl. Earth Obs. Geoinf. 2017, 60, 49–60. [CrossRef] 18. Liberati, A.; Altman, D.G.; Tetzlaff, J.; Mulrow, C.; Gøtzsche, P.C.; Ioannidis, J.P.A.; Clarke, M.; Devereaux, P.J.; Kleijnen, J.; Moher, D. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: Explanation and elaboration. PLoS Med. 2009, 6, e1000100. [CrossRef] 19. Aukema, J.E.; McCullough, D.G.; Von Holle, B.; Liebhold, A.M.; Britton, K.; Frankel, S.J. Historical accumulation of nonindigenous forest pests in the continental United States. BioScience 2010, 60, 886–897. [CrossRef] 20. Cummins, K.W. Trophic relations of aquatic insects. Annu. Rev. Entomol. 1973, 18, 183–206. [CrossRef] 21. Reich, P.B.; Tjoelker, M.G.; Pregitzer, K.S.; Wright, I.J.; Oleksyn, J.; Machado, J.L. Scaling of respiration to nitrogen in leaves, stems and roots of higher land plants. Ecol. Lett. 2008, 11, 793–801. [CrossRef] 22. Katsaruware, R.; Mafongoya, P.; Gubba, A. Responses of insect pests and plant diseases to changing and variable climate: A review. J. Agric. Sci. 2017, 9, 160–168. [CrossRef]

23. Sun, Y.; Yin, J.; Chen, F.; Wu, G.; Ge, F. How does atmospheric elevated CO2 affect crop pests and their natural enemies? Case histories from China. Insect Sci. 2011, 18, 393–400. [CrossRef] 24. Liebhold, A.M.; McCullough, D.G.; Blackburn, L.M.; Frankel, S.J.; Von Holle, B.; Aukema, J.E. A highly aggregated geographical distribution of forest pest invasions in the USA. Divers. Distrib. 2013, 19, 1208–1216. [CrossRef] 25. Hanks, L.M.; Denno, R.F. Natural enemies and plant water relations influence the distribution of an armored scale insect. Ecology 1993, 74, 1081–1091. [CrossRef] 26. Liebhold, A.M.; Yamanaka, T.; Roques, A.; Augustin, S.; Chown, S.L.; Brockerhoff, E.G.; Pyšek, P. Global compositional variation among native and non-native regional insect assemblages emphasizes the importance of pathways. Biol. Invasions 2016, 18, 893–905. [CrossRef] 27. Jamieson, M.A.; Schwartzberg, E.G.; Raffa, K.F.; Reich, P.B.; Lindroth, R.L. Experimental climate warming alters and phytochemistry and performance traits for an outbreak insect herbivore. Glob. Chang. Biol. 2015, 21, 2698–2710. [CrossRef] 28. Youngsteadt, E.; Dale, A.G.; Terando, A.J.; Dunn, R.R.; Frank, S.D. Do cities simulate climate change? A comparison of herbivore response to urban and global warming. Glob. Chang. Biol. 2015, 21, 97–105. [CrossRef] 29. Rand, T.A.; Louda, S.M. Invasive insect abundance varies across the biogeographic distribution of a native host plant. Ecol. Appl. 2006, 16, 877–890. [CrossRef] 30. Edney-Browne, E.; Brockerhoff, E.G.; Ward, D. Establishment patterns of non-native insects in New Zealand. Biol. Invasions 2018, 20, 1657–1669. [CrossRef] 31. Dale, A.G.; Frank, S.D. The effects of urban warming on herbivore abundance and street tree condition. PLoS ONE 2014, 9, e102996. [CrossRef] 32. McKinney, M.L. Effects of urbanization on species richness: A review of plants and animals. Urban Ecosyst. 2008, 11, 161–176. [CrossRef]

33. Sun, Y.; Guo, H.; Ge, F. Plant–aphid interactions under elevated CO2: Some cues from aphid feeding behavior. Front. Plant Sci. 2016, 7, 502. [CrossRef][PubMed] 34. Huberty, A.F.; Denno, R.F. Plant water stress and its consequences for herbivorous insects: A new synthesis. Ecology 2004, 85, 1383–1398. [CrossRef] 35. Zhao, S.; Liu, S.; Zhou, D. Prevalent vegetation growth enhancement in urban environment. Proc. Natl. Acad. Sci. USA 2016, 113, 6313–6318. [CrossRef][PubMed] 36. Pretzsch, H.; Biber, P.; Uhl, E.; Dahlhausen, J.; Schütze, G.; Perkins, D.; Rötzer, T.; Caldentey, J.; Koike, T.; van Con, T.; et al. Climate change accelerates growth of urban trees in metropolises worldwide. Sci. Rep. 2017, 7, 15403. [CrossRef] 37. Meineke, E.K.; Dunn, R.R.; Sexton, J.O.; Frank, S.D. Urban warming drives insect pest abundance on street trees. PLoS ONE 2013, 8, e59687. [CrossRef] 38. Merckx, T.; Van Dyck, H. Urbanization-driven homogenization is more pronounced and happens at wider spatial scales in nocturnal and mobile flying insects. Glob. Ecol. Biogeogr. 2019, 28, 1440–1455. [CrossRef] 39. United Nations, Population Division. World Urbanization Prospects: The 2014 Revision; United Nations, Department of Economic and Social Affairs: New York, NY, USA, 2015. Sustainability 2020, 12, 1582 15 of 15

40. Dai, A. Increasing drought under global warming in observations and models. Nat. Clim. Chang. 2013, 3, 52–58. [CrossRef] 41. Gordo, O.; Sanz, J.J. Temporal trends in phenology of the honey bee Apis mellifera (L.) and the small white Pieris rapae (L.) in the Iberian Peninsula (1952–2004). Ecol. Entomol. 2006, 31, 261–268. [CrossRef] 42. Brooks, D.R.; Bater, J.E.; Clark, S.J.; Monteith, D.T.; Andrews, C.; Corbett, S.J.; Beaumont, D.A.; Chapman, J.W. Large carabid beetle declines in a United Kingdom monitoring network increases evidence for a widespread loss in insect biodiversity. J. Appl. Ecol. 2012, 49, 1009–1019. [CrossRef] 43. Wu, J.; Hobbs, R. Key issues and research priorities in landscape ecology: An idiosyncratic synthesis. Landsc. Ecol. 2002, 17, 355–365. [CrossRef] 44. Bai, X.M.; Shi, P.J.; Liu, Y.S. Society: Realizing China’s urban dream. Nature 2014, 509, 158–160. [CrossRef] [PubMed] 45. Luo, M.; Lau, N. Increasing heat stress in urban areas of eastern China: Acceleration by urbanization. Geophys. Res. Lett. 2018, 45, 13060–13069. [CrossRef] 46. Piao, S.; Ciais, P.; Huang, Y.; Shen, Z.; Peng, S.; Li, J.; Zhou, L.; Liu, H.; Ma, Y.; Ding, Y.; et al. The impacts of climate change on water resources and agriculture in China. Nature 2010, 467, 43–51. [CrossRef][PubMed] 47. Zhai, J.; Huang, J.; Su, B.; Cao, L.; Wang, Y.; Jiang, T.; Fischer, T. Intensity–area–duration analysis of droughts in China 1960–2013. Clim. Dyn. 2017, 48, 151–168. [CrossRef] 48. Han, W.; Liang, C.; Jiang, B.; Ma, W.; Zhang, Y. Major natural disasters in China, 1985–2014: Occurrence and damages. Int. J. Environ. Res. Public Health 2016, 13, 1118. [CrossRef][PubMed] 49. Tu, K.; Yan, Z.; Wang, Y. A spatial cluster analysis of heavy rains in China. Atmos. Ocean. Sci. Lett. 2011, 4, 36–40. 50. Liu, X.; Wang, S.; Zhou, Y.; Wang, F.; Li, W.; Liu, W. Regionalization and spatiotemporal variation of drought in China based on standardized precipitation evapotranspiration index (1961–2013). Adv. Meteorol. 2015, 2015, 950262. [CrossRef] 51. Leite, G.L.D.; Von dos Santos Veloso, R.; Zanuncio, J.C.; Azevedo, A.M.; Silva, J.L.; Wilcken, C.F.; Soares, M.A. Architectural diversity and galling insects on Caryocar brasiliense trees. Sci. Rep. 2017, 7, 16677. [CrossRef] 52. Ziter, C.; Robinson, E.A.; Newman, J.A. Climate change and voltinism in Californian insect pest species: Sensitivity to location, scenario and climate model choice. Glob. Chang. Biol. 2012, 18, 2771–2780. [CrossRef]

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).