Remote sensing for the Spanish forests in the 21st century: a review of advances, needs, and opportunities

[EN] Forest ecosystems provide a host of services and societal benefits, including carbon storage, habitat for fauna, recreation, and provision of wood or non-wood products. In a context of complex demands on forest resources, identifying priorities for biodiversity and carbon budgets require accurate tools with sufficient temporal frequency. Moreover, understanding long term forest dynamics is necessary for sustainable planning and management. Remote sensing (RS) is a powerful means for analysis, synthesis, and report, providing insights and contributing to inform decisions upon forest ecosystems. In this communication we review current applications of RS techniques in Spanish forests, examining possible trends, needs, and opportunities offered by RS in a forestry context. Currently, wall-to-wall optical and LiDAR data are extensively used for a wide range of applications-many times in combination-whilst radar or hyperspectral data are rarely used in the analysis of Spanish forests. Unmanned Aerial Vehicles (UAVs) carrying visible and infrared sensors are gaining ground in acquisition of data locally and at small scale, particularly for health assessments. Forest fire identification and characterization are prevalent applications at the landscape scale, whereas structural assessments are the most widespread analyses carried out at limited extents. Unparalleled opportunities are offered by the availability of diverse RS data like those provided by the European Copernicus programme and recent satellite LiDAR launches, processing capacity, and synergies with other ancillary sources to produce information of our forests. Overall, we live in times of unprecedented opportunities for monitoring forest ecosystems with a growing support from RS technologies.Part of this work was funded by the Spanish Ministry of Science, innovation and University through the project AGL2016-76769-C2-1-R "Influence of natural disturbance regimes and management on forests dynamics. structure and carbon balance (FORESTCHANGE)".Gómez, C.; Alejandro, P.; Hermosilla, T.; Montes, F.; Pascual, C.; Ruiz Fernández, LÁ.; Álvarez-Taboada, F.... (2019). Remote sensing for the Spanish forests in the 21st century: a review of advances, needs, and opportunities. Forest Systems. 28(1):1-33. https://doi.org/10.5424/fs/2019281-14221S133281Ungar S, Pearlman J, Mendenhall J, Reuter D, 2003. Overview of the Earth Observing-1 (EO-1) mission. IEEE T Geosci Remote 41: 1149−1159.Valbuena R, Mauro F, Arjonilla FJ, Manzanera JA, 2011. Comparing Airborne Laser Scanning-Imagery Fusion Methods Based on Geometric Accuracy in Forested Areas. Remote Sens Environ 115(8): 1942-1956.Valbuena R, Mauro F, Rodríguez-Solano R, Manzanera JA, 2012. Partial Least Squares for Discriminating Variance Components in GNSS Accuracy Obtained Under Scots Pine Canopies. Forest Sci 58(2): 139-153.Valbuena R, De Blas A, Martín Fernández S, Maltamo M, Nabuurs GJ, Manzanera JA, 2013a. Within-Species Benefits of Back-projecting Laser Scanner and Multispectral Sensors in Monospecific P. sylvestris Forests. Eur J Remote Sens 46: 401-416.Valbuena R, Maltamo M, Martín-Fernández S, Packalen P, Pascual C, Nabuurs G-J, 2013b. Patterns of covariance between airborne laser scanning metrics and Lorenz curve descriptors of tree size inequality. Can J Remote Sens 39(1): 18-31.Valbuena R, Packalen P, García-Abril A, Mehtätalo L, Maltamo M, 2013c. Characterizing Forest Structural Types and Shelterwood Dynamics from Lorenz-based Indicators Predicted by Airborne Laser Scanning. Can J For Res 43: 1063-1074.Valbuena R, Maltamo M, Packalen P, 2016a. Classification of Multi-Layered Forest Development Classes from Low-Density National Airborne LiDAR Datasets. Forestry 89: 392-341.Valbuena R, Maltamo M, Packalen P, 2016b. Classification of Forest Development Stages from National Low-Density LiDAR Datasets: a Comparison of Machine Learning Methods. Revista de Teledetección 45: 15-25.Valbuena R, Hernando A, Manzanera JA, Martínez-Falero E, García-Abril A, Mola-Yudego B, 2017a. Most Similar Neighbour Imputation of Forest Attributes Using Metrics Derived from Combined Airborne LIDAR and Multispectral Sensors. Int J Digit Earth 11 (12): 1205-1218.Valbuena R, Hernando A, Manzanera JA, Görgens EB, Almeida DRA, Mauro F, García-Abril A, Coomes DA, 2017b. Enhancing of accuracy assessment for forest above-ground biomass estimates obtained from remote sensing via hypothesis testing and overfitting evaluation. Eco Mod 622: 15-26.Valbuena-Rabadán M, Santamaría-Pe-a J, Sanz-Adán F, 2016. Estimation of diameter and height of individual trees for Pinus sylvestris L. based on the individualising of crowns using airborne LiDAR and the National Forest Inventory data. For Sys 25(1): e046Varo-Martínez MA, Navarro-Cerrillo RM, Hernández-Clemente R, Duque-Lazo J, 2017. Semi-automated stand delineation in Mediterranean Pinus sylvestris plantations through segmentation of LiDAR data: The influence of pulse density. Int J Appl Earth Obs 56: 54-64.Vázquez de la Cueva A, 2008. Structural attributes of three forest types in central Spain and Landsat ETM+ information evaluated with redundancy analysis. Int J Remote Sens 29: 5657-5676.Verdú F, Salas J, 2010. Cartografía de áreas quemadas mediante análisis visual de imágenes de satélite en la Espa-a peninsular para el periodo 1991–2005. Geofocus 10: 54–81.Viana-Soto A, Aguado I, Martínez S, 2017. Assessment of post-fire vegetation recovery using fire severity and geographical data in the Mediterranean region (Spain). Environments 4: 90.Vicente-Serrano SG, Pérez-Cabello F, Lasanta T, 2011. Pinus halepensis regeneration after a wildfire in a semiarid environment: assessment using multitemporal Landsat images. Int J Wildland Fire 20Ñ 195-208.Viedma O, Quesada J, Torres I, De Santis A, Moreno JM, 2015. Fire severity in a large fire in a Pinus pinaster forest is highly predictable from burning conditions, stand structure, and topography. Ecosystems 18: 237-250.Yebra M, Chuvieco E, 2009. Generation of a species-specific look-up table for fuel moisture content assessment. IEEE J Selected topics in applied earth observation and RS 2 (1): 21-26.White JC, Wulder MA, Varhola A, Vastaranta M, Coops NC, Cook BD, Pitt D, Woods M, 2013. A best practices guide for generating forest inventory attributes from airborne laser scanning data using an area-based approach. Natural Resources Canada, Canadian Forest Service, Canadian Wood Fibre Centre, Victoria, BC. Information Report FI-X-010, 39 pp.White JC, Wulder MA, Hobart GW, Luther JE, Hermosilla T, Griffiths P, Coops NC, Hall RJ, Hostert P, Dyk A, Guindon L, 2014. Pixel-based image compositing for large-area dense time series applications and science. Can J Remote Sens 40 (3): 192-212.White JC, Coops NC, Wulder MA, Vastaranta M, Hilker T, Tompalski P, 2016. Remote sensing technologies for enhancing forest inventories: a review. Can J Remote Sens 42: 619-641.White JC, Wulder MA, Hermosilla T, Coops NC, Hobart GW, 2017. A nationwide characterization of 25 years of forest disturbance and recovery for Canada using Landsat time series. Remote Sens Environ 194: 303-321.Wulder MA, 1998. Optical remote-sensing techniques for the assessment of forest inventory and biophysical parameters. Progr Phys Geog 22 (4): 449-476.Wulder MA, Dymond CC, 2004. Remote sensing in survey of Mountain Pine impacts: review and recommendations. MPBI Report. Canadian Forest Service. Natural Resources Canada, Victoria, BC, Canada. 89 pp.Wulder MA, Masek JG, Cohen WB, Loveland TR, Woodcock CE, 2012. Opening the archive: how free data has enabled the science and monitoring promise of Landsat. Remote Sens Environ 122: 2-10.Wulder MA, Hilker T, White JC, Coops NC, Masek JG, Pflugmacher D, Crevier Y, 2015. Virtual constellations for global terrestrial monitoring. Remote Sens Environ 170: 62-76.Wulder MA, White JC, Loveland TR, Woodcock CE, Belward AS, Cohen WB, Fosnight EA, Shaw J, Masek JG, Roy DP, 2016. The global Landsat archive: Status, consolidation, and direction. Remote Sens Environ 185: 271-283.Xie Q, Zhu J, Wang Ch, Fu H, López-Sánchez JM, Ballester-Berman JD, 2017. A modified dual-baseline PolInSAR method for forest height estimation. Remote Sens-Basel 9 (8): 819.Xie Y, Sha Z, Yu M, 2008. Remote sensing imagery in vegetation mapping: a review. J Plant Ecol 1 (1): 9-23.Zald HSJ, Wulder MA, White JC, Hilker T, Hermosilla T, Hobart GW, Coops NC, 2016. Integrating Landsat pixel composites and change metrics with LiDAR plots to predictively map forest structure and aboveground biomass in Saskatchewan, Canada. Remote Sens Environ 176: 188-201.Zarco-Tejada PJ, Diaz-Varela R, Angileri V, Loudjani P, 2014. Tree height quantification using very high resolution imagery acquired from an unmanned aerial vehicle (UAV) and automatic 3D photo-reconstruction methods. Eur J Agron 55: 89-99.Zarco-Tejada PJ, Hornero A, Hernández-Clemente R, Beck PSA, 2018. Understanding the temporal dimension of the red-edge spectral region for forest decline detection using high-resolution hyperspectral and Sentinel-2A imagery. ISPRS J Photogramm 137: 134-148

