20220608 12460800 1.0 TRUE TRUE Vertical_Displacement_SJV_DWR_1949_to_2005 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 extract some of the riverine and hilly areas that produced more extreme values and other artifacts purely due to the analysis approach (i.e. not actual real elevation change). * - some of the areas with more than 20 feet of subsidence were omitted from this clipping, because they were in heavily subsided areas and may be “real subsidence.”The resulting subsidence product should be perceived in light of the above. Some of the collar of the San Joaquin Valley shows large changes, but that is simply due to the analysis method. Also, individual grid cells may or may not be comparing the same real features. Errors are baked into both comparison datasets. However, it is important to note that the large areas of subsidence in the primary agriculture area agree fairly well with a cruder USGS subsidence map of the Valley based on extensometer data. We have confidence that the big picture story these results show us is largely correct, and that the magnitudes of subsidence are somewhat reasonable. The contour set can serve as the baseline to support future comparisons using more recent or future data as it becomes available. It should be noted there are two key versions of the data. The “Final Deliverables” from 2/2017 were delivered to support the initial Public Affairs press release. Subsequent improvements were made in coverage and blunder correction as time permitted (it should be noted this occurred in the midst of the Oroville Dam emergency) to produce the final as of 3/12/2017. Further improvements in overall quality and filtering could occur in the future if time and needs demand it. Update (4/3/2019, Ben Brezing): The raster was further "smoothed" to remove artifacts that result from comparing the high resolution NextMAP DEM to the lower resolution DEM that was derived from the 1950’s quad map contours. The smoothing was accomplished by removing raster cells with values that are more than 0.5 feet different than adjacent cells (25 meter cell size), as well as the adjacent cells. The resulting raster was then resampled to a raster with 100 meter cell size using cubic resampling technique and was then converted to a point feature class. The point feature class was then interpolated to a raster with 250 meter cell size using the IDW technique, a fixed search radius of 1250 meters and power=2. The resulting raster was clipped to a smaller extent to remove "noisier" areas around the edges of the Central Valley while retaining coverage for the main area of interest. subsidence San Joaquin Valley elevation California DWR DWR Assessment of subsidence in the San Joaquin Valley between the 1940s/early 1950s and 2005. California Department of Water Resources, Geodetic Branch /9j/4AAQSkZJRgABAQEAAAAAAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0a HBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIy MjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCADIASwDAREA AhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQA AAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3 ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWm p6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEA AwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSEx BhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElK U1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3 uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+gAo AKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgA oAKACgAoAKACgAoAKACgAoAKACgAoAKAIricW8JkZWb5lUKuMkkgAc+5FAEX2ub/AKB1z/31H/8A F0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31 H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1 z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/ AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AH2ub/AKB1z/31H/8AF0AC3redFHJaTxeYxVWYoRnB PZiegNAFqgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKAKmo/wDHsn/XeH/0YtAF ugAoAKACgAoAKACgAoAKACgAoAKACgAoAKAKt3/x82H/AF3P/ot6ALVABQAUAFABQAUAFABQAUAF ABQAUAFABQAUAFABQAUAFABQAUAVNR/49k/67w/+jFoAt0AFABQAUAFABQAUAFABQAUAFABQAUAF ABQBVu/+Pmw/67n/ANFvQBaoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgCpqP/ AB7J/wBd4f8A0YtAFugAoAKACgAoAKACgAoAKACgAoAKACgAoAKAKt3/AMfNh/13P/ot6ALVABQA UAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAVNR/49k/67w/+jFoAt0AFABQAUAFABQA UAFABQAUAFABQAUAFABQBVu/+Pmw/wCu5/8ARb0AWqACgAoAKACgAoAKACgAoAKACgAoAKACgAoA KACgAoAKACgAoAqaj/x7J/13h/8ARi0AW6ACgAoAY8scZw7gH3rkxGPw2Gko1pqLfdlKMnshv2mH /nov51z/ANs5f/z+j94/Zz7EgIIBByD0r0YTjOKnF3T1RD0FqgCgAoAKACgAoAKACgAoAq3f/HzY f9dz/wCi3oAtUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQBU1H/j2T/rvD/wCj FoAt0AFABQBVu4FdGlJO5V4x0r5niLKaOIpTxk2+aEdLWtpd66efc2o1GnymZX5mdpahvWjRUKgq Pzr6nLuJ6uEpQoOCcY+t7GE6Kk27mkGB6EGv0WNWE/hkmclmLWggoAKACgAoAKACgAoAq3f/AB82 H/Xc/wDot6ALVABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAVNR/wCPZP8ArvD/ AOjFoAt0AFABQBHMpeB1UZJHFcOZ0Z1sHUpU1eTTsVBpSTZjdDg1+ONOLs9z0TFvJr6HXIwZVFky Daq9Q2ec+oIIrtpxpSoPT3v0PQowoyw709/9DctrgoBKgUkj61eX5hUy+s6tOKb21vp9zR5tSnf3 WbdfsJ5wUAFABQAUAFABQAUAVbv/AI+bD/ruf/Rb0AWqACgAoAKACgAoAKACgAoAKACgAoAKACgA oAKACgAoAKACgAoAqaj/AMeyf9d4f/Ri0AW6ACgAoAKAM3Ubcj96hxkjPHevgeKcq9nP67TWj+Je ff59fP1OvD1Pssx5JlMyCZUyHIUdcAdD7HOK+SUXyvlPQjB8r5exLCAkpVZiVwcIR0555/Golqrt ETu1do0ra6kWQKzFgxHJ5xX0eTZ9iqVdUqsnKMmld3bXpqt/+DZ7HHUpRaujTr9LOMKACgAoAKAC gAoAq3f/AB82H/Xc/wDot6ALVABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAVNR /wCPZP8ArvD/AOjFoAt0AFABQAUAR3EK3FvJEwyrjFcmOwyxWHnRfVfjuvxLpzcJKS6HiniW+1jR /G4QpL9mdwI42xtkBIJ29Oegr4/DZfh3l0/a6VYtJrto9/6sff5dSw+JwVla9m2+qa2v5f10O9tr hii3BAHmOw2huPXv36mvk6kFfk7HztSmruHY04pVcB43DDPBHNYRlOlNSWjWpxyi1pJFn7Zcf89P 0Fez/rJmn/P3/wAlj/kY+xh2Hx30it853D6V2YPirGU5/wC0PnWnRKy62slr6kyoRexpKwdQynIP Sv0ShXhXpqrTd4vY5GmnZi1qIKACgAoAq3f/AB82H/Xc/wDot6ALVABQAUAFABQAUAFABQAUAFAB QAUAFABQAUAFABQAUAFABQAUAVNR/wCPZP8ArvD/AOjFoAt0AFABQAUAFAFHU9Js9WtmguoVYHo2 0EqfUZrgxmXUcUveVpd1v9/byOnDYqrhp89NnD2v2/TpJbO/ZibdmCyFsfaYunX+8Mcjr371+c4+ hRjWlGl0bT8mn+T/AOAfR1PY10qtL7XT+V/5P/gGlBc2VkivGZGDg5bBbbjnB9O9edOFSo7Pock6 dWq2pdDSiuIZ8+VKrEdcHpXNKEo7o5JU5w+JEtQQa9r/AMesf0r9ayD/AJFtL0/VnBV+Nk1ewZhQ AUAFAFW7/wCPmw/67n/0W9AFqgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKAKmo /wDHsn/XeH/0YtAFugAoAKACgAoAKAOZ8bWRudCuGhUGdULRnphh0OfWvj+IqFKliaOIStzO0vNe Z7GS1lDExUtr6+nU4DwbqKX9rc2UtyXlSPcjNjlSvTHtyOa8LM6LpSjUitG/1PqM4w7ozjVjGyb1 +/8AU426vb/SvEjTszxPHJkBDtDJnsPQgCvqsty7D5jhZUoyXO17qbtdvzatddrpu57ElSnhVG14 tavez/M9X07xVZzaSt7K5VXxjdx1OPw5FfDYvKq9DESw8laUdz46rltT2vJDU7DTb1biNFRRt25D A9a+s4czZz5cA4W5U9b+fa3n3PCxNFwbbNCvsDlCgAoAKAKt3/x82H/Xc/8Aot6ALVABQAUAFABQ AUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAVNR/49k/67w/8AoxaALdABQAUAFABQAUAIyqww ygj3FZ1KVOorVIp+quNNrY46fwTYWOtR6tZbom5WRVPDA9z718RnuExGCwns1LmpN9UuZP17f1Y9 2Gb1a1F0KuvYZqXh60voivlR9MbXXIP+ePyr5ShjKlN3uaYfH1aMr3Z5XrUK6Jqn2O5iJhPPklwf Jw+QV44BH48+1fXYPnxkOenrLvbV6fmfZ4Op9Zoe1g7Nde+nXXoeueE9Ta9EKlVKGLKup4IwOlYc PUPZ5k+9n+aPhs1wypXfmdbX354YUAFABQBVu/8Aj5sP+u5/9FvQBaoAKACgAoAKACgAoAKACgAo AKACgAoAKACgAoAKACgAoAKACgCpqP8Ax7J/13h/9GLQBboAKACgAoAKACgAoAZKgeJlIzkdK5Mf h44jDTpuN7p6efT8SouzTMea3uUhdlj+YA4zX5bUyXHUYupVptRW+3+Z3wqU3JJs4W90SLUJwLyJ bi4kkILynkLnO0YxxkYHt+OeqhjZ4f3qMuVJdO9tz6SljHSg1B2ilt09df6/Tu9HtYbQxwwRLHHG mFVRgDtXVwtUlPMnKTu3F/ofNYurOpeU3ds2q/STzwoAKACgCrd/8fNh/wBdz/6LegC1QAUAFABQ AUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFAFTUf+PZP+u8P/oxaALdABQAUAFABQAUAFABQ BFc820n+7XmZz/yL63+Fl0/jRgtawPKspjG9ehHFfkaqSStc9NVZpct9DT07Zuc5G/oB7V9pwfGg pVJOX7x6JeW7t3/Sxx4i9l2NCvuzlCgAoAKAKt3/AMfNh/13P/ot6ALVABQAUAFABQAUAFABQAUA FABQAUAFABQAUAFABQAUAFABQAUAVNR/49k/67w/+jFoAt0AFABQAUAFABQAUAFADZE8yNkzjIxm ubG4b6zh50L25la5UXyu5S/s3/pr/wCO/wD16+P/ANS/+n//AJL/APbG/wBZ8iAk2NznaZMDjtni vISlkOZarnSWn2b3Xz2enyNV++h2JV1TeuViHXkb/wD61ez/AK5e8k6Nl/i/4CIeFs9WaCncob1G a+zpVFVpxqLqk/vOVqzsLWggoAq3f/HzYf8AXc/+i3oAtUAFABQAUAFABQAUAFABQAUAFABQAUAF ABQAUAFABQAUAFABQBU1H/j2T/rvD/6MWgC3QAUAFABQAUAFABQAUAFABQBy2t3t7a6m2y38yDAG 4AnHy5yfQZr804ghCpmNRSlquX8lse1gqNKpRV5Wf/BKmk291JcNcXEhVJGDDHcV5MvYznGley2b 7G+LqU1FQgtUdmqhECjoBgV+uUKUaNKNKG0UkvlofPt3dxa1EFAFW7/4+bD/AK7n/wBFvQBaoAKA CgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgCpqP/Hsn/XeH/wBGLQBboAKACgAoAKAC gAoAKACgAoAY0MbMSyKSeuRXHVy/CVpupUpxbfVpdClOS0TG/Z4f+eSflUVMqwVSXNOlFv0H7Sfc lAAAAGAK7YxjCKjFWSICqAKAKt3/AMfNh/13P/ot6ALVABQAUAFABQAUAFABQAUAFABQAUAFABQA UAFABQAUAFABQAUAVNR/49k/67w/+jFoAt0AFABQAUAFABQAUAFABQAUAFABQAUAFABQBVu/+Pmw /wCu5/8ARb0AWqACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAqaj/x7J/13h/8A Ri0AW6ACgAoAKACgAoAKACgAoAKACgAoAKACgAoAq3f/AB82H/Xc/wDot6ALVABQAUAFABQAUAFA BQAUAFABQAUAFABQAUAFABQAUAFABQAUAVNR/wCPZP8ArvD/AOjFoAt0AFABQAUAFABQAUAFABQA UAFABQAUAFABQBVu/wDj5sP+u5/9FvQBaoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAo AKACgCpqP/Hsn/XeH/0YtAFugAoAKACgAoAKACgAoAKACgAoAKACgAoAKAKt3/x82H/Xc/8Aot6A LVABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUAFABQAUARXMAuYfLLsnzKwZcZBBBHUEdRQ BD9kn/6CNz/3zH/8RQAfZJ/+gjc/98x//EUAH2Sf/oI3P/fMf/xFAB9kn/6CNz/3zH/8RQAfZJ/+ gjc/98x//EUAH2Sf/oI3P/fMf/xFAB9kn/6CNz/3zH/8RQAfZJ/+gjc/98x//EUAH2Sf/oI3P/fM f/xFAB9kn/6CNz/3zH/8RQAfZJ/+gjc/98x//EUAH2Sf/oI3P/fMf/xFAB9kn/6CNz/3zH/8RQAf ZJ/+gjc/98x//EUAH2Sf/oI3P/fMf/xFAALJvPikkvJ5fKYsqsEAztK84UHoTQBboAKACgAoAKAC gAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKACgAoAKA CgAoAKACgAoAKACgAoAKACgAoA//2Q==