[Shake-dev] Questions about current and future ShakeMap

[Bruce Worden]

Hi Jordi,

Thanks for contacting us. It’s helpful for us to know who is using ShakeMap, and it’s always good to get feedback on our developments. Your questions are good ones, my responses are inline below:

> On Sep 20, 2017, at 3:03 AM, Domingo Ballesta, Jordi (KNMI) <jordi.domingo.ballesta at knmi.nl> wrote:
> 
> Dear all,
> 
> At KNMI (Royal Netherlands Meteorological Institute) we are using ShakeMap for the induced earthquakes in Groningen, with a customised GMPE specific for the region.

For V4, it would be ideal if your custom GMPE were part of the GEM OpenQuake Hazardlib package. It’s not required, but will make things simpler.

> We have some extra requirements that ShakeMap doesn't fulfil, so we are thinking on customising the code. However, before doing that, we would like to know if any of these developments are planned for the version 4 of ShakeMap, or if there is any technical/scientific restriction for doing them.
> 
> Here our questions:
> 
> 1) We would like to use more than the current 3 periods for PSA (0.3, 1.0 and 3.0). Which restrictions should we consider for customising the code? Would it be possible in version 4 to configure a list of periods to compute?

Yes, in V4 you will be able to configure a list of output periods. If the output periods are in the input, the inputs will be used directly in the computation of the outputs. If not, the outputs will be computed by spectral interpolation from the periods available in the inputs. Mixed inputs (i.e., stations with various periods of SA) will be supported.

> 2) Sometimes we would like to obtain the output of the geometric mean of horizontals instead of the maximum value. The same with the unit, cm/s^2 instead of %g for PGA. Again, is there any restriction or will it be configurable in v4?

In V4, you will be able to specify an intensity measure component, as long as conversion equations are available to and from that component. Ideally, of course, the data would be provided with that component, but if not, the conversions can be applied. Currently, the main output file contains grids in units of log(g). I don’t know if it is worth making that configurable, since it is fairly easy to convert them when reading the grids.

Best,
Bruce Worden
USGS


> 
> Many thanks for your help.
> 
> Kind regards,
> 
> Jordi Domingo Ballesta
> Advisor IT and Data Management
> ............................................................
> KNMI
> Ministry of Infrastructure and the Environment
> Utrechtseweg 297 | 3731 GA De Bilt | A.0.06
> Postbus 201 | 3730 AE De Bilt
> ............................................................
> T +31 (0)30 22 06 670
> M +31 (0)6 17 017 546
> jordi.domingo.ballesta at knmi.nl
> http://rdsa.knmi.nl
> ............................................................
> 
> 
> _______________________________________________
> Shake-dev mailing list
> Shake-dev at geohazards.usgs.gov
> https://geohazards.usgs.gov/mailman/listinfo/shake-dev