Measuring urban activities using Foursquare data and network analysis: a case study of Murcia (Spain)

Among social networks, Foursquare is a useful reference for identifying recommendations about local stores, restaurants, malls or other activities in the city. In this article, we consider the question of whether there is a relationship between the data provided by Foursquare regarding users’ tastes and preferences and fieldwork carried out in cities, especially those connected with business and leisure. Murcia was chosen for case study for two reasons: its particular characteristics and the prior knowledge resulting from the fieldwork. Since users of this network establish, what may be called, a ranking of places through their recommendations, we can plot these data with the objective of displaying the characteristics and peculiarities of the network in this city. Fieldwork from the city itself gives us a set of facilities and services observed in the city, which is a physical reality. An analysis of these data using a model based on a network centrality algorithm establishes a classification or ranking of the nodes that form the urban network. We compare the data extracted from the social network with the data collected from the fieldwork, in order to establish the appropriateness in terms of understanding the activity that takes place in this city. Moreover, this comparison allows us to draw conclusions about the degree of similarity between the preferences of Foursquare users and what was obtained through the fieldwork in the city.     

Two-dimensional numerical modeling of tectonic and metamorphic histories at active continental margins

The evolution of an active continental margin is simulated in two dimensions, using a finite difference thermomechanical code with half-staggered grid and marker-in-cell technique. The effect of mechanical properties, changing as a function of P and T, assigned to different crustal layers and mantle materials in the simple starting structure is discussed for a set of numerical models. For each model, representative P–T paths are displayed for selected markers. Both the intensity of subduction erosion and the size of the frontal accretionary wedge are strongly dependent on the rheology chosen for the overriding continental crust. Tectonically eroded upper and lower continental crust is carried down to form a broad orogenic wedge, intermingling with detached oceanic crust and sediments from the subducted plate and hydrated mantle material from the overriding plate. A small portion of the continental crust and trench sediments is carried further down into a narrow subduction channel, intermingling with oceanic crust and hydrated mantle material, and to some extent extruded to the rear of the orogenic wedge underplating the overriding continental crust. The exhumation rates for (ultra)high pressure rocks can exceed subduction and burial rates by a factor of 1.5–3, when forced return flow in the hanging wall portion of the self-organizing subduction channel is focused. The simulations suggest that a minimum rate of subduction is required for the formation of a subduction channel, because buoyancy forces may outweigh drag forces for slow subduction. For a weak upper continental crust, simulated by a high pore pressure coefficient in the brittle regime, the orogenic wedge and megascale melange reach a mid- to upper-crustal position within 10–20 Myr (after 400–600 km of subduction). For a strong upper crust, a continental lid persists over the entire time span covered by the simulation. The structural pattern is similar in all cases, with four zones from trench toward arc: (a) an accretionary complex of low-grade metamorphic sedimentary material; (b) a wedge of mainly continental crust, with medium-grade HP metamorphic overprint, wound up and stretched in a marble cake fashion to appear as nappes with alternating upper and lower crustal provenance, and minor oceanic or hydrated mantle interleaved material; (c) a megascale melange composed of high-pressure and ultrahigh-pressure metamorphic oceanic and continental crust, and hydrated mantle, all extruded from the subduction channel; (d) zone represents the upward tilted frontal part of the remaining upper plate lid in the case of a weak upper crust. The shape of the P–T paths and the time scales correspond to those typically recorded in orogenic belts. Comparison of the numerical results with the European Alps reveals some similarities in their gross structural and metamorphic pattern exposed after collision. A similar structure may be developed at depth beneath the forearc of the Andes, where the importance of subduction erosion is well documented, and where a strong upper crust forms a stable lid.     

三维数值模拟研究俯冲区火山的时空活动性

