Saturday, July 4, 2015

Week 7 - AppsGIS - Coastal Flooding

This week we worked with Coastal Flooding from sea level rise and from storm surge. For the sea level rise, we looked at two scenarios, one with a 3 foot rise, and one with a 6 foot rise. We then mapped the results and how it relates to population density. Later, for storm surge analysis, we compared two different DEM models.

In order to do look at the impact of sea level rise, we started with a DEM raster of the area, which we then wanted to extract only those cells that would be impacted by the associated rise in sea level. To this, the Reclassify tool was used, and the data was reclassified so that only the values of interest (to 3 feet or 6 feet) were included, all others were changed to "NoData". The resulting attribute table was looked at to determine the number of cells that are flooded (value of 1) and not flooded (value of 0). 
We then looked at the Properties of the Layer in order to determine that the cells of the raster are 9 m2. This was multiplied by the number of cells within each floodzone to determine the area.

To analyze the depth of flooding, I used the results from the Reclassify Tool as the input in the Times tool, in order to create a new raster of only the floodzone elevations. In order to get the flooded depth, used the Minus Tool, with either the equivalent of either 3 feet or 6 feet, and the results from the Times Tool as the second input. These results were then mapped against the population density of the census tracts (this map is only for a rise of 6 feet):
Figure 1. Map of District of Honolulu showing impact of 6 foot sea level rise, and
includes population densities of the region.

We then wanted to look at the social vulnerability of the are, with data from the 2010 Census for the area. In order to do this, we first had to determine which census blocks were located within the floodzones. For our analysis, we chose to select those whose centroid was located in the floodzone, using the Select By Location tool. Looking at the Attribute Table, we were able to determine how many blocks were selected, and what population this affects. This can be done simply looking at the Statistics for the column of choice, as only the selected rows are summarized.

Next, we had to add fields in the Census tracts layer for each of the groups of interest, percent white residents, percent owner-occupied homes and percent homes with people over the age of 65. Table joins were then conducted in order to copy over the information. Each join was removed prior to joining another. The Field Calculator was used to fill in the data. 

This was then repeated for the Census Blocks layer, with additional fields added for the population of white residents, owner-occupied homes, and those over the age of 65. For our analysis, the census blocks did not include the make-up of each block, so the census tract data was used. For this, it was assumed that the population composition for the whole was equal for the part. A Table Join was created so that the percentage could be copied over. This percentage was then used to determine the size of the population for each of the three social statuses. Used the Select by Location tool to select the blocks that had their centroid located within the floodzone, as above, for both 3 and 6 feet. Then used the Statistics function to get the sum for each population in order to fill out the table in Deliverable 7. To get the values for the nonflooded areas, simply switched the selection. Did this for all of the variables in the table. Then divided the populations for each variable by the total population for that category (3 feet flooded, 3 feet not-flooded, etc).

The results were as follows:

Variable
Entire District
3 Feet Scenario
6 Feet Scenario


Flooded
Not-flooded
Flooded
Not-flooded
Total Population
1,360,301
8,544
1,351,757
60,005
1,300,296
% White
24.7 %
36.8 %
24.7 %
29.6 %
24.5 %
% Owner-occupied
58.2 %
32.2 %
58.3 %
38.1 %
59.1 %
% 65 and older
14.3 %
17.11 %
14.3 %
17.0 %
14.2 %
Table 1. Percent of population represented by each group in flooded areas compared to not-flooded areas. Results are for 3 and 6 feet of sea level rise.
After this analysis, we looked at storm surge in Collier County, Florida. The purpose for this analysis was to compare the results of two different DEMs, one by USGS created using older methods (and lower resolution), and one created using Lidar. In order to compare the two, percent error in omission (not including data that should be) and percent error in commission (false positives). For our analysis, the Lidar data was treated as accurate for the calculations. The results were quite different, and showed the major benefits of Lidar techniques, with most errors of commission being above 100%.

Being a resident of Florida, this week's analysis was quite intriguing and interesting to do. I especially liked being able to look at how seeming small levels of sea level rise can affect such large areas. I certainly used the NOAA Sea Level Viewer to see if my house was in danger, especially since our yard floods quite a bit in heavy rains. I expect to use this new knowledge quite a bit in my future.

No comments:

Post a Comment