Hestia SWIFL Anthropogenic Heat (AH) Data Product - AZ_urban Version 1.0 - WRF 2km grid
created: Feb 6, 2025

Contact:
     Anna Kato (anna.kato@nau.educ)
     Kevin Gurney (kevin.gurney@nau.edu)
     Northern Arizona University

HESTIA-SWIFL DATA PRODUCT FOR AZ URBAN AREAS
The Hestia-SWIFL AZ_urban v1.0 AH data represent emissions due to the combustion of fossil fuel 
and cement production in Arizona urban areas. The exact spatial extent is defined by the ASU
Weather Research Forecast (WRF) simulation grid. 
The emissions are generated using a bottom-up/engineering approach and are tied to results generated by 
the Vulcan Project version 4, an effort to quantify space/time-resolved FFCO2 & AH emissions for the entire 
United States landscape. A large number of data sources are combined to best estimate FFCO2 and AH
emissions at fine scales such as air quality emissions data, traffic flow data, building information, 
sociodemographic information, and fuel statistics. The native spatial resolution of 
the Hestia FFCO2 & AH emissions data product is a combination of points, lines, and polygons 
dictated primarily by the underlying data sources and the Vulcan emissions output. 
The output made available here places this information into a regularized 2km x 2km grid at hourly and annual 
temporal resolutions.

FAIR USE DISCLAIMER
The Hestia SWIFL Anthropogenic Heat data product is an estimation of AH emissions at very fine scales in 
time and space. It should be considered a "climatology" of emissions rather than the "weather" of 
emissions such that the estimates represent "typical" emissions at a given time and place (average 
conditions). Hence, it is not appropriate to use the data in comparison to short-term "campaign 
style" atmospheric measurements (e.g. 5 days of continuous monitoring at a specific location) 
without consideration of and reference to the mismatch between the measurement and the Hestia 
estimation approach. Users are encouraged to contact Kevin Gurney for updates and consultation 
on such potential use.

SECTORAL COMPOSITION
This product provides AH emissions for two emissions sources:
"transportation" : all mobile souces (Onroad, Nonroad, Airport, Railroad)
"elecp_com_ind_cement": point source emissions (Electricity Production, Commercial, Industrial, cement production)
Note that the AH emissions are Scope 1 only, therefore the emissions driven by electricity consumption are located at 
power production facilities. Also, elecp_com_ind_cement emissions data do NOT include the heat from buildings.

DATA FILE NAMES AND FILE STRUCTURE 
The Hestia-SWIFL AH emissions are represented using two time-resolutions: annual and hourly for 
the 4 year timespan of 2019 to 2022, inclusive.  The 2km x 2km gridded AH emissions are stored in 
netCDF files and represent the emissions resulting from integration of all point, line, and polygon 
elements within a gridcell (using area or length proportions for line and polygon elements that straddle 
gridcells). The annual files contain a 246-row by 231-column matrix with double-precision values 
arranged in row-major order - each row runs longitudinally and each column latitudinally with the first 
grid cell located at the northwestern corner. Each of the hourly emissions files contain 8760 (for the 
years other than 2020) or 8784 (for the year 2020) grid arrays, each of which is a double-precision 
matrix of the same structure as that in the annual files. The time domain reflects local time. The grid 
details are provided in Table 1. All the gridded AH emissions are in units of watt per square meters (W/m2). 
Sector-specific results are available upon request.
 
Table 1. 2km x 2km grid parameters

Parameter
  dimensions:
    X = 231;
    Y = 246;
    T = 8760;
  variables:
    double X(X=231);
      :units = "2km";
      :long_name = "X";

    double Y(Y=246);
      :units = "2km";
      :long_name = "Y";

    int T(T=8760);
      :units = "hour";
      :long_name = "T";
}

Projection: "AZ_North_America_Lambert_Conformal_Conic"

WKT:
PROJCRS["AZ_North_America_Lambert_Conformal_Conic",
    BASEGEOGCRS["NAD83",
        DATUM["North American Datum 1983",
            ELLIPSOID["GRS 1980",6378137,298.257222101,
                LENGTHUNIT["metre",1]],
            ID["EPSG",6269]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["Degree",0.0174532925199433]]],
    CONVERSION["unnamed",
        METHOD["Lambert Conic Conformal (2SP)",
            ID["EPSG",9802]],
        PARAMETER["Latitude of false origin",40,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8821]],
        PARAMETER["Longitude of false origin",-112,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8822]],
        PARAMETER["Latitude of 1st standard parallel",30,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8823]],
        PARAMETER["Latitude of 2nd standard parallel",60,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8824]],
        PARAMETER["Easting at false origin",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8826]],
        PARAMETER["Northing at false origin",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8827]]],
    CS[Cartesian,2],
        AXIS["(E)",east,
            ORDER[1],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]],
        AXIS["(N)",north,
            ORDER[2],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]]]
				

The gridded output files follow a naming convention that indicates the key attributes of the file. For 
example, the file "AZ_urban_transportation.hourly.2km.heat.watts_per_m2.2019.nc" indicates the grid resolution (2km), 
the emissions domain (AZ_urban), the sector (e.g. "transportation", "pt_elecp_com_ind_cement"), 
the time resolution (e.g. "annual", "hourly"), emissions data ("heat.watts_per_m2"), and the data year (2019). 
For questions on file formats and updates, please contact Kevin Gurney and/or Anna Kato.