Rapid microsatellite marker development for African mahogany ( Khaya senegalensis, Meliaceae) using next-generation sequencing and assessment of its intra-specific genetic diversity.

Khaya senegalensis (African mahogany or dry-zone mahogany) is a high-value hardwood timber species with great potential for forest plantations in northern Australia. The species is distributed across the sub-Saharan belt from Senegal to Sudan and Uganda. Because of heavy exploitation and constraints on natural regeneration and sustainable planting, it is now classified as a vulnerable species. Here, we describe the development of microsatellite markers for K. senegalensis using next-generation sequencing to assess its intra-specific diversity across its natural range, which is a key for successful breeding programs and effective conservation management of the species. Next-generation sequencing yielded 93943 sequences with an average read length of 234bp. The assembled sequences contained 1030 simple sequence repeats, with primers designed for 522 microsatellite loci. Twenty-one microsatellite loci were tested with 11 showing reliable amplification and polymorphism in K. senegalensis. The 11 novel microsatellites, together with one previously published, were used to assess 73 accessions belonging to the Australian K. senegalensis domestication program, sampled from across the natural range of the species. STRUCTURE analysis shows two major clusters, one comprising mainly accessions from west Africa (Senegal to Benin) and the second based in the far eastern limits of the range in Sudan and Uganda. Higher levels of genetic diversity were found in material from western Africa. This suggests that new seed collections from this region may yield more diverse genotypes than those originating from Sudan and Uganda in eastern Africa.

Tropical dryland agroforestry on clay soils: : Analysis of systems based on Acacia senegal in the Blue Nile region, Sudan

Acacia senegal, the gum arabic producing tree, is the most important component in traditional dryland agroforestry systems in the Blue Nile region, Sudan. The aim of the present study was to provide new knowledge on the potential use of A. senegal in dryland agroforestry systems on clay soils, as well as information on tree/crop interaction, and on silvicultural and management tools, with consideration on system productivity, nutrient cycling and sustainability. Moreover, the aim was also to clarify the intra-specific variation in the performance of A. senegal and, specifically, the adaptation of trees of different origin to the clay soils of the Blue Nile region. In agroforestry systems established at the beginning of the study, tree and crop growth, water use, gum and crop yields, nutrient cycling and system performance were investigated for a period of four years (1999 to 2002). Trees were grown at 5 x 5 m and 10 x 10 m spacing alone or in mixture with sorghum or sesame; crops were also grown in sole culture. The symbiotic biological N2 fixation by A. senegal was estimated using the 15N natural abundance (δ15N) procedure in eight provenances collected from different environments and soil types of the gum arabic belt and grown in clay soil in the Blue Nile region. Balanites aegyptiaca (a non-legume) was used as a non-N-fixing reference tree species, so as to allow 15N-based estimates of the proportion of the nitrogen in trees derived from the atmosphere. In the planted acacia trees, measurements were made on shoot growth, water-use efficiency (as assessed by the δ13C method) and (starting from the third year) gum production. Carbon isotope ratios were obtained from the leaves and branch wood samples. The agroforestry system design caused no statistically significant variation in water use, but the variation was highly significant between years, and the highest water use occurred in the years with high rainfall. No statistically significant differences were found in sorghum or sesame yields when intercropping and sole crop systems were compared (yield averages were 1.54 and 1.54 ha-1 for sorghum and 0.36 and 0.42 t ha-1 for sesame in the intercropped and mono-crop plots, respectively). Thus, at an early stage of agroforestry system management, A. senegal had no detrimental effect on crop yield, but the pattern of resource capture by trees and crops may change as the system matures. Intercropping resulted in taller trees and larger basal and crown diameters as compared to the development of sole trees. It also resulted in a higher land equivalent ratio. When gum yields were analysed it was found that a significant positive relationship existed between the second gum picking and the total gum yield. The second gum picking seems to be a decisive factor in gum production and could be used as an indicator for the total gum yield in a particular year. In trees, the concentrations of N and P were higher in leaves and roots, whereas the levels of K were higher in stems, branches and roots. Soil organic matter, N, P and K contents were highest in the upper soil stratum. There was some indication that the P content slightly increased in the topsoil as the agroforestry plantations aged. At a stocking of 400 trees ha-1 (5 x 5 m spacing), A. senegal accumulated in the biomass a total of 18, 1.21, 7.8 and 972 kg ha-1of N, P, K and OC, respectively. Trees contributed ca. 217 and 1500 kg ha-1 of K and OC, respectively, to the top 25-cm of soil over the first four years of intercropping. Acacia provenances of clay plain origin showed considerable variation in seed weight. They also had the lowest average seed weight as compared to the sandy soil (western) provenances. At the experimental site in the clay soil region, the clay provenances were distinctly superior to the sand provenances in all traits studied but especially in basal diameter and crown width, thus reflecting their adaptation to the environment. Values of δ13C, indicating water use efficiency, were higher in the sand soil group as compared to the clay one, both in leaves and in branch wood. This suggests that the sand provenances (with an average value of -28.07 ) displayed conservative water use and high drought tolerance. Of the clay provenances, the local one (Bout) displayed a highly negative (-29.31 ) value, which indicates less conservative water use that resulted in high productivity at this particular clay-soil site. Water use thus appeared to correspond to the environmental conditions prevailing at the original locations for these provenances. Results suggest that A. senegal provenances from the clay part of the gum belt are adapted for a faster growth rate and higher biomass and gum productivity as compared to provenances from sand regions. A strong negative relationship was found between the per-tree gum yield and water use efficiency, as indicated by δ13C. The differences in water use and gum production were greater among provenance groups than within them, suggesting that selection among rather than within provenances would result in distinct genetic gain in gum yield. The relative δ15N values ( ) were higher in B. aegyptiaca than in the N2-fixing acacia provenances. The amount of Ndfa increased significantly with age in all provenances, indicating that A. senegal is a potentially efficient nitrogen fixer and has an important role in t agroforestry development. The total above-ground contribution of fixed N to foliage growth in 4-year-old A. senegal trees was highest in the Rahad sand-soil provenance (46.7 kg N ha-1) and lowest in the Mazmoom clay-soil provenance (28.7 kg N ha-1). This study represents the first use of the δ15N method for estimating the N input by A. senegal in the gum belt of Sudan. Key words: Acacia senegal, agroforestry, clay plain, δ13C, δ15N, gum arabic, nutrient cycling, Ndfa, Sorghum bicolor, Sesamum indicum