On residual circulation of the Ria of Vigo, using a 3-D baroclinic model

The present paper studies circulation in the Ria of Vigo (northwest Spain) by means of a 3-D baroclinic model coupled with a Lagrangian model of dispersion. The model is examined to obtain the residual currents' underlying periodic tidal movements. These currents are related to the quality of water because they determine the residence time. The model is forced at the mouth of the estuary with the most important tidal harmonics, and at the innermost zone with the freshwater discharge due to the presence of the Oitaven River. Results show that the area under study has a two-layered circulation, where most of water enters the estuary through the southern mouth and leaves it through the northern one.Se estudia la circulación de la ría de Vigo (noroeste de España) por medio de un modelo tridimensional baroclínico acoplado a un modelo lagrangiano de dispersión. Los resultados provenientes del modelo se utilizan para la obtención de la corriente residual que subyace a los movimientos periódicos de marea. Esta corriente se relaciona con la calidad del agua, pues determina su tiempo de residencia. El modelo es forzado en la boca de la ría a través de las constituyentes armónicas de marea más importantes y con el agua dulce introducida en la parte interior del estuario por el río Oitaven. Los resultados obtenidos muestran una circulación en doble capa, además de una entrada principal del agua por la boca sur y una salida de ésta por la boca norte.Instituto Español de Oceanografí

