Seismic¶
Base class for all seismic models¶
- class burnman.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).
Class for 1D Models¶
- class burnman.seismic.SeismicTable[source]¶
Bases:
burnman.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].
- 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)
- 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.seismic.PREM[source]¶
Bases:
burnman.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.seismic.Slow[source]¶
Bases:
burnman.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.seismic.Fast[source]¶
Bases:
burnman.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.seismic.STW105[source]¶
Bases:
burnman.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.seismic.IASP91[source]¶
Bases:
burnman.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.seismic.AK135[source]¶
Bases:
burnman.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.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].