cimp5_1951_2000_av: Model-Averaging of Decadal Mean Seafloor Climatic Data during...

In chihlinwei/SCC85: Seafloor Climate Change Projections from CIMP5

Model-Averaging of Decadal Mean Seafloor Climatic Data during 1951 to 2000

Description

Ensemble average of the decadal mean seafloor (bottom-most) projections based on historical scenario of CIMP5 experiment.

Format

A RasterBrick object of 4 raster layers:

epc_av_1951_to_2000

Mean export POC flux to seafloor (mg C m^-2 d^-1)

o2_av_1951_to_2000

Mean dissolved oxygen concentration at seafloor (mol m^-3)

ph_av_1951_to_2000

Mean pH at seafloor

thetao_av_1951_to_2000

Mean potential temperature at seafllor (K)

arag_av_1951_to_2000

Aragonite Concentration (mol m-3)

calc_av_1951_to_2000

Calcite Concentration (mol m-3)

co3_av_1951_to_2000

Mole Concentration of Carbonate expressed as Carbon in Sea Water (mol m-3)

co3satarag_av_1951_to_2000

Carbonate ion concentration for seawater in equilibrium with pure aragonite (mol m-3)

co3satcalc_av_1951_to_2000

Carbonate ion concentration for seawater in equilibrium with pure calcite (mol m-3)

aragsat_av_1951_to_2000

Aragonite Saturation State

calcsat_av_1951_to_2000

Calcite Saturation State

Details

Yearly mean projections from 1951 to 2000 were calculated from each of the GFDL-ESM-2G, IPSL-CM5A-MR and MPI-ESM-MR. Decadal mean was calculated from the yearly mean and than averaged across the three models.

Source

https://esgf-node.llnl.gov/search/esgf-llnl/

Examples

# Mask the raster brick by 200 to 2000 m
r0 <- mask(cimp5_1951_2000_av, mask2000)

# Only show color between 1 to 99 percentile
r <- raster()
for(i in 1:4){
  d <- subset(r0, i)
  ma <- quantile(d, 0.99)
  mi <- quantile(d, 0.01)
  d[d>ma] <- ma
  d[d<mi] <- mi
  if(i==2) d <- d*15.9994 # to mg/L
  if(i==3) d <- log10(d)*(-1) # to total scale
  if(i==4) d <- d-273.15 # to degree C
  r <- addLayer(r, d)
}
names(r) <- names(r0)

# plot on google earth
library(plotKML)
for(i in 1:4){
  if(i==1|i==4) cols <- jet.col.log(100) else cols <- jet.colors(100)
  plotKML(subset(r, i), folder.name=names(r)[i],
          colour_scale = cols,
          raster_name = paste(names(r)[i], "png", sep="."))
}

chihlinwei/SCC85 documentation built on Sept. 15, 2023, 4:47 a.m.