Multistation study of nighttime scintillations in low latitudes: Evidence of control by equatorial F region irregularities

VHF nighttime scintillations, recorded during a high solar activity period at a meridian chain of stations covering a magnetic latitude belt of 3°–21°N (420 km subionospheric points) are analyzed to investigate the influence of equatorial spread F irregularities on the occurrence of scintillation at latitudes away from the equator. Observations show that saturated amplitude scintillations start abruptly about one and a half hours after ground sunset and their onset is almost simultaneous at stations whose subionospheric points are within 12°N latitude of the magnetic equator, but is delayed at a station whose subionospheric point is at 21°N magnetic latitude by 15 min to 4 hours. In addition, the occurrence of postsunset scintillations at all the stations is found to be conditional on their prior occurrence at the equatorial station. If no postsunset scintillation activity is seen at the equatorial station, no scintillations are seen at other stations also. The occurrence of scintillations is explained as caused by rising plasma bubbles and associated irregularities over the magnetic equator and the subsequent mapping of these irregularities down the magnetic field lines to the F region of higher latitudes through some instantaneous mechanism; and hence an equatorial control is established on the generation of postsunset scintillation-producing irregularities in the entire low-latitude belt.

Odour, dust and non-methane volatile organic-compound emissions from tunnel-ventilated layer-chicken sheds: a case study of two farms

