Skip to content

Commit

Permalink
port marrakesh branch files over
Browse files Browse the repository at this point in the history
  • Loading branch information
nicolepaul committed Dec 21, 2023
1 parent 0ef3fd5 commit 6649993
Show file tree
Hide file tree
Showing 34 changed files with 12,262 additions and 0 deletions.
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
ID_0,NAME_0,FATALITIES,FATALITIES_GROUND_SHAKING,INJURIES,DISPLACED_POPULATION,AFFECTED_POPULATION,AFFECTED_UNITS,DAMAGED_UNITS,COLLAPSED_UNITS,UNIT_TYPES,ECONOMIC_LOSSES,INSURED_LOSSES,CURRENCY,INDUCED_EFFECTS,FLAG,REFERENCE,COMMENTS
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
ID_1,NAME_1,TAXONOMY,GEM_TAXONOMY,AFFECTED_UNITS,DAMAGED_UNITS,DS1_SLIGHT,DS2_MODERATE,DS3_EXTENSIVE,DS4_COMPLETE,COLLAPSED_UNITS,UNIT_TYPES,CAUSE_OF_DAMAGE,REFERENCE,COMMENTS
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
ID_1,NAME_1,FATALITIES,FATALITIES_GROUND_SHAKING,INJURIES,INJURIES_LIGHT,INJURIES_SEVERE,INJURIES_CRITICAL,HOMELESS,SHELTERED,DISPLACED_POPULATION,AFFECTED_POPULATION,CAUSE_OF_INJURY/DEATH,REFERENCE,COMMENTS
23 changes: 23 additions & 0 deletions Morocco/20230908_M6.9_MarrakeshSafi/Impact/README.md
Original file line number Diff line number Diff line change
@@ -0,0 +1,23 @@
# IMPACT DATA COLLECTION


Information at different administrative levels is collected for the consequences of the event,
which include building damage, economic losses, and human impact due to the earthquake and its induced effects.

The file `Impact_All_ID_0.csv` provides the national summary of consequences reported by multiple sources.

When information is available at finer resolutions, dedicated files for each administrative level
available are generated for the building and human impact.

- `Impact_Buildings_ID_n.csv`: Damaged building data at administrative level n.
- `Impact_Human_ID_n.csv`: Data of human impact information at administrative level n.

_Note: When available, building-level information is included._


## Metadata

