openquake.hazardlib.geo.surface package

Submodules

openquake.hazardlib.geo.surface.base module

Module openquake.hazardlib.geo.surface.base implements BaseSurface and BaseSurface.

class openquake.hazardlib.geo.surface.base.BaseSurface(mesh=None)

Bases: object

Base class for a surface in 3D-space.

get_area()

Compute area as the sum of the mesh cells area values.

get_azimuth(mesh)

This method computes the azimuth of a set of points in a openquake.hazardlib.geo.mesh instance. The reference used for the calculation of azimuth is the middle point and the strike of the rupture. The value of azimuth computed corresponds to the angle measured in a clockwise direction from the strike of the rupture.

Parameters:

mesh – An instance of openquake.hazardlib.geo.mesh

Returns:

An instance of numpy.ndarray

get_azimuth_of_closest_point(mesh)

Compute the azimuth between point in mesh and the corresponding closest point on the rupture surface.

Parameters:

mesh – An instance of openquake.hazardlib.geo.mesh

Returns:

An numpy.ndarray instance with the azimuth values.

get_bounding_box()

Compute surface bounding box from surface mesh representation. That is extract longitudes and latitudes of mesh points and calls: openquake.hazardlib.geo.utils.get_spherical_bounding_box()

Returns:

A tuple of four items. These items represent western, eastern, northern and southern borders of the bounding box respectively. Values are floats in decimal degrees.

get_closest_points(mesh)

For each point from mesh find a closest point belonging to surface.

Parameters:

mesh – Mesh of points to find closest points to.

Returns:

Mesh of the same shape as mesh with closest surface’s points on respective indices.

get_hypo_location(mesh_spacing, hypo_loc=None)

The method determines the location of the hypocentre within the rupture

Parameters:

• mesh – Mesh of points

• mesh_spacing – The desired distance between two adjacent points in source’s ruptures’ mesh, in km. Mainly this parameter allows to balance the trade-off between time needed to compute the distance between the rupture surface and a site and the precision of that computation.

• hypo_loc – Hypocentre location as fraction of rupture plane, as a tuple of (Along Strike, Down Dip), e.g. a hypocentre located in the centroid of the rupture would be input as (0.5, 0.5), whereas a hypocentre located in a position 3/4 along the length, and 1/4 of the way down dip of the rupture plane would be entered as (0.75, 0.25).

Returns:

Hypocentre location as instance of Point

get_joyner_boore_distance(mesh)

Compute and return Joyner-Boore (also known as Rjb) distance to each point of mesh.

Parameters:

mesh – Mesh of points to calculate Joyner-Boore distance to.

Returns:

Numpy array of closest distances between the projections of surface and each point of the mesh to the earth surface.

get_middle_point()

Compute coordinates of surface middle point.

The actual definition of middle point depends on the type of surface geometry.

Returns:

instance of openquake.hazardlib.geo.point.Point representing surface middle point.

get_min_distance(mesh)

Compute and return the minimum distance from the surface to each point of mesh. This distance is sometimes called Rrup.

Parameters:

mesh – Mesh of points to calculate minimum distance to.

Returns:

A numpy array of distances in km.

get_resampled_top_edge(angle_var=3.0)

This methods computes a simplified representation of a fault top edge by removing the points that are not describing a change of direction, provided a certain tolerance angle.

Parameters:

angle_var (float) – Number representing the maximum deviation (in degrees) admitted without the creation of a new segment

Returns:

A Line representing the rupture surface’s top edge.

get_rx_distance(mesh)

Compute distance between each point of mesh and surface’s great circle arc.

Distance is measured perpendicular to the rupture strike, from the surface projection of the updip edge of the rupture, with the down dip direction being positive (this distance is usually called Rx).

In other words, is the horizontal distance to top edge of rupture measured perpendicular to the strike. Values on the hanging wall are positive, values on the footwall are negative.

Parameters:

mesh – Mesh of points to calculate Rx-distance to.

Returns:

Numpy array of distances in km.

get_ry0_distance(mesh)

Compute the minimum distance between each point of a mesh and the great circle arcs perpendicular to the average strike direction of the fault trace and passing through the end-points of the trace.

Parameters:

mesh – Mesh of points to calculate Ry0-distance to.

Returns:

