MSZSI: Multi-Scale Zonal Statistics [AgriClimate] Inventory -------------------------------------------------------------------------------------- MSZSI is a data extraction tool for Google Earth Engine that aggregates time-series remote sensing information to multiple administrative levels using the FAO GAUL data layers. The code at the bottom of this page (metadata) can be pasted into the Google Earth Engine JavaScript code editor and ran at https://code.earthengine.google.com/. Please refer to the associated publication: Peter, B.G., Messina, J.P., Breeze, V., Fung, C.Y., Kapoor, A. and Fan, P., 2024. Perspectives on modifiable spatiotemporal unit problems in remote sensing of agriculture: evaluating rice production in Vietnam and tools for analysis. Frontiers in Remote Sensing, 5, p.1042624. https://www.frontiersin.org/journals/remote-sensing/articles/10.3389/frsen.2024.1042624 Input options: [1] Country of interest [2] Start and end year [3] Start and end month [4] Option to mask data to a specific land-use/land-cover type [5] Land-use/land-cover type code from CGLS LULC [6] Image collection for data aggregation [7] Desired band from the image collection [8] Statistics type for the zonal aggregations [9] Statistic to use for annual aggregation [10] Scaling options [11] Export folder and label suffix Output: Two CSVs containing zonal statistics for each of the FAO GAUL administrative level boundaries Output fields: system:index, 0-ADM0_CODE, 0-ADM0_NAME, 0-ADM1_CODE, 0-ADM1_NAME, 0-ADMN_CODE, 0-ADMN_NAME, 1-AREA_PERCENT_LULC, 1-AREA_SQM_LULC, 1-AREA_SQM_ZONE, 2-X_2001, 2-X_2002, 2-X_2003, ..., 2-X_2020, .geo PREPROCESSED DATA DOWNLOAD The datasets available for download contain zonal statistics at 2 administrative levels (FAO GAUL levels 1 and 2). Select countries from Southeast Asia and Sub-Saharan Africa (Cambodia, Indonesia, Lao PDR, Myanmar, Philippines, Thailand, Vietnam, Burundi, Kenya, Malawi, Mozambique, Rwanda, Tanzania, Uganda, Zambia, Zimbabwe) are included in the current version, with plans to extend the dataset to contain global metrics. Each zip file is described below and two example NDVI tables are available for preview. Key: [source, data, units, temporal range, aggregation, masking, zonal statistic, notes] Currently available: MSZSI-V2_V-NDVI-MEAN.tar: [NASA-MODIS, NDVI, index, 2001–2020, annual mean, agriculture, mean, n/a] MSZSI-V2_T-LST-DAY-MEAN.tar: [NASA-MODIS, LST Day, °C, 2001–2020, annual mean, agriculture, mean, n/a] MSZSI-V2_T-LST-NIGHT-MEAN.tar: [NASA-MODIS, LST Night, °C, 2001–2020, annual mean, agriculture, mean, n/a] MSZSI-V2_R-PRECIP-SUM.tar: [UCSB-CHG-CHIRPS, Precipitation, mm, 2001–2020, annual sum, agriculture, mean, n/a] MSZSI-V2_S-BDENS-MEAN.tar: [OpenLandMap, Bulk density, g/cm3, static, n/a, agriculture, mean, at depths 0-10-30-60-100-200] MSZSI-V2_S-ORGC-MEAN.tar: [OpenLandMap, Organic carbon, g/kg, static, n/a, agriculture, mean, at depths 0-10-30-60-100-200] MSZSI-V2_S-PH-MEAN.tar: [OpenLandMap, pH in H2O, pH, static, n/a, agriculture, mean, at depths 0-10-30-60-100-200] MSZSI-V2_S-WATER-MEAN.tar: [OpenLandMap, Soil water, % at 33kPa, static, n/a, agriculture, mean, at depths 0-10-30-60-100-200] MSZSI-V2_S-SAND-MEAN.tar: [OpenLandMap, Sand, %, static, n/a, agriculture, mean, at depths 0-10-30-60-100-200] MSZSI-V2_S-SILT-MEAN.tar: [OpenLandMap, Silt, %, static, n/a, agriculture, mean, at depths 0-10-30-60-100-200] MSZSI-V2_S-CLAY-MEAN.tar: [OpenLandMap, Clay, %, static, n/a, agriculture, mean, at depths 0-10-30-60-100-200] MSZSI-V2_E-ELEV-MEAN.tar: [MERIT, [elevation, slope, flowacc, HAND], [m, degrees, km2, m], static, n/a, agriculture, mean, n/a] Coming soon MSZSI-V2_C-STAX-MEAN.tar: [OpenLandMap, Soil taxonomy, category, static, n/a, agriculture, area sum, n/a] MSZSI-V2_C-LULC-MEAN.tar: [CGLS-LC100-V3, LULC, category, 2015–2019, mode, none, area sum, n/a] Data sources: https://developers.google.com/earth-engine/datasets/catalog/MODIS_006_MOD13Q1 https://developers.google.com/earth-engine/datasets/catalog/MODIS_006_MOD11A2 https://developers.google.com/earth-engine/datasets/catalog/UCSB-CHG_CHIRPS_PENTAD https://developers.google.com/earth-engine/datasets/catalog/OpenLandMap_SOL_SOL_BULKDENS-FINEEARTH_USDA-4A1H_M_v02 https://developers.google.com/earth-engine/datasets/catalog/OpenLandMap_SOL_SOL_ORGANIC-CARBON_USDA-6A1C_M_v02 https://developers.google.com/earth-engine/datasets/catalog/OpenLandMap_SOL_SOL_PH-H2O_USDA-4C1A2A_M_v02 https://developers.google.com/earth-engine/datasets/catalog/OpenLandMap_SOL_SOL_WATERCONTENT-33KPA_USDA-4B1C_M_v01 https://developers.google.com/earth-engine/datasets/catalog/OpenLandMap_SOL_SOL_CLAY-WFRACTION_USDA-3A1A1A_M_v02 https://developers.google.com/earth-engine/datasets/catalog/OpenLandMap_SOL_SOL_SAND-WFRACTION_USDA-3A1A1A_M_v02 https://developers.google.com/earth-engine/datasets/catalog/OpenLandMap_SOL_SOL_GRTGROUP_USDA-SOILTAX_C_v01...

