File:East sibiria at lgm 2.svg
Original file (SVG file, nominally 1,053 × 744 pixels, file size: 930 KB)
Captions
Summary
[edit]DescriptionEast sibiria at lgm 2.svg |
English: East Siberia at Last Glacial Maximum (LGM). Almost all sites are above 9.5 ° C Tjuly isotherm. There was no sure habitation above 58 °N latitude. This means that tundra north from treeline was unhabitated. |
|||
Date | ||||
Source | Own work | |||
Author | Merikanto | |||
SVG development InfoField | This map was created with an unknown SVG tool.
|
LGM human occupation sites data is from Kuzmin&Keates 2016
Kuzmin, Yaroslav & G. Keates, Susan. (2016). Siberia and neighboring regions in the Last Glacial Maximum: did people occupy northern Eurasia at that time?. Archaeological and Anthropological Sciences.
doi:10.1007/s12520-016-0342-z
Image reproduction scripts in Ubuntu etc.
Init script
- get LGM climate data geotiff
mkdir data
cd data
wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmtn_2-5m.zip
wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmtx_2-5m.zip
wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmpr_2-5m.zip
wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmbi_2-5m.zip
unzip cclgmtn_2-5m.zip
unzip cclgmtx_2-5m.zip
unzip cclgmpr_2-5m.zip
unzip cclgmbi_2-5m.zip
rm *.zip
cd ..
mkdir plot
mkdir teksti
- init mapping utilzzz
sudo apt install gdal-bin python-gdal nco cdo ctioga2 libnetcdf libnetcdf-dev -y
- sumenna.bash
sudo apt install saga
- gcc -o ncblit1 ncblit1.c -lm -lnetcdf
- gcc -o ncprobe1 ncprobe1.c -lm -lnetcdf
Data downloading script, in ccsm dir
sorckanta1="cclgmtn"
sorckanta2="cclgmtx"
sorckanta3="cclgmpr"
sorckanta4="cclgmtk"
psorckanta1="cclgmbi"
cd data
for i in {1..12..1}
do
sorcnimi1=$sorckanta1$i.tif
sorcnimi2=$sorckanta2$i.tif
sorcnimi3=$sorckanta3$i.tif
sorcnimi4=$sorckanta4$i.tif
nesorcnimi1=$sorckanta1$i.nc
nesorcnimi2=$sorckanta2$i.nc
nesorcnimi3=$sorckanta3$i.nc
nesorcnimi4=$sorckanta4$i.nc
gdal_calc.py -A $sorcnimi1 -B $sorcnimi2 --outfile=$sorcnimi4 --calc="(A+B)/20"
gdal_translate -ot Float32 -of netCDF $sorcnimi1 $nesorcnimi1
gdal_translate -ot Float32 -of netCDF $sorcnimi2 $nesorcnimi2
gdal_translate -ot Float32 -of netCDF $sorcnimi3 $nesorcnimi3
gdal_translate -ot Float32 -of netCDF $sorcnimi4 $nesorcnimi4
done
for i in {1..15..1}
do
psorcnimi1=$psorckanta1$i.tif
pesorcnimi1=$psorckanta1$i.nc
gdal_calc.py -A $psorcnimi1 --outfile=temp.tif --calc="A"
gdal_translate -ot Float32 -of netCDF temp.tif $pesorcnimi1
done
cd ..
Gaussian interpolation script, w/saga GIS
ingrid1="../ccsm/data/cclgmtk7.nc"
ingrid2="../ccsm/data/cclgmtk1.nc"
ingrid3="../ccsm/data/cclgmbi1.nc"
outgridbase1="tk7_gauss1"
outgridbase2="tk1_gauss1"
outgridbase3="ta_gauss1"
saga_cmd grid_filter 1 -INPUT $ingrid3 -RESULT $outgridbase3 -SIGMA 6 -RADIUS 6
gdal_translate -of GMT $outgridbase3".sdat" temp.nc
gdal_calc.py -A temp.nc --outfile="temp2.tif" --calc="(A/10)"
gdal_translate -of netCDF temp2.tif $outgridbase3".nc"
exit
saga_cmd grid_filter 1 -INPUT $ingrid1 -RESULT $outgridbase1 -SIGMA 6 -RADIUS 6
gdal_translate -of GMT $outgridbase1".sdat" temp.nc
gdal_calc.py -A temp.nc --outfile="temp2.tif" --calc="(A)"
gdal_translate -of netCDF temp2.tif $outgridbase1".nc"
saga_cmd grid_filter 1 -INPUT $ingrid2 -RESULT $outgridbase2 -SIGMA 6 -RADIUS 6
gdal_translate -of GMT $outgridbase2".sdat" temp.nc
gdal_calc.py -A temp.nc --outfile="temp2.tif" --calc="(A)"
gdal_translate -of netCDF temp2.tif $outgridbase2".nc"
CertainLGM sites loydot.txt, charchoal
091.43 52.97 24 0 14 BL Ui 1
087.95 55.90 24 0 14 BL Shestakovo
091.10 55.96 24 0 14 BL Novoselovo
092.52 55.33 24 0 14 BL Derbina 4
093.23 56.05 24 0 14 BL Gosudarev
103.50 52.90 24 0 14 BL Malta
103.42 52.95 24 0 14 BL Buret
100.33 58.33 24 0 14 BL Ust’-Kova
108.33 57.83 24 0 14 BL Alekseevsk 1
129.58 51.83 24 0 14 BL Ust’-Ulma 1
127.66 50.33 24 0 14 BL Khodulikha
142.40 46.77 24 0 14 BL Ogonki 5
Unvertain sites, bones
083.52 57.75 BL Mogochino
084.60 51.30 BL Anui 2, layer 3
090.95 55.01 BL Sabanikha
091.42 55.13 BL Kashtanka 1
092.05 55.20 BL Shlenka
108.50 50.18 BL Studenoe 2
112.50 57.50 BL Tesa
117.47 59.93 BL Khayrgas
133.75 63.17 BL Ikhine 2
ann takayrat.txt
-15
jan tk1kayrat.txt
-40
jul tk7kayrat.txt
9.5
ta.cpt
-30 100/200/100 10.0 100/200/100 ; arctic
10.0 255/127/50 20 255/127/50 ; temporal, subtropic.
20.0 255/100/100 30 255/100/100 ; subtropic. tropic
B blue
F white
N 127.5
Image gen
- siberia
lona=80
lonb=140
lata=40
latb=80
projektio=-JQ120/10.0i
- paletti="zones.cpt"
paletti="ta.cpt"
tk7="../ccsm/data/cclgmtk7.nc"
tk1="../ccsm/data/cclgmtk1.nc"
tka="../ccsm/data/cclgmtbi12.nc"
tk7s1="tk7_gauss1.nc"
tk1s1="tk1_gauss1.nc"
tas1="ta_gauss1.nc"
ta="tk7kayrat.txt"
tb="tk1kayrat.txt"
tc="takayrat.txt"
gmt makecpt -Crainbow -T-20/20/5 > t.cpt
- gmt makecpt -Crainbow -T9/10/5 > tk.cpt
gmt grdimage $tk7 -R$lona/$lonb/$lata/$latb $projektio -C$paletti -K > out.ps
gmt pscoast -K -O $projektio -R$lona/$lonb/$lata/$latb -B10g5/5g5 -S180/180/255 -W1p,200/200/255 \
-I1/1p,200/200/255 -I2/1p,200/200/255 -I3/0.5p,200/200/255 -I4/0.25p,200/200/255 \
-P -N1/1.5p,127/0/0,.- -Dc >>out.ps
gmt grdcontour $tk7s1 -R$lona/$lonb/$lata/$latb $projektio -W0.75,red@25 -C$ta -K -O >> out.ps
- gmt grdcontour $tk1s1 -R$lona/$lonb/$lata/$latb $projektio -W2,blue@25 -C$tb -K -O >> out.ps
- gmt grdcontour $tas1 -R$lona/$lonb/$lata/$latb $projektio -W3,darkgreen@25 -C$tc -K -O >> out.ps
cat loydot.txt | gmt psxy -R$lona/$lonb/$lata/$latb $projektio -O -K -G100/0/0 -Sc0.75 >> out.ps
cat loydote.txt | gmt psxy -R$lona/$lonb/$lata/$latb $projektio -O -K -G0/0/0 -Sc0.5 >> out.ps
gmt psxy -R$lona/$lonb/$lata/$latb $projektio -O -K -G255/0/0 -St0.75i << EOF >> out.ps
103.5 52.9
EOF
gmt pstext -R -J -O -K -F+f60,Helvetica-Bold,white=thin >> out.ps << END
100 65 Siberia
END
gmt pstext -R -J -O -K -F+f40,Times-Italic,blue >> out.ps << END
109 56 Baikal
END
gmt pstext -R -J -O -K -F+f75,Helvetica-Bold >> out.ps << END
120 75 LGM sites
END
gmt pslegend -R -J -O -K -DjTL+w2.5i+o0.4i -F+pthick+ithinner+gwhite --FONT_ANNOT_PRIMARY=17p,Helvetica-Bold<< EOF >> out.ps
S 0.1i t0.75 0.08i red thinnest 0.3i Baikal site
S 0.1i c0.5 0.08i 100/0/0 thinnest 0.3i Bone sites
S 0.1i c0.5 0.08i black thinnest 0.3i Charcoal sites
S 0.1i - 0.08i red thick,red 0.3i Tjuly 9.5 oC
- S 0.1i f0.1i+l+t 0.25i blue 0.25p 0.3i This is a fault
EOF
- gmt psscale -Dx14c/1.25c+w10c/0.5c+h -B+tcontinuous -Ct.cpt -By+loC -O -Ac >> out.ps
ps2pdf out.ps out.pdf
pdf2svg out.pdf out.svg
- svgo out.svg
Licensing
[edit]- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
- share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 07:14, 7 December 2018 | 1,053 × 744 (930 KB) | Merikanto (talk | contribs) | User created page with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.
File usage on other wikis
The following other wikis use this file:
- Usage on fi.wikipedia.org
Metadata
This file contains additional information such as Exif metadata which may have been added by the digital camera, scanner, or software program used to create or digitize it. If the file has been modified from its original state, some details such as the timestamp may not fully reflect those of the original file. The timestamp is only as accurate as the clock in the camera, and it may be completely wrong.
Width | 842pt |
---|---|
Height | 595pt |