Numpy array of distances in km.

get_surface_boundaries()

Returns the boundaries in the same format as a multiplanar surface, with two lists of lons and lats

get_surface_boundaries_3d()

Returns the boundaries as three lists of lons, lats, depths

get_top_edge_centroid()

Return Point representing the surface’s top edge centroid.

get_top_edge_depth()

Return minimum depth of surface’s top edge.

Returns:

Float value, the vertical distance between the earth surface and the shallowest point in surface’s top edge in km.

abstract get_width()

openquake.hazardlib.geo.surface.base.downsample_mesh(mesh, tol=1.0)

Returns a mesh sampled at a lower resolution - if the difference in azimuth is larger than the specified tolerance a turn is assumed

Returns:

Downsampled mesh as instance of :class: openquake.hazardlib.geo.mesh.RectangularMesh

openquake.hazardlib.geo.surface.base.downsample_trace(mesh, tol=1.0)

Downsamples the upper edge of a fault within a rectangular mesh, retaining node points only if changes in direction on the order of tol are found

Returns:

Downsampled edge as a numpy array of [long, lat, depth]

openquake.hazardlib.geo.surface.base.surface_to_arrays(surface)

Parameters:

surface – a (Multi)Surface object or a list of simple surfaces

Returns:

a list of S arrays of shape (3, Sy, Sz)

openquake.hazardlib.geo.surface.base.to_arrays(geom)

Parameters:

geom – an array [num_surfaces, shape_y, shape_z …, coords]

Returns:

a list of num_surfaces arrays with shape (3, shape_y, shape_z)

openquake.hazardlib.geo.surface.base.to_geom_lons_lats(surface)

Convert a (multi)surface into a float32 array

Returns:

geom, lons, lats

openquake.hazardlib.geo.surface.complex_fault module

Module openquake.hazardlib.geo.surface.complex_fault defines ComplexFaultSurface.

class openquake.hazardlib.geo.surface.complex_fault.ComplexFaultSurface(mesh)

Bases: BaseSurface

Represent a complex fault surface as 3D mesh of points (not necessarily uniformly spaced across the surface area).

Parameters:

mesh – Instance of RectangularMesh representing surface geometry.

Another way to construct the surface object is to call from_fault_data().

classmethod check_aki_richards_convention(edges)

Verify that surface conforms with Aki and Richard convention (i.e. surface dips right of surface strike) Test with 2 adjacent edges to allow for very large, curved surfaces

This method doesn’t have to be called by hands before creating the surface object, because it is called from from_fault_data().

classmethod check_fault_data(edges, mesh_spacing)

Verify the fault data and raise ValueError if anything is wrong.

This method doesn’t have to be called by hands before creating the surface object, because it is called from from_fault_data().

check_proj_polygon()

A common user error is to create a ComplexFaultSourceSurface from invalid fault data (e.g. mixing the order of vertexes for top and bottom edges). Therefore, we want to restrict every complex source to have a projected enclosing polygon that is not a multipolygon.

classmethod check_surface_validity(edges)

Check validity of the surface.

Project edge points to vertical plane anchored to surface upper left edge and with strike equal to top edge strike. Check that resulting polygon is valid.

This method doesn’t have to be called by hands before creating the surface object, because it is called from from_fault_data().

classmethod from_fault_data(edges, mesh_spacing)

Create and return a fault surface using fault source data.

Parameters:

• edges – A list of at least two horizontal edges of the surface as instances of openquake.hazardlib.geo.line.Line. The list should be in top-to-bottom order (the shallowest edge first).

• mesh_spacing – Distance between two subsequent points in a mesh, in km.

Returns:

An instance of ComplexFaultSurface created using that data.

Raises:

ValueError – If requested mesh spacing is too big for the surface geometry (doesn’t allow to put a single mesh cell along length and/or width).

Uses check_fault_data() for checking parameters.

get_dip()

Return the fault dip as the average dip over the mesh.

The average dip is defined as the weighted mean inclination of all the mesh cells. See openquake.hazardlib.geo.mesh.RectangularMesh.get_mean_inclination_and_azimuth()

Returns:

The average dip, in decimal degrees.

get_strike()

Return the fault strike as the average strike over the mesh.

