Skip to contents

Introduction

The herdr package performs environmental impact assessments based on the 2019 Refinement to the 2006 IPCC Guidelines. The core of the model is the Energy Balance, where Gross Energy (GE) intake is derived from the animal’s physiological requirements.


I. Net Energy for Maintenance (NEmNE_m) (Eq. 10.3)

NEmNE_m is the energy required to maintain the animal’s energy balance without gain or loss of weight.

Formula: NEm=CfiWeight0.75NE_m = Cf_i \cdot \text{Weight}^{0.75}

Variable Description Source in the Package
NEmNE_m Net energy for maintenance (MJ/day) Calculated
CfiCf_i Maintenance coefficient based on animal species/status ipcc_coefficients.csv
Weight Live weight of the individual animal (kg) weights.csv

II. Net Energy for Activity (NEaNE_a)

NEaNE_a represents energy used for grazing, walking, or other activities.

  • Cattle: NEa=CaNEmNE_a = C_a \cdot NE_m (Eq. 10.4)
  • Sheep and Goats: NEa=CaWeightNE_a = C_a \cdot \text{Weight} (Eq. 10.5)
Variable Description Source in the Package
NEaNE_a Net energy for activity (MJ/day) Calculated
CaC_a Activity coefficient (varies by feeding situation) ipcc_coefficients.csv
NEmNE_m Maintenance energy Result of Eq. 10.3
Weight Live weight of the animal (kg) weights.csv

III. Net Energy for Growth (NEgNE_g)

Energy required for tissue deposition (weight gain).

  • Cattle and Buffalo: (Eq. 10.6) NEg=22.02(BWCMW)0.75WG1.097NE_{g} = 22.02 \cdot \left( \frac{BW}{C \cdot MW} \right)^{0.75} \cdot WG^{1.097}

  • Sheep and Goats: (Eq. 10.7) NEg=WGlamb/kid(a+0.5b(BWi+BWf))365NE_g = \frac{WG_{lamb/kid} \cdot (a + 0.5b(BW_i + BW_f))}{365}

Variable Description Source in the Package
BWBW Average live body weight of the animal (kg) weights.csv
MWMW Mature body weight of an adult female (kg) weights.csv
WGWG Daily weight gain (kg/day) weights.csv
C,a,bC, a, b IPCC species-specific growth constants ipcc_coefficients.csv
BWi/BWfBW_i / BW_f Initial and final weight for the period (kg) weights.csv

IV. Net Energy for Lactation (NElNE_l)

Energy required for milk production based on volume and fat content.

  • Cattle: NEl=Milk(1.47+0.40Fat)NE_l = \text{Milk} \cdot (1.47 + 0.40 \cdot \text{Fat}) (Eq. 10.8)
  • Sheep and Goats: NEl=MilkEVmilkNE_l = \text{Milk} \cdot EV_{milk} (Eq. 10.9)
Variable Description Source in the Package
NElNE_l Net energy for lactation (MJ/day) Calculated
Milk Average daily milk yield (kg/day) livestock_definitions.csv
Fat Milk fat content percentage (%) livestock_definitions.csv
EVmilkEV_{milk} Energy required to produce 1 kg of milk (MJ/kg) ipcc_coefficients.csv

V. Net Energy for Work, Wool, and Pregnancy

  • Work (NEworkNE_{work}): NEwork=0.10NEmHoursNE_{work} = 0.10 \cdot NE_m \cdot \text{Hours}
  • Wool (NEwoolNE_{wool}): NEwool=PrwoolEVwool365NE_{wool} = \frac{\text{Pr}_{wool} \cdot EV_{wool}}{365}
  • Pregnancy (NEpNE_p): NEp=CpregnancyNEmNE_p = C_{pregnancy} \cdot NE_m
Variable Description Source in the Package
Hours Number of hours worked per day (draft animals) livestock_definitions.csv
Prwool\text{Pr}_{wool} Annual wool production (kg/year) livestock_definitions.csv
EVwoolEV_{wool} Energy per kg of wool (IPCC default 24 MJ/kg) Constant
CpregnancyC_{pregnancy} Pregnancy energy coefficient ipcc_coefficients.csv

VI. Efficiency Ratios (REMREM and REGREG)

