hestia-earth-models 0.65.3__py3-none-any.whl → 0.65.5__py3-none-any.whl

hestia_earth/models/ipcc2019/co2ToAirCarbonStockChange_utils.py

@@ -47,7 +47,6 @@ _NOMINAL_ERROR = 75
 Carbon stock measurements without an associated `sd` should be assigned a nominal error of 75% (2*sd as a percentage of
 the mean).
 """
-_DEFAULT_YEARS_SINCE_LUC_EVENT = 999
 _TRANSITION_PERIOD_YEARS = 20
 _TRANSITION_PERIOD_DAYS = 20 * YEAR  # 20 years in days
 _VALID_DATE_FORMATS = {
@@ -96,6 +95,7 @@ class _InventoryKey(Enum):
     LAND_USE_SUMMARY = "land-use-summary"
     LAND_USE_CHANGE_EVENT = "luc-event"
     YEARS_SINCE_LUC_EVENT = "years-since-luc-event"
+    YEARS_SINCE_INVENTORY_START = "years-since-inventory-start"
 
 
 CarbonStock = NamedTuple("CarbonStock", [
@@ -611,7 +611,6 @@ def _create_compile_inventory_function(
         )
 
         inventory = _squash_inventory(
-            cycle_id,
             cycle_inventory,
             carbon_stock_inventory,
             land_use_inventory,
@@ -1008,14 +1007,20 @@ def _compile_land_use_inventory(
         is_luc_event = detect_land_use_change_func(land_use_summary, prev_land_use_summary)
 
         time_delta = current_year - prev_year
-        prev_years_since_luc_event = result.get(prev_year, {}).get(_InventoryKey.YEARS_SINCE_LUC_EVENT, 0)
+        prev_years_since_luc_event = (
+            result.get(prev_year, {}).get(_InventoryKey.YEARS_SINCE_LUC_EVENT, _TRANSITION_PERIOD_YEARS)
+        )
+        prev_years_since_inventory_start = result.get(prev_year, {}).get(_InventoryKey.YEARS_SINCE_INVENTORY_START, 0)
+
         years_since_luc_event = time_delta if is_luc_event else prev_years_since_luc_event + time_delta
+        years_since_inventory_start = prev_years_since_inventory_start + time_delta
 
         update_dict = {
             current_year: {
                 _InventoryKey.LAND_USE_SUMMARY: land_use_summary,
                 _InventoryKey.LAND_USE_CHANGE_EVENT: is_luc_event,
                 _InventoryKey.YEARS_SINCE_LUC_EVENT: years_since_luc_event,
+                _InventoryKey.YEARS_SINCE_INVENTORY_START: years_since_inventory_start
             }
         }
         return result | update_dict
@@ -1030,9 +1035,7 @@ def _compile_land_use_inventory(
             start_year: {
                 _InventoryKey.LAND_USE_SUMMARY: summarise_land_use_func(
                     land_cover_nodes_by_year.get(start_year, [])
-                ),
-                _InventoryKey.LAND_USE_CHANGE_EVENT: False,
-                _InventoryKey.YEARS_SINCE_LUC_EVENT: _DEFAULT_YEARS_SINCE_LUC_EVENT
+                )
             }
         }
     ) if should_run else {}
@@ -1062,7 +1065,6 @@ def _sorted_merge(*sources: Union[dict, list[dict]]) -> dict:
 
 
 def _squash_inventory(
-    cycle_id: str,
     cycle_inventory: dict,
     carbon_stock_inventory: dict,
     land_use_inventory: dict,
@@ -1075,8 +1077,6 @@ def _squash_inventory(
 
     Parameters
     ----------
-    cycle_id : str
-        The unique identifier of the cycle being processed.
 
     cycle_inventory : dict
         A dictionary representing the share of emissions for each cycle, grouped by year.
@@ -1162,12 +1162,7 @@ def _squash_inventory(
     )))
 
     def should_run_group(method: MeasurementMethodClassification, year: int) -> bool:
-        carbon_stock_inventory_group = carbon_stock_inventory.get(method, {}).get(year, {})
-        share_of_emissions_group = cycle_inventory.get(year, {})
-
-        has_emission = _InventoryKey.CO2_EMISSION in carbon_stock_inventory_group.keys()
-        is_relevant_for_cycle = cycle_id in share_of_emissions_group.get(_InventoryKey.SHARE_OF_EMISSION, {}).keys()
-        return all([has_emission, is_relevant_for_cycle])
+        return _InventoryKey.CO2_EMISSION in carbon_stock_inventory.get(method, {}).get(year, {}).keys()
 
     def squash(result: dict, year: int) -> dict:
         update_dict = next(
@@ -1175,7 +1170,10 @@ def _squash_inventory(
                 {
                     year: {
                         **_get_land_use_change_data(year, land_use_inventory),
-                        **reduce(merge, [carbon_stock_inventory[method][year], cycle_inventory[year]], dict())
+                        **reduce(merge, [
+                            carbon_stock_inventory.get(method, {}).get(year, {}),
+                            cycle_inventory.get(year, {})
+                        ], dict())
                     }
                 } for method in measurement_method_ranking if should_run_group(method, year)
             ),
@@ -1205,15 +1203,22 @@ def _get_land_use_change_data(
         )
     )
 
-    delta_time = closest_inventory_year - year if closest_inventory_year else 0
-    inventory_data = land_use_inventory.get(closest_inventory_year, {})
+    time_delta = closest_inventory_year - year if closest_inventory_year else 0
+    prev_years_since_luc_event = (
+        land_use_inventory.get(closest_inventory_year, {}).get(_InventoryKey.YEARS_SINCE_LUC_EVENT)
+    )
+    prev_years_since_inventory_start = (
+        land_use_inventory.get(closest_inventory_year, {}).get(_InventoryKey.YEARS_SINCE_INVENTORY_START)
+    )
 
-    years_since_luc_event = (
-        inventory_data.get(_InventoryKey.YEARS_SINCE_LUC_EVENT, _DEFAULT_YEARS_SINCE_LUC_EVENT) - delta_time
+    years_since_luc_event = prev_years_since_luc_event - time_delta if prev_years_since_luc_event else None
+    years_since_inventory_start = (
+        prev_years_since_inventory_start - time_delta if prev_years_since_inventory_start else None
     )
 
     return {
-        _InventoryKey.YEARS_SINCE_LUC_EVENT: years_since_luc_event
+        _InventoryKey.YEARS_SINCE_LUC_EVENT: years_since_luc_event,
+        _InventoryKey.YEARS_SINCE_INVENTORY_START: years_since_inventory_start
     }
 
 
@@ -1408,25 +1413,54 @@ def create_run_function(
         """
         data = inventory[year]
         years_since_luc_event = data[_InventoryKey.YEARS_SINCE_LUC_EVENT]
