Hvasbath Beat Versions Save

Bugfix release! Many fixes for GNSS related functionality.

Contributors: Hannes Vasyura-Bathke @hvasbath, Semih Ergintav @sergintav

Contributors: Hannes Vasyura-Bathke @hvasbath

New major version is not backward compatible! Also requires python3.9+!

Added

• new inference mode: "bem" for Boundary Element Modeling
• allow for multi-source type inference: e.g. MTSource, RectangularSource
• added parameter defaults module and config functionality

Changed

• using pymc v5 and pytensor instead of pymc v3 and theano --> makes old setups incompatible
• n_sources in config is now list of integers, previous: integer
• source_type to: list of source_types, variable changed to source_types
• adopted ruff linting
• replace nose with pytest testing

Contributors: Hannes Vasyura-Bathke @hvasbath

Added

• heart.ResultPoint: extend by event object from MAP
• plotting.seismic,geodetic: add standardized residuals to fit plots
• plotting.marginals: add transforms to MTQT parameters plotting marginals
• plotting.marginals: multi-source support for correlation_hist
• plotting.seismic: "station_variance_reductions" plot
• added new source_type: "SFSource"

Changed

• plotting.seismic.seismic_fits: subplots are ordered row wise distance based, channels in cols

Fixed

• sampler.base.iter_parallel_chains: fix chunksize determination
• minor logging messages

Contributors: Hannes Vasyura-Bathke @hvasbath

Added

• covariance.GeodeticNoiseAnalyser: parametrize the residual noise allowing for non-Toeplitz/import
• plotting.geodetic/seismic added standardized residual histograms to residuals
• plotting.geodetic: new plot "geodetic_covariances"
• plotting.geodetic: add "individual" plot_projection to scene_fits to show stdz residuals
• plotting.seismic: fuzzy_bb, lune, hudson and fuzzy_mt_decomp support n_sources > 1
• plotting.seismic: allow plotting of fuzzy_bb for RectangularSource

Fixed

• plotting.marginals: stage_posteriors fixed axis unification and erroneous histogram plotting
• docs: short_installation fix python version to 3.8
• heart: pol_synthetics allow for RectangularSource
• covariance: estimation of variance on amplitude spectra instead of complex spectra

BUGFIX release

Contributors: Hannes Vasyura-Bathke @hvasbath

Fixed

• FFI: do not init wavemaps during SeismicComposite init
• heart: flexible versioning for geodetic GF backend
• docs: correctly state python3.8 in installation instructions

Contributors: Hannes Vasyura-Bathke @hvasbath

Added

• plotting.marginals.traceplot: source_idxs argument can be slice e.g. 1:10 to take mean of patch attributes
• heart.seis_derivative: function to numerically calculate derivatives wrt source attributes for waveforms
• heart.py, config.py: allow for version control for fomosto backends through "version" argument at gf_config
• utility: functions for slice to string conversion and vice versa
• config; NonlinearGFConfig added version attribute to specify backend version to use for GF calculation

Changed

• plotting.marginals.traceplot: CDF plotting of multiple distributions in a single subplot marks quantiles only once

Fixed

• plotting.marginals.traceplot: multipage output repeating variables

Contributors: Hannes Vasyura-Bathke @hvasbath

Added

FFI

---

• add calculation of coupling to derived variables at Fault for the summarize method

plotting

---

• plot_projection argument available for stage_posteriors plot: cdf, pdf, kde
• update Example 8 to showcase the use of this functionality

Fixed

• stage_posterior plot supports multipage output, figsize fixed to fractions of A4
• multievent waveform_fits plot returns separate figures for each sub-event

Contributors: Mahdi Hamidbeygi @mahdihamidbeygi, Hannes Vasyura-Bathke @hvasbath

Added

Polarity

---

• polarity module/composite/dataset-type for inference of wave onset polarities
• joint inference of P, Sh and Sv wave-polarities
• drawing of piercing points on fuzzyBB plot if polarity data exists
• drawing of fuzzyBBs for Sh and Sv waves
• add tutorial Example 8 for polarity inference

Amplitude Spectra

---

• added string-choice domain to WaveMaps to choose time or frequency domain for inference
• plotting: added spectra fits in waveform fits plot

Continuous Integration (CI)

---

• adapted installation scheme to fulfill PEP517 through pip and .toml
• providing beat binary packages on PyPi
• github actions workflow for code formatting and wheel building
• pre-commit hooks for spellchecking, black, formatting

Changed

• split plotting module into plotting directory with submodules
• plotting: stage_posteriors change hist color to grey if only single density is shown
• plotting: fuzzyBB dashed white and black line instead of red for MAP
• docs: updated (short) installation instructions to use the package manager pip instead of setup.py

Added:

FFI

---

• chain initialisation with lsq solution for seismic data
• resolution based discretization:
  • add improved method of Atzori et al. 2019
  • automatic damping estimation
  • add Tutorial 4b to demonstrate the usage and functionality

Plotting

---

• 3d_slip_distribution: allows plotting of selected segments, slip-deficit, coupling, slip_variation
• gnss_fits: added Variance Reduction histograms if nensemble > 1
• slip_distribution: common colorscale and spatial scale accross subfaults
• moment_rate: adjusted size and formatting to be publication ready
• correlation_hist: plot source related correlations for nsources > 1, colorcoding

Changed:

• docs: rename former example 4 to example 4a

Fixed:

• Bin backend: loading of traces with corrupt file header
• beat import: just import sampled variables
• RingFaultSource: allow sampling of slip-direction (sign)

Added:

General

---

• added export method to geodetic data composite
• added functionality for multi-event estimation that are seperated longer time and waveforms are not stacked
• Wavemaps have quantity argument that allow for sampling in displacement, velocity or acceleration space
• added get_variance_reductions method to seismic composite

Tensile Dislocations

---

• added the opening_fraction argument to the RectangularSource to allow to model tensile opening/closing
• for geometry mode, as well as ffi mode to estimate distributed opening / closing

Finite Fault

---

• discretization options for the fault into patches:
  • uniform
  • resolution based (after Atzori & Antonioli 2011 GJI) for geodetic data valid only
• support for multiple subfaults in FFI mode:
  • Laplacian smoothing with different functional forms:
    • distance based: ** Gaussian ** Exponential (after Radiget et al. 2011)
    • only neighbors: ** nearest_neighbor for single subfault only
  • geodetic data
  • seismic data with nucleation point for each subfault
• allow for station corrections
• easier fault geometry setup:
  • beat check --what=geometry allows to pipe fault setup into Talpa for GUI interactive fault editing
  • beat import --results=saved_geometry.yaml --mode=ffi imports the saved_geometry.yaml to the beat config_ffi

Hierarchicals

---

• generalized hierarchical corrections to support various types:
  • geodetic:
    • ramps
    • Euler pole rotation

Plotting

---

• slip_distribution allows for variable patch-size
• plot gnss_fits for horizontal and vertical GNSS components
• waveform misfits:
  • with plot_projection =individual allows to get individual source contributions
  • added variance reduction histograms on top right
  • added time-shift histograms on bottom left (if they were sampled)
• new station map with GMT, allows time shift plotting
• FFI exports geometry object for fast and easy interactive inspection in sparrow
• Lune plot for moment tensors (Tape and Tape 2015)

Changed:

• filter attribute in wavemaps is now a list of filter objects which allows to chain several filters

Fixed:

• extended support for GNSS data from globk
• SMC saving of sampler state does not dump full model graph
• multiprocess FFI performance increase
• MTQTSource magnitude scaling