¿Dónde hay mapaches y dónde intensificar su búsqueda en la Península Ibérica?

Las invasiones biológicas son unas de las principales amenazas para la conservación de la biodiversidad. A nivel global hay una preocupación creciente sobre las especies invasoras, tanto por los daños medioambientales que pueden producir, como por el impacto económico asociado a su gestión y, en algunos casos, el riesgo para la salud pública. El mapache (Procyon lotor) es un mesocarnívoro nativo de Norteamérica y Centroamérica con una gran plasticidad ecológica, lo que le ha permitido adaptarse sin dificultad a diferentes ambientes fuera de su área original. En la actualidad, el mapache ocupa diferentes zonas de Eurasia, siendo la Península Ibérica una de ellas. La presencia de mapaches en los ecosistemas ibéricos puede suponer una amenaza para la biodiversidad, de una dimensión aún desconocida. En este estudio se revisa y actualiza la distribución de la especie en Europa. A partir de la distribución obtenida, se pretende identificar aquellas zonas que suponen nuevas áreas para su establecimiento potencial en la Península Ibérica. Se realizaron modelos de distribución de especies, para lo cual se aplicó la función de favorabilidad a la distribución actualizada del mapache y a un conjunto de variables ambientales. Las variables antrópicas tuvieron un papel importante para explicar la distribución actual de esta especie invasora, lo cual indicaría que el patrón de invasión aún refleja los lugares de suelta. El modelo obtenido muestra zonas muy favorables para la especie que incluye los lugares de presencia estable y reproducción confirmada, en el centro y noroeste de España. Además, las áreas favorables identifican otras zonas con favorabilidad elevada y que aún no han sido ocupadas por el mapache. Estos territorios son aquellos donde el riesgo de invasión puede ser elevado en un futuro cercano, como es el caso del norte de la provincia de León. Esta aproximación permite identificar aquellas zonas en las que focalizar los esfuerzos de monitoreo.Universidad de Málaga. Campus de Excelencia Internacional Andalucía Tec

