138 results found
Collins GS, Lynch E, McAdam R, et al., 2017, A numerical assessment of simple airblast models of impact airbursts, Meteoritics & Planetary Science, Vol: 52, Pages: 1542-1560, ISSN: 1086-9379
Asteroids and comets 10–100 m in size that collide with Earth disrupt dramatically in the atmosphere with an explosive transfer of energy, caused by extreme air drag. Such airbursts produce a strong blastwave that radiates from the meteoroid's trajectory and can cause damage on the surface. An established technique for predicting airburst blastwave damage is to treat the airburst as a static source of energy and to extrapolate empirical results of nuclear explosion tests using an energy-based scaling approach. Here we compare this approach to two more complex models using the iSALE shock physics code. We consider a moving-source airburst model where the meteoroid's energy is partitioned as two-thirds internal energy and one-third kinetic energy at the burst altitude, and a model in which energy is deposited into the atmosphere along the meteoroid's trajectory based on the pancake model of meteoroid disruption. To justify use of the pancake model, we show that it provides a good fit to the inferred energy release of the 2013 Chelyabinsk fireball. Predicted overpressures from all three models are broadly consistent at radial distances from ground zero that exceed three times the burst height. At smaller radial distances, the moving-source model predicts overpressures two times greater than the static-source model, whereas the cylindrical line-source model based on the pancake model predicts overpressures two times lower than the static-source model. Given other uncertainties associated with airblast damage predictions, the static-source approach provides an adequate approximation of the azimuthally averaged airblast for probabilistic hazard assessment.
Forman LV, Bland PA, Timms NE, et al., 2017, Defining the mechanism for compaction of the CV chondrite parent body, Geology, Vol: 45, Pages: 559-562, ISSN: 1943-2682
The Allende meteorite, a relatively unaltered member of the CV carbonaceous chondrite group, contains primitive crystallographic textures that can inform our understanding of early Solar System planetary compaction. To test between models of porosity reduction on the CV parent body, complex microstructures within ~0.5-mm-diameter chondrules and ~10-μm-long matrix olivine grains were analyzed by electron backscatter diffraction (EBSD) techniques. The large area map presented is one of the most extensive EBSD maps to have been collected in application to extraterrestrial materials. Chondrule margins preferentially exhibit limited intragrain crystallographic misorientation due to localized crystal-plastic deformation. Crystallographic preferred orientations (CPOs) preserved by matrix olivine grains are strongly coupled to grain shape, most pronounced in shortest dimension <a>, yet are locally variable in orientation and strength. Lithostatic pressure within plausible chondritic model asteroids is not sufficient to drive compaction or create the observed microstructures if the aggregate was cold. Significant local variability in the orientation and intensity of compaction is also inconsistent with a global process. Detailed microstructures indicative of crystal-plastic deformation are consistent with brief heating events that were small in magnitude. When combined with a lack of sintered grains and the spatially heterogeneous CPO, ubiquitous hot isostatic pressing is unlikely to be responsible. Furthermore, Allende is the most metamorphosed CV chondrite, so if sintering occurred at all on the CV parent body it would be evident here. We conclude that the crystallographic textures observed reflect impact compaction and indicate shock-wave directionality. We therefore present some of the first significant evidence for shock compaction of the CV parent body.
Rae A, Collins GS, Grieve RAF, et al., 2017, Complex crater formation: Insights from combining observations of shock pressure distribution with numerical models at the West Clearwater Lake impact structure, Meteoritics & Planetary Science, Vol: 52, Pages: 1330-1350, ISSN: 1086-9379
Large impact structures have complex morphologies, with zones of structural uplift that can be expressed topographically as central peaks and/or peak rings internal to the crater rim. The formation of these structures requires transient strength reduction in the target material and one of the proposed mechanisms to explain this behavior is acoustic fluidization. Here, samples of shock-metamorphosed quartz-bearing lithologies at the West Clearwater Lake impact structure, Canada, are used to estimate the maximum recorded shock pressures in three dimensions across the crater. These measurements demonstrate that the currently-observed distribution of shock metamorphism is strongly controlled by the formation of the structural uplift. The distribution of peak shock pressures, together with apparent crater morphology and geological observations, is compared with numerical impact simulations to constrain parameters used in the block-model implementation of acoustic fluidization. The numerical simulations produce craters that are consistent with morphological and geological observations. The results show that the regeneration of acoustic energy must be an important feature of acoustic fluidization in crater collapse, and should be included in future implementations. Based on the comparison between observational data and impact simulations we conclude that the West Clearwater Lake structure had an original rim (final crater) diameter of 35–40 km and has since experienced up to ~2 km of differential erosion.
Smith R, 2017, Numerical modelling of tsunami generated by deformable submarine slides
Submarine slides can generate tsunami waves that cause significant damage and loss of life. Numerical modelling of submarine slide generated waves is complex and computationally challenging, but is useful to understand the nature of the waves that are generated, and identify the important factors in determining wave characteristics which in turn are used in risk assessments. In this work, the open-source, finite-element, unstructured mesh fluid dynamics framework Fluidity is used to simulate submarine slide tsunami using a number of different numerical approaches. First, three alternative approaches for simulating submarine slide acceleration, deformation and wave generation with full coupling between the slide and water in two dimensions are compared. Each approach is verified against benchmarks from experimental and other numerical studies, at different scales, for deformable submarine slides. There is good agreement to both laboratory results and other numerical models, both with a fixed mesh and a dynamically adaptive mesh, tracking important features of the slide geometry as the simulation progresses. Second, Fluidity is also used in a single-layer Bousinesq approximation in conjunction with a prescribed velocity boundary condition to model the propagation of slide tsunami in two and three dimensions. A new, efficient approach for submarine slide tsunami that accounts for slide dynamics and deformation is developed by imposing slide dynamics, derived from multi-material simulations. Two submarine slides are simulated in the Atlantic Ocean, and these generate waves up to 10 m high at the coast of the British Isles. Results indicate the largest waves are generated in the direction of slide motion. The lowest waves are generated perpendicular to the slide motion. The slide velocity and acceleration are the most important factors in determining wave height. Slides that deform generate higher waves than rigid slides, although this effect is of secondary importance f
Collins GS, Moon Formation: Punch Combo or Knock-out Blow, Nature Geoscience, ISSN: 1752-0908
Morgan JV, 2016, The formation of peak rings in large impact craters, Science, Vol: 354, Pages: 878-882, ISSN: 0036-8075
Large impacts provide a mechanism for resurfacing planets through mixing near-surface rocks with deeper material. Central peaks are formed from the dynamic uplift of rocks during crater formation. As crater size increases, central peaks transition to peak rings. Without samples, debate surrounds the mechanics of peak-ring formation and their depth of origin. Chicxulub is the only known impact structure on Earth with an unequivocal peak ring, but it is buried and only accessible through drilling. Expedition 364 sampled the Chicxulub peak ring, which we found was formed from uplifted, fractured, shocked, felsic basement rocks. The peak-ring rocks are cross-cut by dikes and shear zones and have an unusually low density and seismic velocity. Large impacts therefore generate vertical fluxes and increase porosity in planetary crust.
Johnson BC, Blair DM, Collins GS, 2016, Formation of the Orientale lunar multiring basin, Science, Vol: 354, Pages: 441-444, ISSN: 0036-8075
Multiring basins, large impact craters characterized by multiple concentric topographic rings, dominate the stratigraphy, tectonics, and crustal structure of the Moon. Using a hydrocode, we simulated the formation of the Orientale multiring basin, producing a subsurface structure consistent with high-resolution gravity data from the Gravity Recovery and Interior Laboratory (GRAIL) spacecraft. The simulated impact produced a transient crater, ~390 kilometers in diameter, that was not maintained because of subsequent gravitational collapse. Our simulations indicate that the flow of warm weak material at depth was crucial to the formation of the basin’s outer rings, which are large normal faults that formed at different times during the collapse stage. The key parameters controlling ring location and spacing are impactor diameter and lunar thermal gradients.
Kring DA, Kramer GY, Collins GS, et al., 2016, Peak-Ring Structure and Kinematics from a Multi-disciplinary Study of the Schrödinger Impact Basin, Nature Communications, Vol: 7, ISSN: 2041-1723
The Schrödinger basin on the lunar farside is ~320 km in diameter and the best-preservedpeak-ring basin of its size in the Earth–Moon system. Spectral and photogeologic analyses ofdata from the Moon Mineralogy Mapper instrument on the Chandrayaan-1 spacecraft and theLunar Reconnaissance Orbiter Camera (LROC) on the LRO spacecraft indicate the peak ring iscomposed of anorthositic, noritic, and troctolitic lithologies that were juxtaposed by severalcross-cutting faults during peak ring formation. Hydrocode simulations indicate the lithologieswere uplifted from depths up to 30 km, representing the crust of the lunar farside. Combining2geological and remote-sensing observations with numerical modeling, here we show a DisplacedStructural Uplift model is best for peak rings, including that in the K-T Chicxulub impact crateron Earth. These results may help guide sample selection in lunar sample return missions that arebeing studied for the multi-agency International Space Exploration Coordination Group.Determining which lunar landing site may yield information about the lunar interior is veryimportant with impact basins usually the best sites. Kring et al. provide a geological map of theSchrödinger basin on the moon via a multidisciplinary approach of remote sensing and numericalmodeling.
Forman LV, Bland PA, Timms NE, et al., 2016, Hidden Secrets Of Deformation: Impact-Induced Compaction Within A CV Chondrite, Earth and Planetary Science Letters, Vol: 452, Pages: 133-145, ISSN: 0012-821X
The CV3 Allende is one of the most extensively studied meteorites in worldwide collections. It is currently classified as S1—essentially unshocked—using the classification scheme of Stöffler et al. (1991), however recent modelling suggests the low porosity observed in Allende indicates the body should have undergone compaction-related deformation. In this study, we detail previously undetected evidence of impact through use of Electron Backscatter Diffraction mapping to identify deformation microstructures in chondrules, AOAs and matrix grains. Our results demonstrate that forsterite-rich chondrules commonly preserve crystal-plastic microstructures (particularly at their margins); that low-angle boundaries in deformed matrix grains of olivine have a preferred orientation; and that disparities in deformation occur between chondrules, surrounding and non-adjacent matrix grains. We find heterogeneous compaction effects present throughout the matrix, consistent with a highly porous initial material. Given the spatial distribution of these crystal-plastic deformation microstructures, we suggest that this is evidence that Allende has undergone impact-induced compaction from an initially heterogeneous and porous parent body. We suggest that current shock classifications (Stöffler et al., 1991) relying upon data from chondrule interiors do not constrain the complete shock history of a sample.
Miljkovic K, Collins GS, Miljkovic K, et al., 2016, Subsurface morphology and scaling of lunar impact basins, Journal of Geophysical Research: Planets, Vol: 121, Pages: 1695-1712, ISSN: 2169-9097
Impact bombardment during the first billion years after the formation of the Moon produced at least several tens of basins. The Gravity Recovery and Interior Laboratory (GRAIL) mission mapped the gravity field of these impact structures at significantly higher spatial resolution than previous missions, allowing for detailed subsurface and morphological analyses to be made across the entire globe. GRAIL-derived crustal thickness maps were used to define the regions of crustal thinning observed in centers of lunar impact basins, which represents a less unambiguous measure of a basin size than those based on topographic features. The formation of lunar impact basins was modeled numerically by using the iSALE-2D hydrocode, with a large range of impact and target conditions typical for the first billion years of lunar evolution. In the investigated range of impactor and target conditions, the target temperature had the dominant effect on the basin subsurface morphology. Model results were also used to update current impact scaling relationships applicable to the lunar setting (based on assumed target temperature). Our new temperature-dependent impact-scaling relationships provide estimates of impact conditions and transient crater diameters for the majority of impact basins mapped by GRAIL. As the formation of lunar impact basins is associated with the first ~700 Myr of the solar system evolution when the impact flux was considerably larger than the present day, our revised impact scaling relationships can aid further analyses and understanding of the extent of impact bombardment on the Moon and terrestrial planets in the early solar system.
Davison M, Collins GS, Bland PA, 2016, MESOSCALE MODELLING OF IMPACT COMPACTION OF PRIMITIVE SOLAR SYSTEM SOLIDS, 79th Annual Meeting of the Meteoritical-Society, Publisher: WILEY-BLACKWELL, Pages: A221-A221, ISSN: 1086-9379
Lyons RJ, Ciesla FJ, Bowling TJ, et al., 2016, THE EFFECT OF EARLY IMPACTS ON IRON METEORITE COOLING RATES, 79th Annual Meeting of the Meteoritical-Society, Publisher: WILEY-BLACKWELL, Pages: A433-A433, ISSN: 1086-9379
Davison TM, Collins GS, Bland PA, et al., 2016, MESOSCALE MODELLING OF THE COMPACTION OF WATER-RICH ASTEROIDS BY LOW-VELOCITY IMPACTS, 79th Annual Meeting of the Meteoritical-Society, Publisher: WILEY-BLACKWELL, Pages: A222-A222, ISSN: 1086-9379
Derrick JG, Rutherford ME, Davison TM, et al., 2016, INTERROGATING HETEROGENEOUS COMPACTION OF METEORITIC MATERIAL AT THE MESOSCALE USING ANALOG EXPERIMENTS AND NUMERICAL MODELS, 79th Annual Meeting of the Meteoritical-Society, Publisher: WILEY-BLACKWELL, Pages: A228-A228, ISSN: 1086-9379
Davison TM, Collins G, Bland P, 2016, Mesoscale Modeling Of Impact Compaction Of Primitive Solar System Solids, Astrophysical Journal, Vol: 821, ISSN: 1538-4357
We have developed a method for simulating the mesoscale compaction of early solar system solids in low velocity impact events, using the iSALE shock physics code. Chondrules are represented by nonporous disks, placed within a porous matrix. By simulating impacts into bimodal mixtures over a wide range of parameter space (including the chondrule-to-matrix ratio, the matrix porosity and composition and the impact velocity), we have shown how each of these parameters influences the shock processing of heterogeneous materials. The temperature after shock processing shows a strong dichotomy: matrix temperatures are elevated much higher than the chondrules, which remain largely cold. Chondrules can protect some matrix from shock compaction, with shadow regions in the lee side of chondrules exhibiting higher porosity that elsewhere in the matrix. Using the results from this mesoscale modelling, we show how the ε − α porous compaction model parameters depend on initial bulk porosity. We also show that the timescale for the temperature dichotomy to equilibrate is highly dependent on the porosity of the matrix after the shock, and will be on the order of seconds for matrix porosities of less than 0.1, and on the order of 10’s to 100’s seconds for matrix porosities of ~ 0.3–0.5. Finally, we have shown that the composition of the post-shock material is able to match the bulk porosity and chondrule-to-matrix ratios of meteorite groups such as carbonaceous chondrites and unequilibrated ordinary chondrites.
Smith R, Hill J, Collins GS, et al., 2016, Comparing approaches for numerical modelling of tsunami generation by deformable submarine slides, Ocean Modelling, Vol: 100, Pages: 125-140, ISSN: 1463-5003
Tsunami generated by submarine slides are arguably an under-consideredrisk in comparison to earthquake-generated tsunami. Numerical simulationsof submarine slide-generated waves can be used to identify the important factorsin determining wave characteristics. Here we use Fluidity, an open sourcefinite element code, to simulate waves generated by deformable submarineslides. Fluidity uses flexible unstructured meshes combined with adaptivitywhich alters the mesh topology and resolution based on the simulationstate, focussing or reducing resolution, when and where it is required. Fluidityalso allows a number of different numerical approaches to be taken tosimulate submarine slide deformation, free-surface representation, and wavegeneration within the same numerical framework. In this work we use amulti-material approach, considering either two materials (slide and waterwith a free surface) or three materials (slide, water and air), as well as asediment model (sediment, water and free surface) approach. In all casesthe slide is treated as a viscous fluid. Our results are shown to be consistentwith laboratory experiments using a deformable submarine slide, anddemonstrate good agreement when compared with other numerical models.The three different approaches for simulating submarine slide dynamics andtsunami wave generation produce similar waveforms and slide deformationgeometries. However, each has its own merits depending on the application.Mesh adaptivity is shown to be able to reduce the computational cost withoutcompromising the accuracy of results.
Johnson BC, Collins GS, Minton DA, et al., 2016, Spherule layers, crater scaling laws, and the population of ancient terrestrial craters, Icarus, Vol: 271, Pages: 350-359, ISSN: 1090-2643
Ancient layers of impact spherules provide a record of Earth's early bombardment history. Here, we compare different bombardment histories to the spherule layer record and show that 3.2-3.5 Ga the flux of large impactors (10-100 km in diameter) was likely 20-40 times higher than today. The E-belt model of early Solar System dynamics suggests that an increased impactor flux during the Archean is the result of the destabilization of an inward extension of the main asteroid belt (Bottke, W.F., Vokrouhlický, D., Minton, D., Nesvorný, D., Morbidelli, A., Brasser, R., Simonson, B., Levison, H.F., 2012. Nature 485, 78–81). Here, we find that the nominal flux predicted by the E-belt model is 7-19 times too low to explain the spherule layer record. Moreover, rather than making most lunar basins younger than 4.1 Gyr old, the nominal E-belt model, coupled with a corrected crater diameter scaling law, only produces two lunar basins larger than 300 km in diameter. We also show that the spherule layer record when coupled with the lunar cratering record and careful consideration of crater scaling laws can constrain the size distribution of ancient terrestrial impactors. The preferred population is main-belt-like up to ∼50 km in diameter transitioning to a steep distribution going to larger sizes.
Asphaug E, Collins GS, Jutzi M, 2015, Global Scale Impacts, Asteroids IV, Editors: Michel, DeMeo, Bottke, Publisher: University of Arizona Press, Pages: 661-678, ISBN: 9780816532131
Global scale impacts modify the physical or thermal state of a substantial fraction of a target asteroid. Specific effects include accretion, family formation, reshaping, mixing and layering, shock and frictional heating, fragmentation, material compaction, dilatation, stripping of mantle and crust, and seismic degradation. Deciphering the complicated record of global scale impacts, in asteroids and meteorites, will lead us to understand the original planet-forming process and its resultant populations, and their evolution in time as collisions became faster and fewer. We provide a brief overview of these ideas, and an introduction to models.
Baker DMH, Head JW, Collins GS, et al., 2015, The formation of peak-ring basins: working hypotheses and path forward in using observations to constrain models of impact-basin formation, Icarus, Vol: 273, Pages: 146-163, ISSN: 0019-1035
Impact basins provide windows into the crustal structure and stratigraphy of planetary bodies; however, interpreting the stratigraphic origin of basin materials requires an understanding of the processes controlling basin formation and morphology. Peak-ring basins (exhibiting a rim crest and single interior ring of peaks) provide important insight into the basin-formation process, as they are transitional between complex craters with central peaks and larger multi-ring basins. New image and altimetry data from the Lunar Reconnaissance Orbiter as well as a suite of remote sensing datasets have permitted a reassessment of the origin of lunar peak-ring basins. We synthesize morphometric, spectroscopic, and gravity observations of lunar peak-ring basins and describe two working hypotheses for the formation of peak rings that involve interactions between inward collapsing walls of the transient cavity and large central uplifts of the crust and mantle. Major facets of our observations are then compared and discussed in the context of numerical simulations of peak-ring basin formation in order to plot a course for future model refinement and development.
Milbury C, Johnson BC, Melosh HJ, et al., 2015, Preimpact porosity controls the gravity signature of lunar craters, Geophysical Research Letters, Vol: 42, Pages: 9711-9716, ISSN: 1944-8007
We model the formation of lunar complex craters and investigate the effect of preimpact porosity on their gravity signatures. We find that while preimpact target porosities less than ~7% produce negative residual Bouguer anomalies (BAs), porosities greater than ~7% produce positive anomalies whose magnitude is greater for impacted surfaces with higher initial porosity. Negative anomalies result from pore space creation due to fracturing and dilatant bulking, and positive anomalies result from destruction of pore space due to shock wave compression. The central BA of craters larger than ~215 km in diameter, however, are invariably positive because of an underlying central mantle uplift. We conclude that the striking differences between the gravity signatures of craters on the Earth and Moon are the result of the higher average porosity and variable porosity of the lunar crust.
Ormo J, Melero-Asensio I, Housen K, et al., Scaling and reproducibility of craters produced at the Experimental Projectile Impact Chamber (EPIC), Centro de Astrobiología, Spain, Meteoritics & Planetary Science, ISSN: 1086-9379
Monteux J, Collins GS, Tobie G, et al., 2015, Consequences of large impacts on Enceladus' core shape, Icarus, Vol: 264, Pages: 300-310, ISSN: 1090-2643
The intense activity on Enceladus suggests a differentiated interior consisting of a rocky core, an internal ocean and an icy mantle. However, topography and gravity data suggests large heterogeneity in the interior, possibly including significant core topography. In the present study, we investigated the consequences of collisions with large impactors on the core shape. We performed impact simulations using the code iSALE2D considering large differentiated impactors with radius ranging between 25 and 100 km and impact velocities ranging between 0.24 and 2.4 km/s. Our simulations showed that the main controlling parameters for the post-impact shape of Enceladus’ rock core are the impactor radius and velocity and to a lesser extent the presence of an internal water ocean and the porosity and strength of the rock core. For low energy impacts, the impactors do not pass completely through the icy mantle. Subsequent sinking and spreading of the impactor rock core lead to a positive core topographic anomaly. For moderately energetic impacts, the impactors completely penetrate through the icy mantle, inducing a negative core topography surrounded by a positive anomaly of smaller amplitude. The depth and lateral extent of the excavated area is mostly determined by the impactor radius and velocity. For highly energetic impacts, the rocky core is strongly deformed, and the full body is likely to be disrupted. Explaining the long-wavelength irregular shape of Enceladus’ core by impacts would imply multiple low velocity (<2.4 km/s) collisions with deca-kilometric differentiated impactors, which is possible only after the LHB period.
Potter RWK, Kring DA, Collins GS, 2015, Scaling of basin-sized impacts and the influence of target temperature, Geological Society of America Special Papers, Vol: 518, Pages: 99-113, ISSN: 0072-1077
We produce a set of scaling laws for basin-sized impacts using data from a suiteof lunar basin numerical models. The results demonstrate the importance of preimpacttarget temperature and thermal gradient, which are shown to greatly infl uencethe modifi cation phase of the impact cratering process. Impacts into targets withcontrasting thermal properties also produce very different crustal and topographicprofi les for impacts of the same energy. Thermal conditions do not, however, signifi -cantly infl uence the excavation stage of the cratering process; results demonstrate,as a consequence of gravity-dominated growth, that transient crater radii are generallywithin 5% of each other over a wide range of thermal gradients. Excavationdepth-to-diameter ratios for the basin models (~0.12) agree well with experimental,geological, and geophysical estimates, suggesting basins follow proportional scaling.This is further demonstrated by an agreement between the basin models andPi- scaling laws based upon fi rst principles and experimental data. The results of thiswork should also be applicable to basin-scale impacts on other silicate bodies, includingthe Hadean Earth.
Forman LV, Bland PA, Timms NE, et al., 2015, RECOVERING THE PRIMORDIAL IMPACT HISTORY OF CHONDRITES IN UNPRECEDENTED DETAIL USING MASSIVE EBSD DATASETS, 78th Annual Meeting of the Meteoritical-Society, Publisher: WILEY-BLACKWELL, ISSN: 1086-9379
Muxworthy AR, Bland PA, Collins G, et al., 2015, MAGNETIC FABRICS IN ALLENDE: IMPLICATIONS FOR MAGNETIC REMANENCE ACQUISITION., 78th Annual Meeting of the Meteoritical-Society, Publisher: WILEY, ISSN: 1086-9379
Miljkovic K, Wieczorek MA, Collins GS, et al., 2015, Excavation of the lunar mantle by basin-forming impact events on the Moon, EARTH AND PLANETARY SCIENCE LETTERS, Vol: 409, Pages: 243-251, ISSN: 0012-821X
Jacobs CT, Goldin TJ, Collins GS, et al., 2015, An improved quantitative measure of the tendency for volcanic ash plumes to form in water: implications for the deposition of marine ash beds, JOURNAL OF VOLCANOLOGY AND GEOTHERMAL RESEARCH, Vol: 290, Pages: 114-124, ISSN: 0377-0273
Collins GS, 2015, Rock Avalanche, Encyclopedia of Planetary Landforms, Publisher: Springer New York, Pages: 1807-1811, ISBN: 9781461431336
Bland PA, Collins GS, Davison TM, et al., 2014, Pressure-temperature evolution of primordial solar system solids during impact-induced compaction, NATURE COMMUNICATIONS, Vol: 5, ISSN: 2041-1723
Collins GS, 2014, Numerical simulations of impact crater formation with dilatancy, JOURNAL OF GEOPHYSICAL RESEARCH-PLANETS, Vol: 119, Pages: 2600-2619, ISSN: 2169-9097
This data is extracted from the Web of Science and reproduced under a licence from Thomson Reuters. You may not copy or re-distribute this data in whole or in part without the written consent of the Science business of Thomson Reuters.