An observational study was undertaken to measure odour and dust (PM10 and PM2.5) emission rates and identify non-methane volatile organic compounds (NMVOCs) and odorants in the exhaust air from two tunnel-ventilated layer-chicken sheds that were configured with multi-tiered cages and manure belts. The study sites were located in south-eastern Queensland and the West Gippsland region of Victoria, Australia. Samples were collected in summer and winter on sequential days across the manure-belt cleaning cycle. Odour emissions ranged from 58 to 512 ou/s per 1000 birds (0.03-0.27 ou/s.kg) and dust emission rates ranged 0.014-0.184 mg/s per 1000 birds for PM10 and 0.001-0.190 mg/s per 1000 birds for PM2.5. Twenty NMVOCs were identified, including three that were also identified as odorants using thermal desorption-gas chromatography-mass spectrometry/olfactometry analysis. Odour emission rates were observed to vary with the amount of manure accumulation on the manure belts, being lowest 2-4 days after removing manure. Odour emission rates were also observed to vary with diurnal and seasonal changes in ventilation rate. Dust emissions were observed to increase with ventilation rate but not with manure accumulation. Some NMVOCs were identified at both farms and in different seasons whereas others were observed only at one farm or in one season, indicating that odorant composition was influenced by farm-specific practices and season.

Spatial impact of projected changes in rainfall and temperature on wheat yields in Australia

Climate projections over the next two to four decades indicate that most of Australia’s wheat-belt is likely to become warmer and drier. Here we used a shire scale, dynamic stress-index model that accounts for the impacts of rainfall and temperature on wheat yield, and a range of climate change projections from global circulation models to spatially estimate yield changes assuming no adaptation and no CO2 fertilisation effects. We modelled five scenarios, a baseline climate (climatology, 1901–2007), and two emission scenarios (“low” and “high” CO2) for two time horizons, namely 2020 and 2050. The potential benefits from CO2 fertilisation were analysed separately using a point level functional simulation model. Irrespective of the emissions scenario, the 2020 projection showed negligible changes in the modelled yield relative to baseline climate, both using the shire or functional point scale models. For the 2050-high emissions scenario, changes in modelled yield relative to the baseline ranged from −5 % to +6 % across most of Western Australia, parts of Victoria and southern New South Wales, and from −5 to −30 % in northern NSW, Queensland and the drier environments of Victoria, South Australia and in-land Western Australia. Taking into account CO2 fertilisation effects across a North–south transect through eastern Australia cancelled most of the yield reductions associated with increased temperatures and reduced rainfall by 2020, and attenuated the expected yield reductions by 2050.

Inventories and significance of the genetic resources of an African mahogany species (Khaya senegalensis (Desr.) A. Juss.) assembled and further developed in Australia.

The forest tree species Khaya senegalensis (Desr.) A. Juss. occurs in a belt across 20 African countries from Senegal-Guinea to Sudan-Uganda where it is a highly important resource. However, it is listed as Vulnerable (IUCN 2015-3). Since introduction in northern Australia around 1959, the species has been planted widely, yielding high-value products. The total area of plantations of the species in Australia exceeds 15,000 ha, mostly planted in the Northern Territory since 2006, and includes substantial areas across 60-70 woodlots and industrial plantations established in north-eastern Queensland since the early-1990s and during 2005-2007 respectively. Collaborative conservation and tree improvement by governments began in the Northern Territory and Queensland in 2001 based on provenance and other trials of the 1960s-1970s. This work has developed a broad base of germplasm in clonal seed orchards, hedge gardens and trials (clone and progeny). Several of the trials were established collaboratively on private land. Since the mid-2000s, commercial growers have introduced large numbers of provenance-bulk and individual-tree seedlots to establish industrial plantations and trials, several of the latter in collaboration with the Queensland Government. Provenance bulks (>140) and families (>400) from 17 African countries are established in Australia, considered the largest genetic base of the species in a single country outside Africa. Recently the annual rate of industrial planting of the species in Australia has declined, and R&D has been suspended by governments and reduced by the private sector. However, new commercial plantings in the Northern Territory and Queensland are proposed. In domesticating a species, the strategic importance of a broad genetic base is well known. The wide range of first- and advanced-generation germplasm of the species established in northern Australia and documented in this paper provides a sound basis for further domestication and industrial plantation and woodlot expansion, when investment conditions are favourable

