RainfallNZ / CASM_Horizons

CASM data processing specific to Horizons
1 stars 0 forks source link

Additional scenarios #3

Open RainfallNZ opened 4 years ago

RainfallNZ commented 4 years ago

Additional scenarios have been requested. See email from Ton to Tim Kerr "FW: PC2 - additional scenarios" Monday 29th June at 7:25 and the attachment "20200626_Scenarios.xlsx" The scenarios are described in the spreadsheet which was prepared by Abby at Environment Southland. Additional details are included in the email exchange between Abby and Ton from 6:42pm June 29th. And then again in an email from Ton to Tim Kerr 29th June at 11:06 am) For the "Point source" scenario, Tim Cox may be able to simply update his model directly with new point source data from

Operative plan (Horizon's Scenario 4)

What if all properties met the operative plan's requirements? For non-intensive properties use MPI Cox-corrected losses. For all intensive land use the Operative Table 14.2 Year 20 values.

Operative plan with Dairy conversions (Horizon's Scenario 5)

What if, within each water management sub zone 11 % expansion of Dairy (from sheep and beef land) and the operative Table 14.2 Year 20 was adhered to? As for Scenario 4 except with landuse changes.

Overseer 2012-2013 rates (Horizon's Scenario 6)

What if loss rates continued as they were in 2012-13? For non intensive properties use MPICoxCorrected. For intensive properties use the properties 2012-13 Overseer-assessed loss rates (if it is available) otherwise use the MPI Cox-corrected.

Natural state scenario (Horizon's Scenario 7)

What if there were no anthropogenic nutrient sources? Implemented by converting all landuse to "Non-productive land (inactive)" and removing all point sources.
The native land use leach rates are taken from the MPI leach rate table and do not account for the variations implemented by Tim Cox to help the initial nutrient budget to balance. The Cox-variations on leach rate were sub-management-class based, not land use based. This means implementing the Cox-variations would lead to spatial variations in leach rates that cannot be attributed to spatial variations in soil, climate or slope. The imprecision of the native scenario (through not considering e.g. change in soil depths, river flows and wetlands) makes this scenario highly uncertain and only really indicative.

Federated Farmers Scenario (Horizon's Scenario 8a and b)

What if non-compliant farms were able to reduce their nutrient losses to 90 % of their 2012-23 values, or to the 75th percentile of all the 2012-13 values for the water management zone?
This is an attempt to model an approach suggested by Federated Farmers see pages 33 and 50 of their submission. Further clarification of the desired scenario was provided by email from Craig Depree (Dairy NZ) and Abbie Matthews (Horizons) dated 28th August with subject "dairy scenario" and forwarded to Ton Snelder, then to Tim Kerr and Tim Cox on the same day. See the very end of the email for the scenario specification.
Non-intensive land is left at the MPI-Cox-Corrected values. Consented land continue with their consented values. Non-consented intensive land use has its N-loss reduced to the greater of:

To be clear, the Water Management Zone baseline data, not the Water Management SUB-zone data has been used to determine the 75th percentiles. Properties without Overseer data are set to PC2 Table 13.2 Year 20.

For variant (b), as per above but in addition:
determine the average leach rate in each water management zone for the unconsented properties that have "Overseer" leach rate estimates. Apply this to all the unconsented dairy areas in a water management zone without "Overseer" estimates.
This generalises the suggested discretionary pathway, at least for the water management zones where "Overseer" values exist.

Potato Scenario (Horizon's Scenario 9)

What if potatoes were grown wherever they could be?
Ton is investigating this with help from Potato NZ (Nic Conland and Chris Keenan). See email from Ton to Chris and Nic titled "representation of more potato growing in scenario modelling" of 1st July 2020 at 12:57 pm, cc'd to Tim Kerr All LUC 1:3 land use is converted to potato growing on a 1 in 5 year crop rotation. The potato N-loss rates determined by Bloomer et al, 2020 report (Section 10) at 46 kg/ha/year. Implementation of this is achieved in two different ways:

  1. adjusting the leach rate at all LUC 1-3 locations to equal:
    46 0.2 + Existing 0.8
    That is, implement a time-weighted leach rate, where 1/5th of the time there are potato leach rates, and 4/5ths of the time it is the existing leach rate.

  2. adjusting the leach rate at all Dairy LUC 1-3 locations as in 1 above.

On 22 nd July, Nic Conland sent Ton a power point document with some alternative potatoe leach rates of 24, 13 and 9 kg/ha for "Bloomer" "good management practice" and "best management practice". The "Bloomer" rate does not match that in the Bloomer et al (2020) report so clarification is still required about these numbers.

_Bloomer, D., O'Brien, G., Posthuma, L. 2020. Modelled Loss Of Nutrients From Vegetable Growing Scenarios In Horowhenua - OverseerFM 6.3.2 Modelling of Representative Rotations within Hokio_1a, Hokio_1b, and Ohau1b WMSZ. Report by Page Bloomer Associates for Pattle Delamore Partners. Hastings.

Consents limited to PC2 Table 14.2 Year 20 (Horizon's Scenario 10)

Non-intensive land is at the MPICoxCorrected values Consented intensive land is at the lesser of the consented rate or the Plan Change 2 Table 14.2 Year 20 rate All other intensive land is at the Plan Change 2 Table 14.2 Year 20 rate

Point Source Scenario (bonus scenario offered by Ton)

Replace the point source data (which is from 2012) with 2017 data. Ton has provided this in a spreadsheet "pointSourcesTN_2012And2017.xlsx" (also in the projects sharepoint"Data" directory) that Ton provided (see email Ton to Tim Kerr) . Tim Cox may be able to implement this directly without any preprocessing.

Consent Impact Scenario (bonus last minute scenario)

What impact has consenting had on water quality?
This scenario explores how the water quality has changed by requiring consenting. Default to Cox-calibrated leach rates everywhere. For consented farms, set the leach rate to the consented value. For the unconsented farms with Overseer estimates from 2012-13, use those rates.
This scenario is intended to be compared with Scenario 6. The hypothesis is that there is an improvement in water quality resulting from consenting.

Horticulture at alternative loss rates (Scenario 13 A and B)

Original version, later called the A variant: This scenario adjusts the horticulture leach rate in the south west climate region from 71 kg/ha/year to 46. 46 is the lower rate from the "Bloomer" report for vegetables.
B Variant: As above but change to ?? kg/ha/year. This is derived from an WSP modelling referenced in section 3 of the Resumed Joint Witness Statement of the Farming and Growing Management Practices from 12th August 2020. This found that under good management practices, mixed crop vegetable rotations led to N-loss reduction of 36 to 94 %. The 36 % reduction was used as it is the more conservative option.
These scenarios were implemented by reclassifying the Cox-calibrated Leachrate raster. The leach rate of 71 kg/ha/year is unique to the SW climate region for vegies. After converting all 71's to the new vegie loss rate and the normal processing was applied. The re-calssifying was done using a function called "LeachRateReclassifier"

Proposed year 20 (pathways A and D) -Scenario 14

This scenario was requested by Christine Foster, Planing Consultant to Horizons. See email from Christine to Abby and Ton of 14th August 2020 subject "PC2 Water Quality Modelling Scenarios" and forwarded on to Tim Kerr and Tim Cox on Saturday 15th August 2020.
Consented Dairy plus all other dairy reduced by 9 kg.ha.year (down to the PC2Table 14.2 year 20 values) plus horticulture at some reduced rates assuming implementation of good management practices. Christine initially assumed the vegetable N-loss reductionw were t be agreed within the resumed joint witness statement of Experts Farming and Growing Management Practices 12th August 2020. This joint witness statement has been provided directly from Christine (see email from Christine to Ton and Abby 17th August 2020, subject "PC2 - Resumed Farm Systems Modelling". (Note that the resumed joint witness statement will eventually be made available from the Horizon's website).
As discussed in Christine's 17th August email, the joint witness statement does NOT agree on commercial vegetable growing N-loss reductions. The nutrient losses for commercial vegetable growing are discussed under issue 3 of the joint witness statement: While no agreement was reached on specified Nitrogen reductions, the statement describes WSP modeling N-loss reduction for commercial vegetable growing. This found that from a 2014 baseline without any good management practice implementation, the N-loss may be reduced by:

This scenario is to include a variant, where the unconsented dairy farms are reduced by 18 kg/ha/year. The unconsented dairy farms included in the Horizons spreadsheet are to be included explicitly.

These scenarios are to be compared with the:

  1. "Operative" scenario (scenario 4)'
  2. "Post regulation" (scenario 12)'
  3. "Pre-regulation" (scenario 6).
RainfallNZ commented 4 years ago

Point Source Scenario details.

Implementing the Point Source Scenario required using new point source data provided by Ton. These data included a new point source site for the 2017 scenario: "Winstone Pulp WWTP at oxpond waste", which is associated with nzsegment 7196647. This segment has not been included in the CASM network before. This is the first scenario that requires a change in the CASM network. This scenario has a changed tributary network, and changed point source sites. There is also a point source site that has been removed: "NZ Pharmaceuticals wastewater". nzsegment 7240047. The removal of this site has no effect on the CASM network as it is not at the end of a tributary.

RainfallNZ commented 4 years ago

Native Land use scenario details.

This was done by loading a modified Land use spatial file where all the anthropogenic land uses had been changed to native.
This was used (in combination with the MPI leach rate look up table) to create a Native Leach Rate raster.
A SubZoneLanduseLUCShapeFile also needed to be prepared, and then aggregated on the unique combinations of water management sub zone, land use and LUC class. From there, the CASM input file was generated without any other interference.

RainfallNZ commented 4 years ago

Dairy Expansion details

Note: Initial implementation was not correct. Issues were corrected with the update on 22nd July 2020.
Dairy expansion was handled by:

  1. Find how much dairy land there was in each water management sub zone
  2. Calculate 11 % of the WMSZ's dairy area total
  3. Assume that conversions to dairy will occur from sheep and beef land prioritising the lower LUC land areas.
  4. Within each WMSZ the total land area in sheep and beef for LUC 1 was found.
  5. If there was enough area available, it was reduced by the required 11 % of the total WMSZ dairy area.
  6. If an equivalent Dairy-LUC node existed, its area was increased by the same amount.
  7. If there was not enough sheep and beef area available in the current LUC class, all the land was converted, and the process (steps 4 to 6 above) was repeated for the next higher LUC class. This continued until enough land had been converted.

In some water management sub zones that were to have dairy expansion, there was not already a Dairy-LUC diffuse node that could be increased in size to match the sheep and beef-LUC diffuse node that was to be reduced in size.
Originally this resulted in the sheep and beef area reducing, but no dairy expansion.
The solution was to create an additional Dairy node with the same attributes as the sheep and beef node being converted, and ammend the area, landuse, location and export coefficient attributes. Area was set to match the amount that the sheep and beef was reduced, land use was changed to "Dairy". location was increased by 0.005
Determining the export coefficient is not straight forward as it is not known where the dairy conversions are to occur other than somewhere in the relevant WMSZ -sheep and beef - LUC areas.
The export coefficient was set by:

  1. finding which grid cells related to the relevant WMSZ- sheep and beef - LUC combination
  2. finding the MPI predictor values for those grid cells
  3. changing the land use predictor from sheep and beef to dairy
  4. then looking up the MPI leach rate for each grid cell (with land use set to dairy)
  5. averaging the leach rate for all the cells.
RainfallNZ commented 4 years ago

Dairy Sector / Federated Farmers Scenario 8 implementation

Scenario 8a was created by analysisng the Horizons master dairy farm data spreadsheet.
This was done in the "Some-scenario-data-preparation" code chunk with the "Scenario.8V2a.Fed.Farmers" option.
For Scenario 8b things got complicated. Firstly, the leach rate was adjusted to the PC2Year20 version on intensive land use. Then the "Some-scenario-data-preparation" code chunk with the "Scenario.8V2b.Fed.Farmers" option was run to generate a "CombinedClasses" data frame. This has all the diffuse nodes with unconsented "Overseer" properties, and their areas and their leach rates and loads. Using this table, the leach rate raster and the ZoneLandUseLUC raster a new leachrate raster was created. This new raster increased the leachrate in any dairy/WMSZ combinations that had higher values determined from the non-consented Overseer rates. This was all done in the function "DairyLeachToWMZRates" function.
Once the leach rate was updated, the same process as used for scenario 8a was applied to update the leach rates for the consented and the unconsented "Overseer" properties.
Any water management zones without non-consented "Overseer" farms were not adjusted, other than to meet PC2 year 20 leach rates.