Climate indicators are used in several statistical models for many research areas and are specially important for modelling Climate Sensitive Diseases (CSD) incidence. Those models usually adopts a lattice structure, where its data is aggregated at administrative boundaries (e.g. disease incidence), but climate indicators are usually presented in a continuous regular grid format.

To make climate indicators compatible with lattice structures, zonal statistics may be adopted. Zonal statistics are descriptive statistics calculated using a set of cells that spatially intersects a given spatial boundary. For each boundary in a map, statistics like average, maximum value, minimum value, standard deviation, and sum are obtained to represent the cell's values that intersect the boundary.

This dataset present zonal statistic of climate indicators computed from Copernicus ERA5-Land daily aggregates for the Brazilian municipalities, from 1950 to 2022.

This dataset contains 14 parquet-format files with monthly data.

    File
    Indicator
    Unit


    aet.parquet
    Actual Evapotranspiration
    mm


    def.parquet
    Climate Water Deficit
    mm


    pdsi.parquet
    Palmer Drought Severity Index (PDSI)
    unitless


    pet.parquet
    Precipitation
    mm


    ppt.parquet
    Potential evapotranspiration
    mm


    q.parquet
    Runoff
    mm


    soil.parquet
    Soil Moisture
    mm


    srad.parquet
    Downward surface shortwave radiation
    W/m2


    swe.parquet
    Snow water equivalent
    mm


    tmax.parquet
    Maximun Temperature
    °C


    tmin.parquet
    Minimum Temperature
    °C


    vap.parquet
    Vapor pressure
    kPa


    vpd.parquet
    Vapor Pressure Deficit
    kpq


    ws.parquet
    Wind speed
    m/s

