Withdrawal to Availability Ratio

TOP ＞ RESULT ＞ Global 0.5deg Grid ＞ Annual ＞ Annual Withdrawal to Availability Ratio

(Original Dir.:data/global/Global(*.asc file))

Index of Water Risk

Withdrwal to availability ratio
This is expressed as (W - S) / Q, where, during a certain period, respectively. W is written in Withdrawal data, and Q is shown in Discharge data.
Larger value of this index means that area(grid cell) is under severer water stress. Value more than 0.4 is especially regarded as index of fatal water stress for human lives.
 

Data (Corrected Version)

  1. Withdrawal to Availability Ratio
  2. W-to-A ratio for each alpha value:
  3. Accumulated Histogram of population under each w-to-a values:

Figures (Corrected Version)

Withdrawal to Availability Ratio, Rws(WRI version)

Change in Rws(WRI version, 1995-2050)

Difference
Ratio

Withdrawal to Availability Ratio Rws(D,I = WRI, Agri = EPIC real)

Withdrawal to Availability Ratio for each alpha value

For river discharge at each alpha, see Annual River Discharge page.

各αに対する渇水指標ごとの人口累積ヒストグラム

(postscript)

Rws for each scenario

データ (Saruhashi Version)

計算方法は下記参照．

  1. 各年の渇水指標
  2. αの値を変えた時の渇水指標(1995). 各αに対する年間河川流量の値は年河川流量の項目参照．ファイル名はAlphaRatio${YEAR}.ascの形．
  3. αを変えた時の，渇水指標別面積の累積ヒストグラム

図 (Saruhashi Version)

αの値を変えた時の渇水指標(1995). 各αに対する年間河川流量の値は年河川流量の項目参照．

各αに対する渇水指標ごとの人口累積ヒストグラム

(postscript)

計算方法およびプログラム

計算「アルゴリズム」

プログラム

（プログラムはまだ公開工事中）

Woking Directory : data/global/Global

Withdrawal to Availability Ratio Rws (no scenario)

foreach YEAR ( 1995 2050 )
	setenv YEAR2 `bin/year2.sh $YEAR`
	setenv BASE ratio${YEAR}-2
	bin/asc_calc XYZ/DIA${YEAR2}-2_106m3.asc '-' XYZ/desal${YEAR2}_106m3.asc \
		'/' ../../discharge/annual/XYZ/Mean${YEAR}_106m3.asc \
	> XYZ/${BASE}.asc
	bin/asc2xyz < XYZ/${BASE}.asc > XYZ/${BASE}.xyz
	bin/globalxyz2image.sh "$BASE"
end

Rws in 2050 (scenario 1: Population only, alpha=1.0)

setenv YEAR_DISC 1995
setenv YEAR_DIA  2050
setenv SCENARIO  1
setenv BASE      "ratio${YEAR_DIA}-2-S${SCENARIO}"
bin/asc_calc XYZ/DIA${YEAR_DIA}-2-S${SCENARIO}_106m3.asc '-' XYZ/desal${YEAR_DIA}_106m3.asc \
	'/' ../../discharge/annual/XYZ/Mean${YEAR_DISC}_106m3.asc \
> XYZ/${BASE}.asc
bin/asc2xyz < XYZ/${BASE}.asc > XYZ/${BASE}.xyz
bin/globalxyz2image.sh "$BASE"

Rws in 2050 (scenario 2: Population + Climatic Change, alpha=1.0)

setenv YEAR_DISC 2050
setenv YEAR_DIA  2050
setenv SCENARIO  2
setenv BASE      "ratio${YEAR_DIA}-2-S${SCENARIO}"
bin/asc_calc XYZ/DIA${YEAR_DIA}-2-S${SCENARIO}_106m3.asc '-' XYZ/desal${YEAR_DIA}_106m3.asc \
	'/' ../../discharge/annual/XYZ/Mean${YEAR_DISC}_106m3.asc \
> XYZ/${BASE}.asc
bin/asc2xyz < XYZ/${BASE}.asc > XYZ/${BASE}.xyz
bin/globalxyz2image.sh "$BASE"

Rws in 2050 (scenario 4: Population + Industrial development + Climatic Change, alpha=1.0)

setenv YEAR_DISC 2050
setenv YEAR_DIA  2050
setenv SCENARIO  4
setenv BASE      "ratio${YEAR_DIA}-2-S${SCENARIO}"
bin/asc_calc XYZ/DIA${YEAR_DIA}-2-S${SCENARIO}_106m3.asc '-' XYZ/desal${YEAR_DIA}_106m3.asc \
	'/' ../../discharge/annual/XYZ/Mean${YEAR_DISC}_106m3.asc \
> XYZ/${BASE}.asc
bin/asc2xyz < XYZ/${BASE}.asc > XYZ/${BASE}.xyz
bin/globalxyz2image.sh "$BASE"

各αの渇水指標図を求めるスクリプト

foreach YEAR (1995 2050)
  setenv YEAR2  `bin/year2.sh $YEAR`
  foreach RATIO ( 0.0 0.1 0.2 0.5 1.0 )
    setenv RUNOFF Alpha${YEAR}_${RATIO}_106m3
    setenv USE    DIA${YEAR2}-2_106m3
    setenv DESAL  desal${YEAR2}_106m3
    setenv RESULT  AlphaRatio${YEAR}-2_${RATIO}
    bin/asc_calc XYZ/${USE}.asc "-" XYZ/${DESAL}.asc \
        '/' XYZ/${RUNOFF}.asc > XYZ/${RESULT}.asc
    bin/asc2xyz --lack 0 XYZ/${RESULT}.asc XYZ/${RESULT}.xyz
    bin/globalxyz2image.sh $RESULT
  end
