This webmap is a subset of Global Landcover 1992 - 2020 Image Layer. You can access the source data from here. This layer is a time series of the annual ESA CCI (Climate Change Initiative) land cover maps of the world. ESA has produced land cover maps for the years 1992-2020. These are available at the European Space Agency Climate Change Initiative website.Time Extent: 1992-2020Cell Size: 300 meterSource Type: ThematicPixel Type: 8 Bit UnsignedData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: GlobalSource: ESA Climate Change InitiativeUpdate Cycle: Annual until 2020, no updates thereafterWhat can you do with this layer?This layer may be added to ArcGIS Online maps and applications and shown in a time series to watch a "time lapse" view of land cover change since 1992 for any part of the world. The same behavior exists when the layer is added to ArcGIS Pro.In addition to displaying all layers in a series, this layer may be queried so that only one year is displayed in a map. This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro with a query set to display just one year. Then, an area count of land cover types may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from other years to show a trend.To sum up area by land cover using this service, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Different Classifications Available to MapFive processing templates are included in this layer. The processing templates may be used to display a smaller set of land cover classes.Cartographic Renderer (Default Template)Displays all ESA CCI land cover classes.*Forested lands TemplateThe forested lands template shows only forested lands (classes 50-90).Urban Lands TemplateThe urban lands template shows only urban areas (class 190).Converted Lands TemplateThe converted lands template shows only urban lands and lands converted to agriculture (classes 10-40 and 190).Simplified RendererDisplays the map in ten simple classes which match the ten simplified classes used in 2050 Land Cover projections from Clark University.Any of these variables can be displayed or analyzed by selecting their processing template. In ArcGIS Online, select the Image Display Options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left hand menu. From the Processing Template pull down menu, select the variable to display.Using TimeBy default, the map will display as a time series animation, one year per frame. A time slider will appear when you add this layer to your map. To see the most current data, move the time slider until you see the most current year.In addition to displaying the past quarter century of land cover maps as an animation, this time series can also display just one year of data by use of a definition query. For a step by step example using ArcGIS Pro on how to display just one year of this layer, as well as to compare one year to another, see the blog called Calculating Impervious Surface Change.Hierarchical ClassificationLand cover types are defined using the land cover classification (LCCS) developed by the United Nations, FAO. It is designed to be as compatible as possible with other products, namely GLCC2000, GlobCover 2005 and 2009.This is a heirarchical classification system. For example, class 60 means "closed to open" canopy broadleaved deciduous tree cover. But in some places a more specific type of broadleaved deciduous tree cover may be available. In that case, a more specific code 61 or 62 may be used which specifies "open" (61) or "closed" (62) cover.Land Cover ProcessingTo provide consistency over time, these maps are produced from baseline land cover maps, and are revised for changes each year depending on the best available satellite data from each period in time. These revisions were made from AVHRR 1km time series from 1992 to 1999, SPOT-VGT time series between 1999 and 2013, and PROBA-V data for years 2013, 2014 and 2015. When MERIS FR or PROBA-V time series are available, changes detected at 1 km are re-mapped at 300 m. The last step consists in back- and up-dating the 10-year baseline LC map to produce the 24 annual LC maps from 1992 to 2015.Source dataThe datasets behind this layer were extracted from NetCDF files and TIFF files produced by ESA. Years 1992-2015 were acquired from ESA CCI LC version 2.0.7 in TIFF format, and years 2016-2018 were acquired from version 2.1.1 in NetCDF format. These are downloadable from ESA with an account, after agreeing to their terms of use. https://maps.elie.ucl.ac.be/CCI/viewer/download.phpCitationESA. Land Cover CCI Product User Guide Version 2. Tech. Rep. (2017). Available at: maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdfMore technical documentation on the source datasets is available here:https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=doc*Index of all classes in this layer:10 Cropland, rainfed11 Herbaceous cover12 Tree or shrub cover20 Cropland, irrigated or post-flooding30 Mosaic cropland (>50%) / natural vegetation (tree, shrub, herbaceous cover) (<50%)40 Mosaic natural vegetation (tree, shrub, herbaceous cover) (>50%) / cropland (<50%)50 Tree cover, broadleaved, evergreen, closed to open (>15%)60 Tree cover, broadleaved, deciduous, closed to open (>15%)61 Tree cover, broadleaved, deciduous, closed (>40%)62 Tree cover, broadleaved, deciduous, open (15-40%)70 Tree cover, needleleaved, evergreen, closed to open (>15%)71 Tree cover, needleleaved, evergreen, closed (>40%)72 Tree cover, needleleaved, evergreen, open (15-40%)80 Tree cover, needleleaved, deciduous, closed to open (>15%)81 Tree cover, needleleaved, deciduous, closed (>40%)82 Tree cover, needleleaved, deciduous, open (15-40%)90 Tree cover, mixed leaf type (broadleaved and needleleaved)100 Mosaic tree and shrub (>50%) / herbaceous cover (<50%)110 Mosaic herbaceous cover (>50%) / tree and shrub (<50%)120 Shrubland121 Shrubland evergreen122 Shrubland deciduous130 Grassland140 Lichens and mosses150 Sparse vegetation (tree, shrub, herbaceous cover) (<15%)151 Sparse tree (<15%)152 Sparse shrub (<15%)153 Sparse herbaceous cover (<15%)160 Tree cover, flooded, fresh or brakish water170 Tree cover, flooded, saline water180 Shrub or herbaceous cover, flooded, fresh/saline/brakish water190 Urban areas200 Bare areas201 Consolidated bare areas202 Unconsolidated bare areas210 Water bodies

