San Joaquin Valley Subsidence Analysis README. Written: Joel Dudas, 3/12/2017. Amended: Ben Brezing, 4/2/2019. DWR’s Division of Engineering Geodetic Branch received a request in 1/2017 from Jeanine Jones to produce a graphic of historic subsidence in the entirety of the San Joaquin Valley. The task was assigned to the Mapping & Photogrammetry Office and the Geospatial Data Support Section to complete by early February. After reviewing the alternatives, the decision was made to produce contours from the oldest available set of quad maps for which there was reasonable certainty about quality and datum, and to compare that to the most current Valley-wide DEM. For the first requirement, research indicated that the 1950’s vintage quad maps for the Valley were the best alternative. Prior quad map editions are uneven in quality and vintage, and the actual control used for the contour lines was extremely suspect. The 1950’s quads, by contrast, were produced primarily on the basis of 1948-1949 aerial photography, along with control corresponding to that period, and referenced to the National Geodetic Vertical Datum of 1929. For the current set, the most recent Valley-wide dataset that was freely available, in the public domain, and of reasonable accuracy was the 2005 NextMap SAR acquisition (referenced to NAVD88). The primary bulk of the work focused on digitizing the 1950’s contours. First, all of the necessary quads were downloaded from the online USGS quad source https://ngmdb.usgs.gov/maps/Topoview/viewer/#4/41.13/-107.51. Then the entire staff of the Mapping & Photogrammetry Lab (including both the Mapping Office and GDDS staff) proceeded to digitize the contours. Given the short turnaround time constraint and limited budget, certain shortcuts occurred in contour development. While efforts were made to digitize accurately, speed really was important. Contours were primarily focused only on agricultural and other lowland areas, and so highlands were by and large skipped. The tight details of contours along rivers, levees, and hillsides was skipped and/or simplified. In some cases, only major contours were digitized. The mapping on the source quads itself varied….in a few cases on spot elevations on benchmarks were available in quads. The contour interval sometimes varied, even within the quad sheet itself. In addition, because 8 different people were creating the contours, variability exists in the style and attention to detail. It should be understood that given the purpose of the project (display regional subsidence patterns), that literal and precise development of the historic contour sets leaves some things to be desired. These caveats being said, the linework is reasonably accurate for what it is (particularly given that the contours of that era themselves were mapped at an unknown and varying actual quality). The digitizers tagged the lines with Z values manually entered after linework that corresponded to the mapped elevation contours. Joel Dudas then did what could be called a “rough” QA/QC of the contours. The individual lines were stitched together into a single contour set, and exported to an elevation raster (using TopoToRaster in ArcGIS 10.4). Gross blunders in Z values were corrected. Gaps in the coverage were filled. The elevation grid was then adjusted to NAVD88 using a single adjustment for the entire coverage area (2.5’, which is a pretty close average of values in this region). The NextMap data was extracted for the area, and converted into feet. The two raster sets were fixed to the same origin point. The subsidence grid was then created by subtracting the old contour-derived grid from the NextMAP DEM. The subsidence grid that includes all of the values has the suffix “ALL”. Then, to improve the display fidelity, some of the extreme values (above +5’ and below -20’*) were filtered out of the dataset, and the subsidence grid was regenerated for these areas and suffixed with “cut.” The purpose of this cut was to