end

Population under each rank of w-to-a ratio

setenv YEAR 1995
setenv TEMPFILE temp.asc
setenv POPFILE  ../../ciesin/XYZ/glp${YEAR}agiC.asc
foreach RISK ( 0.1 0.2 0.4 )

  bin/asc_calc $POPFILE \
    'where' XYZ/ratio${YEAR}-2.asc '<' $RISK > $TEMPFILE
  bin/sumASC $TEMPFILE
end
\rm -f $TEMPFILE
echo -n 'all : '
bin/sumASC $POPFILE

各αに対して，渇水指標ごとの人口累積ヒストグラムを作る

foreach RATIO ( "0.0" "0.1" "0.2" "0.5" "1.0" )
  set NUMCLASS=100
  set THRE_MAX=2.0
  set INTERVAL=`echo $THRE_MAX $NUMCLASS | awk '{print $1 / $2}'`
  set COUNT=1
  set TEMPFILE=temp_$$.asc
  set RESULTFILE=XYZ/AccumAlphaRatio1995-2_${RATIO}.xyz

  echo 0  0  > $RESULTFILE

  while ( $COUNT <= $NUMCLASS )
     set THRE=`echo $INTERVAL $COUNT | awk '{print $1 * $2}'`
     echo -n $THRE ' '  >> $RESULTFILE

     bin/asc_calc ../../ciesin/XYZ/glp1995agiC.asc \
        '/' 100000000 \
        'where' XYZ/AlphaRatio1995-2_${RATIO}.asc '<' $THRE \
         'sum' >> $RESULTFILE

     set COUNT=`expr $COUNT + 1`
  end
end
\rm -f $TEMPFILE

bin/accumpop.gmt 'Accum*'

Population accumation for each scenarios

foreach SCENARIO ( "1" "2" "4" )
  setenv YEAR       2050
  setenv ALPHA      "1.0"
  setenv NUMCLASS   100
  setenv THRE_MAX   2.0
  setenv INTERVAL   `echo $THRE_MAX $NUMCLASS | awk '{print $1 / $2}'`
  set COUNT=1
  setenv RESULTFILE "XYZ/AccumAlphaRatio${YEAR}-2-S${SCENARIO}_${ALPHA}.xyz"
  setenv POPFILE    "../../ciesin/XYZ/glp${YEAR}agiC.asc"

  echo 0  0  > $RESULTFILE

  while ( $COUNT <= $NUMCLASS )
     set THRE=`echo $INTERVAL $COUNT | awk '{print $1 * $2}'`
     echo $SCENARIO '-' $COUNT ' : ' $NUMCLASS

     echo -n $THRE ' '  >> $RESULTFILE

     bin/asc_calc ${POPFILE} \
        '/' 100000000 \
        'where' XYZ/ratio${YEAR}-2-S${SCENARIO}.asc '<' $THRE \
         'sum' >> $RESULTFILE

     set COUNT=`expr $COUNT + 1`
  end
end

Continental Discharge, Water Use and Withdrawal-to-availability ratio

setenv REGION_FILE_BASE ../../BaseMap/region_code4
setenv REGION_MAP_FILE  ${REGION_FILE_BASE}.asc
setenv REGION_LIST_FILE ${REGION_FILE_BASE}.txt
setenv YEAR 1995
setenv YEAR2 `expr $YEAR % 100`

foreach R ( 1 2 4 5 6 8 )
  setenv REGION_NAME `bin/regionname.sh $R  $REGION_LIST_FILE`
  if ( "$REGION_NAME" != "" ) then
    setenv RUNOFF `bin/calctotal_int ../../discharge/annual/XYZ/Ro.Mean1995_106m3.asc\
      $REGION_MAP_FILE $R | cut -f3`
    setenv POP    `bin/calctotal_int ../../ciesin/XYZ/glp${YEAR}agiC.asc \
      $REGION_MAP_FILE $R | awk '{print $3 / 1000}'`
    setenv INDUS  `bin/calctotal_int XYZ/indus${YEAR2}_106m3.asc \
      $REGION_MAP_FILE $R | awk '{print $3 / 1000}'`
    setenv DOMES  `bin/calctotal_int XYZ/domes${YEAR2}_106m3.asc \
      $REGION_MAP_FILE $R | awk '{print $3 / 1000}'`
    setenv IRR    `bin/calctotal_int XYZ/irr${YEAR2}_106m3.asc \
      $REGION_MAP_FILE $R | awk '{print $3 / 1000}'`
    setenv DIA    `bin/calctotal_int XYZ/DIA${YEAR2}-2_106m3.asc \
      $REGION_MAP_FILE $R | awk '{print $3 / 1000}'`
    echo $R $REGION_NAME $POP $RUNOFF $DOMES $INDUS $IRR $DIA  |\
		 awk '{printf("%3d\t%-15s\t%10.0f\t%7.1f\t%5.1f\t%5.1f\t%5.1f\t%5.1f\t%5.1f\t%5.1f\n", \
			$1, $2, $3, $4, $5, $6, $7, $8, $9, $9 / $5)}'
  else
    echo $R : no such region. skipped
  endif
end
AGATASHI