2.1 System and spatial boundaries
We assessed the cradle-to-gate water and climate change impacts of ceasing dairy production in NSW. Due to concerns about the environmental and animal welfare issues around beef production, we assumed dairy and its co-products were replaced with plant-based functional equivalents (Table 1). Consequences were assessed on a regional basis for three regions; the North Coast, South Coast and Inland (Figure 3). These regions were used because they represent distinct agro-climatic zones. The North Coast of NSW has a sub-tropical climate while the South Coast is temperate and the Inland region of NSW has a Mediterranean climate. In addition, the Inland region of NSW is covered by the one river catchment (i.e. the Murray Darling Basin (MDB)). Where necessary, the system boundaries were expanded to include global production (as described in section 4.2)
Table 1
Parameters, their variation from the baseline value and the change in climate change impacts.
|
|
GHG emissions (Mt CO2-e)
|
Parameter
|
Variation (%)
|
Parameter increase
|
Parameter decrease
|
Proportion of land arable
|
50
|
-0.763018813
|
-0.552734866
|
Proportion of raw milk as cream
|
25
|
-0.600267841
|
-0.715485837
|
Soybean yield
|
25
|
-0.785978998
|
-0.52977468
|
C sequestration in vegetation
|
50
|
-0.683936136
|
-0.631817543
|
Water availability
|
50
|
-0.668027684
|
-0.647725995
|
Water requirements for irrigated soy production
|
25
|
-0.653542705
|
-0.665100396
|
2.2 Production and resource use changes
Dairy production occurs on both arable (i.e. cropland) and non-arable land. Cropland vacated by dairy production was first assumed to be used grow soybeans to replace dairy products with soymilk after which soybeans for tofu production were grown. Growing soybean crops in succession is considered poor agronomics due to the build-up of disease and weeds that reduces crop yields so we assumed that a soybeans were grown in rotation with wheat. The timing of rotations differed between each region due to the North Coast having higher rainfall and the rotations are described in Table 2. Water used for dairy production was segregated into irrigation water and water used for livestock. This was done because in NSW extractions for irrigation are licenced differently to water extracted for livestock. Water used for irrigation on dairy farms was assumed to be redirected to the production of soymilk and then to irrigating soybean crops. Water used for livestock was assumed to be made available to other river users under relevant licencing provisions. Where soybean production in a region could not meet the mass required to replace all products and co-products, it was assumed that demand for global soybean production increased. Conversely, where soybean production exceeded demand for soymilk and tofu, the additional production was assumed to avoid global soybean production. The wheat produced in the soybean-wheat rotation was assumed to reduce demand for global wheat production. Dairy cattle are fed supplements that consist primarily of cereal and plant-based protein meal (e.g. canola meal) and it was assumed that the absence of feeding supplements avoided global production of wheat and protein meal. Cream is separated from the milk during processing. We assumed that cream was either used as unprocessed cream or processed to make butter. The reduction in cream production was assumed to increase demand for global vegetable oil, as a part replacement of fresh cream and to make margarine. The reduction in meat meal and offal was assumed to increase global demand for protein meal. A co-product of soymilk and tofu production is okara. Okara has limited uses19 so was assumed to be applied to wheat paddocks avoiding the need for N fertilisers on a 1:1 N mass basis.
Table 2
The function, functional equivalent and substitution basis for each product and co-product from NSW dairy production systems
Product/co-product
|
Function
|
Functional equivalent
|
Substitution basis
|
Fresh milk
|
Human consumption
|
Soy milk
|
Volume
|
Fresh cream
|
Human consumption
|
Soy milk/vegetable oil (2:1 ratio)
|
Volume
|
Cream in butter
|
Human consumption
|
Vegetable oil
|
Volume
|
Red meat
|
Human consumption
|
Tofu
|
Mass
|
Hide
|
Leather
|
PVC
|
Mass
|
Tallow
|
Multiple
|
Vegetable oil
|
Mass
|
Offal
|
Pet meat
|
Protein meal
|
Protein mass
|
Meat meal
|
Chicken/pork feed
|
Protein meal
|
Protein mass
|
Table 3
Descriptions of irrigated and dryland soy rotations for the North Coast, South Coast and Inland regions used in this study.
Region
|
Irrigated soy rotation
|
Dryland soy rotation
|
North Coast
|
Irrigated soybeans and unirrigated wheat annually.
|
Unirrigated soybeans followed by unirrigated wheat over a two-year period.
|
South Coast
|
Irrigated soybeans and unirrigated wheat over a two year period.
|
Unirrigated soybeans followed by unirrigated wheat over a two year period.
|
Inland
|
Irrigated soybeans and unirrigated wheat over a two year period.
|
Unirrigated soybeans followed by unirrigated wheat over a two year period.
|
Non-arable land is primarily used for growing out heifers to join the milking herd or to dry off cows that are due to calve. Under current market conditions this land would be used to graze beef cattle because the soils are too poor the land is too steep to crop. Our objective of replacing dairy and dairy co-products only with plant-based alternatives meant that the land could not be used to produce plant-based functional equivalents. We therefore assumed that non-arable land was allowed to revegetate and sequester atmospheric C. C sequestration was estimated for a central point in each region using the Full Carbon Accounting model following the guidelines20 for the Human Induced Regeneration methodology developed by the Australian Government to issue carbon credits under the Emissions Reduction Fund21. C sequestration in regenerating areas is minimal during the initial years so the average C sequestration over the first 10 years of the simulation was used for each region. Taking existing agricultural land was assumed to increase demand for marginal global agricultural land and this was accounted for using the iLUC framework of Schmidt and Weidema22.
Regional experts suggest that for many farms in the South Coast region, the most reasonable alternative land use if dairy production were to cease would be peri-urban development. Transforming agricultural land to peri-urban development would result in no meaningful agricultural production and would also increase water consumption because more properties would have access to stock and domestic water rights. This option was not included in the study due to a lack of available information on which to base assumptions.
2.3 Sensitivity analysis
Sensitivity analysis was used to test how sensitive the results of the study were to changes in key parameters (Table 1).
2.4 Mitigation strategies
We also tested what the consequences of ceasing NSW dairy production would be if two climate change mitigation strategies, feeding enteric methane inhibitors and capturing methane from manure collected at the dairy, were implemented across the NSW dairy industry as an theoretical alternative baseline. Roque, Salwen 23 showed that dairy cows fed 1% Asparagopsis armata reduced enteric methane production, dry matter intake and milk production per cow by 67%, 38% and 11%, respectively. To test the impact of enteric methane reduction inhibitors on results we assumed that enteric methane emissions and milk production per cow declined by those amounts. The reduction in total dry matter intake per cow would mean that more cows could be milked using the existing resources, so we assumed that the number of cows in the milking herd increased by 50%. This 50% increase is lower than achievable based on changes to dry matter intake alone but increasing the numbers of milking cows will be limited to some extent by infrastructure (e.g., shade, laneways) and limiting the increase to 50% accounts for these limitations. For the methane capture mitigation strategy, it was assumed that effluent ponds were covered, and methane was flared converting it to CO2.
2.5 Life cycle inventory
Inventory for dairy, soybean and wheat production were developed by modifying relevant inventory from AusLCI 24 with values that represented each region and included all farm operations and inputs. Global production of soybean, vegetable oil and protein meal were represented by the appropriate inventory in the ecoinvent v3.5 consequential database 25. Inventory used in the study can be obtained in SimaPro format by contacting the corresponding author of this publication.
2.6 Input data
A list of input data and their source used to develop inventory is available in Table 4. Regional experts were used to source or validate input data and consisted of dairy advisors and processors in each relevant region.
Table 4
Input data for the South Coast, North Coast and Inland regions, and the source of the data, used in the study.
|
Region
|
|
Parameter
|
South Coast
|
North Coast
|
Inland
|
Source
|
Milk production (L)
|
540720000
|
275533000
|
258630000
|
Dairy Australia 46
|
Milk/cow/day (L)
|
18.6
|
17.6
|
20.9
|
Dairy Australia 47
|
Area/farm (ha)
|
355
|
321
|
418
|
Dairy Australia 47
|
Cows/farm (no)
|
418
|
314
|
426
|
Dairy Australia 47
|
Proportion of farm arable (%)
|
80
|
70
|
70
|
Expert opinion
|
Cull age of dairy cows (years)
|
7
|
7
|
7
|
Expert opinion
|
Cow liveweight at culling (kg)
|
566
|
544
|
553
|
Dairy Australia 47
|
Calf liveweight at sale (kg)
|
100
|
100
|
100
|
Expert opinion
|
Irrigated soy yield (kg/ha)
|
3000
|
4000
|
3000
|
Expert opinion
|
Dryland soy yield (kg/ha)
|
2000
|
2500
|
1500
|
Expert opinion
|
Dryland wheat yield (kg/ha)
|
3000
|
3000
|
2500
|
Expert opinion
|
Cream in raw milk (%)
|
5
|
5
|
5
|
Dairy Australia 48
|
Proportion of cream sold fresh (%)
|
50
|
50
|
50
|
Expert opinion
|
Daily dry matter intake as supplement (%)
|
40
|
40
|
40
|
Dairy Australia 47
|
Wheat in supplement (%)
|
80
|
80
|
80
|
Expert opinion
|
Canola meal in supplement (%)
|
20
|
20
|
20
|
Expert opinion
|
Average annual lime application (kg/ha)
|
500
|
500
|
500
|
Expert opinion
|
N applied to irrigated land (kg N/ha
|
201
|
213
|
157
|
Dairy Australia 47
|
N applied to non-irrigated pastures (kg N/ha)
|
113
|
133
|
12
|
Dairy Australia 47
|
Irrigation water available (ML/region))
|
12660
|
14641
|
5417
|
Dairy Australia 47
|
Irrigation to grow one ha of soy (ML)
|
6
|
4
|
6
|
Expert opinion
|
Water for livestock (L/head/day)+
|
160
|
160
|
160
|
Dairy Australia 49
|
C in vegetation (t/ha)
|
0.42
|
0.44
|
0.28
|
Australian Government 50
|
Kangaroo methane (L/head/day)
|
3.05
|
3.05
|
3.05
|
Vendl, Clauss 51
|
Kangaroo density (head/ha)
|
2
|
2
|
0.5
|
Expert opinion
|
Dressing percentage (%)
|
50
|
50
|
50
|
Expert opinion
|
Butcher waste (%)
|
32
|
32
|
32
|
Expert opinion
|
Proportion offal (%)
|
10
|
10
|
10
|
Expert opinion
|
Area of hide (m2)
|
4.5
|
4.5
|
4.5
|
Expert opinion
|
Proportion tallow (%)
|
18
|
18
|
18
|
Expert opinion
|
Proportion meat meal (%)
|
30
|
30
|
30
|
Expert opinion
|
Oil in soy (%)
|
20
|
20
|
20
|
Wikipedia
|
soybeans in soymilk (kg/l)
|
0.1
|
0.1
|
0.1
|
Grant and Hicks 52
|
soybeans in tofu (kg/kg)
|
0.4
|
0.4
|
0.4
|
Mejia, Harwatt 53
|
protein content of tofu (%)
|
12
|
12
|
12
|
Supermarket tofu packaging
|
protein content of red meat (%)
|
23
|
23
|
23
|
Wikipedia
|
protein content of meat meal (%)
|
50
|
50
|
50
|
Expert opinion
|
protein content of soy meal (%)
|
45
|
45
|
45
|
Expert opinion
|
+includes water consumed by heifers |
2.7 GHG emissions calculations
GHG emissions associated with fertiliser use on farm, crop residue breakdown and dissolution of lime were calculated using the relevant methods from the Australian National Inventory report 26. Where GHG emissions were dependent on other biophysical processes (e.g., leaching of nitrates), calculations from the Australian National Inventory report were also used. Although it is highly likely that soil organic carbon (SOC) would lost when converting grassland to cropland, the impacts of crop management practices on SOC are highly uncertain 27 so we assumed that there was no change in SOC. Emissions calculations in background inventory that represented global demand of crops were not modified.
2.8 Impact assessment
Impact assessment was done using SimaPro v8.3.028 with the AusLCI indicator set 29. AR6 GWP100 of 273, 29.8 and 27.2 for N2O, CH4 (fossil origin) and CH4 (non-fossil origin), respectively, were used in the study. Blue water stress impacts were assessed using the AWARE methodology5 and, because the North Coast and South Coast regions include multiple catchments, characterization factors for these regions were the average of all catchments in each region. This gave an AWARE characterization factor of 1.5, 1.2 and 92 for the North Coast, South Coast and Inland regions, respectively.