This layer is a subset of Global Landcover 1992- 2020 Layer. This layer is a time series of the annual ESA CCI (Climate Change Initiative) land cover maps of the world. ESA has produced land cover maps for the years 1992-2020. These are available at the European Space Agency Climate Change Initiative website.Time Extent: 1992-2020Cell Size: 300 meterSource Type: ThematicPixel Type: 8 Bit UnsignedData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: GlobalSource: ESA Climate Change InitiativeUpdate Cycle: Annual until 2020, no updates thereafterWhat can you do with this layer?This layer may be added to ArcGIS Online maps and applications and shown in a time series to watch a "time lapse" view of land cover change since 1992 for any part of the world. The same behavior exists when the layer is added to ArcGIS Pro.In addition to displaying all layers in a series, this layer may be queried so that only one year is displayed in a map. This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro with a query set to display just one year. Then, an area count of land cover types may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from other years to show a trend.To sum up area by land cover using this service, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Different Classifications Available to MapFive processing templates are included in this layer. The processing templates may be used to display a smaller set of land cover classes.Cartographic Renderer (Default Template)Displays all ESA CCI land cover classes.*Forested lands TemplateThe forested lands template shows only forested lands (classes 50-90).Urban Lands TemplateThe urban lands template shows only urban areas (class 190).Converted Lands TemplateThe converted lands template shows only urban lands and lands converted to agriculture (classes 10-40 and 190).Simplified RendererDisplays the map in ten simple classes which match the ten simplified classes used in 2050 Land Cover projections from Clark University.Any of these variables can be displayed or analyzed by selecting their processing template. In ArcGIS Online, select the Image Display Options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left hand menu. From the Processing Template pull down menu, select the variable to display.Using TimeBy default, the map will display as a time series animation, one year per frame. A time slider will appear when you add this layer to your map. To see the most current data, move the time slider until you see the most current year.In addition to displaying the past quarter century of land cover maps as an animation, this time series can also display just one year of data by use of a definition query. For a step by step example using ArcGIS Pro on how to display just one year of this layer, as well as to compare one year to another, see the blog called Calculating Impervious Surface Change.Hierarchical ClassificationLand cover types are defined using the land cover classification (LCCS) developed by the United Nations, FAO. It is designed to be as compatible as possible with other products, namely GLCC2000, GlobCover 2005 and 2009.This is a heirarchical classification system. For example, class 60 means "closed to open" canopy broadleaved deciduous tree cover. But in some places a more specific type of broadleaved deciduous tree cover may be available. In that case, a more specific code 61 or 62 may be used which specifies "open" (61) or "closed" (62) cover.Land Cover ProcessingTo provide consistency over time, these maps are produced from baseline land cover maps, and are revised for changes each year depending on the best available satellite data from each period in time. These revisions were made from AVHRR 1km time series from 1992 to 1999, SPOT-VGT time series between 1999 and 2013, and PROBA-V data for years 2013, 2014 and 2015. When MERIS FR or PROBA-V time series are available, changes detected at 1 km are re-mapped at 300 m. The last step consists in back- and up-dating the 10-year baseline LC map to produce the 24 annual LC maps from 1992 to 2015.Source dataThe datasets behind this layer were extracted from NetCDF files and TIFF files produced by ESA. Years 1992-2015 were acquired from ESA CCI LC version 2.0.7 in TIFF format, and years 2016-2018 were acquired from version 2.1.1 in NetCDF format. These are downloadable from ESA with an account, after agreeing to their terms of use. https://maps.elie.ucl.ac.be/CCI/viewer/download.phpCitationESA. Land Cover CCI Product User Guide Version 2. Tech. Rep. (2017). Available at: maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdfMore technical documentation on the source datasets is available here:https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=doc*Index of all classes in this layer:10 Cropland, rainfed11 Herbaceous cover12 Tree or shrub cover20 Cropland, irrigated or post-flooding30 Mosaic cropland (>50%) / natural vegetation (tree, shrub, herbaceous cover) (<50%)40 Mosaic natural vegetation (tree, shrub, herbaceous cover) (>50%) / cropland (<50%)50 Tree cover, broadleaved, evergreen, closed to open (>15%)60 Tree cover, broadleaved, deciduous, closed to open (>15%)61 Tree cover, broadleaved, deciduous, closed (>40%)62 Tree cover, broadleaved, deciduous, open (15-40%)70 Tree cover, needleleaved, evergreen, closed to open (>15%)71 Tree cover, needleleaved, evergreen, closed (>40%)72 Tree cover, needleleaved, evergreen, open (15-40%)80 Tree cover, needleleaved, deciduous, closed to open (>15%)81 Tree cover, needleleaved, deciduous, closed (>40%)82 Tree cover, needleleaved, deciduous, open (15-40%)90 Tree cover, mixed leaf type (broadleaved and needleleaved)100 Mosaic tree and shrub (>50%) / herbaceous cover (<50%)110 Mosaic herbaceous cover (>50%) / tree and shrub (<50%)120 Shrubland121 Shrubland evergreen122 Shrubland deciduous130 Grassland140 Lichens and mosses150 Sparse vegetation (tree, shrub, herbaceous cover) (<15%)151 Sparse tree (<15%)152 Sparse shrub (<15%)153 Sparse herbaceous cover (<15%)160 Tree cover, flooded, fresh or brakish water170 Tree cover, flooded, saline water180 Shrub or herbaceous cover, flooded, fresh/saline/brakish water190 Urban areas200 Bare areas201 Consolidated bare areas202 Unconsolidated bare areas210 Water bodies