PerCapita_CO2_Footprint_InDioceses_FULLBurhans, Molly A., Cheney, David M., Gerlt, R.. . “PerCapita_CO2_Footprint_InDioceses_FULL”. Scale not given. Version 1.0. MO and CT, USA: GoodLands Inc., Environmental Systems Research Institute, Inc., 2019.MethodologyThis is the first global Carbon footprint of the Catholic population. We will continue to improve and develop these data with our research partners over the coming years. While it is helpful, it should also be viewed and used as a "beta" prototype that we and our research partners will build from and improve. The years of carbon data are (2010) and (2015 - SHOWN). The year of Catholic data is 2018. The year of population data is 2016. Care should be taken during future developments to harmonize the years used for catholic, population, and CO2 data.1. Zonal Statistics: Esri Population Data and Dioceses --> Population per dioceses, non Vatican based numbers2. Zonal Statistics: FFDAS and Dioceses and Population dataset --> Mean CO2 per Diocese3. Field Calculation: Population per Diocese and Mean CO2 per diocese --> CO2 per Capita4. Field Calculation: CO2 per Capita * Catholic Population --> Catholic Carbon FootprintAssumption: PerCapita CO2Deriving per-capita CO2 from mean CO2 in a geography assumes that people's footprint accounts for their personal lifestyle and involvement in local business and industries that are contribute CO2. Catholic CO2Assumes that Catholics and non-Catholic have similar CO2 footprints from their lifestyles.Derived from:A multiyear, global gridded fossil fuel CO2 emission data product: Evaluation and analysis of resultshttp://ffdas.rc.nau.edu/About.htmlRayner et al., JGR, 2010 - The is the first FFDAS paper describing the version 1.0 methods and results published in the Journal of Geophysical Research.Asefi et al., 2014 - This is the paper describing the methods and results of the FFDAS version 2.0 published in the Journal of Geophysical Research.Readme version 2.2 - A simple readme file to assist in using the 10 km x 10 km, hourly gridded Vulcan version 2.2 results.Liu et al., 2017 - A paper exploring the carbon cycle response to the 2015-2016 El Nino through the use of carbon cycle data assimilation with FFDAS as the boundary condition for FFCO2."S. Asefi‐Najafabady P. J. Rayner K. R. Gurney A. McRobert Y. Song K. Coltin J. Huang C. Elvidge K. BaughFirst published: 10 September 2014 https://doi.org/10.1002/2013JD021296 Cited by: 30Link to FFDAS data retrieval and visualization: http://hpcg.purdue.edu/FFDAS/index.phpAbstractHigh‐resolution, global quantification of fossil fuel CO2 emissions is emerging as a critical need in carbon cycle science and climate policy. We build upon a previously developed fossil fuel data assimilation system (FFDAS) for estimating global high‐resolution fossil fuel CO2 emissions. We have improved the underlying observationally based data sources, expanded the approach through treatment of separate emitting sectors including a new pointwise database of global power plants, and extended the results to cover a 1997 to 2010 time series at a spatial resolution of 0.1°. Long‐term trend analysis of the resulting global emissions shows subnational spatial structure in large active economies such as the United States, China, and India. These three countries, in particular, show different long‐term trends and exploration of the trends in nighttime lights, and population reveal a decoupling of population and emissions at the subnational level. Analysis of shorter‐term variations reveals the impact of the 2008–2009 global financial crisis with widespread negative emission anomalies across the U.S. and Europe. We have used a center of mass (CM) calculation as a compact metric to express the time evolution of spatial patterns in fossil fuel CO2 emissions. The global emission CM has moved toward the east and somewhat south between 1997 and 2010, driven by the increase in emissions in China and South Asia over this time period. Analysis at the level of individual countries reveals per capita CO2 emission migration in both Russia and India. The per capita emission CM holds potential as a way to succinctly analyze subnational shifts in carbon intensity over time. Uncertainties are generally lower than the previous version of FFDAS due mainly to an improved nightlight data set."Global Diocesan Boundaries:Burhans, M., Bell, J., Burhans, D., Carmichael, R., Cheney, D., Deaton, M., Emge, T. Gerlt, B., Grayson, J., Herries, J., Keegan, H., Skinner, A., Smith, M., Sousa, C., Trubetskoy, S. “Diocesean Boundaries of the Catholic Church” [Feature Layer]. Scale not given. Version 1.2. Redlands, CA, USA: GoodLands Inc., Environmental Systems Research Institute, Inc., 2016.Using: ArcGIS. 10.4. Version 10.0. Redlands, CA: Environmental Systems Research Institute, Inc., 2016.Boundary ProvenanceStatistics and Leadership DataCheney, D.M. “Catholic Hierarchy of the World” [Database]. Date Updated: August 2019. Catholic Hierarchy. Using: Paradox. Retrieved from Original Source.Catholic HierarchyAnnuario Pontificio per l’Anno .. Città del Vaticano :Tipografia Poliglotta Vaticana, Multiple Years.The data for these maps was extracted from the gold standard of Church data, the Annuario Pontificio, published yearly by the Vatican. The collection and data development of the Vatican Statistics Office are unknown. GoodLands is not responsible for errors within this data. We encourage people to document and report errant information to us at data@good-lands.org or directly to the Vatican.Additional information about regular changes in bishops and sees comes from a variety of public diocesan and news announcements.GoodLands’ polygon data layers, version 2.0 for global ecclesiastical boundaries of the Roman Catholic Church:Although care has been taken to ensure the accuracy, completeness and reliability of the information provided, due to this being the first developed dataset of global ecclesiastical boundaries curated from many sources it may have a higher margin of error than established geopolitical administrative boundary maps. Boundaries need to be verified with appropriate Ecclesiastical Leadership. The current information is subject to change without notice. No parties involved with the creation of this data are liable for indirect, special or incidental damage resulting from, arising out of or in connection with the use of the information. We referenced 1960 sources to build our global datasets of ecclesiastical jurisdictions. Often, they were isolated images of dioceses, historical documents and information about parishes that were cross checked. These sources can be viewed here:https://docs.google.com/spreadsheets/d/11ANlH1S_aYJOyz4TtG0HHgz0OLxnOvXLHMt4FVOS85Q/edit#gid=0To learn more or contact us please visit: https://good-lands.org/Esri Gridded Population Data 2016DescriptionThis layer is a global estimate of human population for 2016. Esri created this estimate by modeling a footprint of where people live as a dasymetric settlement likelihood surface, and then assigned 2016 population estimates stored on polygons of the finest level of geography available onto the settlement surface. Where people live means where their homes are, as in where people sleep most of the time, and this is opposed to where they work. Another way to think of this estimate is a night-time estimate, as opposed to a day-time estimate.Knowledge of population distribution helps us understand how humans affect the natural world and how natural events such as storms and earthquakes, and other phenomena affect humans. This layer represents the footprint of where people live, and how many people live there.Dataset SummaryEach cell in this layer has an integer value with the estimated number of people likely to live in the geographic region represented by that cell. Esri additionally produced several additional layers World Population Estimate Confidence 2016: the confidence level (1-5) per cell for the probability of people being located and estimated correctly. World Population Density Estimate 2016: this layer is represented as population density in units of persons per square kilometer.World Settlement Score 2016: the dasymetric likelihood surface used to create this layer by apportioning population from census polygons to the settlement score raster.To use this layer in analysis, there are several properties or geoprocessing environment settings that should be used:Coordinate system: WGS_1984. This service and its underlying data are WGS_1984. We do this because projecting population count data actually will change the populations due to resampling and either collapsing or splitting cells to fit into another coordinate system. Cell Size: 0.0013474728 degrees (approximately 150-meters) at the equator. No Data: -1Bit Depth: 32-bit signedThis layer has query, identify, pixel, and export image functions enabled, and is restricted to a maximum analysis size of 30,000 x 30,000 pixels - an area about the size of Africa.Frye, C. et al., (2018). Using Classified and Unclassified Land Cover Data to Estimate the Footprint of Human Settlement. Data Science Journal. 17, p.20. DOI: http://doi.org/10.5334/dsj-2018-020.What can you do with this layer?This layer is unsuitable for mapping or cartographic use, and thus it does not include a convenient legend. Instead, this layer is useful for analysis, particularly for estimating counts of people living within watersheds, coastal areas, and other areas that do not have standard boundaries. Esri recommends using the Zonal Statistics tool or the Zonal Statistics to Table tool where you provide input zones as either polygons, or raster data, and the tool will summarize the count of population within those zones. https://www.esri.com/arcgis-blog/products/arcgis-living-atlas/data-management/2016-world-population-estimate-services-are-now-available/