VAMOS: a Pathfinder for the HAWC Gamma-Ray Observatory

VAMOS was a prototype detector built in 2011 at an altitude of 4100m a.s.l. in the state of Puebla, Mexico. The aim of VAMOS was to finalize the design, construction techniques and data acquisition system of the HAWC observatory. HAWC is an air-shower array currently under construction at the same site of VAMOS with the purpose to study the TeV sky. The VAMOS setup included six water Cherenkov detectors and two different data acquisition systems. It was in operation between October 2011 and May 2012 with an average live time of 30%. Besides the scientific verification purposes, the eight months of data were used to obtain the results presented in this paper: the detector response to the Forbush decrease of March 2012, and the analysis of possible emission, at energies above 30 GeV, for long gamma-ray bursts GRB111016B and GRB120328B.Comment: Accepted for pubblication in Astroparticle Physics Journal (20 pages, 10 figures). Corresponding authors: A.Marinelli and D.Zaboro

Constraining the pˉ/p\bar{p}/p Ratio in TeV Cosmic Rays with Observations of the Moon Shadow by HAWC

An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov (HAWC) observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. We present our first upper limits on the $\bar{p}/p$ fraction, which in the absence of any direct measurements, provide the tightest available constraints of $\sim1\%$ on the antiproton fraction for energies between 1 and 10 TeV.Comment: 10 pages, 5 figures. Accepted by Physical Review

All-particle cosmic ray energy spectrum measured by the HAWC experiment from 10 to 500 TeV

We report on the measurement of the all-particle cosmic ray energy spectrum with the High Altitude Water Cherenkov (HAWC) Observatory in the energy range 10 to 500 TeV. HAWC is a ground based air-shower array deployed on the slopes of Volcan Sierra Negra in the state of Puebla, Mexico, and is sensitive to gamma rays and cosmic rays at TeV energies. The data used in this work were taken from 234 days between June 2016 to February 2017. The primary cosmic-ray energy is determined with a maximum likelihood approach using the particle density as a function of distance to the shower core. Introducing quality cuts to isolate events with shower cores landing on the array, the reconstructed energy distribution is unfolded iteratively. The measured all-particle spectrum is consistent with a broken power law with an index of $-2.49\pm0.01$ prior to a break at $(45.7\pm0.1$) TeV, followed by an index of $-2.71\pm0.01$. The spectrum also respresents a single measurement that spans the energy range between direct detection and ground based experiments. As a verification of the detector response, the energy scale and angular resolution are validated by observation of the cosmic ray Moon shadow's dependence on energy.Comment: 16 pages, 11 figures, 4 tables, submission to Physical Review

