Skip to content

MicroInverse is a Python package for inversion of a transport operator from tracer data

License

Notifications You must be signed in to change notification settings

hainegroup/MicroInverse

Repository files navigation

MicroInverse

Documentation Status

MicroInverse is a Python package for inversion of a transport operator from tracer data.

It is based on the simple stochastic climate model approximation

http://latex.codecogs.com/gif.latex?%5Cfrac%7Bd%7D%7Bdt%7D%20%5Cmathbf%7Bx%7D%28t%29%20%3D%20%5Cmathbf%7BB%7D%5Cmathbf%7Bx%7D%28t%29%20+%20%5Cmathbf%7Bf%7D%28t%29

Where x is the vector of tracer anomaly timeseries, B is the transport operator, and f is the forcing of the system. Assuming that the forcing has a shorter decorrelation timescale than the tracer we can solve for the transport operator:

http://latex.codecogs.com/gif.latex?%5Cmathbf%7BB%7D%3D%5Cfrac%7B1%7D%7B%5Ctau%7D%5Clog%20%5Cleft%28%5Cleft%5B%20%5Cmathbf%7Bx%7D%28t+%5Ctau%29%5Cmathbf%7Bx%7D%5ET%28t%29%5Cright%20%5D%20%5C%20%5Cleft%5B%5Cmathbf%7Bx%7D%28t%29%5Cmathbf%7Bx%7D%5ET%28t%29%20%5Cright%5D%5E%7B-1%7D%5Cright%29

Where tau is the chosen decorrelation timescale which should be larger than the forcing decorrelation timescale, but smaller than the decorrelation timescale of the tracer.

In practice tau is hard to choose a priori which is why we suggest first inverting your data at multiple values of tau and combining the results afterwards using MicroInverse.MicroInverse_utils.combine_Taus().

MicroInverse will also relate B to velocity, diffusivity, and decay via advection-diffusion-relaxation equation (see Nummelin et al. (2018) for details)

Features

  • Estimates of velocity, diffusivity, and decay timescale from a timeseries of 2D tracer.

Credits

This package is based on work by Nummelin et al. (2018) and Jeffress and Haine (2014a, 2014b)

Package is created with Cookiecutter using the audreyr/cookiecutter-pypackage project template.

About

MicroInverse is a Python package for inversion of a transport operator from tracer data

Resources

License

Stars

Watchers

Forks

Packages

No packages published