The average strike is defined as the weighted mean azimuth of all the mesh cells. See openquake.hazardlib.geo.mesh.RectangularMesh.get_mean_inclination_and_azimuth()

Returns:

The average strike, in decimal degrees.

get_width()

Return surface’s width (that is surface extension along the dip direction) in km.

The width is computed as the average width along the surface. See openquake.hazardlib.geo.mesh.RectangularMesh.get_mean_width()

classmethod surface_projection_from_fault_data(edges)

Get a surface projection of the complex fault surface.

Parameters:

edges – A list of horizontal edges of the surface as instances of openquake.hazardlib.geo.line.Line.

Returns:

Instance of Polygon describing the surface projection of the complex fault.

property tor

Returns:

top of rupture line

openquake.hazardlib.geo.surface.complex_fault.complex_fault_node(edges)

Parameters:

edges – a list of lists of points

Returns:

a Node of kind complexFaultGeometry

openquake.hazardlib.geo.surface.complex_fault.edge_node(name, points)

Parameters:

• name – ‘faultTopEdge’, ‘intermediateEdge’ or ‘faultBottomEdge’

• points – a list of Point objects

Returns:

a Node of kind faultTopEdge, intermediateEdge or faultBottomEdge

openquake.hazardlib.geo.surface.gridded module

Module openquake.hazardlib.geo.surface.gridded defines GriddedSurface.

class openquake.hazardlib.geo.surface.gridded.GriddedSurface(mesh=None)

Bases: BaseSurface

Gridded surface defined by an unstructured cloud of points. This surface type is required for a proper implementation of some subduction interface surfaces included int the Japan 2012 model.

Note that currently we support only one rupture-site typology i.e. since this the only one that can be unambiguosly computed.

Parameters:

mesh – An unstructured mesh of points ideally representing a rupture surface. Must be an instance of Mesh

classmethod from_points_list(points)

Create a gridded surface from a list of points.

Parameters:

points – A list of Point

Returns:

An instance of GriddedSurface

get_area()

Compute surface’s area in squared km.

Returns:

Float value, the surface area

get_bounding_box()

Compute surface geographical bounding box.

Returns:

A tuple of four items. These items represent western, eastern, northern and southern borders of the bounding box respectively. Values are floats in decimal degrees.

get_dip()

Compute surface’s dip as decimal degrees in a range (0, 90].

The actual definition of the dip might depend on surface geometry.

Returns:

The dip angle in degrees

get_joyner_boore_distance(mesh)

Compute and return Joyner-Boore (also known as Rjb) distance to each point of mesh.

Parameters:

mesh – Mesh of points to calculate Joyner-Boore distance to.

Returns:

Numpy array of closest distances between the projections of surface and each point of the mesh to the earth surface.

get_middle_point()

Compute coordinates of surface middle point.

The actual definition of middle point depends on the type of surface geometry.

Returns:

instance of openquake.hazardlib.geo.point.Point representing surface middle point.

get_rx_distance(mesh)

Compute distance between each point of mesh and surface’s great circle arc.

Distance is measured perpendicular to the rupture strike, from the surface projection of the updip edge of the rupture, with the down dip direction being positive (this distance is usually called Rx).

In other words, is the horizontal distance to top edge of rupture measured perpendicular to the strike. Values on the hanging wall are positive, values on the footwall are negative.

Parameters:

mesh – Mesh of points to calculate Rx-distance to.

Returns:

Numpy array of distances in km.

get_ry0_distance(mesh)

Parameters:

mesh – Mesh of points

get_strike()

Compute surface’s strike as decimal degrees in a range [0, 360).

The actual definition of the strike might depend on surface geometry.

Returns:

The strike angle in degrees

get_surface_boundaries()

Returns:

(lons, lats) for the 5 points of the bounding box

get_surface_boundaries_3d()

Returns:

(lons, lats, depths) for the 5 points of the bounding box

get_top_edge_depth()

Compute minimum depth of surface’s top edge.

Returns:

Float value, the vertical distance between the earth surface and the shallowest point in surface’s top edge in km.

get_width()

Compute surface’s width (that is surface extension along the dip direction) in km.

The actual definition depends on the type of surface geometry.

Returns:

Float value, the surface width

property surface_nodes

Parameters:

points – a list of Point objects

Returns:

a Node of kind ‘griddedSurface’