The global satellite remote sensing software market is experiencing robust growth, driven by increasing demand across diverse sectors. While precise figures for market size and CAGR aren't provided, considering the technological advancements and applications in agriculture (precision farming, crop monitoring), water conservancy (flood management, irrigation optimization), forest management (deforestation monitoring, resource assessment), and the public sector (urban planning, disaster response), a conservative estimate places the 2025 market size at approximately $2 billion. This figure reflects the substantial investments in satellite imagery acquisition and analysis capabilities worldwide. The market is further fueled by the rising adoption of cloud-based solutions, enhancing accessibility and scalability of software platforms. Trends such as the integration of AI and machine learning for automated image processing, the proliferation of high-resolution satellite imagery, and the increasing availability of open-source software are accelerating market expansion. However, factors such as the high cost of specialized software licenses and the need for skilled professionals to operate the sophisticated systems act as restraints. The market is segmented by application (agriculture, water conservancy, forest management, public sector, others) and software type (open-source, non-open-source). The North American and European markets currently hold significant shares, but the Asia-Pacific region is witnessing rapid growth due to increasing infrastructure development and government initiatives promoting geospatial technologies. This dynamic market landscape presents lucrative opportunities for both established players and emerging companies in the years to come. The forecast period (2025-2033) anticipates continued growth, with a projected CAGR of approximately 12%, driven by the aforementioned technological advancements and broadening applications across various industry verticals. The competitive landscape is comprised of both major players like ESRI, Trimble, and PCI Geomatica, offering comprehensive suites of software, and smaller, specialized companies focusing on niche applications or open-source solutions. The market is characterized by both proprietary and open-source software options. Open-source solutions like QGIS and GRASS GIS offer cost-effective alternatives, particularly for research and smaller organizations, while commercial solutions provide advanced functionalities and support. The increasing availability of cloud-based solutions is blurring the lines between these segments, with hybrid models emerging that combine the benefits of both. Future growth will be significantly influenced by collaborations between software providers and satellite imagery providers, fostering a more integrated ecosystem and streamlining the data acquisition and processing workflow. The market will continue to benefit from advancements in satellite technology, producing higher-resolution, more frequent, and more affordable imagery.

This layer is a subset of Global Landcover 1992- 2020 Layer focused on Fiji. This layer is a time series of the annual ESA CCI (Climate Change Initiative) land cover maps of the world. ESA has produced land cover maps for the years 1992-2020. These are available at the European Space Agency Climate Change Initiative website.Time Extent: 1992-2020Cell Size: 300 meterSource Type: ThematicPixel Type: 8 Bit UnsignedData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: GlobalSource: ESA Climate Change InitiativeUpdate Cycle: Annual until 2020, no updates thereafterWhat can you do with this layer?This layer may be added to ArcGIS Online maps and applications and shown in a time series to watch a "time lapse" view of land cover change since 1992 for any part of the world. The same behavior exists when the layer is added to ArcGIS Pro.In addition to displaying all layers in a series, this layer may be queried so that only one year is displayed in a map. This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro with a query set to display just one year. Then, an area count of land cover types may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from other years to show a trend.To sum up area by land cover using this service, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Different Classifications Available to MapFive processing templates are included in this layer. The processing templates may be used to display a smaller set of land cover classes.Cartographic Renderer (Default Template)Displays all ESA CCI land cover classes.*Forested lands TemplateThe forested lands template shows only forested lands (classes 50-90).Urban Lands TemplateThe urban lands template shows only urban areas (class 190).Converted Lands TemplateThe converted lands template shows only urban lands and lands converted to agriculture (classes 10-40 and 190).Simplified RendererDisplays the map in ten simple classes which match the ten simplified classes used in 2050 Land Cover projections from Clark University.Any of these variables can be displayed or analyzed by selecting their processing template. In ArcGIS Online, select the Image Display Options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left hand menu. From the Processing Template pull down menu, select the variable to display.Using TimeBy default, the map will display as a time series animation, one year per frame. A time slider will appear when you add this layer to your map. To see the most current data, move the time slider until you see the most current year.In addition to displaying the past quarter century of land cover maps as an animation, this time series can also display just one year of data by use of a definition query. For a step by step example using ArcGIS Pro on how to display just one year of this layer, as well as to compare one year to another, see the blog called Calculating Impervious Surface Change.Hierarchical ClassificationLand cover types are defined using the land cover classification (LCCS) developed by the United Nations, FAO. It is designed to be as compatible as possible with other products, namely GLCC2000, GlobCover 2005 and 2009.This is a heirarchical classification system. For example, class 60 means "closed to open" canopy broadleaved deciduous tree cover. But in some places a more specific type of broadleaved deciduous tree cover may be available. In that case, a more specific code 61 or 62 may be used which specifies "open" (61) or "closed" (62) cover.Land Cover ProcessingTo provide consistency over time, these maps are produced from baseline land cover maps, and are revised for changes each year depending on the best available satellite data from each period in time. These revisions were made from AVHRR 1km time series from 1992 to 1999, SPOT-VGT time series between 1999 and 2013, and PROBA-V data for years 2013, 2014 and 2015. When MERIS FR or PROBA-V time series are available, changes detected at 1 km are re-mapped at 300 m. The last step consists in back- and up-dating the 10-year baseline LC map to produce the 24 annual LC maps from 1992 to 2015.Source dataThe datasets behind this layer were extracted from NetCDF files and TIFF files produced by ESA. Years 1992-2015 were acquired from ESA CCI LC version 2.0.7 in TIFF format, and years 2016-2018 were acquired from version 2.1.1 in NetCDF format. These are downloadable from ESA with an account, after agreeing to their terms of use. https://maps.elie.ucl.ac.be/CCI/viewer/download.phpCitationESA. Land Cover CCI Product User Guide Version 2. Tech. Rep. (2017). Available at: maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdfMore technical documentation on the source datasets is available here:https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=doc*Index of all classes in this layer:10 Cropland, rainfed11 Herbaceous cover12 Tree or shrub cover20 Cropland, irrigated or post-flooding30 Mosaic cropland (>50%) / natural vegetation (tree, shrub, herbaceous cover) (<50%)40 Mosaic natural vegetation (tree, shrub, herbaceous cover) (>50%) / cropland (<50%)50 Tree cover, broadleaved, evergreen, closed to open (>15%)60 Tree cover, broadleaved, deciduous, closed to open (>15%)61 Tree cover, broadleaved, deciduous, closed (>40%)62 Tree cover, broadleaved, deciduous, open (15-40%)70 Tree cover, needleleaved, evergreen, closed to open (>15%)71 Tree cover, needleleaved, evergreen, closed (>40%)72 Tree cover, needleleaved, evergreen, open (15-40%)80 Tree cover, needleleaved, deciduous, closed to open (>15%)81 Tree cover, needleleaved, deciduous, closed (>40%)82 Tree cover, needleleaved, deciduous, open (15-40%)90 Tree cover, mixed leaf type (broadleaved and needleleaved)100 Mosaic tree and shrub (>50%) / herbaceous cover (<50%)110 Mosaic herbaceous cover (>50%) / tree and shrub (<50%)120 Shrubland121 Shrubland evergreen122 Shrubland deciduous130 Grassland140 Lichens and mosses150 Sparse vegetation (tree, shrub, herbaceous cover) (<15%)151 Sparse tree (<15%)152 Sparse shrub (<15%)153 Sparse herbaceous cover (<15%)160 Tree cover, flooded, fresh or brakish water170 Tree cover, flooded, saline water180 Shrub or herbaceous cover, flooded, fresh/saline/brakish water190 Urban areas200 Bare areas201 Consolidated bare areas202 Unconsolidated bare areas210 Water bodies