Pseudo-time marching schemes for inverse problems in structural health assessment and medical imaging

We propose a self-regularized pseudo-time marching strategy for ill-posed, nonlinear inverse problems involving recovery of system parameters given partial and noisy measurements of system response. While various regularized Newton methods are popularly employed to solve these problems, resulting solutions are known to sensitively depend upon the noise intensity in the data and on regularization parameters, an optimal choice for which remains a tricky issue. Through limited numerical experiments on a couple of parameter re-construction problems, one involving the identification of a truss bridge and the other related to imaging soft-tissue organs for early detection of cancer, we demonstrate the superior features of the pseudo-time marching schemes.

A Kalman filter based strategy for linear structural system identification based on multiple static and dynamic test data

The problem of identification of stiffness, mass and damping properties of linear structural systems, based on multiple sets of measurement data originating from static and dynamic tests is considered. A strategy, within the framework of Kalman filter based dynamic state estimation, is proposed to tackle this problem. The static tests consists of measurement of response of the structure to slowly moving loads, and to static loads whose magnitude are varied incrementally; the dynamic tests involve measurement of a few elements of the frequency response function (FRF) matrix. These measurements are taken to be contaminated by additive Gaussian noise. An artificial independent variable τ, that simultaneously parameterizes the point of application of the moving load, the magnitude of the incrementally varied static load and the driving frequency in the FRFs, is introduced. The state vector is taken to consist of system parameters to be identified. The fact that these parameters are independent of the variable τ is taken to constitute the set of ‘process’ equations. The measurement equations are derived based on the mechanics of the problem and, quantities, such as displacements and/or strains, are taken to be measured. A recursive algorithm that employs a linearization strategy based on Neumann’s expansion of structural static and dynamic stiffness matrices, and, which provides posterior estimates of the mean and covariance of the unknown system parameters, is developed. The satisfactory performance of the proposed approach is illustrated by considering the problem of the identification of the dynamic properties of an inhomogeneous beam and the axial rigidities of members of a truss structure.

Quasi-periodic, global oscillations in sea level pressure on intraseasonal timescales

The sea level pressure (SLP) variability in 30-60 day intraseasonal timescales is investigated using 25 years of reanalysis data addressing two issues. The first concerns the non-zero zonal mean component of SLP near the equator and its meridional connections, and the second concerns the fast eastward propagation (EP) speed of SLP compared to that of zonal wind. It is shown that the entire globe resonates with high amplitude wave activity during some periods which may last for few to several months, followed by lull periods of varying duration. SLP variations in the tropical belt are highly coherent from 25A degrees S to 25A degrees N, uncorrelated with variations in mid latitudes and again significantly correlated but with opposite phase around 60A degrees S and 65A degrees N. Near the equator (8A degrees S-8A degrees N), the zonal mean contributes significantly to the total variance in SLP, and after its removal, SLP shows a dominant zonal wavenumber one structure having a periodicity of 40 days and EP speeds comparable to that of zonal winds in the Indian Ocean. SLP from many of the atmospheric and coupled general circulation models show similar behaviour in the meridional direction although their propagation characteristics in the tropical belt differ widely.

Tactile processing in human somatosensory and auditory cortices

