diff --git a/tudatpy/kernel/expose_numerical_simulation/expose_environment_setup/expose_gravity_field_variation_setup.cpp b/tudatpy/kernel/expose_numerical_simulation/expose_environment_setup/expose_gravity_field_variation_setup.cpp
index 00587f27..9b0b23c0 100644
--- a/tudatpy/kernel/expose_numerical_simulation/expose_environment_setup/expose_gravity_field_variation_setup.cpp
+++ b/tudatpy/kernel/expose_numerical_simulation/expose_environment_setup/expose_gravity_field_variation_setup.cpp
@@ -192,6 +192,12 @@ void expose_gravity_field_variation_setup(py::module &m) {
           py::arg("minimum_order") = 0,
           get_docstring("single_power_polynomial").c_str() );
 
+    m.def("mode_coupled_solid_body_tide",
+          &tss::modeCoupledSolidBodyGravityFieldVariationSettings,
+          py::arg("deforming_bodies"),
+          py::arg("love_numbers"),
+          get_docstring("mode_coupled_solid_body_tide").c_str() );
+
     m.def("polynomial",
           &tss::polynomialGravityFieldVariationsSettings,
           py::arg("cosine_amplitudes_per_power"),
diff --git a/tudatpy/kernel/expose_numerical_simulation/expose_estimation_setup/expose_estimated_parameter_setup.cpp b/tudatpy/kernel/expose_numerical_simulation/expose_estimation_setup/expose_estimated_parameter_setup.cpp
index bb5449ea..25856398 100644
--- a/tudatpy/kernel/expose_numerical_simulation/expose_estimation_setup/expose_estimated_parameter_setup.cpp
+++ b/tudatpy/kernel/expose_numerical_simulation/expose_estimation_setup/expose_estimated_parameter_setup.cpp
@@ -434,6 +434,13 @@ void expose_estimated_parameter_setup(py::module &m) {
           py::arg("use_complex_love_number") = 0,
           get_docstring("order_varying_k_love_number", 2).c_str() );
 
+    m.def("mode_coupled_k_love_numbers",
+          &tep::modeCoupledTidalLoveNumberEstimatableParameterSettings,
+          py::arg("deformed_body"),
+          py::arg("love_number_indices"),
+          py::arg("deforming_bodies"),
+          get_docstring("mode_coupled_k_love_numbers").c_str() );
+
     m.def("polynomial_gravity_field_variation_amplitudes",
           &tep::polynomialGravityFieldVariationParameter,
           py::arg("body_name"),