Observation of the Crab Nebula with the HAWC Gamma-Ray Observatory

The Crab Nebula is the brightest TeV gamma-ray source in the sky and has been used for the past 25 years as a reference source in TeV astronomy, for calibration and verification of new TeV instruments. The High Altitude Water Cherenkov Observatory (HAWC), completed in early 2015, has been used to observe the Crab Nebula at high significance across nearly the full spectrum of energies to which HAWC is sensitive. HAWC is unique for its wide field-of-view, nearly 2 sr at any instant, and its high-energy reach, up to 100 TeV. HAWC's sensitivity improves with the gamma-ray energy. Above $\sim$1 TeV the sensitivity is driven by the best background rejection and angular resolution ever achieved for a wide-field ground array. We present a time-integrated analysis of the Crab using 507 live days of HAWC data from 2014 November to 2016 June. The spectrum of the Crab is fit to a function of the form $\phi(E) = \phi_0 (E/E_{0})^{-\alpha -\beta\cdot{\rm{ln}}(E/E_{0})}$. The data is well-fit with values of $\alpha=2.63\pm0.03$, $\beta=0.15\pm0.03$, and log$_{10}(\phi_0~{\rm{cm}^2}~{\rm{s}}~{\rm{TeV}})=-12.60\pm0.02$ when $E_{0}$ is fixed at 7 TeV and the fit applies between 1 and 37 TeV. Study of the systematic errors in this HAWC measurement is discussed and estimated to be $\pm$50\% in the photon flux between 1 and 37 TeV. Confirmation of the Crab flux serves to establish the HAWC instrument's sensitivity for surveys of the sky. The HAWC survey will exceed sensitivity of current-generation observatories and open a new view of 2/3 of the sky above 10 TeV.Comment: Submitted 2017/01/06 to the Astrophysical Journa

Daily monitoring of TeV gamma-ray emission from Mrk 421, Mrk 501, and the Crab Nebula with HAWC

We present results from daily monitoring of gamma rays in the energy range $\sim0.5$ to $\sim100$ TeV with the first 17 months of data from the High Altitude Water Cherenkov (HAWC) Observatory. Its wide field of view of 2 steradians and duty cycle of $>95$% are unique features compared to other TeV observatories that allow us to observe every source that transits over HAWC for up to $\sim6$ hours each sidereal day. This regular sampling yields unprecedented light curves from unbiased measurements that are independent of seasons or weather conditions. For the Crab Nebula as a reference source we find no variability in the TeV band. Our main focus is the study of the TeV blazars Markarian (Mrk) 421 and Mrk 501. A spectral fit for Mrk 421 yields a power law index $\Gamma=2.21 \pm0.14_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$ and an exponential cut-off $E_0=5.4 \pm 1.1_{\mathrm{stat}}\pm 1.0_{\mathrm{sys}}$ TeV. For Mrk 501, we find an index $\Gamma=1.60\pm 0.30_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$ and exponential cut-off $E_0=5.7\pm 1.6_{\mathrm{stat}} \pm 1.0_{\mathrm{sys}}$ TeV. The light curves for both sources show clear variability and a Bayesian analysis is applied to identify changes between flux states. The highest per-transit fluxes observed from Mrk 421 exceed the Crab Nebula flux by a factor of approximately five. For Mrk 501, several transits show fluxes in excess of three times the Crab Nebula flux. In a comparison to lower energy gamma-ray and X-ray monitoring data with comparable sampling we cannot identify clear counterparts for the most significant flaring features observed by HAWC.Comment: 18 pages, 10 figures, accepted for publication in The Astrophysical Journa