{3f42970d-b2b5-4d3e-9db4-8cddfd226498} eng; USA utf8 dataset EnviroAtlas Coordinator U.S. Environmental Protection Agency, Office of Research and Development-Sustainable and Healthy Communities Research Program, EnviroAtlas Geospatial Data Owner (919) 541-3832 109 T.W. Alexander Drive Research Triangle Park NC 27709 EnviroAtlas@epa.gov https://www.epa.gov/enviroatlas pointOfContact 2025-07-23 ISO 19115 Geographic Information - Metadata ISO 19115:2003/Cor.1:2006(E) EnviroAtlas - Crop phosphorus removal by 12-digit HUC for 2012 for the Conterminous United States 2016-10-28 publication U.S. EPA Office of Research & Development (ORD) - National Exposure Research Laboratory (NERL) Research Triangle Park NC publisher US Environmental Protection Agency (US EPA) originator This EnviroAtlas national map displays the mean crop phosphorus (P) removal from croplands in the conterminous United States (excluding Hawaii and Alaska) for the year 2012 by 12-digit HUC. These data are based on International Plant Nutrition Institute (IPNI) compilations of county-level major crop harvest and P content of these crops, and cropland area from the USGS's U.S. conterminous wall-to-wall anthropogenic land use trends (NWALT) 2012 land cover data. This dataset was produced by the US EPA to support research and online mapping activities related to EnviroAtlas. EnviroAtlas (https://www.epa.gov/enviroatlas) allows the user to interact with a web-based, easy-to-use, mapping application to view and analyze multiple ecosystem services for the contiguous United States. The dataset is available as downloadable data (https://edg.epa.gov/data/Public/ORD/EnviroAtlas) or as an EnviroAtlas map service. Additional descriptive information about each attribute in this dataset can be found in its associated EnviroAtlas Fact Sheet (https://www.epa.gov/enviroatlas/enviroatlas-fact-sheets). The map, Crop Phosphorus Removal, is one of four EnviroAtlas maps that display P inputs and agricultural crop P demand for the conterminous US. These data could be used either alone or in conjunction with other data layers to help identify areas where P is a significant pollutant source (over application) or where there are opportunities for more efficient management or recycling to meet crop demands. These data could also be used in models that examine the transport and cycling of P across terrestrial and aquatic ecosystems. Information on crop P demand and uptake is, or will be, needed for the development of nutrient reduction strategies, nutrient credit exchanges, and payments for ecosystem services. Information on P uptake/removal from crops (like the 2012 data presented here) is important in order to understand whether P sources are being over or under applied compared to crop needs. Because of its essential role in agriculture, but potential as a pollutant, understanding where and how much P is applied as fertilizer or manure versus used by crops is important to inform management strategies that increase food security and water quality across the U.S. The overall goal of EnviroAtlas is to employ and develop the best available science to map indicators of ecosystem services production, demand, and drivers for the nation. completed EnviroAtlas Coordinator U.S. Environmental Protection Agency, Office of Research and Development-Sustainable and Healthy Communities Research Program, EnviroAtlas Geospatial Data Owner (919) 541-3832 109 T.W. Alexander Drive Research Triangle Park NC 27709 EnviroAtlas@epa.gov https://www.epa.gov/enviroatlas pointOfContact asNeeded Ecosystem Environment theme EPA GIS Keyword Thesaurus Agriculture Air Ecosystem Services EnviroAtlas Environment Environmental Atlas Health Water theme EnviroAtlas Alabama Arizona Arkansas California Colorado Connecticut Delaware Florida Georgia Idaho Illinois Indiana Iowa Kansas Kentucky Louisiana Maine Maryland Massachusetts Michigan Minnesota Mississippi Missouri Montana Nebraska Nevada New Hampshire New Jersey New Mexico New York North Carolina North Dakota Ohio Oklahoma Oregon Pennsylvania Rhode Island South Carolina South Dakota Tennessee Texas United States Utah Vermont Virginia Washington Washington DC West Virginia Wisconsin Wyoming place None https://edg.epa.gov/EPA_Data_License.html otherRestrictions Use Constraints: None. Please check sources, scale, accuracy, currency and other available information. Please confirm that you are using the most recent copy of both data and metadata. Acknowledgement of the EPA would be appreciated. otherRestrictions Access Constraints: None PUBLIC https://project-open-data.github.io/schema/#accessLevel Standard Technical Controls eng; USA environment -128.02839229 -65.20336449 22.7348894 51.6773602 2012-01-01 We distributed the county-specific, per area P removal rates from IPNI to agricultural lands from NWALT (60 m x 60 m pixels) within the corresponding county. In order to correct for some pixels with unrealistically high P removal rates (most likely caused by low agricultural land), we capped P removal at 9,000 kg of P per km2, which is higher than potential P removal for high estimates of corn yield (300 bu/acre) and corn P removal. true Entity and Attribute Information Entity and Attribute Overview: Crop Phosphorus Removal in 2012 was estimated using county-level estimates of crop harvest. We acquired county-level data describing total crop P uptake (kg P yr-1) in 2012 from IPNI, which was compiled using 22 major-crop yields and harvested acres from the USDA National Agricultural Statistics Service summaries and US Census of agriculture (see IPNI detailed methodology at http://nugis.ipni.net/Methods/Removal/ for crops considered and removal rates used). We used the U.S. National wall-to-wall land use trends (NWALT) for 2012, acquired from the USGS, at the scale of 60 m x 60 m. These data were converted to per area rates (kg P km-2 yr-1) of crop P removal by dividing the total P crop removal by the land area (km2) of combined cultivated crop and hay/pasture (agricultural) lands within a county as determined from county-level summarization of the 2012 NWALT layer. We distributed the county-specific per area P removal rates from IPNI to agricultural lands from NWALT (60 m x 60 m pixels) within the corresponding county. In order to correct for some pixels with unrealistically high P removal rates (most likely caused by low agricultural land), we capped P removal at 9,000 kg of P per km2, which is higher than potential removal of high estimates of corn yield (300 bu/acre) and corn P removal. The estimates were then averaged across each 12-digit HUC and converted to kg/ha/year by dividing by 100. Entity and Attribute Detail Citation: https://www.epa.gov/enviroatlas/enviroatlas-fact-sheets EnviroAtlas Coordinator U.S. Environmental Protection Agency, Office of Research and Development-Sustainable and Healthy Communities Research Program, EnviroAtlas Geospatial Data Owner (919) 541-3832 109 T.W. Alexander Drive Research Triangle Park NC 27709 EnviroAtlas@epa.gov https://www.epa.gov/enviroatlas distributor https://enviroatlas.epa.gov/download/CONUS_metrics_FGDB.zip download Download The URL providing direct access to the downloadable dataset. download https://enviroatlas.epa.gov/download/CONUS_metrics_CSV.zip ESRI:ArcGIS API The endpoint of a web service to access the dataset (REST endpoint, WMS GetCapabilities URL, or a SOAP WSDL endpoint). webMapService https://www.epa.gov/enviroatlas/enviroatlas-data information Documentation URL for accessing related documents such as data dictionary, technical information about a dataset, developer documentation, etc. information https://www.epa.gov/enviroatlas information Home Page Alternative landing page used to redirect user to a contextual, Agency-hosted "homepage" for the Dataset or API information dataset Horizontal Positional Accuracy Report Data were collected using methods that have unknown accuracy (EPA National Geospatial Data Policy [NGDP] Accuracy Tier 10). For more information, please see EPA's NGDP at https://www.epa.gov/geospatial/geospatial-policies-and-standards Features represented have not been tested for completeness Features represented have not been tested for completeness Tests for integrity have not been performed. This lineage was generated by the FGDC CSDGM to 19139 transformation Download IPNI crop P removal data at the county scale from the NuGIS website (http://nugis.ipni.net/login/). 2016-07-10T00:00:00 Add the FIPS ID to the county ID numbers to avoid losing leading zeros, so the county ID can be used to match the IPNI data with the land use data. 2016-07-10T00:00:00 Convert P2O5 for crop removal to P by multiplying by 0.4364. 2016-07-10T00:00:00 Convert metric tons to kg by multiplying by 1000 2016-07-10T00:00:00 Copy and paste these numbers to a CSV file for later use 2016-07-10T00:00:00 Download NWALT 2012 land use data layer from USGS and open in ArcGIS 2016-07-10T00:00:00 Overlay with U.S. county layer (make sure projection matches land use layer and reproject the county layer to Albers Equal Area NAD 1983 Albers if it does not match) 2016-07-10T00:00:00 Run the zonal histogram tool to get the land use make up of each county 2016-07-10T00:00:00 Copy and paste the output into Excel and transpose so that the row labels are county identifiers and the column labels are the objects (the land use IDs) 2016-07-10T00:00:00 In a new Excel sheet, multiply the number in a cell (= to the number of grid cells for each land use category) by (60*0.001)^2 to get from 60m by 60m grid cell count to area in km2. Save this new sheet a CSV for later use 2016-07-10T00:00:00 Using R, match the county P removal CSV with the county land use CSV data using the FIPS (county identifier) and divide the total P removal by the NWALT agricultural land use area (combined cultivated crop and hay/pasture, which are coded as P43k and P44k) and save this output as a new CSV 2016-07-10T00:00:00 In Excel, remove "FIPS_" from the county ID numbers with the formula =RIGHT(cell, LEN(cell)-5) to join with the land use layer. 2016-07-10T00:00:00 Sort data so that the largest P removal rate is at the top 2016-07-10T00:00:00 Change the Inf to NAs and change NA to 0 2016-07-10T00:00:00 For values of crop P removal that are above 9,000, replace them with 9,000 (justification for that cut off rate is that it is the highest yield and P removal rate for corn). Save the county identifiers and the transformed removal rates for later use. 2016-07-10T00:00:00 In ArcGIS, used the Extract by Attributes tool to create a new raster that just has the 2 desired NWALT agricultural land uses (cultivated crops and hay/pasture) with "VALUE"=43 or "VALUE"=44 formula and save this layer 2016-07-10T00:00:00 Create a binary version of that layer with the reclassify tool (1 if cell is 43 or 44 and 0 if cell is NA) and save it 2016-07-10T00:00:00 Join the previously saved CSV file containing the P crop removal rates with the U.S. county layer using the join attributes tool (base the join on the county identifiers) and export the data. 2016-07-10T00:00:00 Use the polygon to raster tool to create a raster of the P crop removal rates from the polygons created in the previous step and using the NWALT layer's cell size of 60m. 2016-07-10T00:00:00 Multiply the binary agricultural land use layer by the crop P removal rate raster using the raster calculator tool. This distributes the county-specific annual per-area P removal rate to the cultivated crop and hay/pasture lands 2016-07-10T00:00:00 Used ArcGIS tool Zonal Statistics as Table to calculate mean kg/km2/yr of P removed from each 12-digit HUC. Divided this value by 100 to convert to mean kg/ha/yr for each HUC. 2016-10-28T00:00:00 unknown 2029-07-23 This metadata was automatically generated from the FGDC Content Standard for Digital Geospatial Metadata standard (version FGDC-STD-001-1998) using the 2024-02-09T11:11:00 version of the FGDC CSDGM to ISO 19139 transform. EnviroAtlas Coordinator U.S. Environmental Protection Agency, Office of Research and Development-Sustainable and Healthy Communities Research Program, EnviroAtlas Geospatial Data Owner (919) 541-3832 109 T.W. Alexander Drive Research Triangle Park NC 27709 EnviroAtlas@epa.gov https://www.epa.gov/enviroatlas pointOfContact