diff --git a/biogeochem/FatesAllometryMod.F90 b/biogeochem/FatesAllometryMod.F90
index 12934256..8ced4209 100644
--- a/biogeochem/FatesAllometryMod.F90
+++ b/biogeochem/FatesAllometryMod.F90
@@ -2155,7 +2155,8 @@ function CumulativeLayerTVAI(icanlayer,     &
        call endrun(msg=errMsg(sourcefile, __LINE__))
     end if
 
-    cum_tvai = cum_tvai + tvai0 + tvai
+    ! Apply the cumulative TVAI at the current layer mid-point (ie 0.5*)
+    cum_tvai = cum_tvai + tvai0 + 0.5*tvai
     
     return
   end function CumulativeLayerTVAI