Text 2011
Text 2011
Text 2011
+
h
h
1 0
1
(1) 257
where,
h
and
h
represent infinitesimal longitudinal and shear strain components parallel 258
to Y
1
and X
1
axes of a two dimensional horizontal reference frame. To derive the two 259
strain components at different points in Andaman and its adjoining areas we apply the 260
following equations (Jaeger, 1964), 261
( )
2
B A+ =( )
2 2
2
h h
+ + .(2) 262
( )
2 2 2
h h
B A + = .(3) 263
where, A and B represent major and minor axes of the total strain ellipse. Simplifying 264
equations (2) and (3) we get, 265
( ) ( ) { } + = + = + 1 4 1
2 2
B A B A
h
, and,...(4) 266
where, is dilation. 267
Now, replacing
h
in eq.(3) by we get 268
( )
2 2 2
h
B A + = . (5a) 269
Replacing by (AB-1)and rearranging we get, 270
( ) ( ) { }
2 2
1 = AB B A
h
.(5) 271
In case of a subduction zone, however, it would be more appropriate to express 272
the total infinitesimal strain in terms of simple inclined transpression model (Jones et al., 273
2004). This is achieved by considering a three dimensional Cartesian reference frame, 274
whose XY plane is coplanar with the deformation zone dipping at an angle . Let the X- 275
axis in this inclined reference frame coincide with the X
1
axis of the horizontal reference 276
system. 277
According to the simple inclined transpression model, any general deformation in 278
the inclined 3D- reference frame can be factorized into pure shears parallel to Y- and Z- 279
axes (
y
and
z
respectively) and simple shears parallel to X- and Z axes (
x
and
z
280
respectively). The model assumes that there is no stretch parallel to X. Expressing the 281
components of pure and simple shear in terms of
h
and
h
we get, 282
( )
h z y
+ = =
1
1
, .. (6) 283
Sin
h
xy
= , and (7) 284
Cot
Cot
h
yz
= , (8) 285
where is the angle measured on the horizontal plane between the shear zone boundary 286
and the shortening direction. 287
We also know (Jones et.al, 2004) that, 288
.
1
h
h
yz
Cot
+
= .(8a) 289
Equating the r.h.s of eq.(8) and (8a) and rearranging, we get, 290
h
h
h
Cot
+
=
1
. . (9) 291
Since in our case,
h
is already known from eq.(5) we rearrange eq.(9) as, 292
( )
h
h h
Cot
+
=
1
, .. (9a) 293
which can be used to obtain the shortening direction at different points over Andaman 294
Islands. 295
From lithospheric configuration of the Benioff zone in the plate subduction 296
regime of Andaman-Sumatra region it has been found that the subducting Indian plate 297
below Andaman dips between 50
o
to 56
o
(Mukhopadhyay et al., 2009). Taking 52.6 as 298
the average angle of dip and replacing it for , all the four components of strain can be 299
easily evaluated. 300
According to Jones et al.(2004), total strain due to simple inclined transpression is given 301
as, 302
D
1
=
|
|
|
\
|
1
0
0 0
0 1
y y zy
y
y xy
. (10) 303
Since in our case all strain components in D
1
would represent infinitesimal strain, 304
subtracting a unit matrix from D
1
we get, 305
d
1
=
|
|
|
\
|
1 0
0 1 0
0 0
1
y zy y
y
xy y
,..(11) 306
where d
1
represents the components of displacement rate or particle velocity in X,Y and 307
Z directions of our reference frame. 308
We next derive from eq.(11) the three dimensional kinematic vorticity number or 309
W
k
(Trusdel, 1953) as, 310
W
k
=
( ) ( )
2 2 2
2
2
2 2
1
1
1 2
xy zy y
y
y
xy zy y
+ +
(
(
\
|
+
+
. .(12) 311
Replacing the different terms in eq.(12) by their corresponding expressions given in (eq.6 312
to eq.9a) we get W
k
at different points over Andaman Islands and its adjoining areas. 313
314
Spatial variation in
y
,
z
,
xy
and
zy
in and around Andaman Islands are shown 315
graphically in Figure-6. It is observed that all the four components gradually increase 316
towards north and south from South Andaman. 317
W
k
calculated at different points over the Andaman Islands varies from 1 near 318
south Andaman to 0.5 near north Andaman and to 0.6 near Little Andaman (Figure-7). 319
The pattern of variation in W
k
implies that south Andaman is subjected to simple shear 320
alone. As pure shear component gradually increases towards both south and north of 321
South Andaman, the W
k
value is modified accordingly. 322
The calculated values of for all the grid points in Andaman region is 323
graphically represented in Figure-8. It shows that is zero around South Andaman and 324
gradually attains higher values towards North and Little Andaman. Since from GPS data 325
South Andaman is known to be moving WSW (Figure-4), we can fix the origin of our 326
coordinate system on the trench axis at such a position where on orienting the X-axis 327
parallel to WSW displacement vector it closely coincides with the trench axis as well 328
(Figure-8). Apparently, such a position is only available near the easterly bend of the 329
subduction zone, north of North Andaman. It implies that this part of the subduction zone 330
around Andaman, is under active transpression at present. 331
The higher angle of convergence in and around Little Andaman is probably related 332
to the presence of a spreading centre in Andaman Sea (Raju et al., 2004), which is 333
forcing the plate to move northwards. This point, however, needs further evaluation 334
because our conclusion is based on observation at a single point only. 335
336
337
5 Discussion 338
The ITRF2005 velocity field obtained from GPS surveys show trench perpendicular 339
azimuth of velocity vector at Little Andaman, which gradually rotates in a clockwise 340
sense and at North Andaman it becomes almost parallel to the trench axis. This indicates 341
non-rigidity of Burma plate (McCaffrey et al., 2000). Similar along-arc change in GPS 342
velocity azimuths at Sumatra subduction zone was observed before the 2004 earthquake 343
(Prawirodirdjo et al., 1997) and at Andaman islands after the earthquake (Paul et al., 344
2007). 345
The azimuthal variations of velocity vectors represent non-uniform deformation 346
pattern in Andaman Islands. Oblique convergence of Indian plate with South East Asia, 347
decomposes the motion into arc-parallel and arc-normal components (Fitch, 1972) 348
resulting in a combination of simple shear (wrench) component and a simultaneous 349
coaxial shortening component leading to transpressional deformation (Fossen and Tikoff, 350
1998). 351
Geometry of the plate margin and degree of obliquity of velocity vectors have major 352
effects on the resultant deformation (Beck, 1991; McCaffrey, 1992). The present analysis 353
shows Indian plate (HYDE GPS station) is moving towards N 35
0
E at a rate of 54 354
mm/yr. Angle of convergence of this velocity varies from 49
0
near Little Andaman to 28
0
355
around North Andaman. Available seismotectonic sections of Andaman region (Ghosh 356
and Mishra, 2008; Mukhopadhyay et al., 2009; Roy, 1992) show the subducting plate is 357
moderately dipping towards east. The oblique convergence between the easterly dipping 358
Indian plate and the Andaman Islands, resolves the horizontal components along strike 359
and dip of the trench. The resultant combination of oblique simple shear and pure shear 360
components is modeled according to the simple inclined transpression model (Jones et 361
al., 2004). 362
In the inclined transpression model the three dimensional kinematic vorticity number 363
varies between 0.5 to 1 implying variability of simple shear and pure shear components. 364
The increment in pure shear towards north is possibly due to, 365
i) east ward bend of the subduction interface; and 366
ii) the already ruptured portion of the trench offers little resistance to slip. 367
Gradual rise in pure shear towards Little Andaman is probably due to the spreading ridge 368
at Andaman sea (Raju et al., 2004). 369
The convergence between Greater India and Asia was approximately normal at Late 370
Paleocene. Afterwards, clockwise rotation of Greater India plate initiated oblique 371
convergence in Sunda arc, which has gradually increased till present (Curray, 2005). The 372
present study also indicates clockwise variation in velocity vectors in Andaman region. 373
Oblique convergence, clockwise rotation and increased convexity of the trench near 374
North Andaman enhance the transpressional activities more towards north. 375
The 2004 mega earthquake has already ruptured ~1500 km (Banerjee et al., 376
2007)between Sumatra and North Andaman (Ghosh and Mishra, 2008). The portion that 377
ruptured during the earthquake does not allow much accumulation of strain. However, 378
beyond the rupture zone the strain continues to accumulate. The process is likely to 379
continue till failure takes place through a major earthquake. 380
In view of the facts mentioned above it may be concluded that north of North 381
Andaman is a potential zone for high magnitude earthquakes. 382
Figure 7 shows the focal plane solution of earthquakes that have occurred in the 383
region (www.globalcmt.org) after 26
th
December, 2004 and before August, 2009. The 384
distribution of earthquake pattern clearly shows concentration of fault activity around 385
North Andaman and Little Andaman. 386
387
6 Acknowledgements 388
The authors express their deep sense of gratitude to Late Sri A. K. Ghosh Roy, 389
Director, Earthquake Geology Division, GSI, ER, who was the key man for formulating 390
the project and establishment of GPS stations at Andaman Islands. They are also indebted 391
to the Deputy Director General, Geological Survey of India, Eastern Region for 392
permitting to publish the work. Authors thankfully acknowledge Dr Bob King, Principal 393
Scientist, MIT for assistance in GAMIT/GLOBK analysis. Andaman administration, 394
Andaman Public Works Department (APWD) and local residents are thankfully 395
acknowledged for extending their support during GPS data collection. 396
397
7 References 398
399
Acharyya, S.K., Ray, K.K. and Roy, D.K., 1989. Tectonostratigraphy and emplacement 400
history of the ophiolite assemblage from Naga Hills and Andaman Island arc, 401
India, Journal Geological Society of India 33, 4-18. 402
Altamimi, Z., Collilieux, X., Legrand, J., Garayt, B. and Boucher, C., 2007. ITRF2005: A 403
new release of the International Terrestrial Reference Frame based on time series 404
of station positions and Earth Orientation Parameters, J. Geophys. Res., 112. 405
Ammon, C., C. Ji, H.-K. Thio, D. Robinson, S. Ni, V. Hjorleifsdottir, H. Kanamori, T. 406
Lay, S. Das, D. Helmberger, G. Ichinose, J. Polet, D. Wald, 2005. Rupture 407
process of the 2004 Sumatra-Andaman earthquake, Science, 308, 1133-1139. 408
Avouac, J.-P. and Tapponnier, P., 1993. Kinematic Model of Active Deformation in 409
Central Asia, Geophys. Res. Lett., 20. 410
Bandopadhyay, P.C. and Ghosh, M., 1998. Fades, petrology and depositional 411
environment of the Tertiary sedimentary rocks around Port Blair, South 412
Andaman- A reply, Journal of Geological Society of India, 52, 612-613. 413
Banerjee, P., Pollitz, F., Nagarajan, B. and Burgmann, R., 2007. Coseismic Slip 414
Distributions of the 26 December 2004 Sumatra-Andaman and 28 March 2005 415
Nias Earthquakes from GPS Static Offsets, Bulletin of the Seismological Society 416
of America, 97(1A), S86-102. 417
Banerjee, P., Pollitz, F.F. and Burgmann, R., 2005. The Size and Duration of the 418
Sumatra-Andaman Earthquake from Far-Field Static Offsets, Science, 308(5729), 419
1769-1772. 420
Beck, M.E., 1991. Coastwise transport reconsidered: lateral displacements in oblique 421
subduction zones and tectonic consequences., Physics of the Earth Planetary 422
Interiors, 68, 1-8. 423
Bock, Y. et al., 2003. Crustal motion in Indonesia from Global Positioning System 424
measurements, J. Geophys. Res., 108. 425
Chlieh, M. et al., 2007. Coseismic Slip and Afterslip of the Great Mw 9.15 Sumatra- 426
Andaman Earthquake of 2004, Bulletin of the Seismological Society of America, 427
97(1A), S152-173. 428
Curray, J.R., 1989. The Sunda Arc: A model for oblique plate convergence, Netherlands 429
Journal of Sea Research, 24(2-3), 131-140. 430
Curray, J.R., 2005. Tectonics and history of the Andaman Sea region, Journal of Asian 431
Earth Sciences, 25(1), 187-232. 432
Curray, J.R., Emmel, F.J. and Moore, D.G., 2002. The Bengal Fan: morphology, 433
geometry, stratigraphy, history and processes, Marine and Petroleum Geology, 434
19(10), 1191-1223. 435
Curray, J.R. and Moore, D.G., 1974. Sedimentary and tectonic processes in Bengal deep- 436
sea fan and geosyncline. In: C.A. Burk, Drake, C. L. (Editor), The Geology of 437
Continental Margins. Springer, New York. 438
Curray, J.R., Moore, D.G., Lawver, L.A., Emmel, F.J., Raitt, R.W., Henry, M., 439
Kieckhefer, R., 1979. Tectonics of the Andaman Sea and Burma., In: Watkins, J., 440
Montadert, L., Dickerson, P.W. (Eds.), Geological and Geophysical Investigations 441
of Continental Margins. American Association Petroleum Geologists, Memoir, 442
29, 189-198. 443
Dasgupta, S. and Mukhopadhyay, M., 1993. Seismicity and plate deformation below the 444
Andaman Arc, northeast Indian Ocean, Teconophysics, 225, 529-542. 445
Dasgupta, S., Mukhopadhyay, M., Bhattacharya, A. and Jana, T.K., 2003. The geometry 446
of the Burmese-Andaman subducting lithosphere, Journal of Seismology, 7, 155- 447
174. 448
DeMets, C., Gordon, R.G., Argus, D.F. and Stein, S., 1990. Current plate motions, 449
Geophysical Journal International, 101(2), 425-478. 450
Dewey, J.F., Holdsworth, R.E., Strachan, R.A., 1998. Transpression and transtension 451
zones, In: Holdsworth, R.E., Strachan, R.A., Dewey, J.F. (Eds.), Continental 452
Transpressional and Transtensional Tectonics. Special publication of the 453
Geological Society, London, 135, 1-14. 454
Dong, D.N., Herring, T.A. and King, R.W., 1998. Estimating regional deformation from 455
a combination of space and terrestrial geodetic data, Journal of Geodesy, 72, 200- 456
214. 457
Dutton, B.J., 1997. Finite strains in transpression zones with no boundary slip, Journal of 458
structural geology, 19(9), 1189-1200. 459
Fitch, T.J., 1972. Plate Convergence, Transcurrent Faults, and Internal Deformation 460
Adjacent to Southeast Asia and the Western Pacific, J. Geophys. Res., 77 461
Fossen, H. and Tikoff, B., 1993. The deformation matrix for simultaneous simple 462
shearing, pure shearing and volume change, and its application to transpression- 463
transtension tectonics, Journal of structural geology, 15(3-5), 413-422. 464
Fossen, H. and Tikoff, B., 1998. Extended models of transpression and transtension, and 465
application to tectonic settings, In: Holdsworth, R.E., Strachan, R.A., Dewey, J.F. 466
(Eds.), Continental Transpressional and Transtensional Tectonics. Special 467
publication of the Geological Society, London, 135, 15-33. 468
Ghosh, D. and Mishra, O.P., 2008. The 2004 Sumatra-Andaman earthquake sequence and 469
its implications for seismic coupling: Future vulnerability, Indian Minerals, 61(3- 470
4) & 62(1-4), 93-112. 471
Ghosh Roy, A.K., Bardhan, S., Jana, P. and Basir, S.R., 2007. Macroseismic survey in 472
Andaman and Nicobar island in the aftermath of the great earthquake of 26 473
December, 2004., Geological Survey of India; Special publication, 89, 17-50. 474
Ghosh, S.K., 2001. Types of transpressional and transtensional deformation, In: Koyi, 475
H.A., Mancktelow, N.S. (Eds.), Tectonic modelling: A volume in Honour of Hans 476
Ramberg. Geological Society of America Memoir, 193, 1-20. 477
Goodwin, L.B. and Tikoff, B., 2002. Competency contrast, kinematics, and the 478
development of foliations and lineations in the crust, Journal of structural 479
geology, 24(6-7), 1065-1085. 480
Herring, T.A., King, R.W. and McClusky, S.C., 2009. Gamit reference manual; Release 481
10.3, Massachussetts Institute of Technology. 482
Horsman, E. and Tikoff, B., 2007. Constraints on deformation path from finite strain 483
gradients, Journal of structural geology, 29(2), 256-272. 484
Jaeger, J.C., 1964. Elasticity, Fracture and Flow. Methuen, London, pp. 212. 485
Jiang, D. and White, J.C., 1995. Kinematics of rock flow and the interpretation of 486
geological structures, with particular reference to shear zones, Journal of 487
structural geology, 17(9), 1249-1265. 488
Jones, R.R., Holdsworth, R.E., Clegg, P., McCaffrey, K. and Tavarnelli, E., 2004. 489
Inclined transpression, Journal of structural geology, 26(8), 1531-1548. 490
Lay, T. et al., 2005. The Great Sumatra-Andaman Earthquake of 26 December 2004, 491
Science, 308(5725), 1127-1133. 492
Lin, S., Jiang, D., Williams, P.F., 1998. Transpression (or transtension) zones of triclinic 493
symmetry: natural example and theoretical modelling, In: Holdsworth, R.E., 494
Strachan, R.A., Dewey, J.F. (Eds.), Continental Transpressional and 495
Transtensional Tectonics. Special publication of the Geological Society, London, 496
135, 41-57. 497
Mao, A., Harrison, C.G.A. and Dixon, T.H., 1999. Noise in GPS coordinate time series, 498
J. Geophys. Res., 104. 499
McCaffrey, R., 1991. Slip vectors and stretching of the Sumatran fore arc, Geology, 500
19(9), 881-884. 501
McCaffrey, R., 1992. Oblique plate convergence, slip vectors and forearc deformation, 502
Journal of Geophysical Research, 97, 8905-8915. 503
McCaffrey, R. et al., 2000. Strain partitioning during oblique plate convergence in 504
northern Sumatra: Geodetic and seismologic constraints and numerical modeling, 505
J. Geophys. Res., 105. 506
McCaffrey, R. et al., 2007. Strain partitioning during oblique plate convergence in 507
northern Sumatra: Geodetic and seismologic constraints and numerical modeling, 508
J. Geophys. Res., 105 (2000). 509
Michel, G.W. et al., 2001. Crustal motion and block behaviour in SE-Asia from GPS 510
measurements, Earth and Planetary Science Letters, 187(3-4), 239-244. 511
Molnar, P., 1992. Brace-Goetze strength -profiles, the partitioning of strike-slip and 512
thrust faulting at zones of oblique convergence, and the stress-heat flow paradox 513
of the San Andreas fault., In: Evans, B & Wong, T.-F.(eds). Fault Mechanics and 514
Transport Properties of Rocks. Academic Press, New York, 435-459. 515
Mukhopadhyay, B., Acharyya, A., Bhattacharya, A., Dasgupta, S. and Sengupta, S.R., 516
2009. Revisiting the Andaman subduction lithosphere following the 26 December 517
2004 Sumatra earthquake, Indian Journal of Geosciences, 63(1), 1-10. 518
Paul, J., Lowry, A.R., Bilham, R., Sen, S. and Smalley, R., Jr., 2007. Postseismic 519
deformation of the Andaman Islands following the 26 December, 2004 Great 520
SumatraAndaman earthquake, Geophys. Res. Lett., 34. 521
Philip, E. and Peter, M., 1997. The field of crustal velocity in Asia calculated from 522
Quaternary rates of slip on faults, Geophysical Journal International, 130(3), 551- 523
582. 524
Prawirodirdjo, L. et al., 1997. Geodetic Observations of Interseismic Strain Segmentation 525
at the Sumatra Subduction Zone, Geophys. Res. Lett., 24. 526
Raju, K.A.K., Ramprasad, T., Rao, P.S., Rao, B.R. and Varghese, J., 2004. New insights 527
into the tectonic evolution of the Andaman basin, northeast Indian Ocean, Earth 528
and Planetary Science Letters, 221, 145-162. 529
Ray, K.K., 1982. A review of the geology of Andaman and Nicobar islands. . Geological 530
Survey of India. Miscellaneous Publication, pp. 110-125. 531
Ray, S.K. and Acharyya, A., 2007. 26 December 2004 earthquake: Coseismic vertical 532
ground movement in the Andaman Islands, Geological Survey of India, Special 533
publication, 89, 63-81. 534
Robin, P.-Y.F. and Cruden, A.R., 1994. Strain and vorticity patterns in ideally ductile 535
transpression zones, Journal of structural geology, 16(4), 447-466. 536
Roy, S.K., 1992. Accretionary prism in Andaman Forearc, Geological Survey of India, 537
Special publication, 29, 273-278. 538
Sanderson, D.J. and Marchini, W.R.D., 1984. Transpression, Journal of structural 539
geology, 6(5), 449-458. 540
Sella, G.F., Dixon, T.H. and Mao, A., 2002. REVEL: A model for Recent plate velocities 541
from space geodesy, J. Geophys. Res., 107. 542
Sengupta, S., Ray, K.K. and Acharyya, S.K., 1990. Nature of ophiolite occurrences along 543
the eastern margin of the Indian plate and their tectonic significance. , Geology, 544
18, 439-442. 545
Shen, Z.-K. and Jackson, D.D., 2000. Optimal estimation of geodetic strain rates from 546
GPS data, EOS, Transactions of the Americal Geophysical Union, 81(19), S406. 547
Shen, Z.-K., Jackson, D.D. and Ge, B.X., 1996. Crustal deformation across and beyond 548
the Los Angeles basin from geodetic measurements, J. Geophys. Res., 101. 549
Simons, W.J.F. et al., 1999. Observing Plate Motions in S.E. Asia: Geodetic Results of 550
the GEODYSSEA Project, Geophys. Res. Lett., 26. 551
Som, S.K., Shivgotra, Vijay, Saha, Ashim, 2009. Coral microatoll as geodetic tool in 552
North Andaman and Little Andaman, India, Journal of Earth System Science, 553
118(2), 157-162. 554
Subarya, C. et al., 2006. Plate-boundary deformation associated with the great 555
SumatraAndaman earthquake, Nature, 440(7080), 46-51. 556
Teza, G., Pesci, A. and Galgaro, A., 2008. Grid_strain and grid_strain3: Software 557
packages for strain field computation in 2D and 3D environments, Computers & 558
Geosciences, 34(9), 1142-1153. 559
Tikoff, B. and Fossen, H., 1993. Simultaneous pure and simple shear: the unifying 560
deformation matrix, Tectonophysics, 217(3-4), 267-283. 561
Vigny, C. et al., 2005. Insight into the 2004 Sumatra-Andaman earthquake from GPS 562
measurements in southeast Asia, Nature, 436(7048), 201-206. 563
564
Figure captions 565
566
Figure: 1 Velocity of global reference stations at ITRF2005 reference frame derived
from the present analysis.
Figure: 2 ITRF2005 velocity of Andaman sites. Base map shows Andaman geology
compiled after Ghosh Roy et al.,(Ghosh Roy et al., 2007)
Figure: 3 Velocity relative to Eurasia reference frame. Inset for stations at Andaman .
Figure: 4 Velocity relative to India reference frame.
Figure:5 Strain components from GPS velocity vectors. Blue lines for
max
and red
lines for
min
. Bold lines for trench and active faults and dotted lines for
inactive faults.
Figure: 6 Modeled infinitesimal strain components for inclined transpression at
Andaman. (a) Longitudinal strain along y-direction, (b) Longitudinal strain
along z-direction, (c) shear strain in xy-plane, (d) shear strain in yz-plane
Figure: 7 Variation of 3-D infinitesimal kinematical vorticity number with focal plane
solution of earthquakes occurred between 2004 mega-earthquake and
August, 2009.
Figure: 8 Variation of angle in the horizontal plane between the subduction zone
boundary and the direction of overall shortening with defined coordinate
system. Arrow shows India fixed velocity azimuth of BDNB station at South
Andaman.
Table captions 567
568
Table-1 Details of GPS campaign stations used in analysis
Table-2 Station coordinates, observed velocities, residuals and uncertainties
569