------------------------------------ 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.NOAA-21.2025079.231336.diag Volcanic cloud diagnostics path and filename: /output/volcat_diag//geocatHOTOBJ.NOAA-21.2025079.231336.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.NOAA-21.2025079.231336.hdf Hot object path and filename: /output/l2_data/geocatHOTOBJHIRES.NOAA-21.2025079.231336.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-03-20T23:23:30 UTC from and /ancillary/NIFC_DATA/Current-Fire-Location-2025-03-20_232326.geojson/ancillary/NIFC_DATA/Current-Fire-Perimeter-2025-03-20_232326.geojson Number of perimeters: 63 Number of points: 738 Total number of incidents: 801 Perimeter data bounding lon/lat: -115.157, 13.129, 144.974, 48.233 Point data bounding lon/lat: -150.408, 13.132, 144.972, 67.568 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.NOAA-21.2025079.231336.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:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Assessing objects without volcanoes (HOT) 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: 99, Line_Idx: 4987, RP[MW]: 6.42, Elem_Over: 98, Line_Over: 5004, RP_Over[MW]: 4.51, Elem_Traceback: 99, Line_Traceback: 4987 Ref_Obj_Idx: 3, Over_Obj_Idx: 4, Ref_Overlap_Count: 1, Over_Overlap_Count: 1 Elem_Idx: 2081, Line_Idx: 3973, Obj_Idx: 1, RP[MW]: 1.31, RP Modified[MW]: 1.31, Qc: 0, RP Final[MW]: 1.31, Qc Final: 1 Elem_Idx: 2081, Line_Idx: 3974, Obj_Idx: 1, RP[MW]: 1.31, RP Modified[MW]: 1.31, Qc: 0, RP Final[MW]: 1.31, Qc Final: 1 Elem_Idx: 1949, Line_Idx: 4329, Obj_Idx: 2, RP[MW]: 4.65, RP Modified[MW]: 4.65, Qc: 0, RP Final[MW]: 4.65, Qc Final: 1 Elem_Idx: 99, Line_Idx: 4987, Obj_Idx: 3, RP[MW]: 6.42, RP Modified[MW]: 5.38, Qc: 4, RP Final[MW]: 5.38, Qc Final: 4 Elem_Idx: 98, Line_Idx: 5004, Obj_Idx: 4, RP[MW]: 4.51, RP Modified[MW]: .00, Qc: 5, RP Final[MW]: .00, Qc Final: 5 geocat>> VOLCAT_GENERAL_UTILS:RESOLVE_REDUNDANT_VRP(INFORMATION) : Building list of linked objects from spatial overlap Original List (Ref ID=3): Object 1 of 1 --> Id = 4 Unified List (Ref ID=3): Object 1 of 1 --> Id = 4 GROUPS 0 0 0 geocat>> VOLCAT_GROUPING_UTILS:GROUP_SINGULAR_OBJECTS(INFORMATION) : Hot object grouping method: incident based alert 1 1 3 60.72554 -137.3129 4 60.72352 -137.3084 Adding object 4 to group 1 2 2 1 1 1 3 -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: 0003, X_Idx: 0099, Y_Idx: 4987, VRP_Pix: 5.38, VRP_Max: 6.42, VRP_Relative_Diff: 0.16, Rel_Diff_Thresh: 0.21, Local_Max_VRP: 5.38, Internal: F, BT_Pix: 325.0, 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 Group Id: 1, Obj_Idx: 3, Trans_Min: .637, VRP Total: 6.42 Group Id: 1, Obj_Idx: 4, Trans_Min: .622, VRP Total: 4.51 geocat>> VOLCAT_ALERT_UTILS:IDENTIFY_ALERT_OBJECTS(INFORMATION) : Redistributing radiative power to limit energy smearing geocat>> VOLCAT_HISTORY_UTILS:GROUP2HISTORY(INFORMATION) : Could not find an open list slot 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]: 10.93, New Total RP [MW]: 5.38 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.NOAA-21.2025079.231336.hdf geocat>> VOLCAT_ALERT_UTILS:CORRECT_RETRIEVALS(INFORMATION) : Recomputing group stats (so2) VOLCAT_ALERTS>> Updating History, if needed Background= 0 0 60.72554 -137.3129 VOLCAT_ALERTS>> Checking For Alerts To Issue @@**********VOLCAT ALERTS********** STARTING DATE/TIME OF IMAGE: 2025-03-20 23:13:36 [UTC] PRODUCTION DATE/TIME OF ALERT: 2025-03-20 23:46:38 [UTC] PRIMARY INSTRUMENT: NOAA-21 VIIRS PRIMARY DOMAIN: none WMO SPACECRAFT ID: 209 LOCATION/ORBIT: LEO L1 FILE: SVM13_j02_d20250320_t2313355_e2326238_b12225_c20250320233945120090_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: 60.73 [degrees] Longitude of Radiative Center: -137.31 [degrees] Mean Object Date/Time: 2025-03-20 23:18:11 [UTC] Mean Viewing Angle: 68.0 [degrees] Mean Solar Zenith Angle: 64.8 [degrees] Mean Glint Angle: 104.2 [degrees] Water Fraction: .0 [none] Nearby Locations: Skagway Municipality(169.7 km) Yakutat City and Borough(188.1 km) Haines Borough(205.4 km) Hoonah-Angoon Census Area(280.8 km) Juneau City and Borough(315.9 km) Number of background pixels: 1320 [none] Background infrared distance threshold: 5.0 [K] Fraction of water pixels in background: .09 [none] Maximum Observed BT (3.9 um): 324.95 [K] Observed - Background BT (3.9 um): 49.24 [K] BT Anomaly (3.9 um): 25.77 [K] Maximum Observed - Clear RTM BT (3.9 um): 53.79 [K] Clear RTM BT Diff (3.9 um - 10/11/12 um): 51.16 [K] Maximum Observed BTD (3.9-10/11/12 um): 51.41 [K] Observed - Background BTD (3.9-10/11/12 um): 47.88 [K] BTD Anomaly (3.9-10/11/12 um): 30.76 [K] Maximum Observed Radiance (1.6 um): 1.3219 [RU] Observed - Background Radiance (1.6 um): .2913 [RU] Radiance Anomaly (1.6 um): 1.4348 [RU] Maximum Observed Radiance (2.2 um): .4875 [RU] Observed - Background Radiance (2.2 um): .0816 [RU] Radiance Anomaly (2.2 um): .8165 [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: .92 [fraction] Fraction of Surrounding Pixels with a Larger LWIR BTD: .62 [fraction] Fraction of Surrounding Red Channel Pixels that are Brighter: 1.00 [fraction] The red band reflectance: 15.8 [%] Time since last detection: -999.000000 [days] Index of high confidence criteria: 3 [none] Maximum Radiative Power: 5.38 [MW] Maximum Radiative Power Uncertainty: .00 [MW] Total Radiative Power: 5.38 [MW] Total Radiative Power Uncertainty: .00 [MW] Total Radiative Power Anomaly: -999.00 [standard deviations] Total Pixel Area: 0.75 [km^2] Type of Event: MISSING Event Priority Ranking: 0 Country/Countries: USA State/Province(s): Alaska County/Locality(s): Skagway,Yakutat,Haines,Hoonah-Angoon,Juneau 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.NOAA-21.2025079.231336.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 = 9 VOLCAT_ALERTS>> END OF FULL PROGRAM