diff --git a/python/packages/nisar/products/insar/GOFF_writer.py b/python/packages/nisar/products/insar/GOFF_writer.py
index 9c6b033a1..029bbbd8a 100644
--- a/python/packages/nisar/products/insar/GOFF_writer.py
+++ b/python/packages/nisar/products/insar/GOFF_writer.py
@@ -138,8 +138,8 @@ def add_grids_to_hdf5(self):
                 yds=yds,
                 fill_value=255,
             )
-            offset_group['mask'].attrs['valid_min'] = 0
-            offset_group['mask'].attrs['percentage_water'] = 0.0
+            offset_group['mask'].attrs['valid_min'] = np.uint32(0)
+            offset_group['mask'].attrs['percentage_water'] = np.float32(0.0)
             offset_group['mask'].attrs['disclaimer'] = to_bytes(self.water_mask_source)
 
             pixeloffsets_group_name = \
diff --git a/python/packages/nisar/products/insar/GUNW_writer.py b/python/packages/nisar/products/insar/GUNW_writer.py
index cbe3c837b..ff55bf9d0 100644
--- a/python/packages/nisar/products/insar/GUNW_writer.py
+++ b/python/packages/nisar/products/insar/GUNW_writer.py
@@ -330,8 +330,8 @@ def add_grids_to_hdf5(self):
                     yds=yds,
                     fill_value=255,
                 )
-                ds_group['mask'].attrs['valid_min'] = 0
-                ds_group['mask'].attrs['percentage_water'] = 0.0
+                ds_group['mask'].attrs['valid_min'] = np.uint32(0)
+                ds_group['mask'].attrs['percentage_water'] = np.float32(0.0)
                 ds_group['mask'].attrs['disclaimer'] = to_bytes(self.water_mask_source)
 
             for pol in pol_list:
diff --git a/python/packages/nisar/products/insar/InSAR_L1_writer.py b/python/packages/nisar/products/insar/InSAR_L1_writer.py
index 1d122aae2..a12fa0190 100644
--- a/python/packages/nisar/products/insar/InSAR_L1_writer.py
+++ b/python/packages/nisar/products/insar/InSAR_L1_writer.py
@@ -508,7 +508,7 @@ def add_pixel_offsets_to_swaths_group(self):
                                                  " Bits 24-31 are reserved for future use"),
                                     fill_value=255)
             offset_group['mask'].attrs['long_name'] = to_bytes("Valid samples subswath and data anomaly mask")
-            offset_group['mask'].attrs['valid_min'] = 0
+            offset_group['mask'].attrs['valid_min'] = np.uint32(0)
 
             range_offset_path = \
                 os.path.join( self.topo_path,
@@ -738,7 +738,7 @@ def add_interferogram_to_swaths_group(self, is_unwrapped=False):
                                     dtype=np.uint32,
                                     description=mask_description,
                                     fill_value=255)
-            igram_group['mask'].attrs['valid_min'] = 0
+            igram_group['mask'].attrs['valid_min'] = np.uint32(0)
             igram_group['mask'].attrs['long_name'] = to_bytes("Valid samples subswath and data anomaly mask")
 
             range_offset_path = \
diff --git a/python/packages/nisar/workflows/compute_stats.py b/python/packages/nisar/workflows/compute_stats.py
index dfdfa0895..ee8afddf7 100644
--- a/python/packages/nisar/workflows/compute_stats.py
+++ b/python/packages/nisar/workflows/compute_stats.py
@@ -35,10 +35,10 @@ def compute_stats_real_hdf5_dataset(h5_ds):
     h5_ds: h5py.Dataset
         hdf5 dataset object
     """
-    h5_ds.attrs.create('min_value', data=np.nanmin(h5_ds))
-    h5_ds.attrs.create('mean_value', data=np.nanmean(h5_ds))
-    h5_ds.attrs.create('max_value', data=np.nanmax(h5_ds))
-    h5_ds.attrs.create('sample_stddev', data=np.nanstd(h5_ds))
+    h5_ds.attrs.create('min_value', data=h5_ds.dtype.type(np.nanmin(h5_ds)))
+    h5_ds.attrs.create('mean_value', data=h5_ds.dtype.type(np.nanmean(h5_ds)))
+    h5_ds.attrs.create('max_value', data=h5_ds.dtype.type(np.nanmax(h5_ds)))
+    h5_ds.attrs.create('sample_stddev', data=h5_ds.dtype.type(np.nanstd(h5_ds)))
 
 def compute_stats_real_data(raster, h5_ds):
     """
@@ -53,11 +53,11 @@ def compute_stats_real_data(raster, h5_ds):
         stats_obj = isce3.math.compute_raster_stats_float64(raster)[0]
     else:
         stats_obj = isce3.math.compute_raster_stats_float32(raster)[0]
-    h5_ds.attrs.create('min_value', data=stats_obj.min)
-    h5_ds.attrs.create('mean_value', data=stats_obj.mean)
-    h5_ds.attrs.create('max_value', data=stats_obj.max)
+    h5_ds.attrs.create('min_value', data=h5_ds.dtype.type(stats_obj.min))
+    h5_ds.attrs.create('mean_value', data=h5_ds.dtype.type(stats_obj.mean))
+    h5_ds.attrs.create('max_value', data=h5_ds.dtype.type(stats_obj.max))
     h5_ds.attrs.create('sample_stddev',
-                       data=stats_obj.sample_stddev)
+                       data=h5_ds.dtype.type(stats_obj.sample_stddev))
 
 def compute_layover_shadow_water_stats(h5_ds, lines_per_block=1000):
     '''