-        emission_term_id = (
-            land_use_change_emission_term_id if years_since_luc_event <= _TRANSITION_PERIOD_YEARS
-            else management_change_emission_term_id
-        )
+        years_since_inventory_start = data[_InventoryKey.YEARS_SINCE_INVENTORY_START]
+
+        is_luc_emission = bool(years_since_luc_event) and years_since_luc_event <= _TRANSITION_PERIOD_YEARS
+        is_data_complete = bool(years_since_inventory_start) and years_since_inventory_start > _TRANSITION_PERIOD_YEARS
+
+        if is_luc_emission:
+            # If LUC emission allocate emissions to land use change AND add corresponding zero emission to management
+            emission_term_id = land_use_change_emission_term_id
+            zero_emission_term_id = management_change_emission_term_id
+        elif is_data_complete:
+            # If management emission && data complete allocate emissions to management change AND add corresponding
+            # zero emission to management
+            emission_term_id = management_change_emission_term_id
+            zero_emission_term_id = land_use_change_emission_term_id
+        else:
+            # If management emission, but not data complete allocate emissions to management change only
+            emission_term_id = management_change_emission_term_id
+            zero_emission_term_id = None
 
         rescaled_emission = _rescale_carbon_stock_change_emission(
             data[_InventoryKey.CO2_EMISSION], data[_InventoryKey.SHARE_OF_EMISSION][cycle_id]
         )
 