Population density in 2010 within the boundaries of the Narragansett Bay watershed, the Southwest Coastal Ponds watershed, and the Little Narragansett Bay watershed. The methods for analyzing population were developed by the US Environmental Protection Agency ORD Atlantic Coastal Environmental Sciences Division in collaboration with the Narragansett Bay Estuary Program and other partners. Population rasters were generated using the USGS dasymetric mapping tool (see http://geography.wr.usgs.gov/science/dasymetric/index.htm) which uses land use data to distribute population data more accurately than simply within a census mapping unit. The 2010 10m cell population density raster was produced using Rhode Island (2011) state land use data, Massachusetts (2005) state land use, Connecticut (2011) NLCD land use data, and U.S. Census data (2010). To generate a population estimate (number of persons) for any given area within the boundaries of this raster, use the Zonal Statistics as Table tool to sum the 10m cell density values within your zone dataset (e.g., watershed polygon layer). For more information, please reference the 2017 State of Narragansett Bay & Its Watershed Technical Report (nbep.org).

This map features a global estimate of human population for 2016 with a focus on the Caribbean region . Esri created this estimate by modeling a footprint of where people live as a dasymetric settlement likelihood surface, and then assigned 2016 population estimates stored on polygons of the finest level of geography available onto the settlement surface. Where people live means where their homes are, as in where people sleep most of the time, and this is opposed to where they work. Another way to think of this estimate is a night-time estimate, as opposed to a day-time estimate.Knowledge of population distribution helps us understand how humans affect the natural world and how natural events such as storms and earthquakes, and other phenomena affect humans. This layer represents the footprint of where people live, and how many people live there.Dataset SummaryEach cell in this layer has an integer value with the estimated number of people likely to live in the geographic region represented by that cell. Esri additionally produced several additional layers World Population Estimate Confidence 2016: the confidence level (1-5) per cell for the probability of people being located and estimated correctly. World Population Density Estimate 2016: this layer is represented as population density in units of persons per square kilometer.World Settlement Score 2016: the dasymetric likelihood surface used to create this layer by apportioning population from census polygons to the settlement score raster.To use this layer in analysis, there are several properties or geoprocessing environment settings that should be used:Coordinate system: WGS_1984. This service and its underlying data are WGS_1984. We do this because projecting population count data actually will change the populations due to resampling and either collapsing or splitting cells to fit into another coordinate system. Cell Size: 0.0013474728 degrees (approximately 150-meters) at the equator. No Data: -1Bit Depth: 32-bit signedThis layer has query, identify, pixel, and export image functions enabled, and is restricted to a maximum analysis size of 30,000 x 30,000 pixels - an area about the size of Africa.Frye, C. et al., (2018). Using Classified and Unclassified Land Cover Data to Estimate the Footprint of Human Settlement. Data Science Journal. 17, p.20. DOI: https://doi.org/10.5334/dsj-2018-020.What can you do with this layer?This layer is unsuitable for mapping or cartographic use, and thus it does not include a convenient legend. Instead, this layer is useful for analysis, particularly for estimating counts of people living within watersheds, coastal areas, and other areas that do not have standard boundaries. Esri recommends using the Zonal Statistics tool or the Zonal Statistics to Table tool where you provide input zones as either polygons, or raster data, and the tool will summarize the count of population within those zones.

Estimates of malaria cases unseeking care by region in PNG (2011–2019).

File ds01. Global dataset of subsidence rates.

File ds02. Zonal statistics and feature importance.

File ds03. Global raster map of subsidence rates.

This dataset contains zonal-mean atmospheric diagnostics computed from reanalysis datasets on pressure levels. Primary variables include temperature, geopotential height, and the three-dimensional wind field. Advanced diagnostics include zonal covariance terms that can be used to compute, for instance, eddy kinetic energy and eddy fluxes. Terms from the primitive zonal-mean momentum equation and the transformed Eulerian momentum equation are also provided.

This dataset was produced to facilitate the comparison of reanalysis datasets for the collaborators of the SPARC- Reanalysis Intercomparison Project (S-RIP) project. The dataset is substantially smaller in size compared to the full three dimensional reanalysis fields and uses unified numerical methods. The dataset includes all global reanalyses available at the time of its development and will be extended to new reanalysis products in the future.

The tree canopy data was derived from the ArcGIS Living Atlas USA NLCD Tree Canopy Cover raster layer, which represents the canopy cover percentage for each 30-meter size raster cell. The granular raster data was rolled up to the census tract level by averaging the values for all the raster cells within each census tract (Ave Percent Tree Canopy Coverage). This was done with the Zonal Statistics as Table tool, similarly to the Air Quality workflow you saw earlier. The average percent tree canopy values were then used to determine the square kilometer coverage of tree canopy within each census tract (Total Canopy Coverage in Sq Km).Traffic data was also acquired from ArcGIS Living Atlas. The USA Traffic Counts point data was rolled up to the census tract for the city of Los Angeles (with the Summarize Within tool).

This layer shows particulate matter in the air sized 2.5 micrometers of smaller (PM 2.5). The data is aggregated from NASA Socioeconomic Data and Applications Center (SEDAC) gridded data into country boundaries, administrative 1 boundaries, and 50 km hex bins. The unit of measurement is micrograms per cubic meter.The layer shows the annual average PM 2.5 from 1998 to 2016, highlighting if the overall mean for an area meets the World Health Organization guideline of 10 micrograms per cubic meter annually. Areas that don't meet the guideline and are above the threshold are shown in red, and areas that are lower than the guideline are in grey.The data is averaged for each year and over the the 19 years to provide an overall picture of air quality globally. Some of the things we can learn from this layer:What is the average annual PM 2.5 value over 19 years? (1998-2016)What is the annual average PM 2.5 value for each year from 1998 to 2016?What is the statistical trend for PM 2.5 over the 19 years? (downward or upward)Are there hot spots (or cold spots) of PM 2.5 over the 19 years?How many people are impacted by the air quality in an area?What is the death rate caused by the joint effects of air pollution?Choose a different attribute to symbolize in order to reveal any of the patterns above.A space time cube was performed on a multidimensional mosaic version of the data in order to derive an emerging hot spot analysis, trends, and a 19-year average. The country and administrative 1 layers provide a population-weighted PM 2.5 value to emphasize which areas have a higher human impact. Citations:van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2018. Global Annual PM2.5 Grids from MODIS, MISR and SeaWiFS Aerosol Optical Depth (AOD) with GWR, 1998-2016. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC). https://doi.org/10.7927/H4ZK5DQS. Accessed 1 April 2020van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2016. Global Estimates of Fine Particulate Matter Using a Combined Geophysical-Statistical Method with Information from Satellites. Environmental Science & Technology 50 (7): 3762-3772. https://doi.org/10.1021/acs.est.5b05833.Boundaries and population figures:Antarctica is excluded from all maps because it was not included in the original NASA grids.50km hex bins generated using the Generate Tessellation tool - projected to Behrmann Equal Area projection for analysesPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Administrative boundaries from World Administrative Divisions layer from ArcGIS Living Atlas - projected to Behrmann Equal Area projection for analyses and hosted in Web MercatorSources: Garmin, CIA World FactbookPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Country boundaries from Esri 2019 10.8 Data and Maps - projected to Behrmann Equal Area projection for analyses and hosted in Web Mercator. Sources: Garmin, Factbook, CIAPopulation figures attached to the country boundaries come from the World Population Estimate 2016 Sources Living Atlas layer Data processing notes:NASA's GeoTIFF files for 19 years (1998-2016) were first brought into ArcGIS Pro 2.5.0 and put into a multidimensional mosaic dataset.For each geography level, the following was performed: Zonal Statistics were run against the mosaic as a multidimensional layer.A Space Time Cube was created to compare the 19 years of PM 2.5 values and detect hot/cold spot patterns. To learn more about Space Time Cubes, visit this page.The Space Time Cube is processed for Emerging Hot Spots where we gain the trends and hot spot results.The layers are hosted in Web Mercator Auxillary Sphere projection, but were processed using an equal area projection: Behrmann. If using this layer for analysis, it is recommended to start by projecting the data back to Behrmann.The country and administrative layer were dissolved and joined with population figures in order to visualize human impact.The dissolve tool ensures that each geographic area is only symbolized once within the map.Country boundaries were generalized post-analysis for visualization purposes. The tolerance used was 700m. If performing analysis with this layer, find detailed country boundaries in ArcGIS Living Atlas. To create the population-weighted attributes on the country and Admin 1 layers, the hex value population values were used to create the weighting. Within each hex bin, the total population figure and average PM 2.5 were multiplied.The hex bins were converted into centroids and the PM2.5 and population figures were summarized within the country and Admin 1 boundaries.The summation of the PM 2.5 values were then divided by the total population of each geography. This population value was determined by summarizing the population values from the hex bins within each geography.Some artifacts in the hex bin layer as a result of the input NASA rasters. Because the gridded surface is created from multiple satellites, there are strips within some areas that are a result of satellite paths. Some areas also have more of a continuous pattern between hex bins as a result of the input rasters.Within the country layer, an air pollution attributable death rate is included. 2016 figures are offered by the World Health Organization (WHO). Values are offered as a mean, upper value, lower value, and also offered as age standardized. Values are for deaths caused by all possible air pollution related diseases, for both sexes, and all age groups. For more information visit this page, and here for methodology. According to WHO, the world average was 95 deaths per 100,000 people.To learn the techniques used in this analysis, visit the Learn ArcGIS lesson Investigate Pollution Patterns with Space-Time Analysis by Esri's Kevin Bulter and Lynne Buie.

High resolution land cover data set for New York City. This is the 3ft version of the high-resolution land cover dataset for New York City. Seven land cover classes were mapped: (1) tree canopy, (2) grass/shrub, (3) bare earth, (4) water, (5) buildings, (6) roads, and (7) other paved surfaces. The minimum mapping unit for the delineation of features was set at 3 square feet. The primary sources used to derive this land cover layer were the 2010 LiDAR and the 2008 4-band orthoimagery. Ancillary data sources included GIS data (city boundary, building footprints, water, parking lots, roads, railroads, railroad structures, ballfields) provided by New York City (all ancillary datasets except railroads); UVM Spatial Analysis Laboratory manually created railroad polygons from manual interpretation of 2008 4-band orthoimagery. The tree canopy class was considered current as of 2010; the remaining land-cover classes were considered current as of 2008. Object-Based Image Analysis (OBIA) techniques were employed to extract land cover information using the best available remotely sensed and vector GIS datasets. OBIA systems work by grouping pixels into meaningful objects based on their spectral and spatial properties, while taking into account boundaries imposed by existing vector datasets. Within the OBIA environment a rule-based expert system was designed to effectively mimic the process of manual image analysis by incorporating the elements of image interpretation (color/tone, texture, pattern, location, size, and shape) into the classification process. A series of morphological procedures were employed to insure that the end product is both accurate and cartographically pleasing. More than 35,000 corrections were made to the classification. Overall accuracy was 96%. This dataset was developed as part of the Urban Tree Canopy (UTC) Assessment for New York City. As such, it represents a 'top down' mapping perspective in which tree canopy over hanging other features is assigned to the tree canopy class. At the time of its creation this dataset represents the most detailed and accurate land cover dataset for the area. This project was funded by National Urban and Community Forestry Advisory Council (NUCFAC) and the National Science Fundation (NSF), although it is not specifically endorsed by either agency. The methods used were developed by the University of Vermont Spatial Analysis Laboratory, in collaboration with the New York City Urban Field Station, with funding from the USDA Forest Service.

As a part of DVRPC’s long-range planning activities, the Commission is required to maintain forecasts with at least a 20-year horizon, or to the horizon year of the long-range plan. Allocation of growth is forecasted using a land use model, UrbanSim, and working closely with member county planning staffs. DVRPC has prepared regional, county, and municipal-level population and employment forecasts in five-year increments through 2050, using 2015 Census population estimates and 2015 National Establishments Time Series (NETS) employment data as the base. A forthcoming Analytical Data Report will document the forecasting process and methodologies. This data provides the municipal level forecast for our 8 counties outside of Philadelphia, and Planning District-level forecast within Philadelphia. Note: while 2019 land use model results are provided, the forecast was only adopted for 2015, 2020, 2025, 2030, 2035, 2040, 2045, and 2050.

Population density in 1990 within the boundaries of the Narragansett Bay watershed, the Southwest Coastal Ponds watershed, and the Little Narragansett Bay watershed. The methods for analyzing population were developed by the US Environmental Protection Agency ORD Atlantic Coastal Environmental Sciences Division in collaboration with the Narragansett Bay Estuary Program and other partners. Population rasters were generated using the USGS dasymetric mapping tool (see http://geography.wr.usgs.gov/science/dasymetric/index.htm) which uses land use data to distribute population data more accurately than simply within a census mapping unit. The 1990 10m cell population density raster was produced using Rhode Island 1988 state land use data, Massachusetts 1985 state land use, Connecticut 1992 NLCD land use data, and U.S. Census data (1990). To generate a population estimate (number of persons) for any given area within the boundaries of this raster, use the Zonal Statistics as Table tool to sum the 10m cell density values within your zone dataset (e.g., watershed polygon layer). For more information, please reference the 2017 State of Narragansett Bay & Its Watershed Technical Report (nbep.org).

We used a stratified-random sampling approach to estimate the total abundance of Wedge-tailed Shearwater (Ardenna pacifica) nest sites across Kīlauea Point National Wildlife Refuge (KPNWR), Kauaʻi, during 1-7 July 2019. We first identified strata as unique geographic areas of the refuge to account for potential differences in nesting habitat and non-uniform nest site clustering. We then sub-divided strata where we expected high, low, minimal, or no nest site abundance. These distinctions were based on knowledge of shearwater nesting distribution gained while performing extensive ground-searching for tropicbirds across the entire refuge in April and May 2019. We delineated strata boundaries using recent satellite imagery in ArcGIS (version 10.7) and, based on direct observations in the field, refined in order to remove large contiguous areas lacking shearwater presence or nesting habitat. Planar area of each polygon was automatically calculated by ArcGIS. To calculate surface area in each stratum, we obtained a 1/3 arc-second (10-m resolution) digital elevation model from the National Elevation Dataset (USGS 2013). This elevation raster was projected to Universal Transverse Mercator projection (Zone 4 North, North American Datum 1983) and converted to a surface area raster using the raster and sp packages in R (Hijmans 2019; Pebesma & Bivand 2005). We calculated stratum-specific surface area by summing the surface area values of raster cells in each stratum using XToolsPro (version 18; Zonal Statistics tool) in ArcMap (version 10.7). Hijmans RJ. 2019. raster: Geographic Data Analysis and Modeling. R package version 3.0-7. https://CRAN.R-project.org/package=raster Pebesma EJ, Bivand RS. 2005. Classes and methods for spatial data in R. R News 5(2):9-13. U.S. Geological Survey. 2013. USGS NED n23w160 1/3 arc-second 2013 1 x 1 degree ArcGrid: U.S. Geological Survey. Accessed at https://www.sciencebase.gov/catalog/item/581d21dae4b08da350d53be2

This layer contains wildfire hazard potential (WHP) data for the conterminous United States aggregated from states to block groups and 50 km hex bins then enriched with demographic data. The data is from the USDA Forest Service Fire Modeling Institute providing an index of WHP at a 270 meter resolution. Wildfire hazard potential provides information on the relative potential for wildfire that would be difficult for fire crews to contain. "Areas with higher wildfire potential values represent fuels with a higher likelihood of experiencing high-intensity fire with torching, crowning, and other forms of extreme fire behavior." - Fire Modeling Institute. A score of 5 is very high risk and a score between 0-1 is likely non-burnable area such as water or asphalt. "On its own, WHP is not an explicit map of wildfire threat or risk, but when paired with spatial data depicting highly valued resources and assets such as communities, structures, or powerlines, it can approximate relative wildfire risk to those resources and assets. WHP is also not a forecast or wildfire outlook for any particular season, as it does not include any information on current or forecasted weather or fuel moisture conditions. It is instead intended for long-term strategic planning and fuels management."Each layer has been enriched with 2020 Esri demographic attributes to better approximate wildfire hazard risk relating to the human population. This layer is available in a ready to use web map. A hosted imagery layer of this data is available in ArcGIS Living Atlas for additional analysis.Data notes:Zonal Statistics as Table were run against a local copy of the WHP data using US standard geographies as the feature zone input for the analysis. Geographies included are: State, County, Congressional District, ZIP Code, Tract, and Block Group. Statistical tables were joined to geographies. To learn more about zonal statistics, view the documentation here. 50 km hex bins were created using Generate Tessellation and then joined to zonal statistics as described above (step 1).Data was enriched with 2020 Esri Demographics. Attributes include population, households & housing units, growth rate, and calculated variables such as population change over time. To create the population-weighted attributes on the state, congressional district, and county layers, the hex value population values were used to create the weighting. Within each hex bin, the total population figure and average WHP were multiplied.The hex bins were converted into centroids and summarized within the state, congressional district, and county boundaries.The summation of these values were then divided by the total population of each respective geography.

This data is necessary to run the i-SoMPE WP1 R Project available on GitLab and Zenodo.

This NOAA Climate Data Record (CDR) of Zonal Mean Ozone Binary Database of Profiles (BDBP) dataset is a vertically resolved, global, gap-free and zonal mean dataset that was created with a multiple-linear regression model. The dataset has a monthly resolution and spans the period 1979 to 2007. It provides global product in 5 degree zonal bands, and 70 vertical levels of the atmosphere. The regression is based on monthly mean ozone concentrations that were calculated from several different satellite instruments and global ozone soundings. Due to the regression model that was used to create the product, various basis function contributions are provided as unique levels or tiers. To understand the different contributions of basis functions, the data product is provided in five different "Tiers". - Tier 0: raw monthly mean data that was used in the regression model - Tier 1.1: Anthropogenic influences (as determined by the regression model) - Tier 1.2: Natural influences (as determined by the regression model) - Tier 1.3: Natural and volcanic influences (as determined by the regression model) - Tier 1.4: All influences (as determined by the regression model, CDR variable)

Spatial analysis and statistical summaries of the Protected Areas Database of the United States (PAD-US) provide land managers and decision makers with a general assessment of management intent for biodiversity protection, natural resource management, and recreation access across the nation. The PAD-US 3.0 Combined Fee, Designation, Easement feature class in the full geodatabase inventory (with Military Lands and Tribal Areas from the Proclamation and Other Planning Boundaries feature class) was modified to prioritize overlapping designations, avoiding massive overestimation in protected area statistics, and simplified by the following PAD-US attributes to support user needs for raster analysis data: Manager Type, Manager Name, Designation Type, GAP Status Code, Public Access, and State Name. The rasterization process (see processing steps below) prioritized overlapping designations previously identified (GAP_Prity field) in the Vector Analysis File (e.g. Wilderness within a National Forest) based upon their relative biodiversity conservation (e.g. GAP Status Code 1 over 2). The 30-meter Image (IMG) grid Raster Analysis Files area extents were defined by the Census state boundary file used to clip the Vector Analysis File, the data source for rasterization ("PADUS3_0VectorAnalysis_State_Clip_CENSUS2020" feature class from ("PADUS3_0VectorAnalysisFileOtherExtents_Clip_Census.gdb"). Alaska (AK) and Hawaii (HI) raster data are separated from the contiguous U.S. (CONUS) to facilitate analyses at manageable scales. Note, the PAD-US inventory is now considered functionally complete with the vast majority of land protection types (with a legal protection mechanism) represented in some manner, while work continues to maintain updates, improve data quality, and integrate new data as it becomes available (see inventory completeness estimates at: http://www.protectedlands.net/data-stewards/ ). In addition, protection status represents a point-in-time and changes in status between versions of PAD-US may be attributed to improving the completeness and accuracy of the spatial data more than actual management actions or new acquisitions. USGS provides no legal warranty for the use of this data. While PAD-US is the official aggregation of protected areas ( https://www.fgdc.gov/ngda-reports/NGDA_Datasets.html ), agencies are the best source of their lands data.