Additional notes for Regional FishMIP Models [16]
The goal of the FishMIP Model Evaluation Protocol is to understand and reduce uncertainty associated with FishMIP models through model evaluation under historical climate and fishing effort forcings.
This information will allow FishMIP to better target policy initiatives such as IPCC and IPBES by providing more robust uncertainty assessment, as well as advancing the state of FishMIP models for informing vulnerability, impact, and adaptation plans of coastal sea ecosystems and fisheries (requested by the FAO). It will also help move towards a detection and attribution framework.
The focus of studies produced from this simulation round will be a significant advance in at least two ways:
-
forcing models with higher spatial resolution of global climate models and reconstructed historical fishing effort and
-
assessment of models against reconstructed historical catches (as well as biomass and other fisheries and ecological metrics in regions where additional data are available).
This protocol is designed to contribute directly to the ISIMIP 3a Simulation Round on “Evaluation, Detection and Attribution”.
Revised target date for uploading simulation outputs: Jan 30th, 2023
Uploading simulations by this date is essential to ensure enough time for analysis and writing of manuscripts in time for submission to the Special Issue. If you are able to upload results sooner, that would be helpful. If you need more time, please let us know as soon as possible.
Progress towards this deadline will be supported and facilitated through online global and regional modeller workshops to:
-
Ensure correct ESM model inputs and access
-
Ensure fishing drivers work (separate global and regional breakaway groups)
-
Tool sharing & troubleshooting
-
Check model outputs/issues
In this document we describe the general experimental and scenario set-up (Section 3). Further down in Section 4 we include the details of the specific input variables that modellers can use to implement scenarios. In Section 5 we describe the set of outputs to be created. Finally in Sections 6-7 we provide further notes and instructions on how to report and upload model results.
Further information on this protocol can be found here:
https://protocol.isimip.org/#ISIMIP3a/marine-fishery_regional/marine-fishery_global
For this simulation round, we are asking you to run and upload 2 core runs, 2 optional (but preferred) detection and attribution runs, and 2 optional (but helpful) sensitivity test runs, described below.
Each model experiment is a set of model simulations that has a particular goal (e.g. model evaluation). A scenario is a particular setting for forcing drivers that describes how each model run should be set up in the experiment, including both the type of climate forcing (CF) and the type of direct human forcing (DHF), in particular fishing.
The model runs that we are requesting for this simulation round are listed below. Please prioritize the core runs below, and provide the ‘optional’ if possible.
Table 1: Experiment set-up. Each experiment is specified by the climate forcing (CF) and Direct Human Forcing (DHF).
| Experiment | Scenario description | Scenario specifier | Sensitivity specifier |
|---|---|---|---|
Model evaluation 2 core runs |
1.Climate, river inputs, fishing – high res
|
obsclim histsoc |
default |
2.Climate, river inputs, no fishing – high res
|
obsclim nat |
default | |
Riverine influx sensitivity 2 optional (but preferred) runs |
3.Climate, no river input forcing, fishing – high res
|
obsclim histsoc |
1955-riverine-input |
4.Climate, no river input forcing, no fishing – high res
|
obsclim nat |
1955-riverine-input | |
Resolution sensitivity test 2 optional (but helpful )runs |
5.Climate, river input forcing, fishing – low res
|
obsclim histsoc |
60arcmin |
6.Climate, river input forcing, no fishing – low res
|
obsclim nat |
60arcmin | |
Control sensitivity test 2 optional runs |
7.Climate, no river input forcing, fishing – low res
|
obsclim histsoc |
1955-riverine-input+60arcmin |
8.Climate, no river input forcing, no fishing – low res
|
obsclim nat |
1955-riverine-input+60arcmin |
The focal historical period for this model evaluation experiment spans 1961-2010. To capture the transition from a pre-industrial spin-up to 1961 we also provide input for a gradual increase in fishing and environmental variability for the pre-industrial period to 1961.
For fishing effort prior to 1961, we provide input for a nominal spin-up (1841-1860, fishing held constant at 1861 levels) and pre-industrial transition period (1861-1960, reconstructed fishing effort).
To set-up climate-forcing variables for the entire 1841-1960 period, we ask modellers to use the “control run” (ctrlclim) monthly output for the years 1961-1980 (inclusive) on repeat for six cycles. These years have been selected because they correspond with an entire ENSO cycle and because no climate trend is detectable prior to 1980 from the GFDL model.
For models that require longer spin-up prior to 1841, please keep 1841 levels of fishing effort constant and, if needed, repeat the ENSO cycle (e.g. monthly values for 1961-1980 inclusive from ctrlclim) for as many times necessary.
For the ‘no fishing’ runs (nat), the spin-up and pre-industrial transition should not use any fishing effort.
We ask modellers to include all outputs from 1841 onwards for use in our evaluation assessment of model drift. Each output should be saved as two files, the first covering the spin-up and transition period (1841-1960) and the second covering the historical (experiment) period (1961-2010).
Throughout the protocol we use ‘specifiers’ that are shortened names used to denote a particular scenario, variables, or other parameter in the filenames of model inputs and outputs. It is crucial that you also use the same specifiers in your output files.
Correct formatting and naming of output files are essential for model intercomparison and analysis.
Tables 2-4 describe the different scenarios for the model runs described in Table 1. These specifiers are used in the file names of the corresponding input files and should also be used for the names of the output files (see 7. Reporting model results).
| Scenario specifier | Description |
|---|---|
| obsclim | Climate with observed atmospheric forcing and river input forcing used for model evaluation and the detection and attribution task. |
| Scenario specifier | Description |
|---|---|
| histsoc | Varying direct human influences (i.e. historical estimates of fishing effort). |
| nat | No fishing (naturalized run). |
| Scenario specifier | Description |
|---|---|
| default | 0.25 degree resolution climate model inputs. |
| 1955-riverine-input. | No temporal (land-use change) dynamics influencing river influx. |
| 60arcmin | 1 degree resolution climate model inputs. |
| 1955-riverine-input+60arcmin | No temporal (land-use change) dynamics influencing river influx and 1 degree resolution climate model inputs. |
Please remember to use these same specifiers in your output files. More on reporting data can be found at the end of this document.
For modellers new to FishMIP: to access all input data you first need to set up an account with ISIMIP to access the DKRZ server. Please follow the instructions here: https://www.isimip.org/dashboard/accessing-isimip-data-dkrz-server/
| Title | Specifiers | Time period | Reanalysis | Bias adjustment target | Comments | Priority |
|---|---|---|---|---|---|---|
| GFDL-MOM6-COBALT2 | gfdl-mom6-cobalt2_obsclim_<variable>_15arcmin | 1961-2010 | JRA55 | none | River input and JRA-55 driven MOM6-COBALT2 reanalysis. Data on native grid was remapped to a regular 0.25° grid. For further experiment design details see https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2021GL094367 | 1 |
| GFDL-MOM6-COBALT2 | gfdl-mom6-cobalt2_ctrlclim_<variable>_15arcmin | 1961-2010 | JRA55 | none | JRA-55 only driven MOM6-COBALT2 reanalysis (control run). Data on native grid was remapped to a regular 0.25° grid. For further experiment design details see https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2021GL094367 | 2 |
| GFDL-MOM6-COBALT2 | gfdl-mom6-cobalt2_obsclim_<variable>_60arcmin | 1961-2010 | JRA55 | none | River input and JRA-55 driven MOM6-COBALT2 reanalysis. Data on native grid was remapped to a regular 1° grid. For further experiment design details see https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2021GL094367 | 3 |
| GFDL-MOM6-COBALT2 | gfdl-mom6-cobalt2_ctrlclim_<variable>_60arcmin | 1961-2010 | JRA55 | none | JRA-55 only driven MOM6-COBALT2 reanalysis (control run). Data on native grid was remapped to a regular 1° grid. For further experiment design details see https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2021GL094367 | 4 |
Table 6. Climate forcing variables and units for FishMIP 3a simulations. All variables are available on a 0.25 and 1 degree horizontal grid, monthly and annual resolutions. Note: Some variables are available as specific layers extracted from vertically resolved data. Their variable names have been suffixed with -bot (ocean bottom, e.g. o2-bot), -surf (surface values, e.g. pH-surf) or -vint (vertically integrated, e.g. phyc-vint), respectively, or prefixed with int (vertically integrated, e.g. intpp). Temperature is suffixed with b or s for bottom (e.g. tob) or surface (e.g. tos) layers, respectively.
| Variable | Specifier | Unit | Resolution | Datasets |
|---|---|---|---|---|
| Mass Concentration of Total Phytoplankton Expressed as Chlorophyll | chl | kg m-3 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Floor Depth | deptho | m | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Downward Flux of Particulate Organic Carbon | expc-bot | mol m-2 s-1 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Particulate Organic Carbon Content | intpoc | kg m-2 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Primary Organic Carbon Production by All Types of Phytoplankton | intpp | mol m-2 s-1 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Net Primary Organic Carbon Production by Diatoms | intppdiat | mol m-2 s-1 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Net Primary Mole Productivity of Carbon by Diazotrophs | intppdiaz | mol m-2 s-1 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Net Primary Mole Productivity of Carbon by Picophytoplankton | intpppico | mol m-2 s-1 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Maximum Ocean Mixed Layer Thickness Defined by Sigma T | mlotst-0125 | m | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Dissolved Oxygen Concentration | o2, o2-bot o2-surf |
mol m-3 mol m-2 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| pH | ph ph-bot ph-surf |
1 1 1 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Phytoplankton Carbon Concentration | phyc phyc-vint |
mol m-3 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Mole Concentration of Diatoms expressed as Carbon in sea water | phydiat phydiat-vint |
mol m-3 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Mole Concentration of Diazotrophs Expressed as Carbon in Sea Water | phydiaz phydiaz-vint |
mol m-3 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Mole Concentration of Picophytoplankton Expressed as Carbon in Sea Water | phypico phypico-vint |
mol m-3 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Water Salinity | so so-bot so-surf |
‰ ‰ % |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Water Potential Temperature | thetao | °C | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Ocean Model Cell Thickness | thkcello | m | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Water Potential Temperature at Sea Floor | tob | °C | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Surface Temperature | tos | °C | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Water X Velocity | uo | m s-1 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Water Y Velocity | vo | m s-1 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Mole Concentration of Mesozooplankton expressed as Carbon in sea water | zmeso zmeso-vint |
mol m-3 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Mole Concentration of Microzooplankton expressed as Carbon in sea water | zmicro zmicro-vint |
mol m-3 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Zooplankton Carbon Concentration | zooc zooc-vint |
mol m-3 mol m-2 |
0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Net Downward Shortwave Radiation at Sea Water Surface | rsntds | W m-2 | 0.25° , 1° grid | GFDL-MOM6-COBALT2 |
| Sea Ice Area Fraction | siconc | % | 0.25° , 1° grid | Remapped from global JRA-55 re-analysis ice cover |
The monthly climate forcing files for obsclim (60arcmin and 15arcmin) of this simulation protocol can be found on DKRZ here:
levante:/work/bb0820/ISIMIP/ISIMIP3a/InputData/climate/ocean/obsclim/global/monthly/historical/GFDL-MOM6-COBALT2/
The monthly climate forcing files for climate forcing for spin-up and control runs (ctrlclim: onedeg and 15arcmin) of this simulation protocol can be found on DKRZ here:
levante:/work/bb0820/ISIMIP/ISIMIP3a/InputData/climate/ocean/ctrlclim/global/monthly/historical/GFDL-MOM6-COBALT2
The variables deptho and thkcello are fixed through time and can be found in the “fixed/” folder (rather than monthly/).
Production and carbon data for large and small phytoplankton can be derived from the variables in Table 1 by the following:
large = diatoms + diazotrophs
small = picophytoplankton
The GFDL model treats diazotrophs as large phytoplankton as part of their food-web processes.
For regional models, only specific grid cells will be needed from the above global outputs. Please let us know if you require assistance to extract results (e.g. using bounding boxes, masks or shapefiles). This functionality is now partially available (bounding box) through the ISIMIP web-based data portal.
A simple worked example on how to do this for specific regions in R is provided here: https://github.com/Fish-MIP/FishMIP_extracting-data
| Dataset | Included variables (short names) | Time period/ Resolution |
Reference/Source and Comments |
|---|---|---|---|
| File name | socioeconomic/fishing/histsoc/effort_histsoc_1841_2010.csv |
||
| Spatially aggregated fishing effort |
|
|
Sources: Rousseau et al., 2019, PNAS 116 (25) 12238-12243 and Rousseau et al. 2022 in prep. |
| File name | socioeconomic/fishing/histsoc/gridded_industrial_effort_histsoc_1961_2010.csv |
||
| Gridded total industrial fishing effort |
|
|
Rousseau et al. 2022 in prep. |
| File name | socioeconomic/fishing/histsoc/gridded_artisanal_effort_histsoc_1961_2010.csv |
||
| Gridded total artisanal fishing effort |
|
|
Rousseau et al. 2022 in prep. |
| Variable Name | Long name | Unit | Description/notes |
|---|---|---|---|
| Year | (End of the) year when the fishing effort is occurring | Number code | |
| Sector | The fishing sector defined by the law of the country | Name code | I = Industrial and A = artisanal, where artisanal include powered and unpowered artisanal fleets |
| LME | Large Marine Ecosystem Number | Number code | A number code of the Large Marine ecosystem in which the Effort is occurring |
| eez_country_name | Exclusive Economic Zone | Name code | The country-level exclusive economic zone (or high seas) name in which fishing effort is occurring |
| SAUP | A number code for the fishing country, following Sea Around Us numbering | Number code | Ex supranational entities (USSR, Yugoslavia) are disaggregated to their constituent countries. Serbian Fishing Effort included with Montenegro. Crimea included with Ukraine. | |
| Gear | The fishing gear | Name code | Gear names |
| FGroup | The targeted functional group | Name code | Functional groups are in accordance with those used by the Sea Around Us Project |
| NomActive | Nominal fishing effort (i.e., not including the technological creep) of the active fleet | Days at sea X kW | NomActive (of the active fleet; i.e., total) = P (engine power of active the fleet; i.e., total) x DAS (average days at sea of one vessel). Average DAS for one vessel ~ 200 DAS/year. NomActive corresponds to the total (reported, IUU, discards) catch. To find NomActive in DAS do (NomActive/P) X NV |
| Phase | Phase (year ranges) of simulation run | Text | Spin-up, transition, experiment, or validation |
The monthly fishing effort forcing files for the spin-up and experiments (Table 1) of this simulation protocol can be found on DKRZ here:
levante:/work/bb0820/ISIMIP/ISIMIP3a/InputData/socioeconomic/fishing/histsoc/
For global models, the above spatially aggregated fishing effort can be spatially allocated into 0.25 grid cells. This can be achieved using different approaches such as a simple gravity model – e.g. see Coll et al. 2020 but details will depend on model structure.
We are developing a simplified worked example for global modellers to explore and contribute to. This will be made available on github/FishMIP in due course.
While we recommend using the above spatially aggregated effort, for global models that cannot technically carry out spatial allocation of effort, gridded total industrial and artisanal nominal active effort have been provided in the same folder as the file above and are saved as netcdf files. These can be allocated to functional groups (e.g. according to relative biomass) depending on model structure.
Downscaling of the above fishing effort to match regional model inputs is likely to be needed. We request that regional modellers work together in their specific regions to ensure we have clear and common methodologies.
We are developing a worked example for regional modellers to explore and contribute to for their region which will be made available on github/FishMIP in due course.
Modellers are permitted to calibrate or tune their models using historical fisheries catch data (that will also be used for model evaluation) on the condition that only years up to and including 2004 are used in model calibration/tuning.
Modelling groups must keep detailed documentation on how their model was calibrated (e.g. input forcing, calibration data, time domain, spatial domain, fish grouping (size, functional types, total), optimization metric(s), weighting schemes, etc.) to be included in manuscript methods. Written description of sources of calibration data and methods used need to be provided with all simulation outputs. A template will be provided for this documentation in due course.
The fisheries catch data .csv file that can be used for model calibration is here:
levante:/work/bb0820/ISIMIP/ISIMIP3a/InputData/socioeconomic/fishing/histsoc/calibration_catch_histsoc_1850_2004.csv.
The fisheries catch data are already aggregated into the functional groups and spatial zones as the above effort forcing data. The original reference including links to full database is Watson & Tidd, 2018, Marine Policy, 93: 171-177.
Large marine ecosystem (LME) masks in four different spatial resolutions. 0.1°, 0.25°, 0.5° and 1° are available here:
/work/bb0820/ISIMIP/ISIMIP3a/InputData/geo_conditions/fishmip_regions/
Each region has its own variable within each file.
We have also provided conversion tables that can be used to look up LME and SAUP names according to the numeric codes used in the catch and effort files (e.g. LME 22 – North Sea). These files (SAUPnames.csv and LMEnames.csv) are also available here:
/work/bb0820/ISIMIP/ISIMIP3a/InputData/geo_conditions/fishmip_regions/
All spatially gridded outputs should be created as netcdf files. More information on how to prepare these files can be found here. Aspatial regional model results may be saved as .csv files.
In the output files, please label the time variable as “days since 1841-1-1 00:00:00” if the output covers the spin-up and transition period (1841-1960) or “days since 1901-1-1 00:00:00” if the output covers the experiment period (1961-2010).
Table 9: Mandatory output variables for Fisheries and Marine Ecosystem models (global and regional). See notes on additional optional model outputs below. Please use the value 1.e+20 for missing data within your output files. All biomasses are in wet weight (not g C).
| Variable long name | Variable specifier | Unit | Resolution | Comments |
|---|---|---|---|---|
| Total Consumer Biomass Density | tcb | g m-2 | 0.25° grid monthly |
All consumers (trophic level >1, vertebrates and invertebrates) |
| Total Consumer Biomass Density in log10 Weight Bins | tcblog10 | g m-2 | 0.25° grid monthly |
Level dimensions: (time, bins, lat, lon). If the model is size-structured, please provide biomass in equal log 10 g weight bins (1-10g, 10-100g, 100g-1kg, 1-10kg, 10-100kg, >100kg) |
| Total Pelagic Biomass Density | tpb | g m-2 | 0.25° grid monthly |
All pelagic consumers (trophic level >1, vertebrates and invertebrates) |
| Total Demersal Biomass Density | tdb | g m-2 | 0.25° grid monthly |
All demersal consumers (trophic level >1, vertebrates and invertebrates) |
| Total Catch Density (all commercial functional groups / size classes) | tc | g m-2 | 0.25° grid monthly |
Catch at sea (all catch as a result of all effort including reported and IUU) summed for both Industrial and Artisanal sector. |
| Total Industrial Catch Density (all commercial functional groups / size classes) | tic | g m-2 | 0.25° grid monthly |
Catch at sea (all catch as a result of all effort including reported and IUU) for Industrial sector only. |
| Total Catch Density in log10 Weight Bins across both sectors | tclog10 | g m-2 | 0.25° grid monthly |
Level dimensions: (time, bins, lat, lon). If the model is size-structured, please provide biomass in equal log 10 g weight bins (1-10g, 10-100g, 100g-1kg, 1-10kg, 10-100kg, >100kg) |
| Total Pelagic Density Catch across Artisanal and Industrial sectors | tpc | g m-2 | 0.25° grid monthly |
Catch at sea of all pelagic consumers (trophic level >1, vertebrates and invertebrates) |
| Total Demersal Catch Density across Artisanal and Industrial sectors | tdc | g m-2 | 0.25° grid monthly |
Catch at sea of all demersal consumers (trophic level >1, vertebrates and invertebrates) |
| Optional output from global and regional models. All biomasses are in wet weight, not g C. | ||||
| Biomass Density of Small Pelagics <30cm | bp30cm | g m-2 | 0.25° grid monthly |
If a pelagic species and L infinity is <30 cm, include in this variable |
| Biomass Density of Medium Pelagics >=30cm and <90cm | bp30to90cm | g m-2 | 0.25° grid monthly |
If a pelagic species and L infinity is >=30 cm and <90cm, include in this variable |
| Biomass Density of Large Pelagics >=90cm | bp90cm | g m-2 | 0.25° grid monthly |
If a pelagic species and L infinity is >=90cm, include in this variable |
| Biomass Density of Small Demersals <30cm | bd30cm | g m-2 | 0.25° grid monthly |
If a demersal species and L infinity is <30 cm, include in this variable |
| Biomass Density of Medium Demersals >=30cm and <90cm | bd30to90cm | g m-2 | 0.25° grid monthly |
If a demersal species and L infinity is >=30 cm and <90cm, include in this variable |
| Biomass Density of Large Demersals >=90cm | bd90cm | g m-2 | 0.25° grid monthly |
If a demersal species and L infinity is >=90cm, include in this variable |
| Catch Density of Small Pelagics <30cm | cp30cm | g m-2 | 0.25° grid monthly |
Catch at sea of pelagic species with L infinity <30 cm |
| Catch Density of Medium Pelagics >=30cm and <90cm | cp30to90cm | g m-2 | 0.25° grid monthly |
Catch at sea of pelagic species with L infinity >=30 cm and <90 cm |
| Catch Density of Large Pelagics >=90cm | cp90cm | g m-2 | 0.25° grid monthly |
Catch at sea of pelagic species with L infinity >=90 cm |
| Catch Density of Small Demersals <30cm | cd30cm | g m-2 | 0.25° grid monthly |
Catch at sea of demersal species with L infinity <30 cm |
| Catch Density of Medium Demersals >=30cm and <90cm | cd30to90cm | g m-2 | 0.25° grid monthly |
Catch at sea of demersal species with L infinity >=30 cm and <90 cm |
| Catch Density of Large Demersals >=90cm | cd90cm | g m-2 | 0.25° grid monthly |
Catch at sea of demersal species with L infinity >=90 cm |
More specific protocols for each regional model type will be developed through our monthly online regional modeller sessions. Please contact regional FishMIP coordinators for more information.
As a first step, regional modellers will need to provide shapefiles for their respective model domains for us to help with spatial extraction of the above global climate and fishing effort forcing inputs.
Region-specific climate forcing variables will be made available here:
/work/bb0820/ISIMIP/ISIMIP3a/InputData/climate/ocean/<obsclim> or <ctrlclim>/regional/
A .csv file with fishing effort extracted for regional model ecosystems
is also available in the same folder as the global fishing effort data
(../fishing/histsoc), for regional models that have provided
shapefiles.
Regional modellers may wish to make their raw unaggregated output available for more detailed analyses, including for example, a wider range of functional groups/size classes/species and ecosystem indicators. Please discuss this with FishMIP regional coordinators before uploading files.
The specification on how to submit the data, as well as further information and instructions are given on the ISIMIP website at:
https://www.isimip.org/protocol/preparing-simulation-files
It is important that you comply precisely with the formatting specified there, to facilitate the analysis of your simulation results in the ISIMIP framework. Incorrect formatting can seriously delay analyses. The ISIMIP Team will be glad to assist with the preparation of these files if necessary.
File names consist of a series of identifier, separated by underscores. Things to note:
-
Report one variable per file.
-
In filenames, use lowercase letters only.
-
Use underscore (_) to separate identifiers.
-
Variable names consist of a single word without hyphens or underscores.
-
Use hyphens (-) to separate strings within an identifier, e.g. in a model name.
-
Data model is NETCDF4_CLASSIC with minimum compression level of 5.
-
NetCDF file extension is .nc.
-
The relative time axis’ reference date is days since 1841-1-1 00:00:00 if the output covers the spin-up and transition period (1841-1960) or days since 1901-1-1 00:00:00 if the output covers the experiment period (1961-2010). We have provided .csv files to be used for the time dimension in creating NetCDF files based on the 365 days calendar. Please see time_axix_spinup.csv and time_axis_experiment.csv in this repository. The script time_axis.r was used to create these files.
Please name the files in the Fisheries and Marine Ecosystems sector according to the following pattern:
Global models
<model>_<climate-forcing>_<bias-adjustment>_<climate-scenario>_<soc-scenario>_<sens-scenario>_<variable>_<global>_<time-step>_<start-year>_<end-year>.nc
Example:
boats_gfdl-mom6_cobalt2_none_obsclim_histsoc_default_tcb_global_monthly_1961_2010.nc
Regional models
<model>_<climate-forcing>_<bias-adjustment>_<climate-scenario>_<soc-scenario>_<sens-scenario>_<variable>_<region>_<time-step>_<start-year>_<end-year>.nc
Example:
osmose_gfdl-mom6_cobalt2_none_obsclim_histsoc_default_tcb_benguela_monthly_1961_2010.nc
Please see the climate-scenario, soc-scenario, sens-scenario and variable identifiers given in the tables of this document.
Global models
The output files covering the spin-up period (1841-1960) can be saved on DKRZ here:
/work/bb0820/ISIMIP/ISIMIP3a/UploadArea/marine-fishery_global/model_name/temp2
The output files covering the experiment period (1961-2010) can be saved on DKRZ here
/work/bb0820/ISIMIP/ISIMIP3a/UploadArea/marine-fishery_global/model_name/temp
Regional models
The output files covering the spin-up period (1841-1960) can be saved on DKRZ here:
/work/bb0820/ISIMIP/ISIMIP3a/UploadArea/marine-fishery_regional/model_name/temp2
The output files covering the experiment period (1961-2010) can be saved on DKRZ here
/work/bb0820/ISIMIP/ISIMIP3a/UploadArea/marine-fishery_regional/model_name/temp
Please contact FishMIP coordinators or ISIMIP data managers directly (isimip-data@pik‐potsdam.de) if you have any questions or clarifications before submitting files or if you do not find your model’s path on DKRZ as described above.
Please contact FishMIP coordinators if you would like to participate in this simulation round but have encountered issues with any aspect of the protocol.
(For fishing): please provide all assumptions about catchability, technological creep, and model calibration.
Please provide any conversion factors that you have used to convert units.
FishMIP is entirely community-driven, and we appreciate the effort of all involved.