+        zero_emission = CarbonStockChangeEmission(
+            value=array(0),
+            start_date=rescaled_emission.start_date,
+            end_date=rescaled_emission.end_date,
+            method=rescaled_emission.method
+        ) if zero_emission_term_id else None
+
         previous_emission = result.get(emission_term_id)
+        previous_zero_emission = result.get(zero_emission_term_id)
 
-        update_dict = {
+        emission_dict = {
             emission_term_id: (
                 _add_carbon_stock_change_emissions(previous_emission, rescaled_emission) if previous_emission
                 else rescaled_emission
             )
         }
 
-        return result | update_dict
+        zero_emission_dict = {
+            zero_emission_term_id: (
+                _add_carbon_stock_change_emissions(previous_zero_emission, zero_emission) if previous_zero_emission
+                else zero_emission
+            )
+        } if zero_emission_term_id else {}
+
+        return result | emission_dict | zero_emission_dict
 
     def run(cycle_id: str, cycle_start_date: str, cycle_end_date: str, inventory: dict) -> list[dict]:
         """
@@ -1447,9 +1481,13 @@ def create_run_function(
         list[dict]
             A list of [Emission](https://www.hestia.earth/schema/Emission) nodes containing model results.
         """
+
+        def should_run_year(year: int) -> bool:
+            return cycle_id in inventory.get(year, {}).get(_InventoryKey.SHARE_OF_EMISSION, {}).keys()
+
         assigned_emissions = reduce(
             lambda result, year: reduce_emissions(result, year, cycle_id, inventory),
-            inventory.keys(),
+            (year for year in inventory.keys() if should_run_year(year)),
             {}
         )
 

hestia_earth/models/ipcc2019/croppingDuration.py

@@ -1,3 +1,8 @@
+"""
+`croppingDuration` can be added by the uploader for all crops, however the model `IPCC (2019)` will only run for
+`riceGrainInHusk` or `ricePlantFlooded`.
+The model will only run for rice crops as this is the only crop which requires the value for emission recalculations.
+"""
 from hestia_earth.schema import PracticeStatsDefinition
 from hestia_earth.utils.lookup import download_lookup, get_table_value, column_name
 from hestia_earth.utils.tools import safe_parse_float

hestia_earth/models/ipcc2019/organicCarbonPerHa_tier_1_utils.py

@@ -219,18 +219,6 @@ _SOC_REFS = {
 }
 
 _LAND_USE_FACTORS = {
-    IpccLandUseCategory.GRASSLAND: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
-    },
     IpccLandUseCategory.PERENNIAL_CROPS: {
         EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 0.72, "error": 22},
         EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 0.72, "error": 22},
@@ -290,42 +278,6 @@ _LAND_USE_FACTORS = {
         EcoClimateZone.TROPICAL_WET: {"value": 0.82, "error": 17},
         EcoClimateZone.TROPICAL_MOIST: {"value": 0.82, "error": 17},
         EcoClimateZone.TROPICAL_DRY: {"value": 0.93, "error": 11}
-    },
-    IpccLandUseCategory.FOREST: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
-    },
-    IpccLandUseCategory.NATIVE: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
-    },
-    IpccLandUseCategory.OTHER: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
     }
 }
 
@@ -366,30 +318,6 @@ _MANAGEMENT_FACTORS = {
         EcoClimateZone.TROPICAL_MOIST: {"value": 0.9, "error": 8},
         EcoClimateZone.TROPICAL_DRY: {"value": 0.9, "error": 8}
     },
-    IpccManagementCategory.NOMINALLY_MANAGED: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
-    },
-    IpccManagementCategory.FULL_TILLAGE: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
-    },
     IpccManagementCategory.REDUCED_TILLAGE: {
         EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1.05, "error": 4},
         EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 0.99, "error": 3},
@@ -413,18 +341,6 @@ _MANAGEMENT_FACTORS = {
         EcoClimateZone.TROPICAL_WET: {"value": 1.1, "error": 5},
         EcoClimateZone.TROPICAL_MOIST: {"value": 1.1, "error": 5},
         EcoClimateZone.TROPICAL_DRY: {"value": 1.04, "error": 7}
-    },
-    IpccManagementCategory.OTHER: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
     }
 }
 
