Curtis Price 20111120 This data set is a 30-meter resolution raster grid of 2010 population density for the 48 conterminous US states. Data are provided in Arc GRID format (ESRI, Inc.; www.esri.com). Values are integer values, representing persons per sq km. Data are in Albers Conical Equal-Area projection, as were the original source data. Support of analysis for ecological, pesticide nutrient, and other studies to be conducted by NAWQA national teams. Supplemental_Information 2010 Released data from U.S. Census website, November 2011 Complete none planned -125 -66 50 24 none Census 2010 Population density Block group none Conterminous United States Internal use only Internal use only. Curtis Price USGS NAWQA National Synthesis Physical Scientist mailing address

1608 Mountain View Rd

Rapid City SD 57702 USA 605-394-3242 605-355-4520 cprice@usgs.gov Microsoft Windows XP Professional ArcGIS ArcInfo version 10.0 none Data are complete. U.S. Bureau of the Census 2011 http://www.census.gov/geo/www/tiger/shp.html 100000 online zip files, retrieved 11/14/2011 2010 publication date BLKGROUP Polygon boundaries for block groups U.S. Bureau of the Census 2011 ftp://ftp2.census.gov/geo/tiger/TIGER2010BLKPOPHU/ 100000 online zip files, retrieved 11/14/2011 20100101 publication date TABBLOCK 2010 Census block population totals Downloaded data 20111114 Converted polygon shapefiles to a 30-meter raster grid. An ArcGIS 10.0 model builder script was run to convert each set of block group boundaries (by state) to polygons. An "exported" python version of is attached to document the processing steps. > # --------------------------------------------------------------------------- > # script.py > # Created on: 2012-02-23 16:09:53.00000 > # (generated by ArcGIS/ModelBuilder) > # Usage: script <Block_group_shapefile> <Block_group_raster> > # Description: > # --------------------------------------------------------------------------- > > # Import arcpy module > import arcpy > > # Check out any necessary licenses > arcpy.CheckOutExtension("spatial") > > # Load required toolboxes > arcpy.ImportToolbox("E:/work/bg2010/Census To Raster Tools.tbx") > > # Script arguments > Block_group_shapefile = arcpy.GetParameterAsText(0) > if Block_group_shapefile == '#' or not Block_group_shapefile: > Block_group_shapefile = "tl_2010_11_bg10" # provide a default value if unspecified > > Block_group_raster = arcpy.GetParameterAsText(1) > if Block_group_raster == '#' or not Block_group_raster: > Block_group_raster = "E:\\work\\bg10g" # provide a default value if unspecified > > # Local variables: > tabblock2010_11_pophu_Projec1_shp = Block_group_shapefile > Raster_table__2_ = Block_group_shapefile > Raster_table = Block_group_raster > Snap_raster = "%scratchworkspace%\\xxsnap" > > # Process: Project > arcpy.Project_management(Block_group_shapefile, tabblock2010_11_pophu_Projec1_shp, > "PROJCS['USA_Contiguous_Albers_Equal_Area_Conic_USGS_version',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983'," \ > "SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]," \ > "PROJECTION['Albers'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0]," \ > "PARAMETER['Central_Meridian',-96.0],PARAMETER['Standard_Parallel_1',29.5],PARAMETER['Standard_Parallel_2',45.5]," \ > "PARAMETER['Latitude_Of_Origin',23.0],UNIT['Meter',1.0]]", "", > "GEOGCS['GCS_North_American_1983'" \ > "DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]]," \ > "PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]") > > # Process: Create Snap Raster > arcpy.CreateSnapRaster_c2raster("15", "15", "30", Snap_raster, > "PROJCS['USA_Contiguous_Albers_Equal_Area_Conic_USGS_version',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983'," \ > "SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]," \ > "PROJECTION['Albers'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0]," \ > "PARAMETER['Central_Meridian',-96.0],PARAMETER['Standard_Parallel_1',29.5],PARAMETER['Standard_Parallel_2',45.5]," \ > "PARAMETER['Latitude_Of_Origin',23.0],UNIT['Meter',1.0]]") > > # Process: Polygon to Raster > tempEnvironment0 = arcpy.env.snapRaster > arcpy.env.snapRaster = "%scratchworkspace%\\xxsnap" > arcpy.PolygonToRaster_conversion(tabblock2010_11_pophu_Projec1_shp, "FID", Block_group_raster, "MAXIMUM_AREA", "NONE", "30") > arcpy.env.snapRaster = tempEnvironment0 > > # Process: Make Table View > arcpy.MakeTableView_management(Block_group_raster, Raster_table, "", "", "VALUE VALUE VISIBLE NONE;COUNT COUNT VISIBLE NONE") > > # Process: Join Field > arcpy.JoinField_management(Raster_table, "VALUE", Block_group_shapefile, "FID", "GEOID10;ALAND10;AWATER10" 20111115 Tabulated TABBLOCK block data by block group Attributed state-coverage block polygons were downloaded from the FTP site: ftp://ftp2.census.gov/geo/tiger/TIGER2010BLKPOPHU/ These were tabulated by block group using the ArcGIS Frequency tool. 2011115 Joined TABLLOCK summaries to raster attribute table The usbg10 raster table was then joined to the TABBLOCK table results 20111115 Raster Entity_and_Attribute_Overview none Curtis Price USGS NAWQA National Synthesis Physical Scientist mailing address

1608 Mountain View Rd

Rapid City SD 57702 USA 605-394-3242 605-355-4520 cprice@usgs.gov Although these data have been used by the U.S. Geological Survey, U.S. Department of the Interior, no warranty expressed or implied is made by the U.S. Geological Survey as to the accuracy of the data. The act of distribution shall not constitute any such warranty, and no responsibility is assumed by the U.S. Geological Survey in the use of this data, software, or related materials. 20120223 Curtis Price USGS NAWQA National Synthesis Physical Scientist mailing address

1608 Mountain View Rd

Rapid City SD 57702 USA 605-394-3242 605-355-4520 cprice@usgs.gov FGDC Content Standards for Digital Geospatial Metadata FGDC-STD-001-1998 none none