diff --git a/README.md b/README.md
index 89f748c..e1d357c 100644
--- a/README.md
+++ b/README.md
@@ -45,6 +45,11 @@ Copyright © 2022, Battelle Memorial Institute
 
 There are a number of publications that describe the algorithms used in `HexWatershed` in detail. If you make use of `HexWatershed` in your work, please consider including a reference to the following:
 
+* Liao, Chang, Cooper, M. G. PyFlowline a mesh independent river network generator for hydrologic models. Journal of Open Source Software. https://doi.org/10.21105/joss.05446
+
+Liao, Chang, Zhou, T., Xu, D., Tan, Z., Bisht, G., Cooper, M. G., et al. (2023). Topological relationship-based flow direction modeling: Stream burning and depression filling. Journal of Advances in Modeling Earth Systems, 15, e2022MS003487. https://doi.org/10.1029/2022MS003487
+
+* Liao, Chang, Zhou, T., Xu, D., Cooper, M. G., Engwirda, D., Li, H.-Y., Leung, L. R. (2023). Topological relationship-based flow direction modeling: Mesh-independent river networks representation. Journal of Advances in Modeling Earth Systems, 15, e2022MS003089. https://doi.org/10.1029/2022MS003089
 
 * Liao, Chang, Tian Zhou, Donghui Xu, Richard Barnes, Gautam Bisht, Hong-Yi Li, Zeli Tan, et al. (02/2022AD) 2022. “Advances In Hexagon Mesh-Based Flow Direction Modeling”. Advances In Water Resources 160. Elsevier BV: 104099. 
 https://doi.org/10.1016/j.advwatres.2021.104099.
diff --git a/pyhexwatershed/classes/pycase.py b/pyhexwatershed/classes/pycase.py
index 1d0956e..05e3051 100644
--- a/pyhexwatershed/classes/pycase.py
+++ b/pyhexwatershed/classes/pycase.py
@@ -1042,13 +1042,13 @@ def pyhexwatershed_export_all_polygon_variables(self):
                 sFilename_geojson = pBasin.sFilename_variable_polygon  
 
         if self.iMesh_type == 4: #mpas mesh
-            aVariable_json  = ['lSubbasin','dArea','dElevation','dSlope_between', 'dDrainage_area','dDistance_to_watershed_outlet'] #profile not enabled
-            aVariable_geojson = ['subbasin','area','elevation', 'slope', 'drainage_area','travel_distance']
+            aVariable_json  = ['lSubbasin','lHillslope','dArea','dElevation','dSlope_between', 'dDrainage_area','dDistance_to_watershed_outlet'] #profile not enabled
+            aVariable_geojson = ['subbasin','hillslope','area','elevation', 'slope', 'drainage_area','travel_distance']
         else:
-            aVariable_json  = ['lSubbasin','dArea','dElevation','dSlope_between', 'dDrainage_area','dDistance_to_watershed_outlet'] #profile not enabled            
-            aVariable_geojson = ['subbasin','area','elevation', 'slope', 'drainage_area','travel_distance']
+            aVariable_json  = ['lSubbasin','lHillslope','dArea','dElevation','dSlope_between', 'dDrainage_area','dDistance_to_watershed_outlet'] #profile not enabled            
+            aVariable_geojson = ['subbasin','hillslope','area','elevation', 'slope', 'drainage_area','travel_distance']
 
-        aVariable_type= [1,2,2,2,2,2]
+        aVariable_type= [1,1,2,2,2,2,2]
         export_json_to_geojson_polygon(sFilename_json,
                                         sFilename_geojson, 
                                         aVariable_json,