------------------------------------ VOLCAT Alerts Version Information ------------------------------------ $Id: volcat_alerts.f90,v 1.52 2017/01/18 21:13:39 mpav Exp $ $Id: volcat_alert_utils.f90,v 1.565 2017/03/24 00:46:52 mpav Exp $ $Id: volcat_general_utils.f90,v 1.305 2017/03/22 23:44:58 mpav Exp $ $Id: volcat_history_utils.f90,v 1.179 2017/06/06 00:50:12 mpav Exp $ $Id: volcat_structs.f90,v 1.191 2017/03/22 23:44:58 mpav Exp $ $Id: volcat_plume_utils.f90,v 1.31 2017/02/27 21:11:49 mpav Exp $ $Id: volcat_grouping_utils.f90,v 1.128 2017/05/01 22:02:56 mpav Exp $ $Id: volcat_lightning_utils.f90,v 1.9 2017/02/08 19:44:18 mpav Exp $ $Id: volcat_analysis_utils.f90,v 1.9 2017/02/08 19:44:18 mpav Exp $ $Id: volcat_incident_utils.f90,v 1.9 2017/02/08 19:44:18 mpav Exp $ VOLCAT_ALERTS>> Reading measurement platform file ------------------------------------ Measurement platform data version: ------------------------------------ $Id: volcat_sensor_list.txt,v 1.5 2017/01/11 19:39:40 mpav Exp $ VOLCAT_ALERTS>> Reading configuration file ---------------------------------------------------------- Program configuration read from: /output/tmp/volcat/volcat_alert_run_script.config ---------------------------------------------------------- Object Production System: GEOCAT Debug Option: 2 Use IR only: F Volcanic cloud warning log path: /output/tmp/volcat Volcanic cloud warning log filename: Volcanic_Cloud_Alert_Realtime Volcanic cloud warning log path and filename: /output/tmp/volcat/Volcanic_Cloud_Alert_Realtime Volcanic ash time series path: /output/time_series/ash/ Volcanic ice time series path: /output/time_series/ice/ Volcanic hot time series path: none Volcanic so2 time series path: none Volcanic cloud diagnostics path: /output/volcat_diag/ Volcanic cloud diagnostics filename: geocatHOTOBJ.SNPP.2025106.201243.diag Volcanic cloud diagnostics path and filename: /output/volcat_diag//geocatHOTOBJ.SNPP.2025106.201243.diag Volcano information path: /opt/geocat/data_algorithms/volcat_data/ Volcano information filename: Volcat_Alert_FireUSA_Database_Feb21_2022_165111.txt Volcano information path and filename: /opt/geocat/data_algorithms/volcat_data//alerts/Volcat_Alert_FireUSA_Database_Feb21_2022_165111.txt Volcano activity path: /root/git/volcat Volcano activity filename: MISSING Volcano activity path and filename: /root/git/volcat/MISSING Weak plume option for ash: 1 Strong plume option for ash: 0 Weak plume option for SO2: 1 Strong plume option for SO2: 0 Weak plume option for ash: 1 Maximum viewing angle (new ash detect): 80.0 Maximum viewing angle (new ice detect): 90.0 Maximum viewing angle (new hot detect): 90.0 Maximum viewing angle (new so2 detect): 80.0 Noisy ash object threshold (basic criteria): 0.25 Noisy ash object threshold (plume/puff criteria): 0.25 Noisy ash object threshold (neighbor criteria): 0.25 Noisy ash object threshold (temporal criteria): 1.01 Noisy so2 object threshold (basic criteria): 0.25 Noisy so2 object threshold (plume/puff criteria): 0.25 Noisy so2 object threshold (neighbor criteria): 0.25 Noisy so2 object threshold (temporal criteria): 1.01 Min median 11 um tot emissivity for ash object noise screening (basic criteria): 0.50 Min median 11 um tot emissivity for ash object noise screening (plume/puff criteria): 0.50 Min median 11 um tot emissivity for ash object noise screening (neighbor criteria): 0.50 Min median 11 um tot emissivity for ash object noise screening (temporal criteria): 1.00 Min median 11 um tot emissivity for so2 object noise screening (basic criteria): 0.50 Min median 11 um tot emissivity for so2 object noise screening (plume/puff criteria): 0.50 Min median 11 um tot emissivity for so2 object noise screening (neighbor criteria): 0.50 Min median 11 um tot emissivity for so2 object noise screening (temporal criteria): 1.00 Non-uniform ash object threshold (basic criteria): 0.75 Non-uniform ash object threshold (plume/puff criteria): 0.75 Non-uniform ash object threshold (neighbor criteria): 1.01 Non-uniform ash object threshold (temporal criteria): 1.01 Non-uniform so2 object threshold (basic criteria): 0.75 Non-uniform so2 object threshold (plume/puff criteria): 0.75 Non-uniform so2 object threshold (neighbor criteria): 1.01 Non-uniform so2 object threshold (temporal criteria): 1.01 Maximum ash object median tot emiss (basic criteria): 0.90 Maximum ash object median tot emiss (plume/puff criteria): 0.60 Maximum ash object median tot emiss (neighbor criteria): 0.75 Maximum ash object median tot emiss (temporal criteria): 1.01 Maximum so2 object median tot emiss (basic criteria): 0.90 Maximum so2 object median tot emiss (plume/puff criteria): 0.50 Maximum so2 object median tot emiss (neighbor criteria): 0.75 Maximum so2 object median tot emiss (temporal criteria): 1.01 Minimum ash object count (basic/neighbor criteria): 1000000000 Minimum ash object count (plume/puff criteria): 1000000000 Minimum ash object count (temporal criteria): 0 Minimum ice object count (basic/temporal criteria): 3 Minimum ice object count (temporal criteria): 0 Minimum hot object count: 0 Minimum so2 object count (basic/neighbor criteria): 5 Minimum so2 object count (plume/puff criteria): 3 Minimum so2 object count (temporal criteria): 0 Ash object grouping contrast option (basic and neighbor criteria): 1 Ash object grouping contrast option (plume/puff criteria): 1 So2 object grouping contrast option (basic and neighbor criteria): 1 So2 object grouping contrast option (plume/puff criteria): 1 Minimum number of fractional days between ash alerts: 0.50 Minimum number of fractional days between ice alerts: 0.50 Minimum number of fractional days between thermal anomaly alerts: 0.50 Minimum number of fractional days between So2 alerts: 0.50 The number of robust pixels that are needed to be confident that co-registration errors are not to blame (L1 robustness): 50 The number of robust pixels that are needed to be confident that co-registration errors are not to blame (L2 robustness): 50 The number of robust pixels that are needed to be confident that co-registration errors are not to blame (L3 robustness): 20 The number of robust pixels that are needed to be confident that co-registration errors are not to blame (L4 robustness): 10 The fraction of robust area possibly attributed to co-registration errors that gives confidence that robust signal is strictly due to co-registration (L1 robustness): .5000 The fraction of robust area possibly attributed to co-registration errors that gives confidence that robust signal is strictly due to co-registration (L2 robustness): .5000 The fraction of robust area possibly attributed to co-registration errors that gives confidence that robust signal is strictly due to co-registration (L3 robustness): .5000 The fraction of robust area possibly attributed to co-registration errors that gives confidence that robust signal is strictly due to co-registration (L4 robustness): .9500 History grid resolution (detection): .10 History grid resolution (retrievals): 1.00 History file path: /output/history/ History file filename: volcat_history.hdf History file path and filename: /output/history//volcat_history.hdf Ash history usage: 1 Ice history usage: 1 Hot history usage: 1 SO2 history usage: 1 Number of days between valid detections before thermal anomaly is removed from history: 15.00000 Maximum number of volcanoes listed in alerts: 5 Ash object is valid: NO Ash object path: /output/l2_data Ash object filename: none Ash object path and filename: /output/l2_data/none Ice object is valid: NO Ice object path: /output/l2_data Ice object filename: none Ice object path and filename: /output/l2_data/none Hot object is valid: YES Hot object path: /output/l2_data Hot object filename: geocatHOTOBJHIRES.SNPP.2025106.201243.hdf Hot object path and filename: /output/l2_data/geocatHOTOBJHIRES.SNPP.2025106.201243.hdf SO2 object is valid: NO SO2 object path: /output/l2_data SO2 object filename: none SO2 object path and filename: /output/l2_data/none ENI lightning is valid: NO ENI lightning path: none ENI lightning time match buffer (minutes): 15.00 ENI lightning is written to output: NO WWLLN lightning is valid: NO WWLLN lightning path: none WWLLN lightning time match buffer (minutes): 15.00 WWLLN lightning is written to output: NO GLD360 lightning is valid: NO GLD360 lightning path: none GLD360 lightning time match buffer (minutes): 15.00 GLD360 lightning is written to output: NO GLM lightning is valid: NO GLM lightning path: none GLM lightning time match buffer (minutes): 10.00 GLM lightning is written to output: NO Ash group expansion option: 0 Apply ash/ice mask: 1 Interp bad retv ash/ice: 1 Fill ash/ice holes: 1 Identify hot spot artifacts: 3 SO2 group expansion option: 0 Apply SO2 mask: 0 Interp bad retv SO2: 0 Fill SO2 holes: 0 Dust option: -3 Fire option: 1 Valid incidents: 1 Incidents: /ancillary/NIFC_DATA/Current-Fire-Incidents.txt VOLCAT_ALERTS>> Reading Volcano Information and Alert Criteria ------------------------------------ Volcano database version: ------------------------------------ N/A tl_2019_us_county.shp VOLCAT_ALERTS>> Reading Incident Data geocat>> VOLCAT_INCIDENT_UTILS:LOAD_INCIDENTS(INFORMATION) : Reading incident data from: /ancillary/NIFC_DATA/Current-Fire-Incidents.txt Incident database: /ancillary/NIFC_DATA/Current-Fire-Incidents.txt ----------------------------------------------------------------- Creation time: Fire incident data were extracted at 2025-04-16T19:47:17 UTC from and /ancillary/NIFC_DATA/Current-Fire-Location-2025-04-16_194712.geojson/ancillary/NIFC_DATA/Current-Fire-Perimeter-2025-04-16_194712.geojson Number of perimeters: 118 Number of points: 790 Total number of incidents: 908 Perimeter data bounding lon/lat: -113.319, 28.774, -70.293, 48.415 Point data bounding lon/lat: -146.165, 17.714, -66.794, 64.237 VOLCAT_ALERTS>> Reading Object Files VOLCAT_ALERTS>> Reading pixel navigation, if needed geocat>> VOLCAT_ALERT_UTILS:PIXEL_NAVIGATION(INFORMATION) : Reading hires navigation data from : /output//l1_data/geocatL1HIRES.SNPP.2025106.201243.hdf VOLCAT_ALERTS>> Reading pixel L1 data, if needed VOLCAT_ALERTS>> Reading pixel L2 data, if needed VOLCAT_ALERTS>> Ingesting lightning data, if needed VOLCAT_ALERTS>> Reading or Creating History, if needed VOLCAT_ALERTS>> Reading history file VOLCAT_ALERTS>> Identifying Alert Criteria Objects geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Assessing objects with volcanoes (HOT) geocat>> VOLCAT_ALERT_UTILS:HOT_OBJECT_SELECTION_BASIC(INFORMATION) : Eliminating thermal anomaly due to limited immediate contrast Obj Id: 2, Lat: 58.08839, Lon: -116.50, BT MWIR: 313.29 geocat>> VOLCAT_ALERT_UTILS:HOT_OBJECT_SELECTION_BASIC(INFORMATION) : Eliminating thermal anomaly due to limited immediate contrast Obj Id: 2, Lat: 58.08839, Lon: -116.50, BT MWIR: 313.29 geocat>> VOLCAT_ALERT_UTILS:HOT_OBJECT_SELECTION_BASIC(INFORMATION) : Eliminating thermal anomaly due to limited immediate contrast Obj Id: 2, Lat: 58.08839, Lon: -116.50, BT MWIR: 313.29 geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Assessing objects without volcanoes (HOT) geocat>> VOLCAT_ALERT_UTILS:HOT_OBJECT_SELECTION_BASIC(INFORMATION) : Eliminating thermal anomaly due to limited immediate contrast Obj Id: 2, Lat: 58.08839, Lon: -116.50, BT MWIR: 313.29 geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Assessing thermal objects for association with known incidents geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Resolving spatially redundant thermal anomalies geocat>> VOLCAT_GENERAL_UTILS:RESOLVE_REDUNDANT_VRP(INFORMATION) : Modifying radiative power to account for spatial overlap Elem_Idx: 5105, Line_Idx: 1421, Obj_Idx: 1, RP[MW]: 9.97, RP Modified[MW]: 9.97, Qc: 0, RP Final[MW]: 9.97, Qc Final: 1 Elem_Idx: 3226, Line_Idx: 1792, Obj_Idx: 2, RP[MW]: .95, RP Modified[MW]: .95, Qc: 0, RP Final[MW]: .95, Qc Final: 1 Elem_Idx: 1387, Line_Idx: 1795, Obj_Idx: 3, RP[MW]: 3.41, RP Modified[MW]: -999.00, Qc: 0, RP Final[MW]: 3.41, Qc Final: 1 Elem_Idx: 3314, Line_Idx: 2421, Obj_Idx: 4, RP[MW]: 2.70, RP Modified[MW]: 2.70, Qc: 0, RP Final[MW]: 2.70, Qc Final: 1 Elem_Idx: 4951, Line_Idx: 3127, Obj_Idx: 5, RP[MW]: 7.29, RP Modified[MW]: -999.00, Qc: 0, RP Final[MW]: 7.29, Qc Final: 1 Elem_Idx: 2610, Line_Idx: 3266, Obj_Idx: 6, RP[MW]: 2.17, RP Modified[MW]: 2.17, Qc: 0, RP Final[MW]: 2.17, Qc Final: 1 Elem_Idx: 4682, Line_Idx: 3820, Obj_Idx: 7, RP[MW]: 3.82, RP Modified[MW]: 3.82, Qc: 0, RP Final[MW]: 3.82, Qc Final: 1 Elem_Idx: 4164, Line_Idx: 4203, Obj_Idx: 8, RP[MW]: .82, RP Modified[MW]: .82, Qc: 0, RP Final[MW]: .82, Qc Final: 1 geocat>> VOLCAT_GENERAL_UTILS:RESOLVE_REDUNDANT_VRP(INFORMATION) : Building list of linked objects from spatial overlap GROUPS 0 0 0 geocat>> VOLCAT_GROUPING_UTILS:GROUP_SINGULAR_OBJECTS(INFORMATION) : Hot object grouping method: incident based alert 1 1 1 1 1 1 -999 1 geocat>> VOLCAT_ALERT_UTILS:ORGANIZE_GROUP_OUTPUT(INFORMATION) : Fitting polygon using convex hull (no pixel area) geocat>> VOLCAT_ALERT_UTILS:ORGANIZE_GROUP_OUTPUT(INFORMATION) : Fitting polygon using concave hull (no pixel area) geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Identifying hot spot artifacts geocat>> VOLCAT_GENERAL_UTILS:APPLY_VRP_QC(INFORMATION) : Applying standard QC MWIR resolution (m): 375.0, MWIR resolution factor: 5 Group_Id: 001, Object_Id: 0001, X_Idx: 5105, Y_Idx: 1421, VRP_Pix: 9.97, VRP_Max: 9.97, VRP_Relative_Diff: 0.00, Rel_Diff_Thresh: 0.26, Local_Max_VRP: 9.97, Internal: F, BT_Pix: 340.8, BT_Crit: T, Qc: 1 geocat>> VOLCAT_GENERAL_UTILS:APPLY_CLOUDY_VRP_QC(INFORMATION) : Applying cloud scattering QC MWIR resolution (m): 375.0, MWIR resolution factor: 1 geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Redistributing radiative power to limit energy smearing geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Recomputing group stats (hot) VOLCAT_ALERTS>> Correcting Retrievals, if needed geocat>> VOLCAT_ALERT_UTILS:CORRECT_RETRIEVALS(INFORMATION) : Recomputing group stats (ash) geocat>> VOLCAT_ALERT_UTILS:CORRECT_RETRIEVALS(INFORMATION) : Recomputing group stats (ice) geocat>> VOLCAT_ALERT_UTILS:CORRECT_RETRIEVALS(INFORMATION) : Thermal anomaly group radiative power summary: Group_Idx: 1, Orig Total RP [MW]: 9.97, New Total RP [MW]: 9.97 geocat>> VOLCAT_GENERAL_UTILS:FEATURE_RP_TO_SCALED_RASTER(INFORMATION) : Creating feature-based scaled radiative power raster geocat>> VOLCAT_ALERT_UTILS:CORRECT_RETRIEVALS(INFORMATION) : Updating thermal L2: /output//l2_data/geocatL2HIRES.SNPP.2025106.201243.hdf geocat>> VOLCAT_ALERT_UTILS:CORRECT_RETRIEVALS(INFORMATION) : Recomputing group stats (so2) VOLCAT_ALERTS>> Updating History, if needed Background= 0 0 54.01792 -128.6791 VOLCAT_ALERTS>> Checking For Alerts To Issue @@**********VOLCAT ALERTS********** STARTING DATE/TIME OF IMAGE: 2025-04-16 20:12:43 [UTC] PRODUCTION DATE/TIME OF ALERT: 2025-04-16 20:43:02 [UTC] PRIMARY INSTRUMENT: NPP VIIRS PRIMARY DOMAIN: none WMO SPACECRAFT ID: 342 LOCATION/ORBIT: LEO L1 FILE: SVM13_npp_d20250416_t2012431_e2024051_b69796_c20250416203503586837_cspp_dev.h5 L1 FILE(S) - TEMPORAL: NONE VOLCANO DATABASE: /opt/geocat/data_algorithms/volcat_data//alerts/Volcat_Alert_FireUSA_Database_Feb21_2022_165111.txt NUMBER OF ASH CLOUD ALERTS: 0 NUMBER OF Cb ALERTS: 0 NUMBER OF THERMAL ANOMALY ALERTS: 1 NUMBER OF SO2 CLOUD ALERTS: 0 DATABASE VERSION: 3.0 ALGORITHM VERSION: 3.0 --------------------------------------------- POSSIBLE FIRE IDENTIFIED Alert Status: New Feature Alert Confidence: HIGH Identification Method: Enhanced Contextual (Cloud) Latitude of Radiative Center: 54.02 [degrees] Longitude of Radiative Center: -128.68 [degrees] Mean Object Date/Time: 2025-04-16 20:14:02 [UTC] Mean Viewing Angle: 52.5 [degrees] Mean Solar Zenith Angle: 43.8 [degrees] Mean Glint Angle: 51.3 [degrees] Water Fraction: .0 [none] Nearby Locations: Ketchikan Gateway Borough(223.0 km) Prince of Wales-Hyder Census A(341.8 km) Wrangell City and Borough(346.1 km) Petersburg Borough(439.1 km) Number of background pixels: 649 [none] Background infrared distance threshold: 10.0 [K] Fraction of water pixels in background: .12 [none] Maximum Observed BT (3.9 um): 340.80 [K] Observed - Background BT (3.9 um): 55.58 [K] BT Anomaly (3.9 um): 27.62 [K] Maximum Observed - Clear RTM BT (3.9 um): 58.17 [K] Clear RTM BT Diff (3.9 um - 10/11/12 um): 49.61 [K] Maximum Observed BTD (3.9-10/11/12 um): 50.40 [K] Observed - Background BTD (3.9-10/11/12 um): 47.06 [K] BTD Anomaly (3.9-10/11/12 um): 38.35 [K] Maximum Observed Radiance (1.6 um): 4.5032 [RU] Observed - Background Radiance (1.6 um): 3.4217 [RU] Radiance Anomaly (1.6 um): 11.5220 [RU] Maximum Observed Radiance (2.2 um): 1.6973 [RU] Observed - Background Radiance (2.2 um): 1.3718 [RU] Radiance Anomaly (2.2 um): 12.4229 [RU] Likely an instrument artifact (Y or N): N [none] Likely a solar farm reflection (Strong, Moderate, Weak, or No): No [none] BT Time Tendency (3.9 um): Not Available BT Time Tendency Diff (3.9 -10/11/12 um): Not Available Image Interval: Not Available Fraction of Surrounding LWIR Pixels that are Colder: 1.00 [fraction] Fraction of Surrounding Pixels with a Larger LWIR BTD: .94 [fraction] Fraction of Surrounding Red Channel Pixels that are Brighter: .41 [fraction] The red band reflectance: 12.6 [%] Time since last detection: -999.000000 [days] Index of high confidence criteria: 3 [none] Maximum Radiative Power: 9.97 [MW] Maximum Radiative Power Uncertainty: .00 [MW] Total Radiative Power: 9.97 [MW] Total Radiative Power Uncertainty: .00 [MW] Total Radiative Power Anomaly: -999.00 [standard deviations] Total Pixel Area: 0.23 [km^2] Type of Event: MISSING Event Priority Ranking: 0 Country/Countries: USA State/Province(s): Alaska County/Locality(s): Ketchikan Gateway,Prince of Wales-Hyder,Wrangell,Petersburg NWS WFO: NA Interagency Region: NA GACC Region: NA --------------------------------------------- NUMBER OF VALID NON-ALERT ASH GROUPS: 0 NUMBER OF VALID NON-ALERT VOLCANIC CB GROUPS: 0 NUMBER OF VALID NON-ALERT SO2 GROUPS: 0 NUMBER OF VALID NON-ALERT THERMAL ANOMALIES: 0 NUMBER OF VOLCANIC ASH GROUPS MERGED INTO VOLCANIC CB GROUP: 0 NUMBER OF CB GROUPS MERGED INTO ASH GROUP: 0 VOLCAT_ALERTS>> Writing Object Masks geocat>> Attribute Max_Groups reading error from id: 393216 geocat>> VOLCAT_ALERT_UTILS:WRITE_INCIDENT_CSV_LUT(INFORMATION) : Writing incident CSV: /output/l2_data/geocatHOTOBJHIRES.SNPP.2025106.201243.INCIDENT_LUT.csv VOLCAT_ALERTS>> Generating Volcanic Cloud Time Series VOLCAT_ALERTS>> Creating Diagnostic Output VOLCAT_ALERTS>> Cleaning Up geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up volcanoes geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up incidents geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up objects geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up nav geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up hires nav geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up L1 geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up RTM geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up L2 geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up hires L2 geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up Lightning geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up history geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up alerts geocat>> VOLCAT_ALERT_UTILS:CLEAN_UP_ALERTS(INFORMATION) : Cleaning up selected objects VOLCAT_ALERTS>> Total processing time: Days = 0 Hours = 0 Minutes = 0 Seconds = 7 VOLCAT_ALERTS>> END OF FULL PROGRAM