Seismic¶

Base class for all seismic models¶

class burnman.classes.seismic.Seismic1DModel[source]¶

Bases: object

Base class for all the seismological models.

evaluate(vars_list, depth_list=None)[source]¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

internal_depth_list(mindepth=0.0, maxdepth=1e+99, discontinuity_interval=1.0)[source]¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

v_p(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_s(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

v_phi(depth)[source]¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

density(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

G(depth)[source]¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)[source]¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

QK(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

QG(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

depth(pressure)[source]¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

gravity(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

Class for 1D Models¶

class burnman.classes.seismic.SeismicTable[source]¶

Bases: burnman.classes.seismic.Seismic1DModel

This is a base class that gets a 1D seismic model from a table indexed and sorted by radius. Fill the tables in the constructor after deriving from this class. This class uses burnman.seismic.Seismic1DModel

Note: all tables need to be sorted by increasing depth. self.table_depth needs to be defined Alternatively, you can also overwrite the _lookup function if you want to access with something else.

internal_depth_list(mindepth=0.0, maxdepth=10000000000.0, discontinuity_interval=1.0)[source]¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

gravity(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

v_p(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_s(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

QK(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

QG(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

density(depth)[source]¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

bullen(depth)[source]¶: Returns the Bullen parameter only for significant arrays

depth(pressure)[source]¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

radius(pressure)[source]¶

G(depth)¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

evaluate(vars_list, depth_list=None)¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

v_phi(depth)¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

Models currently implemented¶

class burnman.classes.seismic.PREM[source]¶

Bases: burnman.classes.seismic.SeismicTable

Reads PREM (1s) (input_seismic/prem.txt, [DA81]). See also burnman.seismic.SeismicTable.

G(depth)¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

QG(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

QK(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

bullen(depth)¶: Returns the Bullen parameter only for significant arrays

density(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

depth(pressure)¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

evaluate(vars_list, depth_list=None)¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

gravity(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

internal_depth_list(mindepth=0.0, maxdepth=10000000000.0, discontinuity_interval=1.0)¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

radius(pressure)¶

v_p(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_phi(depth)¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

v_s(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

class burnman.classes.seismic.Slow[source]¶

Bases: burnman.classes.seismic.SeismicTable

Inserts the mean profiles for slower regions in the lower mantle (Lekic et al. 2012). We stitch together tables ‘input_seismic/prem_lowermantle.txt’, ‘input_seismic/swave_slow.txt’, ‘input_seismic/pwave_slow.txt’). See also burnman.seismic.SeismicTable.

G(depth)¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

QG(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

QK(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

bullen(depth)¶: Returns the Bullen parameter only for significant arrays

density(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

depth(pressure)¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

evaluate(vars_list, depth_list=None)¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

gravity(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

internal_depth_list(mindepth=0.0, maxdepth=10000000000.0, discontinuity_interval=1.0)¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

radius(pressure)¶

v_p(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_phi(depth)¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

v_s(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

class burnman.classes.seismic.Fast[source]¶

Bases: burnman.classes.seismic.SeismicTable

Inserts the mean profiles for faster regions in the lower mantle (Lekic et al. 2012). We stitch together tables ‘input_seismic/prem_lowermantle.txt’, ‘input_seismic/swave_fast.txt’, ‘input_seismic/pwave_fast.txt’). See also burnman.seismic.Seismic1DModel.

G(depth)¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

QG(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

QK(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

bullen(depth)¶: Returns the Bullen parameter only for significant arrays

density(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

depth(pressure)¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

evaluate(vars_list, depth_list=None)¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

gravity(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

internal_depth_list(mindepth=0.0, maxdepth=10000000000.0, discontinuity_interval=1.0)¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

radius(pressure)¶

v_p(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_phi(depth)¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

v_s(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

class burnman.classes.seismic.STW105[source]¶

Bases: burnman.classes.seismic.SeismicTable

Reads STW05 (a.k.a. REF) (1s) (input_seismic/STW105.txt, [KED08]). See also burnman.seismic.SeismicTable.

G(depth)¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

QG(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

QK(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

bullen(depth)¶: Returns the Bullen parameter only for significant arrays

density(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

depth(pressure)¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

evaluate(vars_list, depth_list=None)¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

gravity(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

internal_depth_list(mindepth=0.0, maxdepth=10000000000.0, discontinuity_interval=1.0)¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

radius(pressure)¶

v_p(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_phi(depth)¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

v_s(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

class burnman.classes.seismic.IASP91[source]¶

Bases: burnman.classes.seismic.SeismicTable

Reads REF/STW05 (input_seismic/STW105.txt, [KED08]). See also burnman.seismic.SeismicTable.

G(depth)¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

QG(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

QK(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

bullen(depth)¶: Returns the Bullen parameter only for significant arrays

density(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

depth(pressure)¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

evaluate(vars_list, depth_list=None)¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

gravity(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

internal_depth_list(mindepth=0.0, maxdepth=10000000000.0, discontinuity_interval=1.0)¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

radius(pressure)¶

v_p(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_phi(depth)¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

v_s(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

class burnman.classes.seismic.AK135[source]¶

Bases: burnman.classes.seismic.SeismicTable

Reads AK135 (input_seismic/ak135.txt, [KEB95]). See also burnman.seismic.SeismicTable.

G(depth)¶

Parameters

depthfloat or array of floats: Shear modulus at given for depth(s) in [Pa].

K(depth)¶

Parameters

depthfloat or array of floats: Bulk modulus at given for depth(s) in [Pa]

QG(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

QGfloat or array of floats: Quality factor (dimensionless) for shear modulus at given depth(s).

QK(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

Qkfloat or array of floats: Quality factor (dimensionless) for bulk modulus at given depth(s).

bullen(depth)¶: Returns the Bullen parameter only for significant arrays

density(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

densityfloat or array of floats: Density at given depth(s) in [kg/m^3].

depth(pressure)¶

Parameters

pressurefloat or array of floats: Pressure(s) [Pa] to evaluate depth at.

Returns

depthfloat or array of floats: Depth(s) [m] for given pressure(s)

evaluate(vars_list, depth_list=None)¶

Returns the lists of data for a Seismic1DModel for the depths provided

Parameters

vars_listarray of str: Available variables depend on the seismic model, and can be chosen from ‘pressure’,’density’,’gravity’,’v_s’,’v_p’,’v_phi’,’G’,’K’,’QG’,’QK’
depth_listarray of floats: Array of depths [m] to evaluate seismic model at.

Returns

Array of values shapes as (len(vars_list),len(depth_list)).

gravity(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate gravity at.

Returns

gravityfloat or array of floats: Gravity for given depths in [m/s^2]

internal_depth_list(mindepth=0.0, maxdepth=10000000000.0, discontinuity_interval=1.0)¶

Returns a sorted list of depths where this seismic data is specified at. This allows you to compare the seismic data without interpolation. The depths can be bounded by the mindepth and maxdepth parameters.

Parameters

mindepthfloat: Minimum depth value to be returned [m]
maxdepth: Maximum depth value to be returned [m]
discontinuity interval: Shift continuities to remove ambigious values for depth, default value = 1 [m]

Returns

depthsarray of floats: Depths [m].

pressure(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

pressurefloat or array of floats: Pressure(s) at given depth(s) in [Pa].

radius(pressure)¶

v_p(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_pfloat or array of floats: P wave velocity at given depth(s) in [m/s].

v_phi(depth)¶

Parameters

depth_listfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_phifloat or array of floats: bulk sound wave velocity at given depth(s) in [m/s].

v_s(depth)¶

Parameters

depthfloat or array of floats: Depth(s) [m] to evaluate seismic model at.

Returns

v_sfloat or array of floats: S wave velocity at given depth(s) in [m/s].

Attenuation Correction¶

burnman.classes.seismic.attenuation_correction(v_p, v_s, v_phi, Qs, Qphi)[source]¶

Applies the attenuation correction following Matas et al. (2007), page 4. This is simplified, and there is also currently no 1D Q model implemented. The correction, however, only slightly reduces the velocities, and can be ignored for our current applications. Arguably, it might not be as relevant when comparing computations to PREM for periods of 1s as is implemented here. Called from burnman.main.apply_attenuation_correction()

Parameters

v_pfloat: P wave velocity in [m/s].
v_sfloat: S wave velocitiy in [m/s].
v_phifloat: Bulk sound velocity in [m/s].
Qsfloat: shear quality factor [dimensionless]
Qphi: float: bulk quality factor [dimensionless]

Returns

v_pfloat: corrected P wave velocity in [m/s].
v_sfloat: corrected S wave velocitiy in [m/s].
v_phifloat: corrected Bulk sound velocity in [m/s].