openquake.hazardlib.geo.surface.multi module

Module openquake.hazardlib.geo.surface.multi defines MultiSurface.

class openquake.hazardlib.geo.surface.multi.MultiSurface(surfaces, msparam=None)

Bases: BaseSurface

Represent a surface as a collection of independent surface elements.

Parameters:

surfaces – List of instances of subclasses of BaseSurface each representing a surface geometry element.

classmethod from_csv(fname: str)

Parameters:

fname – path to a CSV file with header (lon, lat, dep) and 4 x P rows describing planes in terms of corner points in the order topleft, topright, bottomright, bottomleft

Returns:

a MultiSurface made of P planar surfaces

get_area()

Return sum of surface elements areas (in squared km).

get_bounding_box()

Compute bounding box for each surface element, and then return the bounding box of all surface elements’ bounding boxes.

Returns:

A tuple of four items. These items represent western, eastern, northern and southern borders of the bounding box respectively. Values are floats in decimal degrees.

get_dip()

Compute dip of each surface element and return area-weighted average value (in range (0, 90]). Given that dip values are constrained in the range (0, 90], the simple formula for weighted mean is used.

get_joyner_boore_distance(mesh)

For each point in mesh compute the Joyner-Boore distance to all the surface elements and return the smallest value. See superclass method for spec of input and result values.

get_middle_point()

If MultiSurface is defined by a single surface, simply returns surface’s middle point, otherwise find surface element closest to the surface’s bounding box centroid and return corresponding middle point. Note that the concept of middle point for a multi surface is ambiguous and alternative definitions may be possible. However, this method is mostly used to define the hypocenter location for ruptures described by a multi surface (see openquake.hazardlib.source.characteristic.CharacteristicFaultSource.iter_ruptures()). This is needed because when creating fault based sources, the rupture’s hypocenter locations are not explicitly defined, and therefore an automated way to define them is required.

get_min_distance(mesh)

For each point in mesh compute the minimum distance to each surface element and return the smallest value. See superclass method for spec of input and result values.

get_rx_distance(mesh)

Parameters:

mesh – An instance of openquake.hazardlib.geo.mesh.Mesh with the coordinates of the sites.

Returns:

A numpy.ndarray instance with the Rx distance. Note that the Rx distance is directly taken from the GC2 t-coordinate.

get_ry0_distance(mesh)

Parameters:

mesh – An instance of openquake.hazardlib.geo.mesh.Mesh with the coordinates of the sites.

get_strike()

Compute strike of each surface element and return area-weighted average value (in range [0, 360]) using formula from: http://en.wikipedia.org/wiki/Mean_of_circular_quantities Note that the original formula has been adapted to compute a weighted rather than arithmetic mean.

get_surface_boundaries()

Returns the boundaries in the same format as a multiplanar surface, with two lists of lons and lats

get_surface_boundaries_3d()

Returns the boundaries as three lists of lons, lats, depths

get_top_edge_depth()

Compute top edge depth of each surface element and return area-weighted average value (in km).

get_width()

Compute width of each surface element, and return area-weighted average value (in km).

property mesh

Returns:

mesh corresponding to the whole multi surface

property surface_nodes

Returns:

a list of surface nodes from the underlying single node surfaces

openquake.hazardlib.geo.surface.multi.build_msparams(rupture_idxs, secparams, close_sec=None, ry0=False, mon1=<Monitor [runner]>, mon2=<Monitor [runner]>)

Returns:

a structured array of parameters

openquake.hazardlib.geo.surface.multi.build_secparams(sections)

Returns:

an array of section parameters

openquake.hazardlib.geo.surface.planar module

Module openquake.hazardlib.geo.surface.planar contains PlanarSurface.

class openquake.hazardlib.geo.surface.planar.PlanarSurface(strike, dip, top_left, top_right, bottom_right, bottom_left, check=True)

Bases: BaseSurface

Planar rectangular surface with two sides parallel to the Earth surface.

Parameters:

• strike – Strike of the surface is the azimuth from top_left to top_right points.

• dip – Dip is the angle between the surface itself and the earth surface.

Other parameters are points (instances of Point) defining the surface corners in clockwise direction starting from top left corner. Top and bottom edges of the polygon must be parallel to earth surface and to each other.