Based on IPCC Equations 10.14 and 10.15, derived from feed digestibility (DEDE):

  • For Maintenance (REMREM): REM=[1.123(4.092103DE)+(1.126105(DE)2)(25.4DE)]REM = \left[ 1.123 - (4.092 \cdot 10^{-3} \cdot DE) + (1.126 \cdot 10^{-5} \cdot (DE)^2) - \left( \frac{25.4}{DE} \right) \right]
  • For Growth (REGREG): REG=[1.164(5.16103DE)+(1.308105(DE)2)(37.4DE)]REG = \left[ 1.164 - (5.16 \cdot 10^{-3} \cdot DE) + (1.308 \cdot 10^{-5} \cdot (DE)^2) - \left( \frac{37.4}{DE} \right) \right]
Variable Description Source in the Package
REMREM Efficiency of digestible energy for maintenance Calculated
REGREG Efficiency of digestible energy for growth Calculated
DEDE Digestibility of feed as % of gross energy feed_characteristics.csv

VII. Gross Energy (GE)

Gross Energy (GE) is the total daily energy intake required to sustain all physiological functions.

Formula: GE=(NEm+NEa+NEl+NEwork+NEpREM+NEg+NEwoolREG)DE/100GE = \frac{\left( \frac{NE_m + NE_a + NE_l + NE_{work} + NE_p}{REM} + \frac{NE_g + NE_{wool}}{REG} \right)}{DE/100}

Variable Description Source in the Package
GEGE Total daily gross energy requirement (MJ/day) Calculated
NExNE_x Various net energy requirements calculated above Calculated Results
REM/REGREM/REG Efficiency ratios Results of Eq. 10.14/10.15
DEDE Feed digestibility (%) feed_characteristics.csv

VIII. Methane Module: Enteric Fermentation (CH4CH_4)

Methane Emission Factor (EFEF) (Eq. 10.21)

EF=GE(Ym100)36555.65EF = \frac{GE \cdot \left( \frac{Y_m}{100} \right) \cdot 365}{55.65}

Variable Description Source in Package
EFEF Methane emission factor (kg CH4CH_4/head/year) Calculated
GEGE Daily gross energy intake (MJ/day) Result of Section VII
YmY_m Methane conversion factor (%) feed_characteristics.csv
55.65 Energy content of methane (MJ/kg CH4CH_4) Constant

Total Enteric Emissions (ETE_T) (Eq. 10.19)

ET=(P)EF(T,P)(N(T,P)106)E_T = \sum_{(P)} EF_{(T,P)} \cdot \left( \frac{N_{(T,P)}}{10^6} \right)

Variable Description Source in Package
ETE_T Total emissions from category TT (Gg CH4CH_4/year) Calculated
N(T,P)N_{(T,P)} Number of animals in category TT livestock_census.csv

IX. Methane Module: Manure Management (CH4(mm)CH_{4(mm)})

Volatile Solid Excretion Rates (VSVS) (Eq. 10.24)

VS=[GE(1DE100)+(UEGE)][1ASH18.45]VS = \left[ GE \cdot \left( 1 - \frac{DE}{100} \right) + (UE \cdot GE) \right] \cdot \left[ \frac{1 - ASH}{18.45} \right]

Variable Description Source in Package
VSVS Daily volatile solid excretion (kg VS/day) Calculated
UEUE Urinary energy excretion fraction (default 0.04) ipcc_coefficients.csv
ASHASH Ash content of the feed (fraction) feed_characteristics.csv
18.45 Energy content of dietary dry matter (MJ/kg) Constant

Manure Emission Factor (EFEF) (Eq. 10.23)

EF(T)=(VS(T)365)[B0(T)0.67S,kMCFS,k100AWMS(T,S,k)]EF_{(T)} = (VS_{(T)} \cdot 365) \cdot \left[ B_{0(T)} \cdot 0.67 \cdot \sum_{S,k} \frac{MCF_{S,k}}{100} \cdot AWMS_{(T,S,k)} \right]

Variable Description Source in Package
B0(T)B_{0(T)} Maximum methane producing capacity (m3/kgVSm^3/kg VS) ipcc_coefficients.csv
MCFS,kMCF_{S,k} Methane conversion factor for system SS (%) ipcc_mm.csv
AWMSAWMS Fraction of manure handled in system SS manure_management.csv
0.67 Conversion factor (m3CH4m^3 CH_4 to kg CH4CH_4) Constant