@@ -441,18 +357,6 @@ _CARBON_INPUT_FACTORS = {
         EcoClimateZone.TROPICAL_MOIST: {"value": 1.11, "error": 7},
         EcoClimateZone.TROPICAL_DRY: {"value": 1.11, "error": 7}
     },
-    IpccCarbonInputCategory.GRASSLAND_MEDIUM: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
-    },
     IpccCarbonInputCategory.CROPLAND_HIGH_WITH_MANURE: {
         EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1.44, "error": 13},
         EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1.37, "error": 12},
@@ -477,18 +381,6 @@ _CARBON_INPUT_FACTORS = {
         EcoClimateZone.TROPICAL_MOIST: {"value": 1.11, "error": 10},
         EcoClimateZone.TROPICAL_DRY: {"value": 1.04, "error": 13}
     },
-    IpccCarbonInputCategory.CROPLAND_MEDIUM: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
-    },
     IpccCarbonInputCategory.CROPLAND_LOW: {
         EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 0.92, "error": 14},
         EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 0.95, "error": 13},
@@ -500,18 +392,6 @@ _CARBON_INPUT_FACTORS = {
         EcoClimateZone.TROPICAL_WET: {"value": 0.92, "error": 14},
         EcoClimateZone.TROPICAL_MOIST: {"value": 0.92, "error": 14},
         EcoClimateZone.TROPICAL_DRY: {"value": 0.95, "error": 13}
-    },
-    IpccCarbonInputCategory.OTHER: {
-        EcoClimateZone.WARM_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.WARM_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_MOIST: {"value": 1},
-        EcoClimateZone.COOL_TEMPERATE_DRY: {"value": 1},
-        EcoClimateZone.BOREAL_MOIST: {"value": 1},
-        EcoClimateZone.BOREAL_DRY: {"value": 1},
-        EcoClimateZone.TROPICAL_MONTANE: {"value": 1},
-        EcoClimateZone.TROPICAL_WET: {"value": 1},
-        EcoClimateZone.TROPICAL_MOIST: {"value": 1},
-        EcoClimateZone.TROPICAL_DRY: {"value": 1}
     }
 }
 
