• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Spectroscopic detection of forest diseases: a review (1970-2020)

    2022-02-26 10:14:48LorenzoCotrozzi
    Journal of Forestry Research 2022年1期

    Lorenzo Cotrozzi

    Abstract Sustainable forest management is essential to confront the detrimental impacts of diseases on forest ecosystems.This review highlights the potential of vegetation spectroscopy in improving the feasibility of assessing forest disturbances induced by diseases in a timely and cost-effective manner.The basic concepts of vegetation spectroscopy and its application in phytopathology are first outlined then the literature on the topic is discussed.Using several optical sensors from leaf to landscape-level, a number of forest diseases characterized by variable pathogenic processes have been detected, identified and quantified in many country sites worldwide.Overall, these reviewed studies have pointed out the green and red regions of the visible spectrum, the red-edge and the early near-infrared as the spectral regions most sensitive to the disease development as they are mostly related to chlorophyll changes and symptom development.Late disease conditions particularly affect the shortwave-infrared region, mostly related to water content.This review also highlights some major issues to be addressed such as the need to explore other major forest diseases and geographic areas, to further develop hyperspectral sensors for early detection and discrimination of forest disturbances, to improve devices for remote sensing, to implement longterm monitoring, and to advance algorithms for exploitation of spectral data.Achieving of these goals will enhance the capability of vegetation spectroscopy in early detection of forest stress and in managing forest diseases.

    Keywords Forest management · Plant disease detection · Reflectance · Remote sensing · Spectral imaging

    Abbreviations

    ADAR Airborne data acquisition and registration

    AVIRIS Airborne visible/infrared imaging spectrometer

    CAO Carnegie airborne observatory

    CASI Compact airborne spectrographic imager

    CI Chlorophyll index

    DED Dutch elm disease

    DMSC Digital multi-spectral camera

    EVI Enhanced vegetation index

    GA Genetic algorithm

    GEOBIA Geographic object-based image analysis

    GPR Gaussian process regression

    GRSAI Green-red spectral area index

    HS Hyperspectral sensor

    HS-I Hyperspectral imaging sensor

    Lasso Least absolute shrinkage and selection operator

    LIBERTY Leaf incorporating biochemistry exhibiting reflectance and transmittance yields model

    MASTER MODIS (moderate resolution imaging spectroradiometer)/ASTER (advanced spaceborne thermal emission and reflection radiometer)

    MLD Mycosphaerella leaf index

    MNF Minimum noise fraction

    MS Multispectral sensor

    MS-I Multispectral imaging sensor

    NAIP National agriculture imagery program

    NDIwiltNormalized difference wilt index

    NDVI Normalized difference vegetation index

    NDWI Normalized difference water index

    NGRDI Normalized green red difference index

    NIR Near-infrared

    NPV Non-photosynthetic vegetation

    NWI Normalized wilt index

    PC Principal component

    PCA Principal component analysis

    PCR Principal component regression

    PLSDA Partial least squares discriminant analysis

    PLSR Partial least squares regression

    PRI Photochemical reflectance index

    PROSPECT Properties spectra model

    PWD Pine wilt disease

    PWN Pine wilt nematode

    RATIO975Three-band ratio VSI centered at 975 nm

    RATIO1200Three-band ratio VSI centered at 1200 nm

    RF Random forests

    ROD Rapid ‘ōh(huán)i‘a(chǎn) Death

    RVSI Red-edge vegetation stress index

    RWC Relative water content

    SMLR Stepwise multiple linear regression

    SOD Sudden oak death

    SVR Support vector regression

    SWIR Shortwave-infrared

    UAV Unmanned airborne vehicles

    VARI Vegetation atmospherically resistant index

    VIgreen Vegetation index green

    VIS Visible

    VSI Vegetation spectral indices

    WP Wavelength position

    Introduction

    Forests are essential ecosystems, providing numerous ecological, economic, social, and cultural services (Chen and Meentmeyer 2016) and a major portion of the global terrestrial carbon dioxide sink (Pan et al.2011).However, climate change (i.e., changes in global temperatures and rainfall) and global trade have increased forest vulnerability to a number of natural disturbances, including diseases caused by pathogens such as fungi, bacteria, viruses and nematodes (Boyd et al.2013; Wingfield et al.2015).While some pathogens are native to local ecosystems, many recent major disturbances occur from nonindigenous species that may pose dangerous and unpredictable threats to forest health (Boyd et al.2013).Over the past few decades, the frequency and intensity of forest disturbances due to diseases have dramatically increased, leading to wide tree mortality in major forest biomes worldwide (Chen and Meentmeyer 2016).

    Sustainable forest management is fundamental to confront the detrimental consequences of diseases in forest ecosystems.This is especially true when such disturbances have the potential to reduce the dominant native species, causing a permanent change in forest structure.One prerequisite for effective management is to understand the severity and spatial distribution of forest damage.Consequently, mitigation efforts can limit the population and the spread of pathogens on infected or susceptible host trees.Although conventional field mensuration remains the most accurate way to quantify stages of infection, it can be logistically challenging, time-consuming and expensive, especially when pathogen populations reach epidemic levels.

    Remote sensing makes possible a rapid and accurate evaluation of ecosystem functioning and status, and has the potential to be extended to larger spatial scales.Notably, the employment of spectroscopy in remote sensing is an outstanding approach for monitoring vegetation performance under environmental constraints (Cotrozzi et al.2018).Recent developments in space-and airborne sensors have further advanced the ability to collect observation data across multiple spatial, temporal, and spectral scales, making remote sensing feasible to monitor forest disturbances, e.g., variations in biophysical and biochemical parameters, in response to disease outbreaks of varying levels of invasion.Such rapid and accurate delineation of large-area forest damage allows decision makers to take prompt and informed action (Chen and Meentmeyer 2016).

    In order to encourage other research in this area, this paper summarizes the present knowledge on the use of vegetation spectroscopy for the detection of forest damage caused by diseases.First, it briefly reports basic concepts of vegetation spectroscopy, focusing on reflectance spectroscopy.It then reviews the literature concerning the spectroscopic detection of forest diseases.Finally, it highlights knowledge gaps and proposes new perspectives on the topic.The published peer-reviewed literature was searched in the Web of Science (Thompson-ISI, Philadelphia, PA, USA, http://apps.webof knolwe dge.com/) and Scopus (Elsevier, Amsterdam, Netherlands, http://www.scopus.com/) databases, using multiple combinations of “spectroscopy”, “reflectance”, “multispectral”, “hyperspectral”, “forest”, “tree”, “pathology” and “disease” keywords.Database searches were performed in December 2020 back to 1970.Papers that did not match the aim of the present review, e.g., those focused on crop species and/or abiotic factors, were excluded, as well as those not in English.In order to identify other suitable references, the reference lists of articles recorded by this literature search were examined.

    Vegetation spectroscopy

    Basic concepts

    Spectroscopy has been used in plant science over several decades (Gates et al.1965).Radiation striking vegetation, or any other target, may be reflected (diffuse, specular), absorbed, or transmitted (with refraction, Fig.1a).Among the different spectroscopy approaches (Conrad and Bonello 2016), reflectance spectroscopy is the one most used for measuring vegetation in the field.Reflectance spectroscopy relies on the interaction of light with plant structure, chemical configuration and water content reflectance being dependent on variations in absorption as a consequence of vibrational excitation of molecular bonds, e.g., C-H, N-H, and O-H bonds, at specific wavelengths.Reflectance of vegetation can be calculated from the measurement of the vegetation radiance divided by the radiance of white reference panels, and defined as percentage (Kumar et al.2001).

    Fig.1 (a) Incoming light (L) reaching the surface of vegetation may be reflected (R; diffuse, specular), transmitted (T; with refraction) or absorbed (A); (b) during pathogenesis, leaf pathogens, e.g., biotrophic (left) and necrotrophic (right), affect leaf structural and chemical properties, and by this leaf reflectance as well as other leaf optical properties is distinctively altered in the visible spectra (VIS: blue, b; green, g; red, r), near-infrared (NIR) and shortwave-infrared (SWIR); leaf sections re-drawn from Mahlein (2015)

    Within the full range of the solar radiation reaching vegetation (i.e., 350-2500 nm), the typical reflectance profile of a green leaf is characterized by three regions (Fig.2): (1) the visible VIS, 400-700 nm showing low reflectance due to strong absorption by leaf pigments; (2) the near-infrared, NIR, 700-1100 nm characterized by high reflectance and mainly related to leaf cell structure; and, (3) the shortwaveinfrared SWIR, 1100-2500 nm characterized by another reduction of reflectance due to lower incident radiation as well as to the absorption by water and other absorptions by biochemical metabolites, e.g., lignin, starch, cellulose, proteins, phenols (Curran 1989; Asner and Martin 2008; Serbin et al.2015).Another feature of the vegetation reflectance spectrum is the red-edge, i.e., the region in the VIS red-NIR transition zone where reflectance increases drastically, which has been shown to be dependent on chlorophyll content and stress conditions.In light of the above, variations in spectral data can be used for various purposes such as mapping plant functional traits, species biodiversity and distribution, and monitoring plant responses to resource availability and environmental constraints (Kumar et al.2001; Cotrozzi et al.2018).

    Fig.2 Typical full range, i.e., 350-2500 nm, reflectance profile of healthy vegetation; major absorption and reflectance features are indicated.VIS: visible, 400-700 nm; NIR: near-infrared, 700-1100 nm; SWIR: shortwave-infrared, 1100-2500 nm; gray stripes show spectral bands commonly measured by multispectral sensors

    Optical sensors are classified on their spectral resolution and range.Multispectral sensors (MS) measure reflected energy from few to many specific discrete sections, i.e., bands, of the electromagnetic spectrum, e.g., 3-12 bands covering a few tenths of wavelengths each, common only in the VIS-NIR; Fig.2.Hyperspectral sensors (HS) measure energy in narrower and more numerous bands than MS (e.g., hundreds or thousands covering 5-20 nm each).The numerous narrow bands of HS provide a continuous spectral measurement which is more sensitive to subtle variations in reflected light, providing more data complexity and information content.Since the 1990s, several hand-held HS spectroradiometers covering the full VIS-NIR-SWIR range, i.e., 350-2500 nm, have become available and widely used for field vegetation research (Cotrozzi et al.2018).Evidently, MS (or HS covering narrower wavelength ranges) are smaller in volume, lighter in weight, needing fewer internal components for working, and are less expensive than HS in collecting wider wavelength ranges, and for these reasons are more used for remote sensing collections by unmanned airborne vehicles (UAV).

    Reflectance measurements of vegetation are rapid, nondestructive, and relatively economical since the costs to operate such equipment once acquired are minimal compared with the costs to run other standard analytical instrumentation, e.g., spectrophotometry and gas/liquid chromatography systems.This technique also allows to investigate markedly more individual plant traits in situ and in vivo over numerous time periods than standard measurements alone, e.g., those carried out with classic ecophysiological approaches or wet chemistry.In addition, this approach can monitor plant functions over large geographic areas if scaled to remote sensing measurements from space, i.e., satellites, or airborne platforms, e.g., UAV and aircrafts (Cotrozzi et al.2018).

    In a remote sensing context, e.g., forest monitoring, imaging systems are usually preferable compared to nonimaging systems, this being a further categorization of spectral sensors.Non-imaging spectral sensors average the spectral information with spatial resolution over a definite area within their field of view, and are preferable when the user wants to learn more about the radiation/vegetation interaction, and spatial information such as mapped out is not needed e.g., leaf-and ground canopy-level characterization of plant traits.Imaging spectral sensors are different and supply spectral information with spatial resolution of the imaged object; data is generated in the form of a 3D spatial map of spectral variation, where the first two dimensions (x-andy-axes) show the spatial information and the third dimension (z-axis) provides the spectral information (Fig.3).For each pixel in an image, multispectral imaging (MS-I) devices capture a small number of spectral bands, typically three to fifteen, whereas hyperspectral imaging (HS-I) instrumentation capture continuous spectral ranges.Obviously, imaging spectral data needs more pre-processing (e.g., atmospheric calibration, geometric correction and spectral calibration) than non-imaging ones.Nevertheless, being an integration of imaging and conventional point spectroscopy, imaging spectroscopy receives complementary information from both domains (Shaw and Burke 2003).Point spectroscopy provides information to explain the physiology, biochemistry and morphology of plants, while the information from imaging technology is used to design the dynamics of these vegetation characteristics (Asner and Martin 2008; Serbin et al.2015; Singh et al.2015; Mishra et al.2017; Mahlein et al.2018).

    Fig.3 Principle of imaging spectroscopy.Space-and airborne imaging spectral sensors provide spectral information with spatial resolution of the imaged object.Figure has been re-drawn from Shaw and Burke (2003)

    In addition, some physically based models have been successfully used to simulate leaf reflectance and transmittance and to estimate leaf biochemical properties, e.g., radiative transfer models where variation in canopy reflectance can be modeled by relevant plant functional traits that are known to specifically affect light transmission in canopies.The optical PROperties SPECTra (PROSPECT, Féret et al.2017) and the Leaf Incorporating Biochemistry Exhibiting Reflectance and Transmittance Yields (LIBERTY, specifically designed for pine needles; Dawson et al.1998) are among the most widely validated models.

    Advancements in the portability and sensitivity of spectrometers, together with improvements in computational and chemometric efficiency, have allowed the exploitation of reflectance data by using different approaches.First, using simple vegetation spectral indices (VSI) based on the ratio of reflectance at distinct wavelengths that have been developed to estimate several foliar traits, e.g., normalized difference vegetation index, NDVI (Rouse et al.1974); photochemical reflectance index, PRI (Gamon et al.1997); normalized difference water index, NDWI (Gao 1996); chlorophyll index, CI (Merzlyak et al.1999).Second, using multivariate approaches, e.g., partial least squares regression, PLSR (Wold et al.2001), which has become the current standard of chemometric modeling, stepwise multiple linear regression, SMLR (Grossman et al.1996), principal component regression, PCR (Martens 2001), and least absolute shrinkage and selection operator, Lasso (Tibshirani 2011) to directly estimate widely investigated physiological, biochemical and morphological plant traits as a function of the leaf and vegetation hyperspectral profiles (e.g., Serbin et al.2014, 2015; Couture et al.2016; Marchica et al.2019; Cotrozzi et al.2020).Third, testing spectra themselves as a phenotypic expression of the combined signals of the physiological, biochemical and morphological features of leaves under particular environmental conditions (i.e., spectral signatures), analyzable through multivariate approaches such as partial least squares discriminant analysis (PLSDA) or other classification techniques such as support vector machine learning (Cotrozzi and Couture 2020; Calzone et al.2021).These approaches can be applied on untransformed or pre-treated spectra, i.e., derivative spectra is often used for the elimination of background signals and for resolving overlapping spectral features, with remote sensing collections particularly needing data pre-processing, e.g., signal compression and noise reduction (Rasti et al.2018).

    Applications in phytopathology

    Vegetation spectroscopy has shown considerable potential to supply new insights into phytopathology.During host-pathogen interactions, a series of processes occur that induce variations in the plant’s physiology, biochemistry and morphology.As a consequence, these processes during pathogenesis and symptom development affect the optical properties of vegetation (Fig.1b; Mahlein 2015).Therefore, variations in reflectance due to plant pathogens and diseases can be interpreted as highly specific changes in the structure and chemical configuration of the leaf tissue in the course of pathogenesis, e.g., sequence of chlorotic and necrotic tissue or the presentation of characteristic fungal structures.Each host-pathogen interaction has a distinctive spatial and temporal progression and influences specific wavelength ranges during distinct pathogenesis phases (Khaled et al.2018; Mahlein et al.2018).For instance, the occurrence of visible symptoms could be early detected by variations in the VIS spectral region, or the loss of leaf water due to the challenge of a necrotrophic pathogen could be monitored thanks to the water sensitivity of SWIR.As a consequence, reflectance spectroscopy has become one of the most powerful non-invasive approaches used for the detection (deviation from healthy, even at early stage), identification (diagnosis of specific symptoms among others and discrimination of various diseases) and quantification (assessment of disease severity) of plant diseases at multiple scales (Lowe et al.2017; Thomas et al.2018).

    Although plant disease research including such spectral approaches have risen sharply in the last 20 years, most of these studies have been carried out for agricultural applications, specifically within the areas of precision farming, aimed to optimize agricultural procedures and the application of agrochemicals by considering in-field and temporal variability, and plant phenotyping, aimed to assess multiple plant traits to evaluate a specific genotype/environment interaction.Indeed, numerous reviews on the use of vegetation spectroscopy for detecting and monitoring major crop diseases are available (e.g., Bock et al.2010; Mahlein 2015; Lowe et al.2017; Khaled et al.2018; Mahlein et al.2018; Thomas et al.2018; Zhang et al.2020a), and highly suggested to the reader interested in learning the technical and analytical details of the approach.Conversely, the use of vegetation spectroscopy for detecting and monitoring forest diseases has been less investigated.Lower investments and/or greater experimental difficulties might be among the causes of this tendency for forestry, and a review devoted to summarizing the state-of-the-art of this topic was lacking.Chen and Meentmeyer (2016) and Stone and Mohammed (2017) reviewed the application of remote sensing approaches for monitoring forest damage by both insects (mainly) and diseases.

    Spectroscopic detection of forest diseases

    While the idea of applying remote sensing to detect diseaseinduced forest damage was considered many decades ago (e.g., Heller and Bega 1973), only since the late 1990s has the approach began to receive considerable attention for managing disease outbreaks, as shown by reviewed papers reported in Tables 1 and 2.Two reasons possibly explain such slow adoption: (1) an increasing number of studies have showed that the frequency and intensity of forest diseases have significantly increased over the past two decades as a result of globalization, responsible for a worldwide distribution of pests and pathogens, and climate change (Boyd et al.2013), highlighting the need for understanding disease progression to implement effective mitigation strategies; and,(2) the collected spectral data have increased dramatically during the same time period, mainly due to the abatement of data acquisition costs.Geographically, research hotspots have been primarily located in North America and Australia, with several other studies conducted in Europe, East Asia and South Africa.

    Table 1 Spectroscopic assessment of forest diseases (conifers)

    Table 2 Spectroscopic assessment of forest diseases (broadleaf)

    Table 2 (continued)

    In this review, the studies are presented according to the following three criteria: (1) works on conifers (Table 1) are discussed before those on broadleaf species (Table 2); (2) within each tree type, genera mostly studied (i.e., higher number of reviewed papers) are presented first; and, (3) within each genus, diseases mostly studied are listed first.For each disease, a brief description, e.g., causal agent, distribution, symptoms, is given, and then the related studies mentioning the adopted sensors and data analysis approaches, the most disease-sensitive spectral regions, and the major outcomes.

    Pinus spp.

    Sphaeropsis blight

    Sphaeropsis sapinea(orDiplodia pinea) is among the most frequent fungal pathogens affecting conifers worldwide, causing shoot blight, crown wilt, cankers and sapwood stain that have produced heavy losses in pine stands.Usually top and outer branch shoots are first affected, quickly becoming pale green, then yellow, orange and red, and finally being shed, while the lower crown may stay green and of normal density (Stanosz et al.2002).

    Sims et al.(2007) used decision tree analysis to evaluate the crown condition ofPinus radiataD.Don (Monterey pine) stands in southern New South Wales, Australia showing discolored leaves due toS.sapinea.Imagery was collected by the digital multi-spectral camera (DMSC) II system (SpecTerra Services) with four wavelength filters, 680, 720, 740, and 850 nm in September 2002; 680, 700, 720, and 750 nm in September 2003.Imagery was acquired at a pixel resolution of 0.5 × 0.5 m.Spectral VSI and fraction images (see Goodwin et al.2005, reported below) were used to classify crowns into binary condition classes, i.e., proportional discoloration, 0: < 40%, 1: > 40%.The developed decision tree model identified crown discoloration due toS.sapineawith an effective accuracy between 83 (Kappa) and 92% overall accuracy.This study confirmed the outputs reported by Coops et al.(2006) which analyzed the same dataset by evaluating the relation between the crown-condition variables and the spectral VSI through a standard linear regression, and through a classification error matrix, i.e., confusion matrix.The slope of the upper red-edge was the feature mostly related with crown conditions (R2=0.76 and 0.88 for 2002 and 2003, respectively), with an independent classification accuracy of over 90%.In fact, several studies have shown that the shape of the red-edge is dependent on chlorophyll content (Smith et al.2004; Zarco-Tejada et al.2004) and stress conditions (Mutanga and Skidmore 2007).On the same dataset collected in 2002, Goodwin et al.(2005) applied a linear spectral unmixing approach to measure the fractional abundances of the key image endmembers: sunlit canopy, shadow and soil.Applying multiple linear regression and endmember fraction images, high significance rates were reported for the color of the crown leader (a six class variable;R2=0.79) and the overall crown color (proportional discoloration, 0: < 40%, 1: > 40%;R2=0.80).In general, the sunlit canopy image fraction resulted in the most meaningful variable for the estimation of forest condition.

    These outputs are in accordance with Coops et al.(2004) who adopted a similar linear spectral unmixing approach but investigated the potential of the Compact Airborne Spectrographic Imager 2 (CASI-2) to detect the sameS.sapineainfection.Results indicated that four fraction endmember, i.e., sunlit canopy, soil, shadow and non-photosynthetic vegetation (NPV) images could be accurately obtained from the 12 spectral bands of CASI-2 imagery (ranging from 450 to 850 nm, spatial resolution: 0.8 m).Single crowns were detectable on the sunlit crown fraction.Interestingly, NPV showed scattered patches of high fractional abundance where there were aggregates of brown needles.The stepwise regressions for modeling the color of the crown leader and the overall crown color showed highR2(0.75 and 0.72, respectively), with NPV followed by the shadow fractions showing the higher significances.

    Concurrently in the same experimental site, Coops and Stone (2005) collected needle spectra by a field UNISPEC spectrometer (PP Systems; 350-1100 nm), and took needle samples to determine chlorophyllaandb, and water contents.Under increasingS.sapinea-induced disease, needles with a lowering of chlorophyll content showed an increase in reflectance in the VIS region.The shape of the red-edge significantly changed among theS.sapineaseverity classes (i.e., low, moderate, high).Needle chlorophyll concentration was mostly correlated with various VSI, especially those including reflectance at 705 nm (chlorophyll sensitive) and at either 750 or 445 nm (chlorophyll insensitive).A good correlation was also reported between moisture content and the water band index (Pe?uelas et al.1995), suggesting that the inclusion of water sensitive wavelengths (SWIR) may also help to better monitorS.sapineainfection by spectral data.Furthermore, the LIBERTY radiative transfer model was also applied to estimate theoretical needle reflectance (400-2500 nm) by changing chlorophyll and water content input parameters.Such LIBERTY estimated spectra mostly matched both trends and rates of actual spectra, suggesting that the utilization of radiative transfer models, properly parameterized, can produce information in discriminating needle damage severity.

    Pine wilt disease

    Pine wilt disease (PWD) is an important quarantine disease responsible for the overwhelming death ofPinusspecies caused by the pine wood nematode (PWN)Bursaphelenchus xylophilus.This is the only pathogen that is not a microorganism among those discussed in this review.B.xylophilusis transmitted by vector pine sawyer beetles in the genusMonochamus.At present, PWD is primarily distributed in eight countries across North America (United States, Mexico, and Canada), East Asia (China, Japan, South Korea, North Korea) and Europe (Portugal, Spain).However, 52 countries worldwide have listed PWN as a quarantine pathogen as PWN can spread by transport of infested wood.At an early stage, PWD cannot be detected with the naked eye.PWN starts to multiply in xylem cells, obstructing their internal ducts.Trees affected by PWD cease resin secretion and needles gradually become yellowish.Finally, the needles turn brown and the plant dies (Wu et al.2020).

    Spectroscopy assessment of PWD was performed on variousPinusspecies.Ju et al.(2014) and Zhang et al.(2020b) investigated the capability of using hyperspectral data over the 350-1100 nm spectral range, acquired by an ASD Field-Spec spectrometer; Analytical Spectral Devices, to distinguishPinus massonianaLamb.(Chinese red pine) trees infected by PWN from healthy ones.Spectral measurements were collected ca.0.5-1 m above the canopy.Ju et al.(2014) used first derivative spectra and VSI to reduce data dimensionality and to identify the most important wavelengths for detection.The most effective first derivative metric was at 759 nm, and thus selected to classify the health status ofP.massoniana.The NDVI, calculated from reflectance at 810 and 450 nm, in the fully infected stage was strongly related to changes of chlorophyll concentration (r=0.95).The researchers concluded that the combined use of specific spectral features and chlorophyll concentration are an efficient method to detect PWN infection around 30 days after inoculation.Zhang et al.(2020a, b) also observed that samples in different PWD infection states, a five-class variable, generated distinctive reflectance and their firstorder derivative profiles, but they used a different spectra exploitation approach.To indicate such spectral changes, 16 spectral characteristics were extrapolated from the green (510-580 nm), red (620-680 nm), red-edge (680-760 nm) and NIR (780-1100 nm) regions, and were coded as genes constituting the chromosome of a genetic algorithm (GA).Using the optimal spectral features with suitable fitness from the GA, an accurate (R2=0.82) PLSR model was developed to predict PWD severity.

    Kim et al.(2018) focused differently on hyperspectral analysis ofPinus thunbergiiParl.(black pine) inoculated with PWN in Geoje Island, South Korea.Full range leaf reflectance (350-2500 nm) collected with a GER 3700 spectrometer (Geophysical and Environmental Research Corporation) showed PWD-induced changes in the red (ca.600-700 nm) and SWIR (ca.1400-1500 nm) regions within two months from PWN infection.The already existing vegetation atmospherically resistant index (VARI), the vegetation index green (VIgreen) and the normalized wilt index (NWI), as well as the green-red spectral area index (GRSAI) they specifically introduced, constantly reported significant differences between tree infection conditions earlier than the other tested indices.Notably, the GRSAI, which includes the red and green regions, showed less fluctuation over time than the other VSI.

    Recently, Iordache et al.(2020) proposed a machine learning algorithms detection of PWD inPinus pinasterAiton (maritime pine) from multi-(MicaSense Red-Edge M multispectral camera, MicaSense Inc.; five spectral bands centered at 475, 560, 668, 717, and 840 nm) and hyperspectral data (MicroHyperSpec A-series hyperspectral scanner, Headwall Photonics; 380-1100 nm) acquired over two regions of central Portugal affected by PWD.Classification schemes for both multi-and hyperspectral imagery (spatial resolution: 5 and 10 cm, respectively) used sets of 13 selected VSI, and were developed by random forest approaches.These classification schemes achieved high PWD identification overall accuracies (> 91%).Moreover, it was shown that PWD was identifiable by spectral imagery even before symptoms could be detected by the naked eye.The development of PWD was also detectable in the derived maps.

    All these results clearly show the practicability of PWD detection by spectral imaging and non-imaging, as further confirmed by the extensive review by Wu et al.(2020).

    Red band needle blight

    Red band needle blight, caused byDothistroma septosorum, is a fungal disease with a significant impact on pine forests worldwide, affecting over 80 pine species.The development of disease symptoms includes initial chlorosis, rapid loss of cellular integrity, cellular necrosis (bands) and eventual phylloptosis or leaf fall.Infection commonly begins in the lower portions of the crown on older foliage and spreads to higher parts (Bradshaw 2004).

    Smigaj et al.(2019) examined the ability of hyperspectral data for detection of red band needle blight.Airborne VIS-NIR hyperspectral data (450-980 nm, 2 m spatial resolution) were collected in central Scotland using the airborne Aisa FENIX imager (Channel Systems) flown over diseasedPinus contortaDouglas (lodgepole pine) plantations with variable defoliation and infection status.Spectra of foliage affected by red band needle blight were also collected in a youngPinus sylvestrisL.(Scots pine) stand by an ASD FieldSpec Pro spectroradiometer (350-2500 nm).The green, red and early NIR spectral regions were the most sensitive to the development of the disease.First order derivative values of the green to red-edge region performed best at separating healthy and infected leaves with an accuracy around 90%.The best index was the normalized green red difference index (NGRDI), including green and NIR regions, followed by the enhanced vegetation index (EVI).Transferring the results obtained at leaf level to airborne collections carried out on lodgepole pine sites, resulted in a good distinction comparing the more diseased/defoliated plot and the marginally affected pine stand using numerous spectral indices, including PRI, NGRDI and EVI.Smigaj et al.(2019) also investigated LiDAR data and suggested that the combination of data from different sensors could improve the ability to detect red band needle blight.

    In fact, in the early 2000s, leaf and canopy spectral assessments of red band needle blight were, respectively, investigated by Stone et al.(2003) and Coops et al.(2003) overP.radiataplantations in New South Wales, Australia which had been ground evaluated and ranked following the extent ofDothistromaneedle blight (a six-class variable) on an individual tree basis.Stone et al.(2003) used an ASD FieldSpec FR spectroradiometer (350-2500 nm) and found that the largest reflectance difference between healthy and diseased leaves took place on the shoulder of the NIR region at 763 nm.Wavelengths in the 675-691 nm range showed the greatest sensitivity toD.septosporuminfection, followed by those close to 760 and 550 nm (shown by first derivative spectra), while the lowest sensitivity resulted at 434, 493, 506, 709, and 1373 nm.The 709/691 nm reflectance ratio best correlated to needle damage.Among the other tested VSI, an index calculated by the 710-740 nm wavelengths, i.e., the upper red-edge, was also strongly related to disease severity (r=-0.73).Conversely, a weak relation was found between the needle damage category and the slope maximum of the red-edge.

    Coops et al.(2003) acquired airborne remote sensing imagery by the CASI-2 (10 bands ranging from 450 to 850 nm, 0.8 m spatial resolution) and reported thatDothistromaneedle blight was more severe in the lower crown, and significant relationships were found between ground evaluation and crown reflectance features, i.e., lower and upper slope of the red-edge and red-edge vegetation stress index, RVSI, using an ‘halo’ imaging approach which ignored the brightest central pixels of each tree crown.This method showed an accuracy of over 70%.Independent accuracy evaluation of the approach showed that the technique could successfully detect three levels (the original six-class classification system was aggregated into three classes) ofDothistromaneedle blight.

    Blister rust

    Blister rust, due to the invasive fungusCronartium ribicola, is a dominant source of stress and mortality to high-altitudePinus albicaulisEngelm.(whitebark pine) within the Greater Yellowstone Ecosystem, California as well as mountain pine beetle (Dendroctonus ponderosaeHopkins).Prediction and control of blister rust is challenging since the pathogen does not diffuse from tree to tree but by way of wind-borne spores from the population of its alternate hostRibesspp.The development from green to red needles at the outer-most branches is an early symptom of blister rust infection and mountain pine beetle infestation (Geils et al.2010).

    Hatala et al.(2010) analyzed images acquired by the HyMap airborne hyperspectral sensor (HyVista Corporation Pty.Ltd.; 128 bands across the 450-2500 nm spectral range, 2-4 m spatial resolution) to identify the sites of blister rust-induced stress and mortality in the Greater Yellowstone Ecosystem whitebark pine crowns by using sub-pixel mixture-tuned matched-filter analysis.Discrepancies in the spatial trends of blister rust and mountain pine beetle infestation allowed to distinguish areas dominated by mountain pine beetle against blister rust by analyzing variations in the spatial scale of significant stress and mortality clusters estimated by the Ripley’sKalgorithm.

    Pitch canker

    Pitch canker is an episodic disease caused by the fungusFusarium circinatum.It affects onlyPinusforests, withPinus radiatabeing among the most susceptible.Pitch canker was first reported in North Carolina, USA in the 1940s, but the pathogen is now present over the whole southeastern United States and in numerous countries worldwide (Haiti, Chile, South Africa, Japan, Korea, Mexico, Italy, Portugal and Spain).The disease, characterized by strong resin exudation at the infection site on branches, negatively impacts growth and is also able to kill seedlings and mature trees (Wingfield et al.2008).

    Poona and Ismail (2013) explored the capability of transformed high spatial resolution QuickBird satellite imagery (DigitalGlobe; 447-512, 499-594, 620-688, and 755-874 nm, 2.44 spatial resolution) and artificial neural networks to detect and mapP.radiataplantations affected by pitch canker disease in South Africa.Tree crowns were outlined by an automated segmentation and classification procedure within an object-based image analysis environment.Various VSI, including the tasseled cap transformation, were then estimated and included into a neural network model.The feed-forward neural network allowed a disease discrimination with an overall accuracy of 82% and KHAT of 0.65.These outcomes highlight the potential of the approach for future crown-level mapping of the disease at an ecosystem level.

    Pseudotsuga menziesii (Mirbel) Franco (Douglas-fir)

    Laminated root rot

    Phellinus weirii, the fungus causing laminated root rot (this is the only root disease discussed in this review), is among the most significant single natural agent inducing long-term variations in the coastal forests of western North America.Pseudotsuga menziesii(Douglas-fir) is highly susceptible toP.weiriiand is the most economically important host of the fungus.P.weiriiinfection occurs when healthy roots encounter infected stumps or diseased roots of a previous stand,P.weiriigradually induces root decay, leading to weakened uptake of water and nutrients and reduced structural support.Crown symptoms, reduced growth, needle chlorosis and loss, massive production of small cones, usually occur only after the fungus has compromised a significant portion of the root apparatus (Thies 1998).

    Leckie et al.(2004) acquired two sets of CASI imagery (spatial resolution: 0.6 m) over a Douglas-fir site in British Columbia, Canada in 1995 (five spectral bands from 456 to 779 nm) and 1996 (eight spectral bands from 438 to 861 nm).Trees showing different levels of root rot symptoms were evaluated in the field, manually outlined on the imagery, and associated to their spectral properties.The NIR and red bands and band ratio including these two regions resulted in most associated to root rot injury, but a blue region was also effective.Using these four spectral bands, an average classification accuracy of 55-60% was reported for healthy, slightly infected, light to moderately severe, 100% needle loss, snag and shadowed snag conditions, while the accuracy for broader classes was much better, e.g., around 80% if a tolerance of ± 1 class was allowed.Specifically, there was a good relationship between areas of root disease with moderate and severely damage trees, while lightly damaged trees were not reliably identified.

    Quercus spp.

    Sudden oak death

    Sudden oak death (SOD), caused by the pathogenPhytophthora ramorum, was first reported in the mid-1990s on tanoak (Notholithocarpus densiflorus) and oaks (Quercusspp.) in the San Francisco Bay area, then spreading epidemically in several coastal forests in California and Oregon, USA.This disease has also had a devastating impact on larch stands in the United Kingdom.Phytophthora ramorumhas been either recorded as a regulated species or listed in forest legislations in 68 countries worldwide.Symptoms include bark cankers, leaf spots and twig dieback (Grünwald et al.2019).

    Sudden oak death has received the most attention in terms of remote sensing.He et al.(2019) combined multisensor remote sensing with species distribution modeling to map SOD in a California forested region of 80,000 ha from 2005 to 2016.To permit both the fine-scale monitoring of disease development and adequate description at broader scales, they used a spectral unmixing approach to extract sub-pixel disease infection produced by a species distribution model.They calibrated/validated the approach with imagery from high-spatial resolution National Agriculture Imagery Program (NAIP) of the USDA Farm Service Agency, and hyperspectral Airborne Visible/Infrared Imaging Spectrometer (AVIRIS; 360-2500 nm, 3 m spatial resolution) data, Google Earth? imagery, and field investigation.The outcomes revealed an annual SOD infection of 7%, with 76-83% overall mapping accuracy.The combination of multi-sensor remote sensing and species distribution modeling markedly decreased the overestimation of disease impacts as compared to the unique use of remote sensing sensors.

    A similar remote sensing study was performed by Pu et al.(2008) using a CASI-2 dataset (2 m spatial resolution) for detecting SOD-induced tree stress and mortality in another California forest, including costa live oak (Quercus agrifolia), black oak (Q.kelloggii), and valley oak (Q.lobata, no SOD host) living in mixed stands with mature madrone (Arbutus menziesii) and California bay trees (Umbellularia californic).They first developed a multilevel classification (0 to 3) design to enhance classification efficiency.The CASI raw data were then converted to reflectance and rectified for topography, and a principal component (PC) transformation of all 48 bands and the visible and NIR bands were individually carried out to select features from the data.Lastly, the processed CASI imagery was classified by a maximum likelihood approach and tested for classification accuracies.The outcomes showed that classification accuracy of oak tree stress mostly increased as the detailed classification level increased, ranging from 3 to 76%.Furthermore, PCs from VIS and NIR regions led to more accurate classification than the PCs from all 48 CASI bands.

    The work by Pu et al.(2008) was a follow-up of two studies at leaf level (Pu et al.2003, 2004), where reflectance, collected by an ASD FieldSpec ProFR spectroradiometer; 350-2500 nm, and relative water content (RWC, %) were determined from hundreds of coast live oak (Q.agrifolia) leaf samples with three levels of SOD symptoms, healthy, infected and newly dead.Pu et al.(2003) conducted correlation analysis between RWC and absorption features, i.e., wavelength position (WP), absorption feature depth, width, and area at the 975, 1200, and 1750 nm positions, and threeband ratio VSI centered at 975 and 1200 nm (RATIO975and RATIO1200, respectively).Among these features, area showed the highest and more stable relation with RWC.Pu et al.(2004) extracted the maximum or minimum first derivative and its related WP from 10 spectral slopes along each reflectance profile.Using all 306 spectra collected on healthy, infected and newly dead leaves, maximum or minimum derivative values at both sides of the reflectance depressions near 1200, 1400, and 1900 nm, the WPs on the right of the 1200 nm depression, and both sides of the 1400 nm and 1940 depressions were highly correlated with RWC.With a selection of 260 spectra collected only on green and green-yellowish leaves, discarding newly dead leaves, the WPs at the right side of the 1400 nm depression, at the left side of the 1940 nm depression and at the red region showed a constant relation with leaf RWC.As a result, RATIO975and RATIO1200indices, as well as the spectral feature on the right of the 1400 nm spectral depression may be crucial spectral indicators for evaluating water status and for early detection of SOD before the occurrence of symptoms.Overall, this study confirms the importance of including the water sensitive SWIR region for monitoring necrotrophic pathogens.

    Changing perspective, Hultquist et al.(2014) and Chen et al.(2015) examined the effects of fire and SOD on forests, focusing on the evaluation of post-fire burn severity where trees experienced three stages of disease development, pre-fire: early-trees retaining dried foliage and fine twigs, middle-trees losing fine crown fuels, and late-stage, trees falling.The experimental activities were carried out by using Geographic Object-Based Image Analysis (GEOBIA) to MASTER MODIS (Moderate Resolution Imaging Spectroradiometer)/ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) airborne images with high-spatial (4 m) and high-spectral resolutions (50 bands covering VIS, NIR and SWIR spectral regions, as well as other longer infrared wavelengths, i.e., 400-14,000 nm) that were collected immediately after the fire for rapid assessment.Hultquist et al.(2014) investigated the capability of the Gaussian process regression (GPR), random forests (RF) and support vector regression (SVR) machine learning techniques.RF showed the best performance in burn severity assessment, and lower sensitivity to the use of multiple combinations of remote sensing variables as well.Chen et al.(2015) also investigated two commonly used wavelength reduction processes, PCA, principal component analysis, and MNF, minimum noise fraction, for the outline of image objects and the consecutive accuracy of burn severity models using either PCA or MNF derived variables.Both PCA and MNF were promising for balancing computation effi-ciency and the accuracy of burn severity models in forests subjected to the middle and late stages of SOD progression, but not for those in the early stage of the disease.

    A number of studies also focused on mapping SOD using the airborne data acquisition and registration (ADAR) 5500 imager (four spectral bands: blue, 450-550, 520-610, 610-700, and 780-920 nm) in California (Kelly 2002; Kelly and Meentmeyer 2002; Kelly and Liu 2004; Kelly et al.2004; Liu et al.2006, 2007; Guo et al.2007).These studies are not presented in detail since their focus was on the imaging approach instead of the spectroscopy one.

    Japanese oak wilt

    Japanese oak wilt disease, caused by the fungusRaffaelea quercivoraand transmitted by the bark beetlePlatypus quercivorus, is responsible for considerable mortality of Fagaceae species in Japan.Following infection,R.quercivoraimpedes the rising of sap through the vessels, finally leading to tree death (Kubono and Ito 2002).

    Since the detection of affected trees is the only strategy to avoid tree mortality, Uto et al.(2008, 2009) proposed an automatic detection approach using remotely sensed hyper-or multi-spectral data to overcome the problem.After the collection of leaf spectra by an Imspector V10 imaging spectrograph (Specim, Spectral imaging; 400-1000 nm), they first introduced a normalized wilt index (NWI) which was then validated on hyperspectral remote sensing data and a multispectral satellite imagery acquired by the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER, NASA Jet Propulsion Laboratory).Since the performance of this index on leaves showing various autumnal tints was not initially verified, this gap was covered by Uto et al.(2011), who proposed another index, i.e., normalized difference wilt index, NDIwiltto be applied for wilt identification in leaves collected in October.

    Oak wilt (and bur oak blight)

    Oak wilt is a severe disease of oaks (Quercusspp.) over large regions of the eastern and central United States, and is caused by the invasive fungal pathogenBretziella fagacearum, previously known asCeratocystis fagacearumBretz, which spreads in the outer sapwood.Oak wilt is transmitted by bark beetles and sap-feeding beetles, but can also be transmitted over short distances via root grafts between adjacent trees.After infection, trees respond by producing tyloses and gum which consequently reduce the flow of water and nutrients in infected vascular tissues, leading to tree wilt and finally death (Koch et al.2010).

    The ability of reflectance measurements to detect oak wilt disease was firstly shown by Everitt et al.(1999).They performed a study in south-central Texas, USA where airborne reflectance collections in the VIS green (520-600 nm), VIS red (630-690 nm) and NIR (760-900 nm) spectral bands were made with a Barnes modular multispectral radiometer (Barnes Engineering Company) on branches ofQuercus fusiformis(escarpment live oak) showing different stages of disease.Reflectance data collected in the field showed that dead, diseased and healthy trees had divergent VIS and NIR levels; oak wilt resulted in increasing green and red values and in decreasing NIR.

    More recently, Fallon et al.(2020) evaluated the capability of spectral data to discriminate oak wilt from other processes of decline.They exposed greenhouse-grown seedlings (Quercus ellipsoidalisandQ.macrocarpa) to drought, or to oak wilt or bur oak blight fungi (Tubakia iowensis, a common leafspot fungus).They collected leaf and canopy reflectance by a Spectra Vista HR 1024i spectroradiometer (Spectra Vista Corporation; 400-2400 nm) and instantaneous photosynthetic and stomatal conductance levels, then ran PLSDA to estimate the exposure type from hyperspectral data.They also detected oak wilt before the onset of symptoms, accurately classified the disease in symptomatic samples, and showed that classification accuracy from spectral data increased with oak wilt-induced reductions in photosynthetic performance.Wavelengths found to be mostly important for oak wilt diagnosis occurred only in NIR and SWIR spectral regions, and were related to water status, non-structural carbohydrates and photosynthetic processes.They showed that hyperspectral models could discriminate oak wilt from other factors of tree damage, and that detection was related to biological processes linked to oak wilt infection and progression.They also found that within the canopy, symptom variety could negatively affect detection, but that symptomatic leaves and canopies were accurately diagnosed by spectral information.

    Eucalyptus globulus Labill.

    Mycosphaerella leaf disease

    NumerousMycosphaerellaspp.have been related to leaf diseases of manyEucalyptusspp., generally specified as Mycosphaerella Leaf Disease (MLD).Juvenile leaves and shoots ofEucalyptustrees are mostly affected by MLD, with symptoms including leaf spots, phylloptosis or premature leaf drop, twig cankers and reduced tree growth.However, manyMycosphaerellaspp.are also able to infect adult leaves owing to their ability to develop a proto-appressorium that allows direct cuticle penetration.As a consequence, sometimes trees are infected by an array ofMycosphaerellaspp.(Hunter et al.2011).

    Early investigations carried out spectral sensitivity analysis at the leaf level to select wavelengths capable of detecting increasing MLD severity onEucalyptus globulusLabill.Refelctance measurements were performed by Pietrzykowski et al.(2006) on healthy and infected juvenileE.globulusleaves collected in a plantation in Tasmania, Australia using an ASD FieldSpec FR spectroradiometer (400-1000 nm).Reflectance at 760 nm gave the largest difference between the minimum and maximum of infection severity (a six-class variable).The wavelength most sensitive to disease severity was at 680 nm, followed by those between 480 and 500 nm.Conversely, reflectance near 708 nm and 550 nm was relatively insensitive to disease severity.The VSI most related to leaf infection was the ratio between refelctance at 678 and 550 nm (r=0.84), which was also used in a linear regression model for accurately estimating leaf infection severity from an independent data set (r=0.95).This index was then tested on the same Tasmanian plantation site with airborne multispectral imagery (digital multi-spectral camera; four spectral bands at 550, 680, 740, and 780 nm, 0.5 m × 0.5 m pixel size) to detect necrosis at a site characterized by strong MLD infection (Pietrzykowski et al.2007).Models to detect crown-level necrosis using this index reported a classification accuracy of 67% (Kappa=0.54) according to field evaluations, and a similar accuracy for a model developed to detect crown-level defoliation which included reflectance at 780 nm.Although maps drawn by these linear models showed precisely the pattern and gradient of MLD severity recorded in the field, this accuracy could likely be improved by utilizing spectral data at a higher resolution, allowing for the calculation of more complex VSI and the exploitation of whole spectral regions by multivariate approaches.

    Therefore, Barry et al.(2011) carried out two studies focused on necrosis detection in pot-grownE.globulusseedlings, one in which necrosis was due to artificial infection ofTetratosphaeriaspp.and one in which necrosis was due to injuries superimposed on plants already showing reddening or chlorosis.Leaf refelctance profiles were acquired with the same ASD FieldSpec FR spectroradiometer (400-1000 nm) reported in Pietrzykowski et al.(2006).By spectral sensitivity analysis, they reported that 679-695 nm wavelengths were most sensitive to necrosis occurrence, while 706-726 nm ones were least sensitive.Anad-hocdeveloped VSI (i.e., ratio of reflectance at 688 and 721 nm), was capable of discriminating necrotic and non-necrotic treatments, while other important VSI could not.However, multivariate approaches, e.g., penalized discriminant analysis and Lasso, applied on the whole spectrum enabled much greater discrimination of all necrotic treatments compared to the VSI.Wavelengths between 679 and 695 nm were only sporadically selected as key wavelengths; therefore, results were not similar to the spectral sensitivity data.

    Metrosideros polymorpha Gaudich.

    Rapid ‘ōh(huán)i‘a(chǎn) Death

    Rapid ‘ōh(huán)i‘a(chǎn) Death (ROD) is a disease greatly killing a multitude ofMetrosideros polymorphatrees (‘ōh(huán)i ‘a(chǎn), Myrtaceae), a native keystone species in Hawaii.In 2010, many ‘ōh(huán)i’a trees began showing browning foliage, rapidly followed by tree mortality.In following years, the extent of ROD spread substantially.Two previously unknown fungal pathogens in the genusCeratocystiswere isolated as the cause of infection.First symptoms include leaf and/or branch yellowing, wilting, and then browning within several days to weeks.This browning stage is quite apparent some extent specific to ROD infection.Leaves finally drop from the infected tree over several weeks to months until only bare branches remain (Keith et al.2015).

    Two concomitant studies (Asner et al.2018; Vaughn et al.2018) were carried out to spectroscopically assess ROD.Vaughn et al.(2018) used the Carnegie Airborne Observatory (CAO) aircraft laboratory to acquire laserguided imaging spectroscopy data (350-2485 nm) and high-resolution imagery (24 cm) over more than 500,000 ha of Hawaii in summer 2017.They established a method to map individual tree crowns matching the symptoms of both active (brown; desiccated ‘ōh(huán)i’a crowns) and past (leafless tree crowns) ROD infection by the combination of two machine learning methods, SVM and a gradient boosting machine.Using a conservative classification design to curtail false-positives, model sensitivity levels were 87% and 83%, and precision levels 97% and 95% for browning and leafless crowns, respectively.Over the Hawaii Island, they recorded 43,134 individual crowns displaying the browning phase of ROD, while leafless crowns 547,666 in total and more scattered over the island.

    The associated study by Asner et al.(2018) used field collections of leaf spectral data with laboratory chemical investigations and airborne remote sensing to develop a spectral signature for ROD.Leaf spectra were collected with an ASD FieldSpec 3 spectroradiometer (350-2500 nm), while remote data were collected with the same CAO reported above.PLSR-models were developed to predict several leaf traits.Reflectance spectra predicted leaf mass per area, chlorophyll, phenol, carbon and water concentrations with high accuracies (R2: 0.85-0.99), while lower prediction accuracies were found for nitrogen, non-structural carbohydrates, cellulose and tannin concentrations.They reported that around 80% of ROD-infected plants showed distinct decreases in chlorophyll, moisture and non-structural carbohydrates, which overall led to significant and consistent variations in leaf reflectance over the VIS and SWIR regions.Leaf-level outcomes were scaled at the canopy level using airborne laser-guided imaging spectroscopy, with quantitative spectral separability of normal green-leaf canopies from suspected ROD-infected brown-leaf canopies in the VIS and SWIR spectrum.These achievements provided the spectral-chemical basis for detection, mapping and monitoring of the development of ROD in Hawaii.

    Ulmus americana L.

    Dutch elm disease

    Dutch elm disease (DED), caused by some ascomycete fungi of the genusOphiostoma(Ophiostoma ulmis.l.) affecting elm trees (Ulmusspp.), is one of the most destructive diseases in forest pathology.DED spread and infection of suitable hosts is due to a synchrony between life cycles of the host tree, the pathogen and its vectors, elm bark beetles.The first DED symptom is commonly an upper branch with leaves starting to wither and yellowing in summer, months before autumn leaf shedding.This gradually extents to the rest of the tree, with further dieback of branches (Santini and Faccoli 2015).

    Reflectance measurements of DED-infected trees were collected by Wilson et al.(1998) under laboratory conditions in the 350-2500 nm wavelength range using an ASD Field-Spec FR spectroradiometer.Different reflectance profiles were reported on leaf samples from trees showing variable DED infection (with early and late symptoms).Early disease stage resulted in a rapid increase in green and red regions and a small decrease in NIR (pigment variations), while the late disease stage induced a decrease in green and NIR and an increase in SWIR reflectance (water loss).These outcomes were useful for selection of suitable remote sensors to early detect and monitor DED in forests.

    Conclusions and future perspectives

    This review highlights the capability of vegetation spectroscopy in improving the assessment of forest disturbances induced by diseases in a timely and cost-effective manner.Using several optical sensors from leaf to landscape-level, a number of diseases characterized by variable pathogenic processes have been detected, identified and quantified in numerous sites worldwide.Overall, these studies have pointed to the green and red regions of the VIS, the red-edge and the early NIR as the spectral regions, i.e.those related to chlorophyll changes and symptom development, as being the most sensitive to disease development.Early disease stages usually produce a rapid increase in green and red reflectance and a small decrease in early NIR, whereas late disease stages result in a decrease of green and NIR regions and an increase in SWIR, mostly related to water content.

    This review also highlights some major issues to be addressed in the future.First, many major forest diseases and geographic areas remain unexplored by spectroscopic studies.Second, further development of high-resolution, cost-effective and portable spectral sensors are required for improving the assessment of forest diseases, since hyperspectral data have mostly shown their capability to assist with early detection of tree stress, e.g.slight decline in chlorophyll and leaf water content, providing forest managers an opportunity to perform efficient disease control.Third, since it has been demonstrated that leaf-level spectroscopy can be effectively used as ground reference or training input for airborne-based platforms and scaled to field and landscape levels, more studies should focus on linking spectral data collected at different scales.Fourth, since sensors mounted on airborne platforms make data acquisition an expensive process, and only a few satellite sensors are operational to date (their application has been restrained because of limited spatial coverage and high spectral noise), a potential solution is to assemble hyperspectral sensors on UAV which, although largely used in agriculture, have not been operated for forest disease assessment.One limitation, however, is the obligation to meet UAV regulations and policies that may vary considerably from region to region.Fifth, data archives, especially those of hyperspectral sensors, have to increase both temporal and spatial coverage since long-term monitoring is required to assess historical impacts of diseases on forests.Sixth, since the disturbances affecting the same forested regions may come from a range of sources, and it is also possible that one disturbance may influence forest responses to another, further studies should focus on differentiating between disease-caused forest damage and other types of damage.Since one major challenge is that single sensors are typically not suitable to complete this task, together with improvements in sensitivity and resolution of instrumentation, a likely solution may be the development of a multisensory approach.Seventh, accordingly to previous issues, algorithms for modelling spectra-disturbance relationships will need to be continually refined or redeveloped to take advantage of new data and novel landscape changes caused by pathogens.

    AcknowledgementsI gratefully acknowledge Cristina Nali and Giacomo Lorenzini for their expert advices and helpful suggestions.

    FundingOpen access funding provided by Università di Pisa within the CRUI-CARE Agreement.

    Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creat iveco mmons.org/l icen ses/ by/4.0/.

    国产有黄有色有爽视频| 欧美色视频一区免费| 19禁男女啪啪无遮挡网站| 一区二区日韩欧美中文字幕| 亚洲成人久久性| 国产亚洲av高清不卡| 搡老熟女国产l中国老女人| 国产在线精品亚洲第一网站| 国产激情久久老熟女| 久久精品国产99精品国产亚洲性色 | 国产亚洲精品久久久久久毛片| 制服人妻中文乱码| 91国产中文字幕| 高清黄色对白视频在线免费看| 在线天堂中文资源库| 日韩欧美一区二区三区在线观看| 巨乳人妻的诱惑在线观看| 一区二区三区激情视频| 亚洲av成人一区二区三| 中文字幕av电影在线播放| 69精品国产乱码久久久| 国内毛片毛片毛片毛片毛片| 久久人妻福利社区极品人妻图片| 韩国精品一区二区三区| 免费在线观看视频国产中文字幕亚洲| 热99re8久久精品国产| www.www免费av| 中亚洲国语对白在线视频| 狠狠狠狠99中文字幕| 露出奶头的视频| 日韩欧美三级三区| 亚洲自拍偷在线| 久久欧美精品欧美久久欧美| 亚洲精品中文字幕在线视频| 亚洲中文字幕日韩| av免费在线观看网站| 中文字幕精品免费在线观看视频| 亚洲欧美一区二区三区久久| 国产精品电影一区二区三区| 精品久久久久久久久久免费视频 | 搡老熟女国产l中国老女人| 高潮久久久久久久久久久不卡| 校园春色视频在线观看| 亚洲免费av在线视频| 国产极品粉嫩免费观看在线| 一区在线观看完整版| 国产一区二区三区在线臀色熟女 | 狂野欧美激情性xxxx| 国产亚洲欧美98| 亚洲人成伊人成综合网2020| av福利片在线| 欧美乱妇无乱码| 久久久久国内视频| 久久亚洲真实| 啪啪无遮挡十八禁网站| 日本免费a在线| 男女之事视频高清在线观看| av在线天堂中文字幕 | 男女之事视频高清在线观看| 叶爱在线成人免费视频播放| 亚洲中文av在线| 18禁国产床啪视频网站| 日韩欧美三级三区| 岛国视频午夜一区免费看| 99国产精品免费福利视频| 欧美日本中文国产一区发布| 老汉色av国产亚洲站长工具| 久久香蕉国产精品| 欧美激情极品国产一区二区三区| 欧美激情极品国产一区二区三区| 岛国视频午夜一区免费看| 久久中文看片网| 亚洲视频免费观看视频| 欧美老熟妇乱子伦牲交| 国产精品国产av在线观看| 欧美老熟妇乱子伦牲交| 午夜福利,免费看| 宅男免费午夜| 欧美黄色片欧美黄色片| 久久中文看片网| 99国产精品免费福利视频| 成人免费观看视频高清| 91精品三级在线观看| 久久国产精品男人的天堂亚洲| 99re在线观看精品视频| 亚洲国产精品sss在线观看 | 国产成年人精品一区二区 | 免费av中文字幕在线| 午夜成年电影在线免费观看| 免费在线观看黄色视频的| 18禁黄网站禁片午夜丰满| 国产成人一区二区三区免费视频网站| 操出白浆在线播放| 色婷婷av一区二区三区视频| 人人妻人人澡人人看| 一级作爱视频免费观看| 无人区码免费观看不卡| x7x7x7水蜜桃| 中文字幕色久视频| 国产精品久久电影中文字幕| 高清在线国产一区| 国产精品一区二区在线不卡| 精品午夜福利视频在线观看一区| 免费久久久久久久精品成人欧美视频| 国产亚洲欧美精品永久| 一级片'在线观看视频| 午夜成年电影在线免费观看| 一级毛片女人18水好多| 夜夜躁狠狠躁天天躁| 91大片在线观看| 在线观看免费高清a一片| 一本综合久久免费| av电影中文网址| 777久久人妻少妇嫩草av网站| 美女国产高潮福利片在线看| 多毛熟女@视频| 在线国产一区二区在线| 香蕉久久夜色| 淫妇啪啪啪对白视频| 日韩视频一区二区在线观看| 18禁国产床啪视频网站| 国产午夜精品久久久久久| 亚洲七黄色美女视频| 又大又爽又粗| 久久香蕉激情| 99久久精品国产亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看日本一区| 久久精品亚洲精品国产色婷小说| 老司机在亚洲福利影院| 亚洲 欧美 日韩 在线 免费| 成人国产一区最新在线观看| 亚洲国产精品sss在线观看 | 久久中文字幕人妻熟女| 人成视频在线观看免费观看| 亚洲五月天丁香| 成人黄色视频免费在线看| 午夜老司机福利片| e午夜精品久久久久久久| 成人亚洲精品av一区二区 | 欧美av亚洲av综合av国产av| 99re在线观看精品视频| 99精品欧美一区二区三区四区| 一二三四社区在线视频社区8| 精品久久久久久久毛片微露脸| 亚洲一码二码三码区别大吗| 满18在线观看网站| 在线观看午夜福利视频| 超色免费av| 日韩一卡2卡3卡4卡2021年| 熟女少妇亚洲综合色aaa.| 国产精品电影一区二区三区| 91九色精品人成在线观看| 不卡av一区二区三区| 动漫黄色视频在线观看| 亚洲美女黄片视频| 日本 av在线| 国产一卡二卡三卡精品| 91精品三级在线观看| 成在线人永久免费视频| 最近最新免费中文字幕在线| 99久久国产精品久久久| 级片在线观看| 在线观看舔阴道视频| 一二三四在线观看免费中文在| 亚洲av成人不卡在线观看播放网| 午夜a级毛片| 亚洲精品粉嫩美女一区| 这个男人来自地球电影免费观看| 日本撒尿小便嘘嘘汇集6| 精品久久久久久成人av| 国产精品电影一区二区三区| 色老头精品视频在线观看| 日本免费一区二区三区高清不卡 | 99久久久亚洲精品蜜臀av| 国产精品香港三级国产av潘金莲| 女性生殖器流出的白浆| 高清黄色对白视频在线免费看| 麻豆国产av国片精品| 香蕉久久夜色| 老熟妇乱子伦视频在线观看| 久久婷婷成人综合色麻豆| 国产又色又爽无遮挡免费看| 91精品三级在线观看| 一夜夜www| 后天国语完整版免费观看| 可以免费在线观看a视频的电影网站| 精品一区二区三区av网在线观看| 老司机在亚洲福利影院| 天天躁夜夜躁狠狠躁躁| 欧美在线黄色| 人人妻人人爽人人添夜夜欢视频| 国产有黄有色有爽视频| 亚洲人成电影观看| 久久这里只有精品19| 亚洲欧美精品综合久久99| 怎么达到女性高潮| 亚洲专区国产一区二区| 99国产极品粉嫩在线观看| 欧美黑人欧美精品刺激| 日日干狠狠操夜夜爽| 99精品久久久久人妻精品| 亚洲第一欧美日韩一区二区三区| 精品国产乱子伦一区二区三区| 女性被躁到高潮视频| 999久久久国产精品视频| 婷婷六月久久综合丁香| 国产三级黄色录像| 亚洲精华国产精华精| 免费不卡黄色视频| 欧美最黄视频在线播放免费 | 午夜精品国产一区二区电影| 巨乳人妻的诱惑在线观看| 午夜福利免费观看在线| 18美女黄网站色大片免费观看| 国产亚洲欧美98| 亚洲欧美日韩高清在线视频| 手机成人av网站| 村上凉子中文字幕在线| 亚洲精品美女久久av网站| 51午夜福利影视在线观看| 久99久视频精品免费| 国产高清videossex| 色婷婷久久久亚洲欧美| 12—13女人毛片做爰片一| 国产精品国产av在线观看| 叶爱在线成人免费视频播放| 成人黄色视频免费在线看| 成人亚洲精品av一区二区 | aaaaa片日本免费| 午夜免费激情av| 美女 人体艺术 gogo| 国产激情久久老熟女| 热99国产精品久久久久久7| 麻豆久久精品国产亚洲av | 国产亚洲欧美精品永久| 91九色精品人成在线观看| 国产一区二区三区视频了| 国产免费av片在线观看野外av| 极品教师在线免费播放| 男人操女人黄网站| 久久精品影院6| 久久精品国产亚洲av高清一级| 亚洲欧美激情综合另类| 我的亚洲天堂| 自线自在国产av| 亚洲激情在线av| 亚洲av第一区精品v没综合| 日本一区二区免费在线视频| 制服人妻中文乱码| 水蜜桃什么品种好| 久久国产精品男人的天堂亚洲| 露出奶头的视频| 精品熟女少妇八av免费久了| 多毛熟女@视频| 水蜜桃什么品种好| 99re在线观看精品视频| 国内毛片毛片毛片毛片毛片| 三级毛片av免费| 国产精品野战在线观看 | cao死你这个sao货| 首页视频小说图片口味搜索| 男男h啪啪无遮挡| 超碰97精品在线观看| 午夜免费观看网址| 亚洲aⅴ乱码一区二区在线播放 | 伦理电影免费视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线黄色| 久久人妻福利社区极品人妻图片| 免费在线观看亚洲国产| 久久青草综合色| 亚洲成人免费电影在线观看| 国产高清激情床上av| 亚洲七黄色美女视频| 操美女的视频在线观看| 看片在线看免费视频| 亚洲精品国产区一区二| 中文字幕人妻丝袜一区二区| 母亲3免费完整高清在线观看| 一级片'在线观看视频| 中文字幕精品免费在线观看视频| 老汉色∧v一级毛片| 美女高潮喷水抽搐中文字幕| 国产黄色免费在线视频| 国产高清videossex| 亚洲专区国产一区二区| 久热这里只有精品99| 这个男人来自地球电影免费观看| 精品国产一区二区三区四区第35| 亚洲视频免费观看视频| 国产成人精品久久二区二区91| 国产麻豆69| 真人做人爱边吃奶动态| 亚洲成人免费电影在线观看| 亚洲欧美一区二区三区久久| 亚洲精品成人av观看孕妇| 国产xxxxx性猛交| 99国产精品一区二区蜜桃av| 精品久久蜜臀av无| 成人永久免费在线观看视频| 国产三级在线视频| 国产有黄有色有爽视频| aaaaa片日本免费| 91九色精品人成在线观看| av中文乱码字幕在线| 丝袜美腿诱惑在线| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美网| 人人妻,人人澡人人爽秒播| 国产精品自产拍在线观看55亚洲| 亚洲欧美激情综合另类| 国产熟女xx| 波多野结衣av一区二区av| 日本免费a在线| av网站在线播放免费| 九色亚洲精品在线播放| 久久精品国产清高在天天线| 日韩av在线大香蕉| 多毛熟女@视频| 国产亚洲精品久久久久久毛片| 久热这里只有精品99| av网站免费在线观看视频| 国产精品永久免费网站| 亚洲三区欧美一区| 国产成年人精品一区二区 | 亚洲va日本ⅴa欧美va伊人久久| 午夜影院日韩av| 一夜夜www| www.www免费av| 熟女少妇亚洲综合色aaa.| 水蜜桃什么品种好| 精品国产美女av久久久久小说| 亚洲人成电影观看| 亚洲精品av麻豆狂野| 免费看a级黄色片| 满18在线观看网站| 熟女少妇亚洲综合色aaa.| 老司机在亚洲福利影院| av视频免费观看在线观看| 最近最新免费中文字幕在线| 黄片播放在线免费| 好男人电影高清在线观看| 亚洲精品中文字幕一二三四区| 日本三级黄在线观看| 欧美成人午夜精品| 国产一区二区三区综合在线观看| 久9热在线精品视频| 黄色a级毛片大全视频| 午夜精品国产一区二区电影| 国产精品日韩av在线免费观看 | 极品人妻少妇av视频| 成人精品一区二区免费| 悠悠久久av| 午夜福利在线免费观看网站| aaaaa片日本免费| 免费日韩欧美在线观看| av欧美777| 国产精品成人在线| 黄色怎么调成土黄色| 国产av在哪里看| 51午夜福利影视在线观看| 男女高潮啪啪啪动态图| 精品国产亚洲在线| 国产一区二区在线av高清观看| 日韩人妻精品一区2区三区| 我的亚洲天堂| 久久欧美精品欧美久久欧美| 满18在线观看网站| 中文字幕最新亚洲高清| 不卡一级毛片| 国产免费现黄频在线看| 成年女人毛片免费观看观看9| 变态另类成人亚洲欧美熟女 | 国产xxxxx性猛交| 丝袜人妻中文字幕| 久久国产亚洲av麻豆专区| 亚洲成人免费av在线播放| 69av精品久久久久久| 我的亚洲天堂| 国产三级在线视频| 桃红色精品国产亚洲av| 又黄又粗又硬又大视频| 日韩三级视频一区二区三区| 成人三级做爰电影| 免费在线观看影片大全网站| 久久中文字幕人妻熟女| 成人亚洲精品av一区二区 | 午夜福利在线观看吧| 精品国内亚洲2022精品成人| 欧美日本亚洲视频在线播放| 日韩人妻精品一区2区三区| 精品第一国产精品| 国产欧美日韩一区二区三区在线| 久久中文字幕人妻熟女| 午夜久久久在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲av电影在线进入| 国产精品 欧美亚洲| 精品国产亚洲在线| 嫩草影院精品99| 自拍欧美九色日韩亚洲蝌蚪91| 久热爱精品视频在线9| 成人免费观看视频高清| 欧美一级毛片孕妇| 亚洲精品成人av观看孕妇| 五月开心婷婷网| 亚洲一区高清亚洲精品| 国产精品亚洲一级av第二区| 亚洲中文av在线| 老司机福利观看| 热re99久久精品国产66热6| 一个人免费在线观看的高清视频| 男女高潮啪啪啪动态图| 亚洲 国产 在线| 51午夜福利影视在线观看| 级片在线观看| 涩涩av久久男人的天堂| 超色免费av| 美国免费a级毛片| 欧美老熟妇乱子伦牲交| 免费在线观看黄色视频的| 水蜜桃什么品种好| 国产精品一区二区三区四区久久 | 国产精品一区二区在线不卡| 大型黄色视频在线免费观看| 老鸭窝网址在线观看| 脱女人内裤的视频| 在线永久观看黄色视频| 国产成人一区二区三区免费视频网站| 757午夜福利合集在线观看| 美女高潮喷水抽搐中文字幕| 女人被狂操c到高潮| 激情视频va一区二区三区| 中文字幕人妻丝袜制服| 亚洲男人的天堂狠狠| 老汉色∧v一级毛片| 亚洲aⅴ乱码一区二区在线播放 | 欧美国产精品va在线观看不卡| 91精品三级在线观看| 欧美中文日本在线观看视频| 国产99白浆流出| 天堂中文最新版在线下载| xxx96com| 欧美av亚洲av综合av国产av| 99在线视频只有这里精品首页| 亚洲精品一区av在线观看| 色哟哟哟哟哟哟| 狠狠狠狠99中文字幕| 亚洲精品av麻豆狂野| 午夜影院日韩av| 91国产中文字幕| 成人永久免费在线观看视频| 91精品三级在线观看| 亚洲国产中文字幕在线视频| 欧美一区二区精品小视频在线| 18美女黄网站色大片免费观看| 亚洲av片天天在线观看| 亚洲av美国av| 大码成人一级视频| 女人被躁到高潮嗷嗷叫费观| 母亲3免费完整高清在线观看| 欧美乱码精品一区二区三区| 国产精品二区激情视频| 国产精品永久免费网站| 91成人精品电影| 一级作爱视频免费观看| 国产精品98久久久久久宅男小说| 叶爱在线成人免费视频播放| 午夜免费鲁丝| 久久精品aⅴ一区二区三区四区| 日本wwww免费看| 精品一区二区三区视频在线观看免费 | 99久久99久久久精品蜜桃| 男女做爰动态图高潮gif福利片 | 免费观看人在逋| 麻豆一二三区av精品| 黄片小视频在线播放| 欧美另类亚洲清纯唯美| 91成人精品电影| 精品国产一区二区久久| 亚洲成人免费电影在线观看| 欧美在线一区亚洲| 色综合婷婷激情| 亚洲精品美女久久久久99蜜臀| 美女扒开内裤让男人捅视频| 色在线成人网| 久9热在线精品视频| 国产亚洲精品久久久久久毛片| а√天堂www在线а√下载| 看黄色毛片网站| 亚洲男人天堂网一区| 伦理电影免费视频| 一a级毛片在线观看| 国产精品秋霞免费鲁丝片| 欧美日韩亚洲高清精品| 一区二区三区精品91| 午夜激情av网站| 免费观看人在逋| 麻豆一二三区av精品| 日韩三级视频一区二区三区| 91大片在线观看| 日韩欧美一区视频在线观看| 精品午夜福利视频在线观看一区| 亚洲人成网站在线播放欧美日韩| 日日干狠狠操夜夜爽| 中文字幕av电影在线播放| www.999成人在线观看| 国产片内射在线| 久久狼人影院| 久久久久国产一级毛片高清牌| 18禁观看日本| 欧美色视频一区免费| 国产精品亚洲一级av第二区| 免费搜索国产男女视频| 欧美精品啪啪一区二区三区| 51午夜福利影视在线观看| 香蕉久久夜色| 亚洲性夜色夜夜综合| 窝窝影院91人妻| 午夜免费观看网址| 校园春色视频在线观看| 国产精品久久久av美女十八| 欧美黄色淫秽网站| 色在线成人网| 性少妇av在线| 国产精品影院久久| 亚洲人成伊人成综合网2020| 少妇粗大呻吟视频| 久久久国产成人免费| 精品一区二区三区视频在线观看免费 | 国产午夜精品久久久久久| 午夜免费鲁丝| 天堂动漫精品| 91大片在线观看| 国产91精品成人一区二区三区| 久久久国产一区二区| 天天影视国产精品| 成熟少妇高潮喷水视频| avwww免费| 日韩欧美在线二视频| 日韩 欧美 亚洲 中文字幕| 久久久久国产一级毛片高清牌| 在线观看一区二区三区激情| 日本撒尿小便嘘嘘汇集6| av天堂在线播放| 亚洲人成伊人成综合网2020| 波多野结衣高清无吗| e午夜精品久久久久久久| 国产av一区在线观看免费| av视频免费观看在线观看| svipshipincom国产片| 日本 av在线| 国产精品永久免费网站| 日本 av在线| 国产精品永久免费网站| 国产av一区二区精品久久| 欧美丝袜亚洲另类 | 嫩草影院精品99| 免费观看精品视频网站| 亚洲激情在线av| 国产成人啪精品午夜网站| 日本欧美视频一区| 精品久久蜜臀av无| 搡老岳熟女国产| 99久久综合精品五月天人人| 欧美在线黄色| 多毛熟女@视频| 丰满的人妻完整版| 欧美大码av| 国产精品一区二区三区四区久久 | 国产又爽黄色视频| 日本撒尿小便嘘嘘汇集6| 午夜免费成人在线视频| 日本一区二区免费在线视频| 女人高潮潮喷娇喘18禁视频| 久久精品国产99精品国产亚洲性色 | 在线永久观看黄色视频| 国产乱人伦免费视频| av网站在线播放免费| 天堂√8在线中文| 成人国产一区最新在线观看| 亚洲情色 制服丝袜| 午夜精品在线福利| 亚洲欧美日韩无卡精品| 久久精品亚洲精品国产色婷小说| 脱女人内裤的视频| 成人影院久久| 黑人欧美特级aaaaaa片| 久久久精品国产亚洲av高清涩受| 免费搜索国产男女视频| x7x7x7水蜜桃| 精品第一国产精品| 色在线成人网| 操出白浆在线播放| 亚洲 欧美 日韩 在线 免费| 日韩欧美免费精品| 亚洲欧美日韩高清在线视频| 日韩成人在线观看一区二区三区| 在线观看免费日韩欧美大片| 亚洲一区二区三区不卡视频| 国产欧美日韩一区二区三| 成人18禁在线播放| 精品国产亚洲在线| 丝袜人妻中文字幕| 在线国产一区二区在线| 90打野战视频偷拍视频| 女人被狂操c到高潮| 国产在线观看jvid| 亚洲一区高清亚洲精品| 99久久久亚洲精品蜜臀av| 高清毛片免费观看视频网站 | 精品福利永久在线观看| 怎么达到女性高潮| 国产成年人精品一区二区 |