diff --git a/doc/dynare.texi b/doc/dynare.texi
index 4a6c642ac1f7fcf7696c6848a9bc30b49a385f6f..d661d35b6e42dccd2c5f3c717d26ae6446410f13 100644
--- a/doc/dynare.texi
+++ b/doc/dynare.texi
@@ -4905,7 +4905,11 @@ simulation that crashed prematurely. Shouldn't be used together with
 Name of the file containing previous value for the mode. When
 computing the mode, Dynare stores the mode (@code{xparam1}) and the
 hessian (@code{hh}, only if @code{cova_compute=1}) in a file called
-@file{@var{MODEL_FILENAME}_mode.mat}
+@file{@var{MODEL_FILENAME}_mode.mat}. After a successful run of the estimation
+command, the @code{mode_file} will be disabled to prevent other function calls
+from implicitly using an updated mode-file. Thus, if the mod-file contains subsequent
+@code{estimation} commands, the @code{mode_file} option, if desired, needs to be 
+specified again.
 
 @item mode_compute = @var{INTEGER} | @var{FUNCTION_NAME}
 @anchor{mode_compute}
diff --git a/matlab/dynare_estimation.m b/matlab/dynare_estimation.m
index 7bb49f38038222a859a2cecf17015140a96e4d82..06a220e459644a98ca1a534a8f276913c27533a0 100644
--- a/matlab/dynare_estimation.m
+++ b/matlab/dynare_estimation.m
@@ -199,3 +199,4 @@ if nnobs > 1 && horizon > 0
         end
     end
 end
+options_.mode_file = ''; %delete stored mode-file so that it is not reaccessed in later calls (and in case it was only set by the recursive estimation)