diff --git a/.gitmodules b/.gitmodules
index 75876f6fa..e1ad81bb1 100644
--- a/.gitmodules
+++ b/.gitmodules
@@ -8,8 +8,10 @@
   branch = production/RRFS.v1
 [submodule "ccpp/physics"]
   path = ccpp/physics
-  url = https://github.com/ufs-community/ccpp-physics
-  branch = production/RRFS.v1
+#  url = https://github.com/ufs-community/ccpp-physics
+#  branch = production/RRFS.v1
+  url = https://github.com/JiliDong-NOAA/ccpp-physics    
+  branch = sigmab_restart_fix
 [submodule "upp"]
   path = upp
   url = https://github.com/NOAA-EMC/UPP
diff --git a/ccpp/driver/GFS_restart.F90 b/ccpp/driver/GFS_restart.F90
index 890aeb9fb..f00784679 100644
--- a/ccpp/driver/GFS_restart.F90
+++ b/ccpp/driver/GFS_restart.F90
@@ -179,7 +179,7 @@ subroutine GFS_restart_populate (Restart, Model, Statein, Stateout, Sfcprop, &
       Restart%num3d = Restart%num3d + 9
     endif
     if (Model%rrfs_sd) then
-      Restart%num3d = Restart%num3d + 4
+      Restart%num3d = Restart%num3d + 5
     endif
     !Prognostic area fraction
     if (Model%progsigma) then
@@ -690,6 +690,11 @@ subroutine GFS_restart_populate (Restart, Model, Statein, Stateout, Sfcprop, &
         Restart%data(nb,num)%var3p => Coupling(nb)%chem3d(:,:,3)
       enddo
       num = num + 1
+      Restart%name3d(num) = 'ebu_smoke'
+      do nb = 1,nblks
+        Restart%data(nb,num)%var3p => Coupling(nb)%ebu_smoke(:,:)
+      enddo
+      num = num + 1
       Restart%name3d(num) = 'ext550'
       do nb = 1,nblks
         Restart%data(nb,num)%var3p => Radtend(nb)%ext550(:,:)
diff --git a/ccpp/physics b/ccpp/physics
index c19759964..852f8060d 160000
--- a/ccpp/physics
+++ b/ccpp/physics
@@ -1 +1 @@
-Subproject commit c197599640f63a35867c0cee993fee66a67e1b68
+Subproject commit 852f8060d20642272d707aa82a6222f5a04007fa
diff --git a/io/fv3atm_restart_io.F90 b/io/fv3atm_restart_io.F90
index 7dc8460ba..55bfb6c9f 100644
--- a/io/fv3atm_restart_io.F90
+++ b/io/fv3atm_restart_io.F90
@@ -405,6 +405,8 @@ subroutine fv3atm_checksum (Model, GFS_Data, Atm_block)
         call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%tsnow_ice)
         call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%snowfallac_land)
         call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%snowfallac_ice)
+        call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%acsnow_land)
+        call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%acsnow_ice)
         call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%sfalb_lnd)
         call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%sfalb_lnd_bck)
         call copy_from_GFS_Data(ii1,jj1,isc,jsc,nt,temp2d,GFS_Data(nb)%Sfcprop%sfalb_ice)
diff --git a/io/fv3atm_sfc_io.F90 b/io/fv3atm_sfc_io.F90
index c0bfcf6d9..bfd5d8c0f 100644
--- a/io/fv3atm_sfc_io.F90
+++ b/io/fv3atm_sfc_io.F90
@@ -123,9 +123,9 @@ function Sfc_io_calculate_indices(sfc, Model, reading, warm_start)
     endif
     if (Model%lsm == Model%lsm_ruc .and. warm_start) then
       if (Model%rdlai) then
-        nvar2r = 13
+        nvar2r = 15
       else
-        nvar2r = 12
+        nvar2r = 14
       endif
       nvar3  = 5
     else
@@ -534,6 +534,8 @@ subroutine Sfc_io_fill_2d_names(sfc,Model,warm_start)
       nt=nt+1 ; sfc%name2(nt) = 'tsnow_ice'
       nt=nt+1 ; sfc%name2(nt) = 'snowfall_acc_land'
       nt=nt+1 ; sfc%name2(nt) = 'snowfall_acc_ice'
+      nt=nt+1 ; sfc%name2(nt) = 'acsnow_land'
+      nt=nt+1 ; sfc%name2(nt) = 'acsnow_ice'
       nt=nt+1 ; sfc%name2(nt) = 'sfalb_lnd'
       nt=nt+1 ; sfc%name2(nt) = 'sfalb_lnd_bck'
       nt=nt+1 ; sfc%name2(nt) = 'sfalb_ice'
@@ -1051,6 +1053,8 @@ subroutine Sfc_io_transfer(sfc, reading, Model, Atm_block, Sfcprop, warm_start,
         call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%tsnow_ice)
         call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%snowfallac_land)
         call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%snowfallac_ice)
+        call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%acsnow_land)
+        call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%acsnow_ice)
         call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%sfalb_lnd)
         call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%sfalb_lnd_bck)
         call GFS_Data_transfer(reading,ii1,jj1,isc,jsc,nt,sfc%var2,Sfcprop(nb)%sfalb_ice)