This layer is a time series of the annual ESA CCI (Climate Change Initiative) land cover maps of the world. ESA has produced land cover maps for the years since 1992. These are available at the European Space Agency Climate Change Initiative website.Time Extent: 1992-2019Cell Size: 300 meterSource Type: ThematicPixel Type: 8 Bit UnsignedData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: GlobalSource: ESA Climate Change InitiativeUpdate Cycle: AnnualWhat can you do with this layer?This layer may be added to ArcGIS Online maps and applications and shown in a time series to watch a "time lapse" view of land cover change since 1992 for any part of the world. The same behavior exists when the layer is added to ArcGIS Pro.In addition to displaying all layers in a series, this layer may be queried so that only one year is displayed in a map. This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro with a query set to display just one year. Then, an area count of land cover types may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from other years to show a trend.To sum up area by land cover using this service, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Different Classifications Available to MapFive processing templates are included in this layer. The processing templates may be used to display a smaller set of land cover classes.Cartographic Renderer (Default Template)Displays all ESA CCI land cover classes.*Forested lands TemplateThe forested lands template shows only forested lands (classes 50-90).Urban Lands TemplateThe urban lands template shows only urban areas (class 190).Converted Lands TemplateThe converted lands template shows only urban lands and lands converted to agriculture (classes 10-40 and 190).Simplified RendererDisplays the map in ten simple classes which match the ten simplified classes used in 2050 Land Cover projections from Clark University.Any of these variables can be displayed or analyzed by selecting their processing template. In ArcGIS Online, select the Image Display Options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left hand menu. From the Processing Template pull down menu, select the variable to display.Using TimeBy default, the map will display as a time series animation, one year per frame. A time slider will appear when you add this layer to your map. To see the most current data, move the time slider until you see the most current year.In addition to displaying the past quarter century of land cover maps as an animation, this time series can also display just one year of data by use of a definition query. For a step by step example using ArcGIS Pro on how to display just one year of this layer, as well as to compare one year to another, see the blog called Calculating Impervious Surface Change.Hierarchical ClassificationLand cover types are defined using the land cover classification (LCCS) developed by the United Nations, FAO. It is designed to be as compatible as possible with other products, namely GLCC2000, GlobCover 2005 and 2009.This is a heirarchical classification system. For example, class 60 means "closed to open" canopy broadleaved deciduous tree cover. But in some places a more specific type of broadleaved deciduous tree cover may be available. In that case, a more specific code 61 or 62 may be used which specifies "open" (61) or "closed" (62) cover.Land Cover ProcessingTo provide consistency over time, these maps are produced from baseline land cover maps, and are revised for changes each year depending on the best available satellite data from each period in time. These revisions were made from AVHRR 1km time series from 1992 to 1999, SPOT-VGT time series between 1999 and 2013, and PROBA-V data for years 2013, 2014 and 2015. When MERIS FR or PROBA-V time series are available, changes detected at 1 km are re-mapped at 300 m. The last step consists in back- and up-dating the 10-year baseline LC map to produce the 24 annual LC maps from 1992 to 2015.Source dataThe datasets behind this layer were extracted from NetCDF files and TIFF files produced by ESA. Years 1992-2015 were acquired from ESA CCI LC version 2.0.7 in TIFF format, and years 2016-2018 were acquired from version 2.1.1 in NetCDF format. These are downloadable from ESA with an account, after agreeing to their terms of use. https://maps.elie.ucl.ac.be/CCI/viewer/download.phpCitationESA. Land Cover CCI Product User Guide Version 2. Tech. Rep. (2017). Available at: maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdfMore technical documentation on the source datasets is available here:https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=doc*Index of all classes in this layer:10 Cropland, rainfed11 Herbaceous cover12 Tree or shrub cover20 Cropland, irrigated or post-flooding30 Mosaic cropland (>50%) / natural vegetation (tree, shrub, herbaceous cover) (<50%)40 Mosaic natural vegetation (tree, shrub, herbaceous cover) (>50%) / cropland (<50%)50 Tree cover, broadleaved, evergreen, closed to open (>15%)60 Tree cover, broadleaved, deciduous, closed to open (>15%)61 Tree cover, broadleaved, deciduous, closed (>40%)62 Tree cover, broadleaved, deciduous, open (15-40%)70 Tree cover, needleleaved, evergreen, closed to open (>15%)71 Tree cover, needleleaved, evergreen, closed (>40%)72 Tree cover, needleleaved, evergreen, open (15-40%)80 Tree cover, needleleaved, deciduous, closed to open (>15%)81 Tree cover, needleleaved, deciduous, closed (>40%)82 Tree cover, needleleaved, deciduous, open (15-40%)90 Tree cover, mixed leaf type (broadleaved and needleleaved)100 Mosaic tree and shrub (>50%) / herbaceous cover (<50%)110 Mosaic herbaceous cover (>50%) / tree and shrub (<50%)120 Shrubland121 Shrubland evergreen122 Shrubland deciduous130 Grassland140 Lichens and mosses150 Sparse vegetation (tree, shrub, herbaceous cover) (<15%)151 Sparse tree (<15%)152 Sparse shrub (<15%)153 Sparse herbaceous cover (<15%)160 Tree cover, flooded, fresh or brakish water170 Tree cover, flooded, saline water180 Shrub or herbaceous cover, flooded, fresh/saline/brakish water190 Urban areas200 Bare areas201 Consolidated bare areas202 Unconsolidated bare areas210 Water bodies

