This appendix provides full and detailed description of the database and the methodology used in this study.
Almost the entire GIS analyses were carried out using a Sun Spark Station 10.
The weather generator and BAMnut programs were run on a Pentium III DELL OptiPlex GX1 PC, 1.2 GB, 130,488 KB RAM 17 inch Video Super VGA.
Data transfers between the Sun Spark Station and the Pentium PC were carried out under Microsoft Windows NT.
A Geographic Information System (ARC, Version 7.03, ESRI, Redlands, CA, USA) which has both raster and vector capabilities was used. This system is able to efficiently store geographical referenced database which includes both digital maps and their attribute files.
The weather generator and BAMnut were written in Visual Pascal. BAMnut should be run on at least a 486PC with 8MB of RAM with Windows 95 as the operating system.
Location of climate data for downloads:
http://ipcc-ddc.cru.uea.ac.uk/cru_data/datadownload/observed/climatology_download.html
The following data were selected and downloaded; Maximum Temperature (ctmx6190.zip), Minimum Temperature (ctmn6190.zip), Precipitation (cpre6190.zip), Wet Days (cwet6190.zip), Vapour pressure (cvap6190.zip), Radiation (crad6190.zip), and Windspeed (cwnd6190.zip).
An extract of the ASCII data file ctmx6190.dat contained in ctmx6190.zp is illustrated below:
grd_sz |
xmin |
ymin |
xmax |
ymax |
n_cols |
n_rows |
n_months |
missing |
||||||||||||||
0.50 |
0.25 |
-89.75 |
359.75 |
89.75 |
720 |
360 |
12 |
-9999 |
||||||||||||||
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
|||||||
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
|||||||
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
|||||||
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
|||||||
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
|||||||
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
-9999 |
Each zipped file (e.g. ctmn6190.zip) contains 12 monthly ASCII files, so it was necessary to separate these files by month. Due to the very large size of each file (i.e. 259,200 records) the separation of these files was done using software. Once files were separated they were saved into individual files.
The header for each monthly ascii file
ncols |
720 |
nrows |
360 |
xllcorner |
-180.0 |
yllcorner |
-90.0 |
cellsize |
0.5 |
NODATA_value |
-9999.0 |
The data downloaded presented a shift at zero latitude, so a shift adjustment had to be made to all the files and this was done automatically using Java software.
The observed climate data from the University of East Anglia Climate Research Unit (CRU) data was found to be the most comprehensive and complete climate data available, however, while displaying this data using Java it was realised that a large number of islands not displayed in the CRU data visualisation page were present (i.e. http://ipcc-ddc.cru.uea.ac.uk/cru_data/visualisation/visual_index.html). To solve this problem, adjustments were made according to PET data from IIASA for the following two reasons; (a) PET data could be used as an input to the model of the present study because it was based on CRU data, and (b) the islands displayed in the PET data were identical to the CRU data visualisation page. The adjustments made are described below:
The following example illustrates how one of the text files was converted from DOS to UNIX using Arc/Info software:
Grid: asciigrid ctmx6190_0.dat ctmx6190_0.txt
Where:
asciigrid = Arc/Info command
ctmx6190_0.dat = ASCII file for maximum temperature for the month of January in DOS format
ctmx6190_0.txt = ASCII file for maximum temperature for the month of January in UNIX format.
The following example illustrates how one of the files was copied and converted into a grid using Arc/Info software:
Grid: asciigrid ctmx6190_0.txt ctmx6190_0
Where:
asciigrid = Arc/Info ASCII to grid conversion command
ctmx6190_0.txt = ASCII file for maximum temperature for the month of January
ctmx6190_0 resulting grid.
A grid was used as a template to standardise grid extensions. The cell size and geographical project of PET for the month of January was chosen as the base grid. A description of the base grid is illustrated below:
The base grid
Description of GRID /sun1 disk4/faogis4/jippe/dem/petjan |
|
Cell Size = 0.500 |
Data Type: Integer |
Number of Rows = 360 |
Number of Values = 289 |
Number of Columns = 720 |
Attribute Data (bytes) = 8 |
Boundary |
Statistics |
Xmin = -180.000 |
Minimum Value = 0.000 |
Xmax = 180.000 |
Maximum Value = 288.000 |
Ymin = -90.000 |
Mean = 63.706 |
Ymax = 90.000 |
Standard Deviation = 70.197 |
Coordinate System Description
Projection GEOGRAPHIC
Units DD Spheriod CLARKE1866
Based on the number of rows and columns defined above, the base grid had a total of grid cells (i.e. 360 x 720 = 259,200). The number of pixels for the global `land' areas alone was 62,482. Data originally stored in ASCII format or having a different resolution or projection were converted to the base grid. Example:
Grid: setwindow petjan
Grid: setmask petjan
Grid: setcell petjan
Grid: tmx6190_0 = ctmx6190_0
Where:
setwindow, setmask and setcell are Arc/Info commands
petjan = PET grid for the month of January
ctmx6190_0 = maximum temperature for the month of January
tmx6190_0 = output grid.
Tables 4, 5, 6 and 7 show statistics of the mean monthly climate datasets for global land areas used in this study as inputs to the model BAMnut developed.
A comparison of climate values between the original CRU data and the adjusted CRU data (i.e. CRU data set to the exact PET format from IIASA) showed no significant differences amongst data, the largest difference was of 0.4 when comparing the mean value of minimum temperature for the month of January.
Statistics of mean monthly air temperature values for global land areas
Month |
Maximum temperature |
Minimum temperature |
||||||
MIN |
MAX |
MEAN |
SD |
MIN |
MAX |
MEAN |
SD |
|
January |
-48.6 |
40.3 |
3.58 |
22.81 |
-56.0 |
26.5 |
-7.66 |
22.00 |
February |
-44.3 |
39.2 |
5.00 |
22.33 |
-54.2 |
26.0 |
-6.65 |
21.79 |
March |
-39.9 |
40.5 |
8.83 |
20.11 |
-56.3 |
25.9 |
-3.33 |
20.10 |
April |
-32.9 |
41.9 |
13.80 |
16.79 |
-45.2 |
26.8 |
1.66 |
16.50 |
May |
-18.6 |
43.3 |
18.40 |
13.39 |
-34.7 |
28.5 |
6.57 |
12.35 |
June |
-7.8 |
44.7 |
21.97 |
10.70 |
-24.6 |
29.2 |
10.19 |
9.68 |
July |
-5.1 |
45.6 |
23.68 |
9.39 |
-21.6 |
31.0 |
12.05 |
8.67 |
August |
-7.0 |
44.8 |
22.76 |
9.93 |
-28.5 |
30.2 |
11.19 |
9.21 |
September |
-15.4 |
42.3 |
19.53 |
12.14 |
-41.4 |
28.5 |
8.18 |
11.23 |
October |
-26.7 |
39.3 |
14.42 |
15.93 |
-44.4 |
26.9 |
3.77 |
14.37 |
November |
-38.5 |
40.1 |
8.57 |
19.85 |
-53.0 |
26.0 |
-2.29 |
18.76 |
December |
-45.7 |
40.9 |
4.74 |
21.96 |
-53.6 |
26.4 |
-6.01 |
20.78 |
Statistics of mean monthly precipitation and wet days values for global land areas
Month |
Rainfall |
Wet days |
||||||
MIN |
MAX |
MEAN |
SD |
MIN |
MAX |
MEAN |
SD |
|
January |
0 |
21.9 |
1.71 |
2.42 |
0 |
30.5 |
10.51 |
7.15 |
February |
0 |
19.5 |
1.74 |
2.48 |
0 |
27.8 |
9.36 |
6.14 |
March |
0 |
19.1 |
1.75 |
2.39 |
0 |
30.6 |
9.72 |
6.19 |
April |
0 |
19.8 |
1.75 |
2.22 |
0 |
29.8 |
8.99 |
5.39 |
May |
0 |
21.4 |
1.80 |
2.22 |
0 |
29.8 |
9.38 |
5.44 |
June |
0 |
37 |
2.12 |
2.57 |
0 |
29 |
9.70 |
5.64 |
July |
0 |
40.1 |
2.35 |
2.69 |
0 |
30.9 |
10.34 |
5.82 |
August |
0 |
39.2 |
2.31 |
2.57 |
0 |
30.9 |
10.42 |
5.83 |
September |
0 |
38.7 |
2.05 |
2.28 |
0 |
28.8 |
9.97 |
5.68 |
October |
0 |
32.7 |
1.80 |
2.08 |
0 |
27.8 |
10.25 |
6.13 |
November |
0 |
26 |
1.76 |
2.15 |
0 |
28.9 |
10.37 |
6.58 |
December |
0 |
20.8 |
1.71 |
2.31 |
0 |
29.8 |
10.62 |
7.03 |
Statistics of mean monthly vapour pressure and radiation values for global land areas
Month |
Vapour pressure |
Radiation |
||||||
MIN |
MAX |
MEAN |
SD |
MIN |
MAX |
MEAN |
SD |
|
January |
0.1 |
31.1 |
7.77 |
8.96 |
0 |
294 |
89.49 |
77.79 |
February |
0.1 |
32 |
8.01 |
9.06 |
0 |
268 |
104.78 |
72.34 |
March |
0.1 |
31.5 |
8.61 |
8.90 |
16 |
262 |
129.33 |
58.09 |
April |
0.1 |
32 |
9.53 |
8.60 |
59 |
286 |
159.31 |
41.47 |
May |
0.1 |
34.6 |
10.86 |
8.05 |
33 |
307 |
181.00 |
42.11 |
June |
0.1 |
35.8 |
12.62 |
7.40 |
22 |
335 |
190.09 |
48.61 |
July |
0.1 |
35.6 |
14.14 |
7.16 |
25 |
332 |
186.70 |
46.79 |
August |
0.1 |
34.2 |
13.83 |
7.39 |
42 |
308 |
166.79 |
47.95 |
September |
0.1 |
35.4 |
12.08 |
7.83 |
18 |
268 |
140.56 |
60.07 |
October |
0.1 |
32.3 |
10.16 |
8.2 |
0 |
282 |
114.61 |
72.79 |
November |
0.1 |
31.8 |
8.72 |
8.55 |
0 |
290 |
95.04 |
79.18 |
December |
0.1 |
31.5 |
7.96 |
8.77 |
0 |
301 |
85.78 |
79.78 |
Statistics of mean monthly wind speed values for global land areas
Month |
Wind speed |
|||
MIN |
MAX |
MEAN |
SD |
|
January |
0.1 |
14.3 |
3.16 |
1.52 |
February |
0.2 |
14.3 |
3.19 |
1.46 |
March |
-2.8 |
14.2 |
3.28 |
1.36 |
April |
-6.1 |
14.1 |
3.38 |
1.34 |
May |
-7.4 |
14.2 |
3.31 |
1.26 |
June |
-7.2 |
24.4 |
3.21 |
1.29 |
July |
-8.6 |
28.9 |
3.15 |
1.35 |
August |
-7.4 |
29 |
3.13 |
1.39 |
September |
-7.1 |
24.4 |
3.18 |
1.41 |
October |
-4.7 |
19.8 |
3.26 |
1.47 |
November |
-2.2 |
14.2 |
3.23 |
1.52 |
December |
-1 |
14.3 |
3.15 |
1.49 |
Table 8 shows the settings that were made for the weather generator. Each folder corresponds to one area of the world, the number of files created from each area and their size in Mb is also shown. Computer memory problems were encountered due to the vast amount of files that the generator had to create, so the information in Table 8 is based on trial-and-error meaning that it is information that should be used to recreate the data or is helpful in creating daily data for another year or even better for enhancing the generator used.
Settings for weather generator
Top left |
Bottom right |
Technical notes |
||||
Folder no. |
Lat. |
Long. |
Lat. |
Long. |
No. Files |
Size Mb |
1 |
90 |
0 |
75 |
180 |
505 |
9.27 |
2 |
75 |
0 |
65 |
180 |
4410 |
84.2 |
3 |
65 |
0 |
60 |
180 |
3581 |
68.3 |
4 |
60 |
0 |
50 |
180 |
5706 |
9.26 |
5 |
50 |
0 |
40 |
180 |
5283 |
100 |
6 |
40 |
0 |
30 |
180 |
4453 |
85 |
7 |
30 |
0 |
20 |
180 |
4560 |
87.1 |
8 |
20 |
0 |
10 |
180 |
2968 |
56.5 |
9 |
10 |
0 |
0 |
180 |
2250 |
42.7 |
10 |
0 |
0 |
-10 |
180 |
1881 |
35.6 |
11 |
-10 |
0 |
-30 |
180 |
4210 |
80.3 |
12 |
-30 |
0 |
-60 |
180 |
1225 |
23 |
13 |
-60 |
0 |
-90 |
180 |
0 |
0 |
14 |
90 |
-180 |
85 |
0 |
0 |
0 |
15 |
85 |
-180 |
83 |
0 |
28 |
0.55 |
16 |
83 |
-180 |
75 |
-90 |
305 |
5.43 |
17 |
83 |
-90 |
75 |
0 |
0 |
0 |
18 |
75 |
-180 |
65 |
-90 |
1945 |
36.9 |
19 |
75 |
-90 |
65 |
0 |
0 |
0 |
20 |
65 |
-180 |
60 |
0 |
2100 |
39.8 |
21 |
60 |
-180 |
50 |
0 |
3054 |
58.2 |
22 |
50 |
-180 |
40 |
0 |
2599 |
49.4 |
23 |
40 |
-180 |
30 |
0 |
2125 |
40.3 |
24 |
30 |
-180 |
20 |
0 |
1302 |
24.5 |
25 |
20 |
-180 |
10 |
0 |
1173 |
22 |
26 |
10 |
-180 |
0 |
0 |
1292 |
24.3 |
27 |
0 |
-180 |
-10 |
0 |
1738 |
32.9 |
28 |
-10 |
-180 |
-30 |
0 |
2534 |
48.2 |
29 |
-30 |
-180 |
-60 |
0 |
1255 |
23.6 |
30 |
-60 |
-180 |
-90 |
0 |
0 |
0 |
62,482 |
Output data format for one of the files created for the United Republic of Tanzania, only the first 10 days of the simulated year (1990) are shown:
@Date |
SRAD |
Tmax |
Tmin |
Rain |
Evap |
MNVP |
Wind |
ACO2 |
90001 |
6.2 |
25.3 |
17.1 |
4.9 |
-99. |
-99. |
-99. |
340.0 |
90002 |
11.8 |
29 |
20.9 |
15 |
-99. |
-99. |
-99. |
340.0 |
90003 |
10.8 |
29.6 |
19 |
6.1 |
-99. |
-99. |
-99. |
340.0 |
90004 |
4.6 |
30.3 |
20.3 |
0.2 |
-99. |
-99. |
-99. |
340.0 |
90005 |
1.6 |
29.4 |
23.6 |
8.4 |
-99. |
-99. |
-99. |
340.0 |
90006 |
9 |
31.5 |
22.2 |
10 |
-99. |
-99. |
-99. |
340.0 |
90007 |
16.5 |
28.5 |
21.2 |
3.5 |
-99. |
-99. |
-99. |
340.0 |
90008 |
9.5 |
25.8 |
19.1 |
1.4 |
-99. |
-99. |
-99. |
340.0 |
90009 |
14.1 |
23.3 |
17.5 |
5.3 |
-99. |
-99. |
-99. |
340.0 |
90010 |
16.7 |
26.6 |
17.4 |
5.4 |
-99. |
-99. |
-99. |
340.0 |
The weather generator takes approximately 160 hours to generate the required data, however, it can now serve as inputs to current and future versions of BAMnut, thus, there is no need to recreate daily data. The present GIS version of BAMnut takes 80 minutes to create the biomass and pod yield files, and it is expected that the Web version of BAMnut will the take about the same amount of time to generate the two files. However, the user will be able to work in other computer applications and leave the model running, and the Web model will have a report progress function, so that the user will not have to wait for model outputs.
The result grids were converted to an equal-area projection (Flat Polar Quartic) to calculate the areas covered by each class in each country of Africa.
Arc: project grid bio11feb_5 biomass_pq llpolq.prj
Arc: project grid pod11feb_5 podyield_pq llpolq.prj
To convert grids to Flat Polar Quartic projection the parameters specified in the file llpolq.prj illustrated below were used:
INPUT
PROJECTION GEOGRAPHIC
UNITS DD
PARAMETERS
OUTPUT
PROJECTION FLAT POLAR QUARTIC
UNITS METERS
PARAMETERS
00 00 00
end
After the biomass and pod yield maps had been converted to a polar quartic projection, they presented a distortion in the projection when displayed. To solve this problem a polar quartic projection mask (named templpolq) was used:
Arc: describe templpolq
Description of Grid templpolq |
|
Cell Size = 9250.000 |
Data Type: Integer |
Number of Rows = 1827 |
Number of Values = 1 |
Number of Columns = 4057 |
Attribute Data (bytes) = 8 |
Boundary |
Statistics |
Xmin = -18761636.000 |
Minimum Value = 1.000 |
Xmax = 18765614.000 |
Maximum Value = 1.000 |
Ymin = -8445739.000 |
Mean = 1.000 |
Ymax = 8454011.000 |
Standard Deviation = 0.000 |
Coordinate System Description
Projection FLAT_POLAR_QUARTIC |
|||
Units |
METERS |
Spheroid |
SPHERE |
Parameters: |
|||
Longitude of projection centre |
0 |
0 |
0.000 |
Procedures for projection adjustment using a mask:
Grid: setmask templpolq
Grid: biomass_pq2 = biomass_pq
Grid: setmask templpolq
Grid: podyield_pq2 = podyield_pq
The coverage WFSCOV, which contains the country boundaries (i.e. areas) was converted to a GRID:
Arc: polygrid wfscov wfscov_grd cntcode
Where:
polygrid = GRID command
wfscov = coverage
wfscov_grd = output grid
cntcode, coverage item containing to country boundary names.
The cell size was set to an equal area projection cell size of 59311.282.
The GRID wfscov_grd was overlaid on the biomass and pod yield maps in order to produce the statistics by country.
Grid: biomass_stat = combine(wfscov_grd, biomass_pq2)
Grid: podyield_stat = combine(wfscov_grd, podyield_pq2)
Where:
biomass_stat and podyield_stat = output grids
combine = GRID command
wfscov_grd = world grid in flat polar quartic projection
biomass_pq2 and podyield_pq2 = biomass and pod yield grids in flat polar quartic projection.
Biomass
ENTER COMMAND >SEL BIOMASS_STAT.VAT
ENTER COMMAND >OUTPUT ../BIOMASS_STAT.TXT
ENTER COMMAND >DISPLAY COUNT,',',WFSCOV_GRD,',',BIOMASS_PQ2 PRINT
ENTER COMMAND > Q STOP
Pod Yield
ENTER COMMAND >SEL PODYIELD_STAT.VAT
ENTER COMMAND >OUTPUT ../PODYIELD_STAT.TXT
ENTER COMMAND >DISPLAY COUNT,',',WFSCOV_GRD,',',PODYIELD_PQ2 PRINT
ENTER COMMAND > Q STOP
The following example illustrates how one of the text files was converted from UNIX to DOS using Arc/Info software:
Grid: gridascii biomass_stat.txt biomass_stat.dat
Where:
gridascii = Arc/Info command
biomass_stat.txt = biomass file in UNIX
biomass_stat.dat = resulting biomass file in DOS format.
For example: 59.30642 x 59.30642 = 3517.251
Suitability class |
Count |
Cell size |
Count x cell size |
Sum |
Percentage |
2 |
155 |
3517.251 |
545174 |
622553.5 |
87.57 |
2 |
22 |
3517.251 |
77379.53 |
622553.5 |
12.42 |
The procedures used to create the map compositions for the maps presented in this study are illustrated below using the biomass grid as an example:
For purposes of analysis and/or illustration, the biomass grid had to be reclassified to a common colour range.
biomass_5 = reclass (biomass, bio11feb_5.rem)
Where:
biomass_5 = reclassified grid file
reclass = GRID function to reclassify (or change) integer values of the input cells using a remap table on a cell-by-cell basis within the analysis window
biomass = original water requirement grid
biomass_5.rem = remap table.
The remap table biomass_5.rem is shown below:
0 - 0 : 1
1-1,500: 2
1,500-4,500: 3
4,500-8,500: 4
>8,500 : 5
The AML, biomass.aml illustrated below was used in ARCPLOT to generate a map composition for the biomass grid:
killmap biomass_mc
mape biomass
pagesize 11.7 8.3
map biomass_mc
linesymbol 5
box 0.1 0.1 11.6 8.2
mapposition cen cen
shadeset BAMnut.shd
gridnodatasymbol white
grids biomass
linecolor black
arcs/sun1disk1/faogistemp/ctry25m
textset font
textsymbol 1
textsize 0.18
keyarea 0.6 1 10.7 1.23
keybox 0.15 0.1
keyseparation .13 .14
keyshade biomass.leg
The maps which are presented in this study were converted to eps (i.e. postscript format). The following example illustrates how the biomass grid (map) was converted into and eps file:
gissw2-faogis>> setenv CANVASCOLOR WHITE
Arcplot: &r biomass.aml
Arcplot: display 1040
Arcplot: biomass.gra
Arcplot: plot biomass_mc box 0.1 0.1 11.6 8.2
Arc: postscript biomass.gra biomass.eps
where:
setenv CANVASCOLOR WHITE = sets the background colour to white
&r biomas.aml = automates the map composition of the grid
display 1040 = command used to save the map composition
biomass.gra = graphics file
plot biomass_mc box 0.1 0.1 11.6 8.2 = map composition which was reduced to suit the size required for publication
postscript = ARC command used to create the postscript file.