俯冲板块的深部脱水使得上覆地幔含水, 从而降低含水地幔的熔点, 导致上覆地幔部分熔融。 部分熔融的地幔柱一旦喷发到地表就是俯冲带火山, 也形成新的地壳。 相对于周围的地幔来讲, 具有较小密度和黏度的部分熔融地幔的时空活动性就控制着俯冲带火山的时空分布特征。 本文主要回顾近年来运用三维热力学岩石力学模型数值模拟研究与板片脱水相关的俯冲带火山活动的时空分布特性。 结果表明, 部分熔融地幔的有效黏度和密度是影响俯冲板片之上的三维地幔柱横向分布特征的主要因素。 高黏度的部分熔融地幔(10²⁰~10²¹ Pa·s )易于形成近平行于海沟的、 长波长(70~100 km)的、 薄的波状地幔柱; 低黏度(10¹⁸~10¹⁹ Pa·s )的熔融地幔易于形成平行于海沟的, 短波长(30~50 km)的蘑菇状地幔柱和垂直于海沟的山脊状地幔柱。 当部分熔融地幔和周围地幔的密度相差小于50 kg/m³时, 在俯冲板片之上只能形成长波长低幅度(宽50~100 km, 高10~15 km)的地幔山丘。 岩浆产率随着时间的变化反映了火山活动的生命周期性。 板块俯冲速度会影响地幔柱形成的深度和范围大小。 高效率熔融提取会增加新地壳增长总量。 低的板块俯冲速度和低的熔融提取效率会增加上地壳(花岗岩质)和中地壳(英安岩质)化学成分的比例。 数值模拟结果可以很好地解释如日本东北、 新西兰、 南阿拉斯加俯冲区火山的横向分布特征。     

Non‐uniform splitting of a single mantle plume by double cratonic roots: Insight into the origin of the central and southern East African Rift System

Using numerical thermo‐mechanical experiments we analyse the role of an active mantle plume and pre‐existing lithospheric thickness differences in the structural development of the central and southern East African Rift system. The plume‐lithosphere interaction model setup captures the essential features of the studied area: two cratonic bodies embedded into surrounding lithosphere of normal thickness. The results of the numerical experiments suggest that localization of rift branches in the crust is mainly defined by the initial position of the mantle plume relative to the cratons. We demonstrate that development of the Eastern branch, the Western branch and the Malawi rift can be the result of non‐uniform splitting of the Kenyan plume, which has been rising underneath the southern part of the Tanzanian craton. Major features associated with Cenozoic rifting can thus be reproduced in a relatively simple model of the interaction between a single mantle plume and pre‐stressed continental lithosphere with double cratonic roots.     

Nanochemo-mechanical signature of organic-rich shales: a coupled indentation–EDX analysis

The organic–inorganic nature of organic-rich source rocks poses several challenges for the development of functional relations that link mechanical properties with geochemical composition. With this focus in mind, we herein propose a method that enables chemo-mechanical characterization of this highly heterogeneous source rock at the micron and submicron length scale through a statistical analysis of a large array of energy-dispersive X-ray spectroscopy (EDX) data coupled with nanoindentation data. The ability to include elemental composition to the indentation probe via EDX is shown to provide a means to identify pure material phases, mixture phases, and interfaces between different phases. Employed over a large array, the statistical clustering of this set of chemo-mechanical data provides access to the properties of the fundamental building blocks of clay-dominated organic-rich source rocks. The versatility of the approach is illustrated through the application to a large number of source rocks of different origin, chemical composition, and organic content. We find that the identified properties exhibit a unique scaling relation between stiffness and hardness. This suggests that organic-rich shale properties can be reduced to their elementary constituents, with several implications for the development of predictive functional relations between chemical composition and mechanical properties of organic-rich source rocks such as the intimate interplay between clay-packing, organic maturity, and mechanical properties of porous clay/organic phase.     

Stagnant lid tectonics:Perspectives from silicate planets,dwarf planets,large moons,and large asteroids

To better understand Earth's present tectonic style-plate tectonics—and how it may have evolved from single plate(stagnant lid) tectonics, it is instructive to consider how common it is among similar bodies in the Solar System. Plate tectonics is a style of convection for an active planetoid where lid fragment(plate) motions reflect sinking of dense lithosphere in subduction zones, causing upwelling of asthenosphere at divergent plate boundaries and accompanied by focused upwellings, or mantle plumes;any other tectonic style is usefully called "stagnant lid" or "fragmented lid". In 2015 humanity completed a 50+ year effort to survey the 30 largest planets, asteroids, satellites, and inner Kuiper Belt objects,which we informally call "planetoids" and use especially images of these bodies to infer their tectonic activity. The four largest planetoids are enveloped in gas and ice(Jupiter, Saturn, Uranus, and Neptune)and are not considered. The other 26 planetoids range in mass over 5 orders of magnitude and in diameter over 2 orders of magnitude, from massive Earth down to tiny Proteus; these bodies also range widely in density, from 1000 to 5500 kg/m~3. A gap separates 8 silicate planetoids with ρ = 3000 kg/m~3 or greater from 20 icy planetoids(including the gaseous and icy giant planets) with ρ = 2200 kg/m~3 or less. We define the "Tectonic Activity Index"(TAI), scoring each body from 0 to 3 based on evidence for recent volcanism, deformation, and resurfacing(inferred from impact crater density). Nine planetoids with TAI = 2 or greater are interpreted to be tectonically and convectively active whereas 17 with TAI 2 are inferred to be tectonically dead. We further infer that active planetoids have lithospheres or icy shells overlying asthenosphere or water/weak ice. TAI of silicate(rocky) planetoids positively correlates with their inferred Rayleigh number. We conclude that some type of stagnant lid tectonics is the dominant mode of heat loss and that plate tectonics is unusual. To make progress understanding Earth's tectonic history and the tectonic style of active exoplanets, we need to better understand the range and controls of active stagnant lid tectonics.     