Precipitation is water released from clouds in the form of rain, sleet, snow, or hail. It is the primary source of recharge to the planet's fresh water supplies. This map contains a historical record showing the volume of precipitation that fell during each month from March 2000 to the present. Snow and hail are reported in terms of snow water equivalent - the amount of water that will be produced when they melt. Dataset SummaryThe GLDAS Precipitation layer is a time-enabled image service that shows average monthly precipitation from 2000 to the present, measured in millimeters. It is calculated by NASA using the Noah land surface model, run at 0.25 degree spatial resolution using satellite and ground-based observational data from the Global Land Data Assimilation System (GLDAS-2.1). The model is run with 3-hourly time steps and aggregated into monthly averages. Review the complete list of model inputs, explore the output data (in GRIB format), and see the full Hydrology Catalog for all related data and information!Phenomenon Mapped: PrecipitationUnits: MillimetersTime Interval: MonthlyTime Extent: 2000/01/01 to presentCell Size: 28 kmSource Type: ScientificPixel Type: Signed IntegerData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: Global Land SurfaceSource: NASAUpdate Cycle: SporadicWhat can you do with this layer?This layer is suitable for both visualization and analysis. It can be used in ArcGIS Online in web maps and applications and can be used in ArcGIS for Desktop. It is useful for scientific modeling, but only at global scales.By applying the "Calculate Anomaly" processing template, it is also possible to view these data in terms of deviation from the mean. Mean precipitation for a given month is calculated over the entire period of record - 2000 to present. Time: This is a time-enabled layer. By default, it will show the first month from the map's time extent. Or, if time animation is disabled, a time range can be set using the layer's multidimensional settings. If you wish to calculate the average, sum, or min/max over the time extent, change the mosaic operator used to resolve overlapping pixels. In ArcGIS Online, you do this in the "Image Display Order" tab. In ArcGIS Pro, use the "Data" ribbon. In ArcMap, it is in the 'Mosaic' tab of the layer properties window. If you do this, make sure to also select a specific variable. The minimum time extent is one month, and the maximum is 8 years. Variables: This layer has three variables: total precipitation, rainfall and snowfall. By default total is shown, but you can select a different variable using the multidimensional filter, or by applying the relevant raster function. Important: You must switch from the cartographic renderer to the analytic renderer in the processing template tab in the layer properties window before using this layer as an input to geoprocessing tool.