See NodalPlane for more detailed definition of strike and dip. Note that these parameters are supposed to match the factual surface geometry (defined by corner points), but this is not enforced or even checked.

Raises:

ValueError – If either top or bottom points differ in depth or if top edge is not parallel to the bottom edge, if top edge differs in length from the bottom one, or if mesh spacing is not positive.

property bottom_left

property bottom_right

property corner_depths

property corner_lats

property corner_lons

classmethod from_(planar_array)

classmethod from_array(array34)

Parameters:

array34 – an array of shape (3, 4) in order tl, tr, bl, br

Returns:

a PlanarSurface instance

classmethod from_corner_points(top_left, top_right, bottom_right, bottom_left)

Create and return a planar surface from four corner points.

The azimuth of the line connecting the top left and the top right corners define the surface strike, while the angle between the line connecting the top left and bottom left corners and a line parallel to the earth surface defines the surface dip.

Parameters:

• top_left (openquake.hazardlib.geo.point.Point) – Upper left corner

• top_right (openquake.hazardlib.geo.point.Point) – Upper right corner

• bottom_right (openquake.hazardlib.geo.point.Point) – Lower right corner

• bottom_left (openquake.hazardlib.geo.point.Point) – Lower left corner

Returns:

An instance of PlanarSurface.

classmethod from_hypocenter(hypoc, msr, mag, aratio, strike, dip, rake, ztor=None)

Create and return a planar surface given the hypocenter location and other rupture properties.

Parameters:

• hypoc – An instance of :class: openquake.hazardlib.geo.point.Point

• msr – The magnitude scaling relationship e.g. an instance of :class: openquake.hazardlib.scalerel.WC1994

• mag – The magnitude

• aratio – The rupture aspect ratio

• strike – The rupture strike

• dip – The rupture dip

• rake – The rupture rake

• ztor – If not None it doesn’t consider the hypocentral depth constraint

classmethod from_ucerf(array43)

Parameters:

array43 – an array of shape (4, 3) in order tl, tr, br, bl

Returns:

a PlanarSurface instance

get_area()

Return surface’s area value (in squared km) obtained as the product of surface length and width.

get_bounding_box()

Compute surface bounding box from plane’s corners coordinates. Calls openquake.hazardlib.geo.utils.get_spherical_bounding_box()

Returns:

A tuple of four items. These items represent western, eastern, northern and southern borders of the bounding box respectively. Values are floats in decimal degrees.

get_closest_points(mesh)

See superclass' method.

get_dip()

Return dip value that was provided to the constructor.

get_joyner_boore_distance(mesh)

See superclass' method.

This is an optimized version specific to planar surface that doesn’t make use of the mesh.

get_middle_point()

Compute middle point from surface’s corners coordinates. Calls openquake.hazardlib.geo.utils.get_middle_point()

get_min_distance(mesh)

See superclass' method.

get_rx_distance(mesh)

See superclass method for spec of input and result values.

This is an optimized version specific to planar surface that doesn’t make use of the mesh.

get_ry0_distance(mesh)

Parameters:

mesh – Mesh of points to calculate Ry0-distance to.

Returns:

Numpy array of distances in km.

See also superclass method for spec of input and result values.

This is version specific to the planar surface doesn’t make use of the mesh

get_strike()

Return strike value that was provided to the constructor.

get_surface_boundaries()

The corners lons/lats in WKT-friendly order (clockwise)

get_surface_boundaries_3d()

The corners lons/lats/depths in WKT-friendly order (clockwise)

get_top_edge_centroid()

Overrides superclass' method in order to avoid creating a mesh.

get_top_edge_depth()

Overrides superclass' method in order to avoid creating a mesh.

get_width()

Return surface’s width value (in km) as computed in the constructor (that is mean value of left and right surface sides).

property length

Return length of the rupture

property mesh

Returns:

a mesh with the 4 corner points tl, tr, bl, br

property surface_nodes

A single element list containing a planarSurface node

property top_left

property top_right

property tor

Returns:

top of rupture line

translate(p1, p2)

Translate the surface for a specific distance along a specific azimuth direction.

Parameters are two points (instances of openquake.hazardlib.geo.point.Point) representing the direction and an azimuth for translation. The resulting surface corner points will be that far along that azimuth from respective corner points of this surface as p2 is located with respect to p1.

