Additional notes for Regional FishMIP Models $$16$$
Reporting model results $$16$$
The goal of the FishMIP "Ocean Futures" Protocol is to extend 3b CMIP climate projections to include:
-
exploring fishing impacts in addition to those of climate, including future projections with fishing under a range of socio-economic scenarios (aligned with SSPs, Maury et al. IN PREP);
-
comparison to a no-climate change preindustrial control baseline (which was not available for the last round of simulations), and
-
if possible, building on the 3a model evaluation simulation round, model benchmarking against observed catches.
Note that this FishMIP Phase 2 protocol builds on the Phase 1 2020 protocol which focussed primarily on results without fishing for inclusion in the IPCC 5th Assessment.
Modellers are welcome to redo all previous 3b runs with version of their models that have had improvmets following the 3a model evaluation round.
Timelines for simulations
| Scenarios | Models | Date |
|---|---|---|
| histsoc,soc_2015 | global, regional | July 1st 2024 |
| soc_OSP1, soc_OSP2 | global, regional | Oct 1st, 2024 |
To aid with progress we will hold specific technical workshops to:
-
Ensure correct OSP integration 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/#ISIMIP3b/marine-fishery_regional/marine-fishery_global
For this simulation round, we are asking you to run XXXXXXXX
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).
Below we summarise the simulation experiments for both Phase 1 and Phase 2 of the 3b simulation round. Modellers that have already completed Phase 1 can skip to Phase 2. 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 | Short description | Historical | Future |
+===+==============================+============================================================================+=============+===========+
| 1 | Pre-industrial control | Climate: No climate change, fixed 1850s CO~2~ levels | picontrol | picontrol |
| | | | | |
| | nat | F ishing: No fishing | nat | nat |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
| 2 | Pr e-i ndustrial control | C limate: No climate change, fixed 1850s CO~2~ levels | picontrol | picontrol |
| | | | | |
| | histsoc | F ishing: H istorical fishing effort, then future fixed at 2015 levels | histsoc | 2015soc |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
| 3 | S SP1 -RCP2.6 | C limate: Simulated historical climate, then SS P1-RCP2.6 climate | h istorical | ssp126 |
| | | | | |
| | nat | Fis hing:No fishing | nat | nat |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
| 4 | SSP1 -RCP2.6 | Climate: Simulated h istorical climate, then SS P1-RCP2.6 climate | h istorical | ssp126 |
| | | | | |
| | histsoc | F ishing: H istorical fishing effort, then future fixed at 2015 levels | histsoc | 2015soc |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
| 5 | S SP1 -RCP2.6 | C limate: Simulated h istorical climate, then SS P1-RCP2.6 climate | h istorical | ssp126 |
| | | | | |
| | OSP1 | F ishing: H istorical fishing effort, then change driven by OSP1 | histsoc | OSP1soc |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
| 6 | S SP5 -RCP8.5 | C limate: Simulated h istorical climate, then SS P5-RCP2.6 climate | h istorical | ssp585 |
| | | | | |
| | nat | F ishing: No fishing | nat | nat |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
| 7 | S SP5 -RCP8.5 | C limate: Simulated h istorical climate, then SS P5-RCP8.5 climate | h istorical | ssp585 |
| | | | | |
| | histsoc | F ishing: H istorical fishing effort, then held fixed at 2015 levels | histsoc | 2015soc |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
| 8 | S SP5 -RCP8.5 | C limate: Simulated h istorical climate, then SS P5-RCP8.5 climate | h istorical | ssp585 |
| | | | | |
| | OSP2soc? | F ishing: H istorical fishing effort, then change driven by OSP2 | histsoc | OSP2soc |
+---+------------------------------+----------------------------------------------------------------------------+-------------+-----------+
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 histirical (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 |
|---|---|
| picontrol | Pre-industrial climate as simulated by the Earth System Models (ESMs) |
| historical | Historical climate as simulated by the ESMs, starting in 1950. |
| ssp126 | SSP1-RCP2.6 climate as simulated by the ESMs |
| ssp585 | SSP5-RCP8.5 climate as simulated by the ESMs. |
| sspXXX | ANY OTHERS? |
| Scenario specifier | Description |
|---|---|
| histsoc | Varying direct human influences in the historical period (1850-2014) (i.e. historical estimates of fishing effort). |
| 2015soc | Fixed year-2015 direct human influences (i.e. fishing effort). |
| OSPX | Future fishing determined by SSP and OSP driver forcings for OSPX |
| OSPX | Future fishing determined by SSP and OSP driver forcings for OSPX |
| nat | No fishing (naturalized run). |
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 | Institution | Original resolution |
|---|---|---|---|
| GFDLESM4 | gfdlesm4 | National Oceanic and Atmospheric Administration, Geophysical Fluid Dynamics Laboratory, Princeton, NJ 08540, USA | 288x180 |
| IPSLCM6ALR | ipslcm6a-lr | Institut Pierre Simon Laplace, Paris 75252, France | 144x143 |
| OTHERS?? |
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 | Sp ecifier | Unit | Res olution | ESM datasets | +=============================================================================+===================+=============+=================+==============+ | Mass Con centration of Total Phy toplankton Expressed as C hlorophyll | chl | kg m-3 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Floor Depth | - * deptho** | m | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Downward Flux of P articulate Organic Carbon | exp c-bot | mol m-2 s-1 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | P articulate Organic Carbon Content | - * intpoc** | kg m-2 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Primary Organic Carbon Production by All Types of Phy toplankton | intpp | mol m-2 s-1 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Net Primary Organic Carbon Production by Diatoms | intp pdiat | mol m-2 s-1 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Net Primary Mole Pr oductivity of Carbon by D iazotrophs | intp pdiaz | mol m-2 s-1 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Net Primary Mole Pr oductivity of Carbon by Pico phy toplankton | intp ppico | mol m-2 s-1 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Maximum Ocean Mixed Layer Thickness Defined by Sigma T | - *mlotst | m | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | -0125** | m | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Dissolved Oxygen Con centration | o2, | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | o 2-bot | mol m-2 | | | | | | | | | | | o2 -surf | mol m-2 | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | pH | ph | 1 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | p h-bot | 1 | | | | | | | | | | | ph -surf | 1 | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Phy toplankton Carbon Con centration | - \ | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | *phyc** | mol m-2 | | | | | | | | | | | phyc -vint | | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Mole Con centration of Diatoms expressed as Carbon in sea water | ph ydiat | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | phydiat -vint | mol m-2 | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Mole Con centration of D iazotrophs Expressed as Carbon in Sea Water | ph ydiaz | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | phydiaz -vint | mol m-2 | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Mole Con centration of Pico phy toplankton Expressed as Carbon in Sea Water | ph ypico | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | phypico -vint | mol m-2 | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Water Salinity | so | ‰ | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | s o-bot | ‰\ | | | | | | ‰ | | | | | so -surf | | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Water Potential T emperature | t hetao | °C | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Ocean Model Cell Thickness | thk cello | m | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Water Potential T emperature at Sea Floor | tob | °C | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Surface T emperature | tos | °C | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Water X Velocity | uo | m s-1 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Water Y Velocity | vo | m s-1 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Mole Con centration of Me soz ooplankton expressed as Carbon in sea water | zmeso | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | zmeso -vint | mol m-2 | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Mole Con centration of Mic roz ooplankton expressed as Carbon in sea water | z micro | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | - *zmicro | mol m-2 | | | | | | | | | | | -vint** | | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Z ooplankton Carbon Con centration | - \ | mol m-3 | 0.25° , 1° grid | GFDL, IPSL | | | | | | | | | *zooc** | mol m-2 | | | | | | | | | | | zooc -vint | | | | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Net Downward Shortwave Radiation at Sea Water Surface | r sntds | W m-2 | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+ | Sea Ice Area Fraction | s iconc | % | 0.25° , 1° grid | GFDL, IPSL | +-----------------------------------------------------------------------------+-------------------+-------------+-----------------+--------------+
The climate forcing input files can be found using the following pattern:
levante:/work/bb0820/ISIMIP/ISIMIP3b/InputData/climate/ocean/uncorrected/<glob
al or regional>/monthly/<climate-scenario>/<climate-forcing>/<climateforcing>_<ensemble-member>_<climate-scenario>_<climatevariable>_global_monthly_<start-year>_<end-year>.nc
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
+-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | Specifier | Included variables (short names) and definitions | Time pe riod /Resolution | Filename | +=========================+================================================================================================================================================+==========================+==========+ | histsoc | NomActive = Nominal fishing effort of the active fleet dis-aggregated by: | - 1850-2015 | | | | | | | | | - eez_country_name = The exclusive economic zone/high seas name in which fishing effort is occurring | - Annual | | | | | | | | | - LME = A number code of the Large Marine ecosystem in which the Effort is occurring | | | | | | | | | | - SAUP = A number code for the fishing country, following Sea Around Us numbering | | | | | | | | | | - Gear = the fishing gear | | | | | | | | | | - FGroup = the targeted functional group | | | | | | | | | | - Sector = the fishing sector defined by the law of the country | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | 2015soc | Final year of values from histsoc repeated until 2100 | - 2015-2100 | | | | | | | | | | - Annual | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | OSP1soc | Which variables? Determined from % change relative to 2015 in SSP1Population, SSP1GDP and relative change to 2015 (drivers of?) fishing effort | - 2015-2100 | | | | | | | | | | - Annual | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | OSP2soc | Which variables? Determined from % change relative to 2015 in SSP5Population, SSP5GDP and relative change to 2015 (drivers of?) fishing effort | - 2015-2100 | | | | | | | | | | - Annual | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | H istorical P opulation | | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | H istorical GDP | | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | OTHERS? | | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | SSP1 P opulation | | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | SSP1 GDP | | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | SSP5 P opulation | | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+ | SSP5 GDP | | | | +-------------------------+------------------------------------------------------------------------------------------------------------------------------------------------+--------------------------+----------+
| Va riable Name | Long name | Unit | Description/notes |
|---|---|---|---|
| Year | (End of the) year when the f ishing effort is occ urring | Number code | |
| Sector | The f ishing sector d efined by the law of the c ountry | Name code | I = Industrial and A = artisanal, where artisanal include powered and unpowered artisanal fleets |
| LME | Large Marine Eco system Number | Number code | A number code of the Large Marine ecosystem in which the Effort is occurring |
| eez_ countr y_name | Exc lusive Ec onomic 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 f ishing co untry, fol lowing Sea Around Us num bering | Number code | Ex supranational entities (USSR, Yugoslavia) are disaggregated to their constituent countries. Serbian Fishing Effort included with Montenegro. Crimea included with Ukrainian. |
| Gear | The f ishing gear | Name code | Gear names |
| FGroup | The ta rgeted func tional group | Name code | Functional groups are in accordance with those used by the Sea Around Us Project |
| Nom Active | N ominal f ishing effort (i.e., not inc luding the t echnol ogical 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 |
| OSP | Variable | Change relative to 2015 |
|---|---|---|
| OSP1 | ???? | ???? |
| OSP2 | ???? | ???? |
TO DO: We provide code examples showing how to implement the OSPs... NEED TO ADD TO GITHUB REPO
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/ISIMIP3b/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 | All consumers (trophic level >1, vertebrates and invertebrates) | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Total Consumer Biomass Density in log10 Weight Bins | tcblog10 | g m-2 | 0.25° grid | Level dimensions: (time, bins, lat, lon). | | | | | | | | | | | monthly | 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 | All pelagic consumers (trophic level >1, vertebrates and invertebrates) | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Total Demersal Biomass Density | tdb | g m-2 | 0.25° grid | All demersal consumers (trophic level >1, vertebrates and invertebrates) | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Total Catch Density (all commercial functional groups / size classes) | tc | g m-2 | 0.25° grid | Catch at sea (all catch as a result of all effort including reported and IUU) summed for both Industrial and Artisanal sector. | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Total Industrial Catch Density (all commercial functional groups / size classes) | tic | g m-2 | 0.25° grid | Catch at sea (all catch as a result of all effort including reported and IUU) for Industrial sector only. | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Total Catch Density in log10 Weight Bins across both sectors | tclog10 | g m-2 | 0.25° grid | Level dimensions: (time, bins, lat, lon). | | | | | | | | | | | monthly | 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 | Catch at sea of all pelagic consumers (trophic level >1, vertebrates and invertebrates) | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Total Demersal Catch Density across Artisanal and Industrial sectors | tdc | g m-2 | 0.25° grid | Catch at sea of all demersal consumers (trophic level >1, vertebrates and invertebrates) | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | 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 | If a pelagic species and L infinity is <30 cm, include in this variable | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Biomass Density of Medium Pelagics >=30cm and <90cm | bp30to90cm | g m-2 | 0.25° grid | If a pelagic species and L infinity is >=30 cm and <90cm, include in this variable | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Biomass Density of Large Pelagics >=90cm | bp90cm | g m-2 | 0.25° grid | If a pelagic species and L infinity is >=90cm, include in this variable | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Biomass Density of Small Demersals <30cm | bd30cm | g m-2 | 0.25° grid | If a demersal species and L infinity is <30 cm, include in this variable | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Biomass Density of Medium Demersals >=30cm and <90cm | bd30to90cm | g m-2 | 0.25° grid | If a demersal species and L infinity is >=30 cm and <90cm, include in this variable | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Biomass Density of Large Demersals >=90cm | bd90cm | g m-2 | 0.25° grid | If a demersal species and L infinity is >=90cm, include in this variable | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Catch Density of Small Pelagics <30cm | cp30cm | g m-2 | 0.25° grid | Catch at sea of pelagic species with L infinity <30 cm | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Catch Density of Medium Pelagics >=30cm and <90cm | cp30to90cm | g m-2 | 0.25° grid | Catch at sea of pelagic species with L infinity >=30 cm and <90 cm | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Catch Density of Large Pelagics >=90cm | cp90cm | g m-2 | 0.25° grid | Catch at sea of pelagic species with L infinity >=90 cm | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Catch Density of Small Demersals <30cm | cd30cm | g m-2 | 0.25° grid | Catch at sea of demersal species with L infinity <30 cm | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Catch Density of Medium Demersals >=30cm and <90cm | cd30to90cm | g m-2 | 0.25° grid | Catch at sea of demersal species with L infinity >=30 cm and <90 cm | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+ | Catch Density of Large Demersals >=90cm | cd90cm | g m-2 | 0.25° grid | Catch at sea of demersal species with L infinity >=90 cm | | | | | | | | | | | monthly | | +------------------------------------------------------------------------------------------------+--------------------+-------+------------+---------------------------------------------------------------------------------------------------------------------------------------------+
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.