Calculating the total volume of water stored in a landscape can be challenging. In addition to lakes and reservoirs, water can be stored in soil, snowpack, or even inside plants and animals, and tracking the all these different mediums is not generally possible. However, calculating the change in storage is easy - just subtract the water output from the water input. Using the GLDAS layers we can do this calculation for every month from January 2000 to the present day. The precipitation layer tells us the input to each cell and runoff plus evapotranspiration is the output. When the input is higher than the output during a given month, it means water was stored. When output is higher than input, storage is being depleted. Generally the change in storage should be close to the change in soil moisture content plus the change in snowpack, but it will not match up exactly because of the other storage mediums discussed above.Dataset SummaryThe GLDAS Change in Storage layer is a time-enabled image service that shows net monthly change in storage from 2000 to the present, measured in millimeters of water. It is calculated by NASA using the Noah land surface model, run at 0.25 degree spatial resolution using satellite and ground-based observational data from the Global Land Data Assimilation System (GLDAS-2.1). The model is run with 3-hourly time steps and aggregated into monthly averages. Review the complete list of model inputs, explore the output data (in GRIB format), and see the full Hydrology Catalog for all related data and information!Phenomenon Mapped: Change in Water StorageUnits: MillimetersTime Interval: MonthlyTime Extent: 2000/01/01 to presentCell Size: 28 kmSource Type: ScientificPixel Type: Signed IntegerData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: Global Land SurfaceSource: NASAUpdate Cycle: SporadicWhat can you do with this layer?This layer is suitable for both visualization and analysis. It can be used in ArcGIS Online in web maps and applications and can be used in ArcGIS for Desktop. It is useful for scientific modeling, but only at global scales.In ArcGIS Pro you can use the built-in raster functions or create your own to create custom extracts of the data. Imagery layers provide fast, powerful inputs to geoprocessing tools, models, or Python scripts in Pro.Online you can filter the layer to show subsets of the data using the filter button and the layer's built-in raster functions.By applying the "Calculate Anomaly" raster function, it is possible to view these data in terms of deviation from the mean, instead of total change in storage. Mean change in storage for a given month is calculated over the entire period of record - 2000 to present.Time: This is a time-enabled layer. By default, it will show the first month from the map's time extent. Or, if time animation is disabled, a time range can be set using the layer's multidimensional settings. If you wish to calculate the average, sum, or min/max change in storage over the time extent, change the mosaic operator used to resolve overlapping pixels. In ArcGIS Online, you do this in the "Image Display Order" tab. In ArcGIS Pro, use the "Data" ribbon. In ArcMap, it is in the 'Mosaic' tab of the layer properties window. The minimum time extent is one month, and the maximum is 8 years. Important: You must switch from the cartographic renderer to the analytic renderer in the processing template tab in the layer properties window before using this layer as an input to geoprocessing tools.

Hurricane Maria brought intense rainfall and caused widespread landsliding throughout Puerto Rico during September 2017. Previous detailed landslide inventories following the hurricane include Bessette-Kirton et al. (2017, 2019). Here we continue that work with an in-depth look at two areas in San Lorenzo, which is a municipality in the east-central part of the main island. To study a characteristic sample of landslides in San Lorenzo, we mapped all visible landslides in two physiographically diverse areas, but all within the San Lorenzo Formation. We used aerial imagery collected between 9-15 October 2017 (Quantum Spatial, Inc., 2017) to map landslide source and runout areas, and 1-m-resolution pre-event and post-event lidar (U.S. Geological Survey, 2018, 2020) as a digital base map for mapping. Difficulties with using these tools arose when aerial imagery was not correctly georeferenced to the lidar, when cloud cover was present in all images of an area, and in interpreting failure modes using only two-dimensional aerial photos. These difficulties with aerial imagery were partially resolved using the lidar. The map data comprises headscarp points, travel distance lines, source area polygons, and affected area polygons that are provided as point, line, and polygon shapefiles that may be viewed using common geographic information systems. Various characteristics of the landslides and their geomorphic settings are included in attribute tables of the mapped features, and this information is described in the "Attribute Summary" document in the accompanying files. Quantitative attributes (e.g., failure travel distance, failure fall height, watershed contributing area, etc.) were determined using tools available with the ESRI ArcGIS Pro v. 3.0.36056 geographic information system. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. References Bessette-Kirton, E.K., Cerovski-Darrian, C., Schulz, W.H., Coe, J.A., Kean, J.W., Godt, J.W., Thomas, M.A. and Hughes, K.S., 2019, Landslides triggered by Hurricane Maria: Assessment of an extreme event in Puerto Rico: GSA Today, v. 29, no. 6. Bessette-Kirton, E.K., Coe, J.A., Godt, J.W., Kean, J.W., Rengers, F.K., Schulz, W.H., Baum, R.L., Jones, E.S., and Staley, D.M., 2017, Map data showing concentration of landslides caused by Hurricane Maria in Puerto Rico: U.S. Geological Survey data release, https://doi.org/10.5066/F7JD4VRF. Quantum Spatial, Inc., 2017 FEMA PR Imagery: https://s3amazonaws.com/fema-cap-imagery/Others/Maria (accessed October 2017). U.S. Geological Survey, 2018, USGS NED Original Product Resolution PR Puerto Rico 2015: http://nationalmap.gov/elevation.html (accessed October 2018). U.S. Geological Survey, 2020, USGS NED Original Product Resolution PR Puerto Rico 2015: http://nationalmap.gov/elevation.html (accessed October 2018).