Returns:

A new PlanarSurface object with the same mesh spacing, dip, strike, width, length and depth but with corners longitudes and latitudes translated.

property width

Return length of the rupture

openquake.hazardlib.geo.surface.planar.build_corners(usd, lsd, rar, area, mag, strike, dip, rake, hdd, lon, lat)

openquake.hazardlib.geo.surface.planar.build_planar(planin, hdd, lon, lat, usd, lsd, rar, shift_hypo=False)

Parameters:

• planin – Surface input parameters as an array of shape (M, N)

• hdd – Hypocenter depths

• lat (lon,) – Longitude and latitude of the hypocenters (scalars)

Returns:

an array of shape (M, N, D, 3)

openquake.hazardlib.geo.surface.planar.build_planar_array(corners, sdr=None, hypo=None, check=False)

Parameters:

• corners – array of shape (4, M, N, D, 3)

• hypo – None or array of shape (M, N, D, 3)

Returns:

a planar_array array of length (M, N, D, 3)

openquake.hazardlib.geo.surface.planar.dot(a, b)

openquake.hazardlib.geo.surface.planar.get_azimuth(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of of shape (N, 3)

Returns:

(U, N) distances

openquake.hazardlib.geo.surface.planar.get_distances_planar(planar, sites, dist_type)

Parameters:

• planar – a planar array of shape (U, 3)

• sites – a filtered site collection with N sites

• dist_type – kind of distance to compute

Returns:

an array of distances of shape (U, N)

openquake.hazardlib.geo.surface.planar.get_repi(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of of shape (N, 3)

Returns:

(U, N) distances

openquake.hazardlib.geo.surface.planar.get_rhypo(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of of shape (N, 3)

Returns:

(U, N) distances

openquake.hazardlib.geo.surface.planar.get_rjb(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of of shape (N, 3)

Returns:

(U, N) values

openquake.hazardlib.geo.surface.planar.get_rupdims(usd, lsd, rar, area, dip)

Parameters:

• usd – upper seismogenic depth

• lsd – lower seismogenic depth

• rar – rupture aspect ratio

• area – area of the surface

• dip – dip angle

Returns:

array of shape 3with rupture length, width and height

The rupture area is calculated using the method get_median_area(). If the calculated rupture width, inclined by the nodal plane’s dip angle, would not fit in between upper and lower seismogenic depth, the rupture width is shrunken to the maximum possible and the rupture length is extended to preserve the same area.

openquake.hazardlib.geo.surface.planar.get_rvolc(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of of shape (N, 3)

Returns:

(U, N) distances

openquake.hazardlib.geo.surface.planar.get_rx(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of of shape (N, 3)

Returns:

(U, N) distances

openquake.hazardlib.geo.surface.planar.get_ry0(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of of shape (N, 3)

Returns:

(U, N) distances

openquake.hazardlib.geo.surface.planar.project(planar, points)

Parameters:

• planar – a planar recarray of shape (U, 3)

• points – an array of euclidean coordinates of shape (N, 3)

Returns:

(3, U, N) values

openquake.hazardlib.geo.surface.planar.project_back(planar, xx, yy)

Parameters:

• planar – a planar recarray of shape (U, 3)

• xx – an array of of shape (U, N)

• yy – an array of of shape (U, N)

Returns:

(3, U, N) values

openquake.hazardlib.geo.surface.simple_fault module

Module openquake.hazardlib.geo.surface.simple_fault defines SimpleFaultSurface.

class openquake.hazardlib.geo.surface.simple_fault.SimpleFaultSurface(mesh)

Bases: BaseSurface

Represent a fault surface as regular (uniformly spaced) 3D mesh of points.

Parameters:

mesh – Instance of RectangularMesh representing surface geometry.

Another way to construct the surface object is to call from_fault_data().

classmethod check_fault_data(fault_trace, upper_seismogenic_depth, lower_seismogenic_depth, dip, mesh_spacing)

Verify the fault data and raise ValueError if anything is wrong.

This method doesn’t have to be called by hands before creating the surface object, because it is called from from_fault_data().

classmethod from_fault_data(fault_trace, upper_seismogenic_depth, lower_seismogenic_depth, dip, mesh_spacing)

Create and return a fault surface using fault source data.

Parameters:

• fault_trace (openquake.hazardlib.geo.line.Line) – Geographical line representing the intersection between the fault surface and the earth surface. The line must be horizontal (i.e. all depth values must be equal). If the depths are not given, they are assumed to be zero, meaning the trace intersects the surface at sea level, e.g. fault_trace = Line([Point(1, 1), Point(1, 2)]).

• upper_seismo_depth – Minimum depth ruptures can reach, in km (i.e. depth to fault’s top edge).

• lower_seismo_depth – Maximum depth ruptures can reach, in km (i.e. depth to fault’s bottom edge).

• dip – Dip angle (i.e. angle between fault surface and earth surface), in degrees.

• mesh_spacing – Distance between two subsequent points in a mesh, in km.

Returns:

An instance of SimpleFaultSurface created using that data.

Uses check_fault_data() for checking parameters.

get_dip()

Return the fault dip as the average dip over the fault surface mesh.

The average dip is defined as the weighted mean inclination of top row of mesh cells. See openquake.hazardlib.geo.mesh.RectangularMesh.get_mean_inclination_and_azimuth()

Returns:

The average dip, in decimal degrees.

classmethod get_fault_patch_vertices(rupture_top_edge, upper_seismogenic_depth, lower_seismogenic_depth, dip, index_patch=1)

Get surface main vertices. Parameters are the same as for from_fault_data(), excluding fault_trace, and mesh spacing.

Parameters:

• rupture_top_edge – A instances of openquake.hazardlib.geo.line.Line representing the rupture surface’s top edge.

• index_patch – Indicate the patch of the fault in order to output the vertices. The fault patch numbering follows the same logic of the right-hand rule i.e. patch with index 1 is the first patch along the trace.

Returns:

Four :class:~openquake.hazardlib.geo.point.Point objects representing the four vertices of the target patch.

get_strike()

Return the fault strike as the average strike along the fault trace.

The average strike is defined as the weighted mean azimuth of top row of mesh cells. See openquake.hazardlib.geo.mesh.RectangularMesh.get_mean_inclination_and_azimuth()

Returns:

The average strike, in decimal degrees.

classmethod get_surface_vertexes(fault_trace, upper_seismogenic_depth, lower_seismogenic_depth, dip)

Get surface main vertexes.

Parameters are the same as for from_fault_data(), excluding mesh spacing.

Returns:

Instance of Polygon describing the surface projection of the simple fault with specified parameters.

get_width()

Return surface’s width (that is surface extension along the dip direction) in km.

The width is computed as the average width along the surface. See openquake.hazardlib.geo.mesh.RectangularMesh.get_mean_width()

classmethod hypocentre_patch_index(hypocentre, rupture_top_edge, upper_seismogenic_depth, lower_seismogenic_depth, dip)

This methods finds the index of the fault patch including the hypocentre.

Parameters:

• hypocentre – Point object representing the location of hypocentre.

• rupture_top_edge – A instances of openquake.hazardlib.geo.line.Line representing the rupture surface’s top edge.

• upper_seismo_depth – Minimum depth ruptures can reach, in km (i.e. depth to fault’s top edge).

• lower_seismo_depth – Maximum depth ruptures can reach, in km (i.e. depth to fault’s bottom edge).

• dip – Dip angle (i.e. angle between fault surface and earth surface), in degrees.

Returns:

An integer corresponding to the index of the fault patch which contains the hypocentre.

classmethod surface_projection_from_fault_data(fault_trace, upper_seismogenic_depth, lower_seismogenic_depth, dip)

Get a surface projection of the simple fault surface.

Parameters are the same as for from_fault_data(), excluding mesh spacing.

Returns:

Instance of Polygon describing the surface projection of the simple fault with specified parameters.

property tor

Returns:

top of rupture line

openquake.hazardlib.geo.surface.simple_fault.simple_fault_node(fault_trace, dip, upper_depth, lower_depth)

Parameters:

• fault_trace – an object with an attribute .points

• dip – dip parameter

• upper_depth – upper seismogenic depth

• lower_depth – lower seismogenic depth

Returns:

a Node of kind simpleFaultGeometry

Module contents

Package openquake.hazardlib.geo.surface contains base and specific implementations of earthquake rupture surfaces.