Tactile sensation plays an important role in everyday life. While the somatosensory system has been studied extensively, the majority of information has come from studies using animal models. Recent development of high-resolution anatomical and functional imaging techniques has enabled the non-invasive study of human somatosensory cortex and thalamus. This thesis provides new insights into the functional organization of the human brain areas involved in tactile processing using magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI). The thesis also demonstrates certain optimizations of MEG and fMRI methods. Tactile digit stimulation elicited stimulus-specific responses in a number of brain areas. Contralateral activation was observed in somatosensory thalamus (Study II), primary somatosensory cortex (SI; I, III, IV), and post-auditory belt area (III). Bilateral activation was observed in secondary somatosensory cortex (SII; II, III, IV). Ipsilateral activation was found in the post-central gyrus (area 2 of SI cortex; IV). In addition, phasic deactivation was observed within ipsilateral SI cortex and bilateral primary motor cortex (IV). Detailed investigation of the tactile responses demonstrated that the arrangement of distal-proximal finger representations in area 3b of SI in humans is similar to that found in monkeys (I). An optimized MEG approach was sufficient to resolve such fine detail in functional organization. The SII region appeared to contain double representations for fingers and toes (II). The detection of activations in the SII region and thalamus improved at the individual and group levels when cardiac-gated fMRI was used (II). Better detection of body part representations at the individual level is an important improvement, because identification of individual representations is crucial for studying brain plasticity in somatosensory areas. The posterior auditory belt area demonstrated responses to both auditory and tactile stimuli (III), implicating this area as a physiological substrate for the auditory-tactile interaction observed in earlier psychophysical studies. Comparison of different smoothing parameters (III) demonstrated that proper evaluation of co-activation should be based on individual subject analysis with minimal or no smoothing. Tactile input consistently influenced area 3b of the human ipsilateral SI cortex (IV). The observed phasic negative fMRI response is proposed to result from interhemispheric inhibition via trans-callosal connections. This thesis contributes to a growing body of human data suggesting that processing of tactile stimuli involves multiple brain areas, with different spatial patterns of cortical activation for different stimuli.

The efficiency of bounding vampires

Vampire bats, Desmodus rotundus, must maximize their feeding cycle of one blood meal per day by being efficient in the stalking and acquisition of their food. Riskin and Hermanson documented the running gait of the common vampire bat and observed they were efficient at running speeds, using longer stride lengths and thus decreased stride frequency. We obtained preliminary data on gait maintained for up to 10 minutes on a moving treadmill belt at speeds ranging from 0.23 to 0.74 m/s, which spanned a range from walking to running gaits. Bats tended to transition between gaits at about 0.40 m/s. Fourteen bats were studied and included four that were able to walk or run for 10 minutes. There was no significant change in either stride duration or frequency associated with an increase in speed. We estimated O2 consumption and CO2 production both before and 5 minutes after exercise, and found that O2 consumption increased 1 minute and 5 minutes after exercise. CO2 levels increased significantly 1 minute after exercise, but tended back towards pre-exercise level 5 minutes after exercise. Two bats were tested for blood O2, CO2 and pH levels. Interestingly, pH levels fell from 7.3 to about 7.0, indicating lactate accumulation.

Simultaneous material selection and geometry design of statically determinate trusses using continuous optimization

In this work, we explore simultaneous geometry design and material selection for statically determinate trusses by posing it as a continuous optimization problem. The underlying principles of our approach are structural optimization and Ashby’s procedure for material selection from a database. For simplicity and ease of initial implementation, only static loads are considered in this work with the intent of maximum stiffness, minimum weight/cost, and safety against failure. Safety of tensile and compression members in the truss is treated differently to prevent yield and buckling failures, respectively. Geometry variables such as lengths and orientations of members are taken to be the design variables in an assumed layout. Areas of cross-section of the members are determined to satisfy the failure constraints in each member. Along the lines of Ashby’s material indices, a new design index is derived for trusses. The design index helps in choosing the most suitable material for any geometry of the truss. Using the design index, both the design space and the material database are searched simultaneously using gradient-based optimization algorithms. The important feature of our approach is that the formulated optimization problem is continuous, although the material selection from a database is an inherently discrete problem. A few illustrative examples are included. It is observed that the method is capable of determining the optimal topology in addition to optimal geometry when the assumed layout contains more links than are necessary for optimality.

Asteroid identification using statistical orbital inversion methods

