diff --git a/jobs/JGLOBAL_ARCHIVE b/jobs/JGLOBAL_ARCHIVE
index b60b73c6df..1ea8fdfeef 100755
--- a/jobs/JGLOBAL_ARCHIVE
+++ b/jobs/JGLOBAL_ARCHIVE
@@ -12,7 +12,7 @@ export CDUMP=${RUN/enkf}
 YMD=${PDY} HH=${cyc} generate_com -rx COM_ATMOS_ANALYSIS COM_ATMOS_BUFR COM_ATMOS_GEMPAK \
   COM_ATMOS_GENESIS COM_ATMOS_HISTORY COM_ATMOS_INPUT COM_ATMOS_MASTER COM_ATMOS_RESTART \
   COM_ATMOS_TRACK COM_ATMOS_WMO \
-  COM_CHEM_HISTORY \
+  COM_CHEM_HISTORY COM_CHEM_ANALYSIS\
   COM_MED_RESTART \
   COM_ICE_HISTORY COM_ICE_INPUT COM_ICE_RESTART \
   COM_OBS COM_TOP \
diff --git a/scripts/exglobal_archive.sh b/scripts/exglobal_archive.sh
index 2695799613..5fea07f4ed 100755
--- a/scripts/exglobal_archive.sh
+++ b/scripts/exglobal_archive.sh
@@ -33,6 +33,7 @@ source "${HOMEgfs}/ush/file_utils.sh"
 
 [[ ! -d ${ARCDIR} ]] && mkdir -p "${ARCDIR}"
 nb_copy "${COM_ATMOS_ANALYSIS}/${APREFIX}gsistat" "${ARCDIR}/gsistat.${RUN}.${PDY}${cyc}"
+nb_copy "${COM_CHEM_ANALYSIS}/${APREFIX}aerostat" "${ARCDIR}/aerostat.${RUN}.${PDY}${cyc}"
 nb_copy "${COM_ATMOS_GRIB_1p00}/${APREFIX}pgrb2.1p00.anl" "${ARCDIR}/pgbanl.${RUN}.${PDY}${cyc}.grib2"
 
 # Archive 1 degree forecast GRIB2 files for verification
diff --git a/ush/hpssarch_gen.sh b/ush/hpssarch_gen.sh
index 0086fd9c19..8707552ca1 100755
--- a/ush/hpssarch_gen.sh
+++ b/ush/hpssarch_gen.sh
@@ -344,6 +344,9 @@ if [[ ${type} == "gdas" ]]; then
     if [[ -s "${COM_ATMOS_ANALYSIS}/${head}oznstat" ]]; then
        echo "${COM_ATMOS_ANALYSIS/${ROTDIR}\//}/${head}oznstat"
     fi
+    if [[ -s "${COM_CHEM_ANALYSIS}/${head}aerostat" ]]; then
+       echo "${COM_CHEM_ANALYSIS/${ROTDIR}\//}/${head}aerostat"
+    fi
     if [[ -s "${COM_ATMOS_ANALYSIS}/${head}radstat" ]]; then
        echo "${COM_ATMOS_ANALYSIS/${ROTDIR}\//}/${head}radstat"
     fi
diff --git a/ush/python/pygfs/task/aero_analysis.py b/ush/python/pygfs/task/aero_analysis.py
index 722ebf0760..8bfcfd9ed3 100644
--- a/ush/python/pygfs/task/aero_analysis.py
+++ b/ush/python/pygfs/task/aero_analysis.py
@@ -279,7 +279,7 @@ def get_berror_dict(self, config: Dict[str, Any]) -> Dict[str, List[str]]:
         b_datestr = to_fv3time(config.BERROR_DATE)
         berror_list = []
 
-        for ftype in ['cor_rh', 'cor_rv', 'stddev']:
+        for ftype in ['stddev']:
             coupler = f'{b_datestr}.{ftype}.coupler.res'
             berror_list.append([
                 os.path.join(b_dir, coupler), os.path.join(config.DATA, 'berror', coupler)
@@ -290,7 +290,10 @@ def get_berror_dict(self, config: Dict[str, Any]) -> Dict[str, List[str]]:
                 berror_list.append([
                     os.path.join(b_dir, tracer), os.path.join(config.DATA, 'berror', tracer)
                 ])
-
+        radius = 'cor_aero_universe_radius'
+        berror_list.append([
+            os.path.join(b_dir, radius), os.path.join(config.DATA, 'berror', radius)
+        ])
         nproc = config.ntiles * config.layout_x * config.layout_y
         for nn in range(1, nproc + 1):
             berror_list.append([