Soils and soil moisture greatly influence the water cycle and have impacts on runoff, flooding and agriculture. Soil type and soil particle composition (sand, clay, silt) affect soil moisture and the ability of the soil to retain water. Soil moisture is also affected by levels of evaporation and plant transpiration, potentially leading to near dryness and eventual drought.Measuring and monitoring soil moisture can ensure the fitness of your crops and help predict or prepare for flash floods and drought. The GLDAS soil moisture data is useful for modeling these scenarios and others, but only at global scales. Dataset SummaryThe GLDAS Soil Moisture layer is a time-enabled image service that shows average monthly soil moisture from 2000 to the present, measured as the millimeters of water contained within four different depth levels. It is calculated by NASA using the Noah land surface model, run at 0.25 degree spatial resolution using satellite and ground-based observational data from the Global Land Data Assimilation System (GLDAS-2.1). The model is run with 3-hourly time steps and aggregated into monthly averages. Review the complete list of model inputs, explore the output data (in GRIB format), and see the full Hydrology Catalog for all related data and information!Phenomenon Mapped: Soil MoistureUnits: MillimetersTime Interval: MonthlyTime Extent: 2000/01/01 to presentCell Size: 28 kmSource Type: ScientificPixel Type: Signed IntegerData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: Global Land SurfaceSource: NASAUpdate Cycle: SporadicWhat can you do with this layer?This layer is suitable for both visualization and analysis. It can be used in ArcGIS Online in web maps and applications and can be used in ArcGIS Desktop. The GLDAS soil moisture data is useful for modeling, but only at global scales. By applying the "Calculate Anomaly" processing template, it is also possible to view these data in terms of deviation from the mean. Mean soil moisture for a given month is calculated over the entire period of record - 2000 to present.Time: This is a time-enabled layer. By default, it will show the first month from the map's time extent. Or, if time animation is disabled, a time range can be set using the layer's multidimensional settings. If you wish to calculate the average, sum, or min/max over the time extent, change the mosaic operator used to resolve overlapping pixels. In ArcGIS Online, you do this in the "Image Display Order" tab. In ArcGIS Pro, use the "Data" ribbon. In ArcMap, it is in the 'Mosaic' tab of the layer properties window. If you do this, make sure to also select a specific variable. The minimum time extent is one month, and the maximum is 8 years. Variables: This layer has five variables, corresponding to different depth levels. By default total is shown, but you can view an individual depth level using the multidimensional filter, or by applying the relevant raster function. Important: You must switch from the cartographic renderer to the analytic renderer in the processing template tab in the layer properties window before using this layer as an input to geoprocessing tools.

When precipitation falls on the surface of the Earth, much of it is captured in storage (e.g. lakes, aquifers, soil moisture, snowpack, and vegetation). Precipitation that exceeds the storage capacity of the landscape becomes runoff, which flows into river systems. Overland flow is the most visible form of runoff, causing erosion and flash floods, but subsurface flow is the larger contributor in many watersheds. Subsurface flow can emerge on the surface through springs, or more commonly, seep into rivers and lakes through their banks. In urban areas, impervious land cover drastically increases the amount of surface runoff generated, which sweeps trash and urban debris into waterways and increases the likelihood and severity of flash floods. In agricultural areas, surface or subsurface runoff can carry excess salts and nutrients, especially nitrogen and phosphorus. This map contains a historical record showing the amount of runoff generated each month from March 200 to present. It is reported in millimeters, so multiply by a surface area to calculate the total volume of runoff.Dataset SummaryThe GLDAS Runoff layer is a time-enabled image service that shows average monthly runoff from 2000 to the present measured in millimeters. It is calculated by NASA using the Noah land surface model, run at 0.25 degree spatial resolution using satellite and ground-based observational data from the Global Land Data Assimilation System (GLDAS-2.1). The model is run with 3-hourly time steps and aggregated into monthly averages. Review the complete list of model inputs, explore the output data (in GRIB format), and see the full Hydrology Catalog for all related data and information!Phenomenon Mapped: RunoffUnits: MillimetersTime Interval: MonthlyTime Extent: 2000/01/01 to presentCell Size: 28 kmSource Type: ScientificPixel Type: Signed IntegerData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: Global Land SurfaceSource: NASAUpdate Cycle: SporadicWhat can you do with this layer?This layer is suitable for both visualization and analysis. It can be used in ArcGIS Online in web maps and applications and can be used in ArcGIS Desktop. t is useful for scientific modeling, but only at global scales.By applying the "Calculate Anomaly" processing template, it is also possible to view these data in terms of deviation from the mean. Mean runoff for a given month is calculated over the entire period of record - 2000 to present.Time: This is a time-enabled layer. By default, it will show the first month from the map's time extent. Or, if time animation is disabled, a time range can be set using the layer's multidimensional settings. If you wish to calculate the average, sum, or min/max over the time extent, change the mosaic operator used to resolve overlapping pixels. In ArcGIS Online, you do this in the "Image Display Order" tab. In ArcGIS Pro, use the "Data" ribbon. In ArcMap, it is in the 'Mosaic' tab of the layer properties window. If you do this, make sure to also select a specific variable. The minimum time extent is one month, and the maximum is 8 years. Variables: This layer has three variables: total runoff, surface flow and subsurface flow. By default total is shown, but you can select a different variable using the multidimensional filter, or by applying the relevant raster function. Important: You must switch from the cartographic renderer to the analytic renderer in the processing template tab in the layer properties window before using this layer as an input to geoprocessing tools.

Ground control points used to ensure Structure from Motion (SfM) terrain models were georectified (Westoby et al. 2012; Wolf 2021) using Emlid R2 RTK (real-time kinematic) GNSS (global navigation satellite system) system consisting of a base station set up over an established known point (established with Canadian Geodetic Survey of Natural Resources Canada (NRCAN) service Canadian Spatial Reference System Precise Point Positioning (CSRS- PPP)) and a rover.

Once the ground control points were surveyed, aerial drone images were acquired. We created flight polygons in Drone Deploy. Pictures were captured with a DJI Mavic II drone with minimum 80 % overlap of photos. Drone Deploy was chosen because it has an option to account for the doming error commonly found in models created from drone imagery and structure from motion (SfM). The doming effect is a systematic error that impacts the DEMs vertical component and can provide errors larger than the usual centimeter level (Sanz-Ablanedo et al. 2020). Generally, each site was flown once in fall of 2020 and once in spring of 2021.