An efficient and statistically robust solution for the identification of asteroids among numerous sets of astrometry is presented. In particular, numerical methods have been developed for the short-term identification of asteroids at discovery, and for the long-term identification of scarcely observed asteroids over apparitions, a task which has been lacking a robust method until now. The methods are based on the solid foundation of statistical orbital inversion properly taking into account the observational uncertainties, which allows for the detection of practically all correct identifications. Through the use of dimensionality-reduction techniques and efficient data structures, the exact methods have a loglinear, that is, O(nlog(n)), computational complexity, where n is the number of included observation sets. The methods developed are thus suitable for future large-scale surveys which anticipate a substantial increase in the astrometric data rate. Due to the discontinuous nature of asteroid astrometry, separate sets of astrometry must be linked to a common asteroid from the very first discovery detections onwards. The reason for the discontinuity in the observed positions is the rotation of the observer with the Earth as well as the motion of the asteroid and the observer about the Sun. Therefore, the aim of identification is to find a set of orbital elements that reproduce the observed positions with residuals similar to the inevitable observational uncertainty. Unless the astrometric observation sets are linked, the corresponding asteroid is eventually lost as the uncertainty of the predicted positions grows too large to allow successful follow-up. Whereas the presented identification theory and the numerical comparison algorithm are generally applicable, that is, also in fields other than astronomy (e.g., in the identification of space debris), the numerical methods developed for asteroid identification can immediately be applied to all objects on heliocentric orbits with negligible effects due to non-gravitational forces in the time frame of the analysis. The methods developed have been successfully applied to various identification problems. Simulations have shown that the methods developed are able to find virtually all correct linkages despite challenges such as numerous scarce observation sets, astrometric uncertainty, numerous objects confined to a limited region on the celestial sphere, long linking intervals, and substantial parallaxes. Tens of previously unknown main-belt asteroids have been identified with the short-term method in a preliminary study to locate asteroids among numerous unidentified sets of single-night astrometry of moving objects, and scarce astrometry obtained nearly simultaneously with Earth-based and space-based telescopes has been successfully linked despite a substantial parallax. Using the long-term method, thousands of realistic 3-linkages typically spanning several apparitions have so far been found among designated observation sets each spanning less than 48 hours.

Seismic Tomography and Earthquake Mechanism Beneath the Central Fennoscandian Shield

The aim of this thesis was to study the seismic tomography structure of the earth s crust together with earthquake distribution and mechanism beneath the central Fennoscandian Shield, mainly in southern and central Finland. The earthquake foci and some fault plane solutions are correlated with 3-D images of the velocity tomography. The results are discussed in relation to the stress field of the Shield and with other geophysical, e.g. geomagnetic, gravimetric, tectonic, and anisotropy studies of the Shield. The earthquake data of the Fennoscandian Shield has been extracted from the Nordic earthquake parameter data base which was founded at the time of inception of the earthquake catalogue for northern Europe. Eight earlier earthquake source mechanisms are included in a pilot study on creating a novel technique for calculating an earthquake fault plane solution. Altogether, eleven source mechanisms of shallow, weak earthquakes are related in the 3-D tomography model to trace stresses of the crust in southern and central Finland. The earthquakes in the eastern part of the Fennoscandian Shield represent low-active, intraplate seismicity. Earthquake mechanisms with NW-SE oriented horizontal compression confirm that the dominant stress field originates from the ridge-push force in the North Atlantic Ocean. Earthquakes accumulate in coastal areas, in intersections of tectonic lineaments, in main fault zones or are bordered by fault lines. The majority of Fennoscandian earthquakes concentrate on the south-western Shield in southern Norway and Sweden. Onwards, epicentres spread via the ridge of the Shield along the west-coast of the Gulf of Bothnia northwards along the Tornio River - Finnmark fault system to the Barents Sea, and branch out north-eastwards via the Kuusamo region to the White Sea Kola Peninsula faults. The local seismic tomographic method was applied to find the terrane distribution within the central parts of the Shield the Svecofennian Orogen. From 300 local explosions a total of 19765 crustal Pg- and Sg-wave arrival times were inverted to create independent 3-D Vp and Vs tomographic models, from which the Vp/Vs ratio was calculated. The 3-D structure of the crust is presented as a P-wave and for the first time as an S-wave velocity model, and also as a Vp/Vs-ratio model of the SVEKALAPKO area that covers 700x800 km2 in southern and central Finland. Also, some P-wave Moho-reflection data was interpolated to image the relief of the crust-mantle boundary (i.e. Moho). In the tomography model, the seismic velocities vary smoothly. The lateral variations are larger for Vp (dVp =0.7 km/s) than for Vs (dVs =0.4 km/s). The Vp/Vs ratio varies spatially more distinctly than P- and S-wave velocities, usually from 1.70 to 1.74 in the upper crust and from 1.72 to 1.78 in the lower crust. Schist belts and their continuations at depth are associated with lower velocities and lower Vp/Vs ratios than in the granitoid areas. The tomography modelling suggests that the Svecofennian Orogen was accreted from crustal blocks ranging in size from 100x100 km2 to 200x200 km2 in cross-sectional area. The intervening sedimentary belts have ca. 0.2 km/s lower P- and S-wave velocities and ca. 0.04 lower Vp/Vs ratios. Thus, the tomographic model supports the concept that the thick Svecofennian crust was accreted from several crustal terranes, some hidden, and that the crust was later modified by intra- and underplating. In conclusion, as a novel approach the earthquake focal mechanism and focal depth distribution is discussed in relation to the 3-D tomography model. The schist belts and the transformation zones between the high- and low-velocity anomaly blocks are characterized by deeper earthquakes than the granitoid areas where shallow events dominate. Although only a few focal mechanisms were solved for southern Finland, there is a trend towards strike-slip and oblique strike-slip movements inside schist areas. The normal dip-slip type earthquakes are typical in the seismically active Kuusamo district in the NE edge of the SVEKALAPKO area, where the Archean crust is ca. 15-20 km thinner than the Proterozoic Svecofennian crust. Two near vertical dip-slip mechanism earthquakes occurred in the NE-SW junction between the Central Finland Granitoid Complex and the Vyborg rapakivi batholith, where high Vp/Vs-ratio deep-set intrusion splits the southern Finland schist belt into two parts in the tomography model.