@@ -556,7 +436,7 @@ def _sample_parameter(
         The sampled parameter as a numpy array with shape `(1, iterations)`.
     """
     parameter_dict = _IPCC_CATEGORY_TO_FACTOR_DICT.get(type(parameter))
-    kwargs = parameter_dict.get(parameter, {}).get(eco_climate_zone, {})
+    kwargs = parameter_dict.get(parameter, {}).get(eco_climate_zone, {"value": 1})
     func = _get_sample_func(kwargs)
     return func(iterations=iterations, seed=seed, **kwargs)
 
@@ -1008,7 +888,7 @@ def _compile_inventory(
 
     # If no `landCover` nodes in `site.management` use `site.siteType` to assign static `IpccLandUseCategory`.
     run_with_site_type = len(filter_list_term_type(management_nodes, [TermTermType.LANDCOVER])) == 0
-    site_type_ipcc_land_use_category = SITE_TYPE_TO_IPCC_LAND_USE_CATEGORY.get(site_type, IpccLandUseCategory.OTHER)
+    site_type_ipcc_land_use_category = SITE_TYPE_TO_IPCC_LAND_USE_CATEGORY.get(site_type, IpccLandUseCategory.UNKNOWN)
 
     grouped_management = group_nodes_by_year(management_nodes)
 
@@ -1213,7 +1093,7 @@ def _assign_ipcc_land_use_category(
         Assigned IPCC land use category.
     """
     DECISION_TREE = _LAND_USE_CATEGORY_DECISION_TREE
-    DEFAULT = IpccLandUseCategory.OTHER
+    DEFAULT = IpccLandUseCategory.UNKNOWN
 
     cover_crop_nodes, land_cover_nodes = split_on_condition(
         filter_list_term_type(management_nodes, [TermTermType.LANDCOVER]),
@@ -1389,7 +1269,7 @@ def _assign_ipcc_management_category(
     """
     decision_tree = _IPCC_LAND_USE_CATEGORY_TO_DECISION_TREE.get(ipcc_land_use_category, {})
     default = _IPCC_LAND_USE_CATEGORY_TO_DEFAULT_IPCC_MANAGEMENT_CATEGORY.get(
-        ipcc_land_use_category, IpccManagementCategory.OTHER
+        ipcc_land_use_category, IpccManagementCategory.NOT_RELEVANT
     )
 
     land_cover_nodes = filter_list_term_type(management_nodes, [TermTermType.LANDCOVER])
@@ -1431,10 +1311,10 @@ def _check_grassland_ipcc_management_category(
     bool
         `True` if the conditions match the specified management category, `False` otherwise.
     """
-    target_term_id = IPCC_MANAGEMENT_CATEGORY_TO_GRASSLAND_MANAGEMENT_TERM_ID.get(key, None)
+    target_term_ids = IPCC_MANAGEMENT_CATEGORY_TO_GRASSLAND_MANAGEMENT_TERM_ID.get(key, None)
     return cumulative_nodes_term_match(
         land_cover_nodes,
-        target_term_ids=target_term_id,
+        target_term_ids=target_term_ids,
         cumulative_threshold=MIN_AREA_THRESHOLD
     )
 
@@ -1443,8 +1323,7 @@ _GRASSLAND_IPCC_MANAGEMENT_CATEGORY_DECISION_TREE = {
     IpccManagementCategory.SEVERELY_DEGRADED: _check_grassland_ipcc_management_category,
     IpccManagementCategory.IMPROVED_GRASSLAND: _check_grassland_ipcc_management_category,
     IpccManagementCategory.HIGH_INTENSITY_GRAZING: _check_grassland_ipcc_management_category,
-    IpccManagementCategory.NOMINALLY_MANAGED: _check_grassland_ipcc_management_category,
-    IpccManagementCategory.OTHER: _check_grassland_ipcc_management_category
+    IpccManagementCategory.NOMINALLY_MANAGED: _check_grassland_ipcc_management_category
 }
 """
 Decision tree mapping IPCC management categories to corresponding check functions for grassland.
@@ -1508,9 +1387,9 @@ Value: Corresponding decision tree for IPCC management categories based on land
 """
 
 _IPCC_LAND_USE_CATEGORY_TO_DEFAULT_IPCC_MANAGEMENT_CATEGORY = {
-    IpccLandUseCategory.GRASSLAND: IpccManagementCategory.NOMINALLY_MANAGED,
-    IpccLandUseCategory.ANNUAL_CROPS_WET: IpccManagementCategory.FULL_TILLAGE,
-    IpccLandUseCategory.ANNUAL_CROPS: IpccManagementCategory.FULL_TILLAGE
+    IpccLandUseCategory.GRASSLAND: IpccManagementCategory.UNKNOWN,
+    IpccLandUseCategory.ANNUAL_CROPS_WET: IpccManagementCategory.UNKNOWN,
+    IpccLandUseCategory.ANNUAL_CROPS: IpccManagementCategory.UNKNOWN
 }
 """
 Mapping of default IPCC management categories for each IPCC land use category.
@@ -1540,7 +1419,9 @@ def _assign_ipcc_carbon_input_category(
         Assigned IPCC Carbon Input Category.
     """
     decision_tree = _DECISION_TREE_FROM_IPCC_MANAGEMENT_CATEGORY.get(ipcc_management_category, {})
-    default = _DEFAULT_CARBON_INPUT_CATEGORY.get(ipcc_management_category, IpccCarbonInputCategory.OTHER)
+    default = _DEFAULT_CARBON_INPUT_CATEGORY.get(
+        ipcc_management_category, IpccCarbonInputCategory.NOT_RELEVANT
+    )
 
     should_run_ = len(management_nodes) > 0
 
@@ -1577,7 +1458,7 @@ def _check_grassland_ipcc_carbon_input_category(
 
 _GRASSLAND_IPCC_CARBON_INPUT_CATEGORY_TO_MIN_NUM_IMPROVEMENTS = {
     IpccCarbonInputCategory.GRASSLAND_HIGH: 2,
-    IpccCarbonInputCategory.GRASSLAND_MEDIUM: 1
+    IpccCarbonInputCategory.GRASSLAND_MEDIUM: 0
 }
 """
 A mapping from IPCC Grassland Carbon Input Categories to the minimum number of improvements required.
@@ -1898,10 +1779,10 @@ Value: Decision tree for Carbon Input Categories corresponding to the management
 """
 
 _DEFAULT_CARBON_INPUT_CATEGORY = {
-    IpccManagementCategory.IMPROVED_GRASSLAND: IpccCarbonInputCategory.GRASSLAND_MEDIUM,
-    IpccManagementCategory.FULL_TILLAGE: IpccCarbonInputCategory.CROPLAND_LOW,
-    IpccManagementCategory.REDUCED_TILLAGE: IpccCarbonInputCategory.CROPLAND_LOW,
-    IpccManagementCategory.NO_TILLAGE: IpccCarbonInputCategory.CROPLAND_LOW
+    IpccManagementCategory.IMPROVED_GRASSLAND: IpccCarbonInputCategory.UNKNOWN,
+    IpccManagementCategory.FULL_TILLAGE: IpccCarbonInputCategory.UNKNOWN,
+    IpccManagementCategory.REDUCED_TILLAGE: IpccCarbonInputCategory.UNKNOWN,
+    IpccManagementCategory.NO_TILLAGE: IpccCarbonInputCategory.UNKNOWN
 }
 """
 A mapping from IPCC Management Categories to default Carbon Input Categories.
@@ -2030,8 +1911,10 @@ def _should_run_inventory_year(group: dict) -> bool:
     """
     Determines whether there is sufficient data in an inventory year to run the tier 1 model.
 
-    1. Check if the land use category is not "OTHER"
-    2. Check if all required keys are present.
+    1. Check if all required keys are present.
+    2. Check if the land use category is not "OTHER" or "UNKNOWN"
+    3. Check if the management category is not "UNKNOWN"
+    4. Check if the carbon input category is not "UNKNOWN"
 
     Parameters
     ----------
@@ -2043,7 +1926,8 @@ def _should_run_inventory_year(group: dict) -> bool:
     bool
         True if the inventory year is valid, False otherwise.
     """
-    return all([
-        group.get(_InventoryKey.LU_CATEGORY) != IpccLandUseCategory.OTHER,
-        all(key in group.keys() for key in _REQUIRED_KEYS),
+    return all(key in group.keys() for key in _REQUIRED_KEYS) and all([
+        group.get(_InventoryKey.LU_CATEGORY) not in [IpccLandUseCategory.OTHER, IpccLandUseCategory.UNKNOWN],
+        group.get(_InventoryKey.MG_CATEGORY) not in [IpccManagementCategory.UNKNOWN],
+        group.get(_InventoryKey.CI_CATEGORY) not in [IpccCarbonInputCategory.UNKNOWN]
     ])

hestia_earth/models/ipcc2019/organicCarbonPerHa_tier_2_utils.py

@@ -765,7 +765,7 @@ def _get_f_2_annual(
     Get the value of `f_2` (the stabilisation efficiencies for structural decay products entering the active pool)
     based on the tillage `IpccManagementCategory`.
 
-    If tillage regime is unknown, `IpccManagementCategory.OTHER` should be assumed.
+    If tillage regime is unknown, `IpccManagementCategory.UNKNOWN` should be assumed.
 
     Parameters
     ----------
@@ -793,7 +793,7 @@ def _get_f_2_annual(
         IpccManagementCategory.FULL_TILLAGE: f_2_full_tillage,
         IpccManagementCategory.REDUCED_TILLAGE: f_2_reduced_tillage,
         IpccManagementCategory.NO_TILLAGE: f_2_no_tillage,
-        IpccManagementCategory.OTHER: f_2_unknown_tillage
+        IpccManagementCategory.UNKNOWN: f_2_unknown_tillage
     }
     default = f_2_unknown_tillage
     return vstack([ipcc_tillage_management_category_to_f_2s.get(till, default) for till in tillage_category_annual])
@@ -1611,7 +1611,7 @@ def _get_grouped_tillage_categories(grouped_cycles):
 
 def _assign_tier_2_ipcc_tillage_management_category(
     cycles: list[dict],
-    default: IpccManagementCategory = IpccManagementCategory.OTHER
+    default: IpccManagementCategory = IpccManagementCategory.UNKNOWN
 ) -> IpccManagementCategory:
     """
     Assigns a tillage `IpccManagementCategory` to a list of HESTIA `Cycle`s.

hestia_earth/models/ipcc2019/organicCarbonPerHa_utils.py

@@ -65,6 +65,7 @@ class IpccLandUseCategory(Enum):
     FOREST = "forest"
     NATIVE = "native"
     OTHER = "other"
+    UNKNOWN = "unknown"
 
 
 SITE_TYPE_TO_IPCC_LAND_USE_CATEGORY = {
@@ -84,7 +85,8 @@ IPCC_LAND_USE_CATEGORY_TO_LAND_COVER_LOOKUP_VALUE = {
     IpccLandUseCategory.ANNUAL_CROPS: "Annual crops",
     IpccLandUseCategory.SET_ASIDE: "Set aside",
     IpccLandUseCategory.FOREST: "Forest",
-    IpccLandUseCategory.NATIVE: "Native"
+    IpccLandUseCategory.NATIVE: "Native",
+    IpccLandUseCategory.OTHER: "Other"
 }
 """
 A dictionary mapping IPCC land use categories to corresponding land cover lookup values in the
@@ -106,15 +108,15 @@ class IpccManagementCategory(Enum):
     FULL_TILLAGE = "full tillage"
     REDUCED_TILLAGE = "reduced tillage"
     NO_TILLAGE = "no tillage"
-    OTHER = "other"
+    UNKNOWN = "unknown"
+    NOT_RELEVANT = "not relevant"
 
 
 IPCC_MANAGEMENT_CATEGORY_TO_GRASSLAND_MANAGEMENT_TERM_ID = {
     IpccManagementCategory.SEVERELY_DEGRADED: "severelyDegradedPasture",
     IpccManagementCategory.IMPROVED_GRASSLAND: "improvedPasture",
     IpccManagementCategory.HIGH_INTENSITY_GRAZING: "highIntensityGrazingPasture",
-    IpccManagementCategory.NOMINALLY_MANAGED: "nominallyManagedPasture",
-    IpccManagementCategory.OTHER: "nativePasture"
+    IpccManagementCategory.NOMINALLY_MANAGED: ["nominallyManagedPasture", "nativePasture"]
 }
 """
 A dictionary mapping IPCC management categories to corresponding grassland management term IDs from the land cover
@@ -146,7 +148,8 @@ class IpccCarbonInputCategory(Enum):
     CROPLAND_HIGH_WITHOUT_MANURE = "cropland high (without manure)"
     CROPLAND_MEDIUM = "cropland medium"
     CROPLAND_LOW = "cropland low"
-    OTHER = "other"
+    UNKNOWN = "unknown"
+    NOT_RELEVANT = "not relevant"
 
 
 CarbonSource = NamedTuple(

hestia_earth/models/linkedImpactAssessment/freshwaterWithdrawalsInputsProduction.py

@@ -30,7 +30,8 @@ REQUIREMENTS = {
 }
 RETURNS = {
     "Indicator": [{
-        "value": ""
+        "value": "",
+        "inputs": ""
     }]
 }
 TERM_ID = 'freshwaterWithdrawalsInputsProduction'

hestia_earth/models/linkedImpactAssessment/landOccupationInputsProduction.py

@@ -31,6 +31,7 @@ REQUIREMENTS = {
 RETURNS = {
     "Indicator": [{
         "value": "",
+        "inputs": "",
         "landCover": ""
     }]
 }

hestia_earth/models/linkedImpactAssessment/landTransformation100YearAverageInputsProduction.py

@@ -31,6 +31,7 @@ REQUIREMENTS = {
 RETURNS = {
     "Indicator": [{
         "value": "",
+        "inputs": "",
         "landCover": "",
         "previousLandCover": ""
     }]

hestia_earth/models/linkedImpactAssessment/landTransformation20YearAverageInputsProduction.py

@@ -31,6 +31,7 @@ REQUIREMENTS = {
 RETURNS = {
     "Indicator": [{
         "value": "",
+        "inputs": "",
         "landCover": "",
         "previousLandCover": ""
     }]