Skip to content

Latest commit

 

History

History
60 lines (36 loc) · 2.6 KB

README.rst

File metadata and controls

60 lines (36 loc) · 2.6 KB

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.