X. N₂O Module: Nitrogen-Based Emissions

I. Nitrogen Balance & Excretion

1. Nitrogen Intake (NintakeN_{intake}) (Eq. 10.31/10.32)

Nintake=(GE18.45)(CP/1006.25)N_{intake} = \left( \frac{GE}{18.45} \right) \cdot \left( \frac{CP/100}{6.25} \right)

Variable Description Source in Package
NintakeN_{intake} Total nitrogen ingested (kg N/day) Calculated
CPCP Crude Protein content of the diet (%) feed_characteristics.csv
6.25 Conversion factor from N to Protein Constant

2. Nitrogen Retention (NretentionN_{retention})

Reflects the logic used in the herdr codebase: * Cattle: Milkprotein=1.9+0.4FatcontentMilk_{protein} = 1.9 + 0.4 \cdot Fat_{content}Nretention=[MilkMilkprotein6.38]+[WG(268(7.03NEgWG))/10006.25]N_{retention} = \left[ \frac{Milk \cdot Milk_{protein}}{6.38} \right] + \left[ \frac{WG \cdot (268 - (7.03 \cdot \frac{NE_g}{WG})) / 1000}{6.25} \right]

Variable Description Source in Package
FatcontentFat_{content} Milk fat percentage (%) livestock_definitions.csv
MilkMilk Daily milk yield (kg/day) livestock_definitions.csv
WGWG Daily weight gain (kg/day) weights.csv
6.38 / 6.25 Specific N-to-Protein conversion factors Constants

3. Annual Excretion (NexcretedN_{excreted})

  • Sheep/Goats: Nexcreted=(Nintake(10.1))365N_{excreted} = (N_{intake} \cdot (1 - 0.1)) \cdot 365
  • Cattle: Nexcreted=(NintakeNretention)365N_{excreted} = (N_{intake} - N_{retention}) \cdot 365

II. Direct N2ON_2O Emissions (Eq. 10.25)

N2Odirect=N(T,P)NexcretedAllocationEF3(4428)N_2O_{direct} = N_{(T,P)} \cdot N_{excreted} \cdot Allocation \cdot EF_3 \cdot \left( \frac{44}{28} \right)

Variable Description Source in Package
AllocationAllocation Fraction of N managed in a specific system manure_management.csv
EF3EF_3 Direct N2ON_2O emission factor for the system ipcc_mm.csv
44/28 Stoichiometric ratio of N2ON_2O to N2N_2 Constant

XI. N₂O Module: Indirect Emissions

I. Nitrogen Losses from Manure

1. Volatilization (NvolatilizationMMSN_{volatilization-MMS}) (Eq. 10.26)

NvolMMS=[(NNexAWMS)FracGasMS]N_{vol-MMS} = \sum \left[ (N \cdot Nex \cdot AWMS) \cdot Frac_{GasMS} \right]

2. Leaching and Runoff (NleachingMMSN_{leaching-MMS}) (Eq. 10.27)

NleachMMS=[(NNexAWMS)FracLeachMS]N_{leach-MMS} = \sum \left[ (N \cdot Nex \cdot AWMS) \cdot Frac_{LeachMS} \right]

Variable Description Source in Package
FracGasMSFrac_{GasMS} Fraction of N lost as volatilized gas (NH3+NOxNH_3 + NO_x) ipcc_mm.csv
FracLeachMSFrac_{LeachMS} Fraction of N lost via leaching into water ipcc_mm.csv

II. Indirect N2ON_2O Calculation (Eq. 10.28/10.29)

  • From Volatilization: N2OG(mm)=(NvolMMSEF4)4428N_2O_{G(mm)} = (N_{vol-MMS} \cdot EF_4) \cdot \frac{44}{28}
  • From Leaching: N2OL(mm)=(NleachMMSEF5)4428N_2O_{L(mm)} = (N_{leach-MMS} \cdot EF_5) \cdot \frac{44}{28}
Variable Description Source in Package
EF4EF_4 Indirect N2ON_2O factor for volatilization ipcc_coefficients.csv
EF5EF_5 Indirect N2ON_2O factor for leaching ipcc_coefficients.csv

📚 References

  • IPCC (2019). Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Volume 4: AFOLU, Chapter 10.