We created orthorectified images and digital terrain models using Agisoft Metashape, a photogrammetric processing software application that uses SfM. We followed the workflow outlined in Bywater-Reyes and Pratt-Sitaula (2022). Once processed, orthorectified imagery and Digital Elevation Models (DEMs) were exported to ArcGIS Pro for additional analysis. Data collection metadata and postprocessing outcomes can be found in this Zenodo repository.

Melting snowpack is a key part of the spring water budget in many parts of the world. Like a natural reservoir, snowpack stores winter precipitation and releases it as runoff over the course of many months. Where summer rains are scarce snowpack provides crucial base flow without which rivers might go dry. Where summer rains are torrential, this exacerbates the flooding and can lead to the loss of lives. This map contains a historical record showing the water stored in snowpack during each month from March 2000 to the present. It is not a map of snow depth, but of snow water equivalent, which is the amount of water that would be produced if all the snow melted. For fresh snow, this can be anywhere from 5% to 20% the depth of the snow, depending on temperature (snow tends to be fluffier at lower temperatures). As the snow settles and melts, it becomes more dense, up to 40% or 50% in the spring. Temperature, albedo (the reflective property of the snow), density, and volume all affect the melting rate of the snowpack. Additionally, melting rate is influenced by wind, relative humidity, air temperature and solar radiation.Dataset SummaryThe GLDAS Snowpack layer is a time-enabled image service that shows average monthly snowpack from 2000 to present, measured in millimeters of snow water equivalent. It is calculated by NASA using the Noah land surface model, run at 0.25 degree spatial resolution using satellite and ground-based observational data from the Global Land Data Assimilation System (GLDAS-2.1). The model is run with 3-hourly time steps and aggregated into monthly averages. Review the complete list of model inputs, explore the output data (in GRIB format), and see the full Hydrology Catalog for all related data and information!Phenomenon Mapped: SnowpackUnits: MillimetersTime Interval: MonthlyTime Extent: 2000/01/01 to presentCell Size: 28 kmSource Type: ScientificPixel Type: Signed IntegerData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: Global Land SurfaceSource: NASAUpdate Cycle: SporadicWhat can you do with this layer?This layer is suitable for both visualization and analysis. It can be used in ArcGIS Online in web maps and applications and can be used in ArcGIS Desktop. Is useful for scientific modeling, but only at global scales. The GLDAS snowpack data is useful for modeling, but only at global scales. By applying the "Calculate Anomaly" processing template, it is also possible to view these data in terms of deviation from the mean, instead of total snowpack. Mean snowpack for a given month is calculated over the entire period of record - 2000 to present.Time: This is a time-enabled layer. By default, it will show the first month from the map's time extent. Or, if time animation is disabled, a time range can be set using the layer's multidimensional settings. If you wish to calculate the average, sum, or min/max over the time extent, change the mosaic operator used to resolve overlapping pixels. In ArcGIS Online, you do this in the "Image Display Order" tab. In ArcGIS Pro, use the "Data" ribbon. In ArcMap, it is in the 'Mosaic' tab of the layer properties window. The minimum time extent is one month, and the maximum is 8 years. Important: You must switch from the cartographic renderer to the analytic renderer in the processing template tab in the layer properties window before using this layer as an input to geoprocessing tools.

Most of us understand the hydrologic cycle in terms of the visible paths that water can take such as rainstorms, rivers, waterfalls and lakes. However, an even larger volume of water flows through the air all around us in two invisible paths: evaporation and transpiration. These two paths together are referred to as evapotranpsiration (ET), and claim 61% of all terrestrial precipitation. Solar radiation, air temperature, wind speed, soil moisture, and land cover all affect the rate of evapotranspiration, which is a major driver of the global water cycle, and key component of most catchments' water budget. This map contains a historical record showing the volume of water lost to evapotranspiration globally during each month from March 2000 to the present.Dataset SummaryThe GLDAS Evapotranspiration layer is a time-enabled image service that shows total actual evapotranspiration monthly from 2000 to the present, measured in millimeters of water loss. It is calculated by NASA using the Noah land surface model, run at 0.25 degree spatial resolution using satellite and ground-based observational data from the Global Land Data Assimilation System (GLDAS-2.1). The model is run with 3-hourly time steps and aggregated into monthly averages. Review the complete list of model inputs, explore the output data (in GRIB format), and see the full Hydrology Catalog for all related data and information!Phenomenon Mapped: EvapotranspirationUnits: MillimetersTime Interval: MonthlyTime Extent: 2000/01/01 to presentCell Size: 28 kmSource Type: ScientificPixel Type: Signed IntegerData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: Global Land SurfaceSource: NASAUpdate Cycle: SporadicWhat can you do with this layer?This layer is suitable for both visualization and analysis. It can be used in ArcGIS Online in web maps and applications and can be used in ArcGIS for Desktop. It is useful for scientific modeling, but only at global scales. By applying the "Calculate Anomaly" processing template, it is also possible to view these data in terms of deviation from the mean. Mean evapotranspiration for a given month is calculated over the entire period of record - 2000 to present.Time: This is a time-enabled layer. By default, it will show the first month from the map's time extent. Or, if time animation is disabled, a time range can be set using the layer's multidimensional settings. If you wish to calculate the average, sum, or min/max over the time extent, change the mosaic operator used to resolve overlapping pixels. In ArcGIS Online, you do this in the "Image Display Order" tab. In ArcGIS Pro, use the "Data" ribbon. In ArcMap, it is in the 'Mosaic' tab of the layer properties window. The minimum time extent is one month, and the maximum is 8 years. Important: You must switch from the cartographic renderer to the analytic renderer in the processing template tab in the layer properties window before using this layer as an input to geoprocessing tools.