Assessing the damages of La Palma Volcano
Last updated
Last updated
Since 11 September 2021, a swarm of seismic activity had been ongoing in the southern part of the Spanish Canary Island of La Palma (Cumbre Vieja region). The increasing frequency, magnitude, and shallowness of the seismic events were an indication of a pending volcanic eruption; which occurred on 16th September, leading to evacuation of people living in the vicinity.
In this tutorial we are going to assess the number of buildings and population that may get affected by the lava flow and its deposits. We’ll also estimate the value of damaged residential properties affected by the volcano eruption.
Access the Data Explorer section from your CARTO Workspace using the navigation menu.
In the Data Explorer page, navigate to CARTO Data Warehouse > demo_data > demo_table.
In this tutorial, we are going to use the following 3 tables:
lapalma_buildings: it contains the buildings in La Palma as obtained from the Spanish cadaster website;
lapalma_sociodemo_parcels: it contains a sample from Unica360’s dataset in the Data Observatory “Cadaster and Sociodemographics (Parcel)”;
lapalma_volcano_lavaflow: it includes the lava flow from the Volcano eruption in La Palma, Spain as measured by the Copernicus satellite on 10/04/2021.
Spend some time exploring the three tables in the Data Explorer.
Select lapalma_buildings and click on "Create map" button on the top.
This will open CARTO Builder with this table added as a layer to a map.
Rename the layer to “La Palma Buildings” and the map title to "Assessing the damages of La Palma Volcano"
Click on the layer to access the layer panel. In this section, you can style the layer according to your preferences. We have set the Fill Color to purple, reduce the opacity to 0,1. Then, we have set the Stroke Color to dark blue.
Let's add the lapalma_sociodemo_parcels source. To do so, follow the below steps:
Select the Add source from button at the bottom left on the page.
Click on the Data Explorer option.
Navigate to CARTO Data Warehouse > demo_data > demo_tables. Search for lapalma_sociodemo_parcels. Once you find it, select it and click on "Add Source".
Once added, a new layer appears on the map. Rename it to "La Palma demographics".
We'll now change the style of La Palma demographics layer. Access the layer panel and set the Fill Color to green and the Outline color to black. Also reduce the Stroke width to 1. Then, style the size of the points based on the population living in the parcel. To do so, select p_t column in the Radius section and set the range from 2 to 25.
Now, we are looking to analyse the number of buildings, their estimated values for residential properties and total population affected by the volcano lava extent. To perform this analysis, we'll use Workflows.
Go back to the Workspace tab in your browser and access Workflows.
In Workflows page, use the "New workflow" button to start a new Workflow. Select CARTO Data warehouse as the connection you want to work with.
From the Sources panel located on the left side, navigate to CARTO Data Warehouse > demo_data > demo_tables and locate lapalma_volcano_lavaflow. Drag and drop the source table into the canvas.
Repeat Step 13 to add lapalma_buildings into the canvas.
Now, use Enrich Polygons component to obtain the total of estimated property value of those residential properties affected by the lava flow as well as the total number of buildings affected. Connect lapalma_volcano_lavaflow as the target polygon and lapalma_buildings as the source. In the Variables section, in the node, add SUM for estimated_prop_value column and COUNT aggregation for numberOfBuildingUnits column. The output result is the lava flow source with the addition of the two new properties.
Add lapalma_sociodemo_parcels source to the canvas.
To obtain the total population affected by the lava flow extent, we will add the Enrich Polygons again. This time, we'll link lapalma_volcano_lavaflow as the target and lapalma_sociodemo_parcels as the source. Then, in the Variables section add SUM of p_t column.
Using the Join component, we'll join both Enriched Polygons output in a single table using the geoid as the common column. To achieve that, add the Join component to the canvas, use geoid as the common column for both sources and select Inner as the join type.
Save the output result as a new table using the Save as Table component. Set the destination to Organization > Private of your CARTO Data Warehouse and rename the output table to lapalma_volcano_lavaflow_enriched. Then, click on "Run".
Now, in the same Workflow, let's perform another analysis. This time, we are going to create a 500 meter buffer around the lava flow, and perform the same aggregations as we have done on Step 14 and Step 15 to compute the total number of buildings and the estimated damaged value of the residential properties within this larger region. To do so, add the Buffer component and link it to lapalma_volcano_lavaflow source. Set the distance to 500 meters. Then, click on "Run".
Afterwards, we'll add Enrich Polygons component, this time connecting the Buffer output as the target source. In the source input we'll connect lapalma_buildings source. Add the same aggregated variables: SUM for estimated_prop_values and COUNT for numberOfBuildingUnits. You can review the output in the Data Preview.
Let's add Enrich Polygons component again, this time to enrich the buffered output of La Palma lava flow with La Palma sociodemographics. In the Variable section of the Enrich Polygons component, add SUM for p_t to obtain the population affected by this buffered extent.
We'll add the Join component to join the output from both Enrich Polygons components. In the Join node, select geoid as the common column from both inputs and set the Join type to Inner.
Use the Select component to keep just the necessary columns using the below statement:
Finally, save the results as a table using the Save as Table component. Navigate to CARTO Data Warehouse > organization > private and save your table as lapalma_volcano_lavaflow_enriched_buffer.
Now let's go back to Builder. We'll first add lapalma_volcano_lavaflow_enriched as a table data source following the below steps:
Access Add source from..
Click on the Data Explorer option.
Navigate to CARTO Data Warehouse > organization > private. Search for lapalma_volcano_lavaflow_enriched. Once you find it, select it and click on "Add Source".
A new layer is added to the map. Rename it to "Lava flow" and move it to the bottom, just below La palma buildings layer.
Access Lava flow layer panel and set the Fill Color in the layer styling to light red.
Now let's add the enriched lava flow which was buffered by 500 meters. To do so, follow
these steps:
Access Add source from..
Click on the Data Explorer option.
Navigate to CARTO Data Warehouse > organization > private. Search for lapalma_volcano_lavaflow_enriched_buffer. Once you find it, select it and click on "Add Source".
Rename the recently added layer to 'Lava flow buffer' and move it to the bottom, just below Lava flow layer.
Set the layer style for Lava flow buffer to very light red. To do so, access the Layer panel and pick the color in the Fill Color section. Also, set the opacity in this section to 0.3 and disable the Stroke Color using the toggle button.
In the Interactions tab, enable interactions for both Lava flow and Lava flow buffer layers. For each column, set the right formatting and rename it to a user-friendly label.
Change the basemap to Google Terrain by navigating to the Basemap tag and selecting Terrain type.
Now, we can add a map description to provide further context about this map to our viewer users. You can use the below markdown description or add your own.
Finally we can make the map public and share the link to anybody in the organization. For that you should go to “Share” on the top right corner and set the map as Public. For more details, see Publishing and sharing maps.
Finally, we can visualize the result!