Samovar: a thermomechanical code for modeling of geodynamic processes in the lithosphere—application to basin evolution

We present a new 2D finite difference code, Samovar, for high-resolution numerical modeling of complex geodynamic processes. Examples are collision of lithospheric plates (including mountain building and subduction) and lithosphere extension (including formation of sedimentary basins, regions of extended crust, and rift zones). The code models deformation of the lithosphere with viscoelastoplastic rheology, including erosion/sedimentation processes and formation of shear zones in areas of high stresses. It also models steady-state and transient conductive and advective thermal processes including partial melting and magma transport in the lithosphere. The thermal and mechanical parts of the code are tested for a series of physical problems with analytical solutions. We apply the code to geodynamic modeling by examining numerically the processes of lithosphere extension and basin formation. The results are directly applicable to the Basin and Range province, western USA, and demonstrate the roles of crust–mantle coupling, preexisting weakness zones, and erosion rate on the evolutionary trends of extending continental regions. Modeling of basin evolution indicates a critical role of syn-rift sedimentation on the basin depth and a governing role of Peierls deformation in cold lithospheric mantle. While the former may increase basin depth by 50%, the latter limits the depth of rift basins by preventing faulting in the subcrustal lithosphere.     

Numerical modeling of protocore destabilization during planetary accretion: Methodology and results

We developed and tested an efficient 2D numerical methodology for modeling gravitational redistribution processes in a quasi spherical planetary body based on a simple Cartesian grid. This methodology allows one to implement large viscosity contrasts and to handle properly a free surface and self-gravitation. With this novel method we investigated in a simplified way the evolution of gravitationally unstable global three-layer structures in the interiors of large metal-silicate planetary bodies like those suggested by previous models of cold accretion [Sasaki, S., Nakazawa, K., 1986. J. Geophys. Res. 91, 9231-9238; Karato, S., Murthy, V.R., 1997. Phys. Earth Planet Interios 100, 61-79; Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819]: an innermost solid protocore (either undifferentiated or partly differentiated), an intermediate metal-rich layer (either continuous or disrupted), and an outermost silicate-rich layer. Long-wavelength (degree-one) instability of this three-layer structure may strongly contribute to core formation dynamics by triggering planetary-scale gravitational redistribution processes. We studied possible geometrical modes of the resulting planetary reshaping using scaled 2D numerical experiments for self-gravitating planetary bodies with Mercury-, Mars- and Earth-size. In our simplified model the viscosity of each material remains constant during the experiment and rheological effects of gravitational energy dissipation are not taken into account. However, in contrast to a previously conducted numerical study [Honda, R., Mizutani, H., Yamamoto, T., 1993. J. Geophys. Res. 98, 2075-2089] we explored a freely deformable planetary surface and a broad range of viscosity ratios between the metallic layer and the protocore (0.001-1000) as well as between the silicate layer and the protocore (0.001-1000). An important new prediction from our study is that realistic modes of planetary reshaping characterized by a high viscosity protocore and low viscosity molten silicate and metal [Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819] may result in the transient exposure of the protocore to the planetary surface and a strongly (up to 8% of the planetary diameter) aspherical deviation of the planetary shape during the early stages of core formation. Exposure of the protocore might happen in the early stages of iron core formation. This process may conceivably convert a large amount of potential energy into temperature increase and a transient strongly non-uniform depth of the magma ocean around the protoplanet. Our simplified model also predicts that the time for metallic core formation out of the metal-rich layer depends mainly on the dynamics of the deformation of the solid strong protocore. In nature this dynamics will be strongly dependent on the effective viscosity of the protocore, which should generally have non-Newtonian pressure-, temperature-, and stress-dependent rheology with strong thermomechanical feedbacks from gravitational energy dissipation.     

Visualization of multi-scale dynamics of hydrous cold plumes at subduction zones

