LPMLEn - A code for estimating heat transport parameters in 1D

LPMLEn - A code for estimating heat transport parameters in 1D

The LPMLEn combines the local polynomial method (LP method) with a maximum likelihood estimator (MLE) to estimate 1D vertical streambed fluxes and thermal diffusivities using time-series from n temperature sensors. It operates in the frequency domain and can use multiple frequencies and sensors simultaneously for the parameter estimation. The LPMLEn is provided here with two models, (i) the semi-infinite domain model where only an upper temperature boundary condition is used to estimate the parameters and (ii) a bounded (finite) domain model where an additional lower local temperature boundary condition is assigned to estimate the parameters for a distinct section of the streambed.

Contents
The MATLAB scirpts that are used to create the figures in the paper are:
- Estimation_with_synthetic_dataset1_and_2_SI_vs_BD.m for Table 1, Fig. S1 and S2.
- Estimation_with_synthetic_dataset3_change_in_D.m for Fig. 1.
- Estimation_with_synthetic_dataset4_change_in_D_from_low_to_high.m for Fig. S3.
- Estimation_with_experimental_dataset.m for Fig. 2b, 2c, 3, 4, 5, S4 and S5.
- ML1_90.txt contains the measurement data of the experimental dataset.
The analysis performed on the dataset in Estimation_with_experimental_dataset.m is resource demanding. For this reason the computational results are saved in Estimation_experimental_dataset_workspace.mat, which can be loaded into MATLAB to bypass the computations.

To start using the LPMLEn, please check the simplified example Example_simplified_LPMLEn.m that uses the function MLEn_hydrology_time.m that only requires the time-series, measurement depths, and model choice as input.
The more advance user may want to use the LP-method (LocalPolyAnal.m) and MLEn (MLEn.m) sepperatly for more control and advanced settings. For this, the Estimation_with_experimental_dataset.m can be used as an example.