Check the definitions of the parameters used in the impact files in the [impact metadata section](https://gitlab.openquake.org/risk/ecd/-/blob/main/metadata.md#impact-data).


## Notes
96 changes: 96 additions & 0 deletions Morocco/20230908_M6.9_MarrakeshSafi/OpenQuake_gmfs/README.md
Original file line number Diff line number Diff line change
@@ -0,0 +1,96 @@
# GROUND MOTION FIELDS USING OPENQUAKE

This folder includes the OpenQuake input files to generate scenario ground motion
fields (gmfs) considering different sources of information:

1. From rupture.
1. From rupture and conditioned to recording stations (with only seismic stations or with seismic and macroseismic stations).
1. From ShakeMap.

The following files are available:

- `stationlist_seismic.csv`: OpenQuake input file with past observations, considering only seismic stations. If no seismic stations are available for the event, the file is not reported.
- `stationlist_all.csv`: OpenQuake input file with past observations, considering seismic and macroseismic stations. If no macroseismic stations are available for the event, the file is not reported.
- `site_model_stations.csv`: OpenQuake input file with the site model (Vs30 values) at the location of the observations.
- `site_model.csv`: OpenQuake input file with the site model for which the ground shaking will be estimated.
- `site_model.png`: Figure of the site model (Vs30 based on USGS data).
- `gmpe_logic_tree_*.xml`: OpenQuake input file with a ground motion model logic tree (gmlt).
- `job_stations_none.ini`: A configuration file (`job.ini`) for running OpenQuake scenario calculations without conditioning the ground shaking to given observations.
- `job_stations.ini`: A configuration file (`job.ini`) for running OpenQuake calculations when conditioning the ground shaking to given observations.

_NOTE: the configuration files are examples for a given rupture geometry and a given gmlt. These parameters can be adjusted to the preferred models._
## Median ground motion fields

The following images present the median ground motion field for the combination of rupture and GMPE that generates the lowest `nominal bias` when conditioning the ground shaking to observations. The folder [Sensitivity](./Sensitivity/) provides the summary of the sensitivity analysis carried out with the available ruptures and GMPEs.

<img src="median_gmf_stations_none.png" height="350">
<img src="median_gmf_stations_all.png" height="350">

## Sensitivity analysis

The following information is available in the `Sensitivity` folder :

- `calculation_summary.csv`: Details on the calculations carried out in the sensitivity analysis.
- `log_calc_*.txt`: File with the OpenQuake running information in a _.txt_ format. This file includes the `nominal bias` estimates for each GMPE present in the logic tree model.
- `gmf_median_PGA_*.csv`: Generated ground motion fields for PGA for a given calculation ID.

### Metadata in file `calculation_summary.csv`

- `calc_id`: unique identifier for the calculation.
- `description`: details on the calculation
- `cal_time`: time of the calculation.
- `recording_stations`: "All" stations (seismic + macroseismic), or only "seismic" stations.
- `gmlt`: Ground Motion Logic Tree used in the calculation.
- `rupture`: Rupture model used in the calculation.
- `gmpe`: Ground Motion Prediction Equation (ground motion model).
- `imt`: Intensity measure type.
- `max_gmv`: Maximum ground motion value.
- `nominal_bias_mean`: A single value useful as a measure of the event bias for a given gmpe. It refers to the mean of the between-event residuals after conditioning the shaking to the recording stations.
- `nominal_bias_stdev`: Standard deviation for the `nominal_bias_mean`.
- `abs_bias`: Absolute value of the `nominal_bias_mean`.


## To run an OpenQuake calculation
### 1. From Rupture

Start a calculation from a given rupture in OpenQuake `xml` format and combine it with a
logic tree for the applicable GMPEs (ground motion prediction equations):

```
rupture + gmpe_logic_tree + job_stations_none.ini
```


### 2. From Rupture and conditional to recording stations

Condition the ground motion fields based on the recording stations available for the event
(`stationlist.csv` files). In this case, the calculation uses a given rupture in
OpenQuake `xml` format, the observed values in the recording stations, and a given ground
motion model.

The conditioning of the ground motion can consider:
- The recording stations with only seismic stations.
- All recording stations (with seismic and macroseismic values).

```
rupture + gmpe_logic_tree + stationlist_seismic.csv + job_stations_seismic.ini
rupture + gmpe_logic_tree + stationlist_all.csv + job_stations_all.ini
```


### 3. From Shakemap

The gmfs computed from the USGS ShakeMap can be generated with the command

```
oq shakemap2gmfs {shakemap_id} {sites.csv} -n 0 -s 0 -t 0 --spatialcorr yes -c yes
example of shakemap_id: us2000ar20
```

It is possible to specify the number of ground motion fields, the truncation level,
the spatial correlation, and the cross-correlation.

In OpenQuake, it is also possible to run scenario risk calculations starting from the USGS
ShakeMap.

Original file line number Diff line number Diff line change
@@ -0,0 +1,18 @@
calc_id,description,cal_time,recording_stations,gmlt,rupture,gmpe,imt,max_gmv,nominal_bias_mean,nominal_bias_stdev,abs_bias
454,"2023_M6.9_Marrakesh median gmf, Stations:All, gmlt:GEM, Rupture:USGS
",21, All, GEM, USGS,AtkinsonBoore2006Modified2011,PGA,0.5556112,0.311,0.101,0.311
454,"2023_M6.9_Marrakesh median gmf, Stations:All, gmlt:GEM, Rupture:USGS
",21, All, GEM, USGS,ChiouYoungs2014,PGA,0.5502717,0.395,0.085,0.395
454,"2023_M6.9_Marrakesh median gmf, Stations:All, gmlt:GEM, Rupture:USGS
",21, All, GEM, USGS,PezeshkEtAl2011NEHRPBC,PGA,0.651501,0.533,0.084,0.533
454,"2023_M6.9_Marrakesh median gmf, Stations:All, gmlt:GEM, Rupture:USGS
",21, All, GEM, USGS,AkkarEtAlRjb2014,PGA,0.8183377,0.714,0.103,0.714
447,"2023_M6.9_Marrakesh median gmf, Stations:None, gmlt:GEM, Rupture:USGS
",10, None, GEM, USGS,AkkarEtAlRjb2014,PGA,0.4993889,,,
447,"2023_M6.9_Marrakesh median gmf, Stations:None, gmlt:GEM, Rupture:USGS
",10, None, GEM, USGS,AtkinsonBoore2006Modified2011,PGA,0.37780565,,,
447,"2023_M6.9_Marrakesh median gmf, Stations:None, gmlt:GEM, Rupture:USGS
",10, None, GEM, USGS,ChiouYoungs2014,PGA,0.3595791,,,
447,"2023_M6.9_Marrakesh median gmf, Stations:None, gmlt:GEM, Rupture:USGS
",10, None, GEM, USGS,"NRCan15SiteTerm
gmpe_name = 'PezeshkEtAl2011NEHRPBC'",PGA,0.38041118,,,
Loading

0 comments on commit 6649993

Please sign in to comment.