Recent increases in the computational power of high-performance computing systems have led to a large gap between the high-resolution runs of numerical simulations—typically approaching 50–100 million tracers and 1–5 million grid points in two dimensions—and the modest resolution of 1–2 million pixels for conventional display devices. This technical problem is further compounded by the variety of fields produced by numerical simulations and the limited bandwidth available through the Internet in the course of collaborative ventures. We have developed a visualization system using the paradigm of web-based inquiry to address these mounting problems. We have employed, as a case study, a problem involving two-dimensional multi-scale dynamics of hydrous cold plumes at subduction zones. A Lagrangian marker method, in which the number of markers varies dynamically, is used to delineate the many different fields, such as temperature, viscosity, strain, and chemical composition. We found commercially available software to be insufficient for our visualization needs and so we were driven to develop a new set of tools tailored to high-resolution, multi-aspect, multi-scale simulations, and adaptable to many other applications in which large datasets involving tens of millions of tracers with many different fields are prevalent. In order to address this gap in visualization techniques, we have developed solutions for remote visualization and for local visualization. Our remote visualization solution is a web-based, zoomable image service (WEB-IS) that requires minimal bandwidth while allowing the user to explore our data through time, across many thermo–physical properties, and through different spatial scales. For local visualization, we found it optimal to use bandwidth-intensive, high-resolution display walls for performing parallel visualization in order to best comprehend the causal and temporal relationships between the multiple physical and chemical properties in a simulation.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s10069-004-0017-2     

Late orogenic, large-scale rotations in the Tien Shan and adjacent mobile belts in Kyrgyzstan and Kazakhstan

Most of Kazakhstan belongs to the central part of the Eurasian Paleozoic mobile belts for which previously proposed tectonic scenarios have been rather disparate. Of particular interest is the origin of strongly curved Middle and Late Paleozoic volcanic belts of island-arc and Andean-arc affinities that dominate the structure of Kazakhstan. We undertook a paleomagnetic study of Carboniferous to Upper Permian volcanics and sediments from several localities in the Ili River basin between the Tien Shan and the Junggar–Alatau ranges in southeast Kazakhstan. Our main goal was to investigate the Permian kinematic evolution of these belts, particularly in terms of rotations about vertical axes, in the hope of deciphering the dynamics that played a role during the latest Paleozoic deformation in this area. This deformation, in turn, can then be related to the amalgamation of this area with Baltica, Siberia, and Tarim in the expanding Eurasian supercontinent. Thermal demagnetization revealed that most Permian rocks retained a pretilting and likely primary component, which is of reversed polarity at three localities and normal at the fourth. In contrast, most Carboniferous rocks are dominated by postfolding reversed overprints of probably “mid-Permian” age, whereas presumably primary components are isolated from a few sites at two localities. Mean inclinations of primary components generally agree with coeval reference values extrapolated from Baltica, whereas declinations from primary as well as secondary components are deflected counterclockwise (ccw) by up to  90°. Such ccw rotated directions have previously also been observed in other Tien Shan sampling areas and in the adjacent Tarim Block to the south. However, two other areas in Kazakhstan show clockwise (cw) rotations of Permian magnetization directions. One area is located in the Kendyktas block about 300 km to the west of the Ili River valley, and the other is found in the Chingiz Range, to the north of Lake Balkhash and about 400 km to the north of the Ili River valley. The timing of the ccw as well as cw rotations is clearly later than the disappearance of any marine basins from northern Tarim, the Tien Shan and eastern Kazakhstan, so that the rotations cannot be attributed to island-arc or Andean-margin plate settings — instead we attribute the rotations to large-scale, east–west (present-day coordinates), sinistral wrenching in an intracontinental setting, related to convergence between Siberia and Baltica, as recently proposed by Natal'in and Şengör [Natal'in, B.A., and Şengör, A.M.C., 2005. Late Palaeozoic to Triassic evolution of the Turan and Scythian platforms: the pre-history of the palaeo-Tethyan closure, Tectonophysics, 404, 175–202.]. Our previous work in the Chingiz and North Tien Shan areas on Ordovician and Silurian rocks suggested relative rotations of  180°, whereas the Permian declination differences are of the order of 90° between the two areas. Thus, we assume that about 50% of the total post-Ordovician rotations are of pre-Late Permian age, with the other half of Late Permian–earliest Mesozoic age. The pre-Late Permian rotations are likely related to oroclinal bending during plate boundary evolution in a supra-subduction setting, given the calc-alkaline character of nearly all of the pre-Late Permian volcanics in the strongly curved belts.