The Palu Metamorphic Complex, NW Sulawesi, Indonesia: Origin and evolution of a young metamorphic terrane with links to Gondwana and Sundaland

The Palu Metamorphic Complex (PMC) is exposed in a late Cenozoic orogenic belt in NW Sulawesi, Indonesia. It is a composite terrane comprising a gneiss unit of Gondwana origin, a schist unit composed of meta-sediments deposited along the SE Sundaland margin in the Late Cretaceous and Early Tertiary, and one or more slivers of amphibolite with oceanic crust characteristics. The gneiss unit forms part of the West Sulawesi block underlying the northern and central sections of the Western Sulawesi Province. The presence of Late Triassic granitoids and recycled Proterozoic zircons in this unit combined with its isotopic signature suggests that the West Sulawesi block has its origin in the New Guinea margin from which it rifted in the late Mesozoic. It docked with Sundaland sometime during the Late Cretaceous. U–Th–Pb dating results for monazite suggest that another continental fragment may have collided with the Sundaland margin in the earliest Miocene. High-pressure (HP) and ultrahigh-pressure (UHP) rocks (granulite, peridotite, eclogite) are found as tectonic slices within the PMC, mostly along the Palu–Koro Fault Zone, a major strike-slip fault that cuts the complex. Mineralogical and textural features suggest that some of these rocks resided at depths of 60–120 km during a part of their histories. Thermochronological data (U–Th–Pb zircon and 40Ar/39Ar) from the metamorphic rocks indicate a latest Miocene to mid-Pliocene metamorphic event, which was accompanied by widespread granitoid magmatism and took place in an extensional tectonic setting. It caused recrystallization of, and new overgrowths on, pre-existing zircon crystals, and produced andalusite–cordierite–sillimanite–staurolite assemblages in pelitic protoliths, indicating HT–LP (Buchan-type) metamorphism. The PMC was exhumed as a core complex at moderate rates (c. 0.7–1.0 mm/yr) accompanied by rapid cooling in the Plio-Pleistocene. Some of the UHP rocks were transported to the surface at significantly higher rates (⩾16 mm/yr). The results of our study do not support recent plate tectonic reconstructions that propose a NW Australia margin origin for the West Sulawesi block (e.g. Hall et al., 2009).