Exploring the associations between cooling centre accessibility and marginalization in Montreal,Toronto, and Vancouver,Canada

Cooling centres provide respite, safety, and social support during extreme heat events for populations that do not have the resources to own or operate in-home air conditioning. The objective of this study was to measure the spatial accessibility of cooling centres and analyze the associations between cooling centre access and marginalization in Montreal, Toronto, and Vancouver, Canada. The potential spatial accessibility of cooling centres within a 15-minute walk was measured at the dissemination area scale using the two-step floating catchment area method. A two-stage modelling approach was used to analyze the associations between cooling centre access and marginalization. Approximately 62%, 58%, and 54% of the populations in Montreal, Toronto, and Vancouver had access to at least one cooling centre. In Montreal and Vancouver, high marginalization areas were more likely to have cooling centre access than low marginalization areas. Of the areas with cooling centre access, smaller access scores were observed in areas with high residential instability. Approximately one-fifth of the areas in each city had no cooling centre access and high marginalization, and may be considered for future cooling centres or programs that improve accessibility to existing centres.     

First Insights into the Technique Used for Heat Treatment of Chert at the Solutrean Site of Laugerie‐Haute,France

The earliest evidence of flint and chert heat treatment was found in the ~21.5–17 ka old European Solutrean culture. The appearance of pyrotechnology as part of the production of stone tools has important implications for our understanding of Upper Palaeolithic technological evolution and the specific adaptations during the last glacial maximum in Europe. However, the techniques and procedures used to heat‐treat rocks during the Solutrean remain poorly understood. No direct archaeological evidence has so far been found and the most promising approach is to understand these techniques by determining the parameters with which flint and chert were heated at that time. In this study, we investigate the heating temperature of 44 heat‐treated laurel‐leaf points from Laugerie‐Haute, using a non‐destructive technique based on infrared spectroscopy. Our results document that most of the artefacts were heated to a narrow interval of temperatures between 250 °C and 300 °C. This indicates a standardized technique that allowed to created similar conditions during successive heating cycles. The implications of these results for our understanding of the technical complexity during the Solutrean must be discussed in the light of different heating techniques used at different places and periods.     

人流联系和经济联系视角下区域城市关联比较——基于手机信令数据和企业关联数据的研究

手机信令数据所代表的人流联系与国内外研究常用的企业关联数据在测度城市关联存在何种差异,尚缺乏明确的研究。本文以江西省北部地区为研究对象,分别以手机信令数据、企业关联数据测度城市关联,比较了两种城市关联与城市之间空间距离、人口规模的关系,比较了两种城市关联网络的层级、结构、腹地。主要研究发现：人流联系不仅集中在城市周边地区,而且集中于较高等级城市之间;经济联系则主要集中在较高等级城市之间。人流联系、经济联系受到城市人口规模影响均较小。人流联系显著地符合幂函数空间衰减规律。手机信令数据和企业关联数据是两种不同的城市关联数据,但可互为补充,较为全面地反映城市关联。     

Multiphase flow simulation through porous media with explicitly resolved fractures

Accurate simulation of multiphase flow in fractured porous media remains a challenge. An important problem is the representation of the discontinuous or near discontinuous behaviour of saturation in real geological formations. In the classical continuum approach, a refined mesh is required at the interface between fracture and porous media to capture the steep gradients in saturation and saturation‐dependent transport properties. This dramatically increases the computational load when large numbers of fractures are present in the numerical model. A discontinuous finite element method is reported here to model flow in fractured porous media. The governing multiphase porous media flow equations are solved in the adaptive mesh computational fluid dynamics code IC‐FERST on unstructured meshes. The method is based on a mixed control volume – discontinuous finite element formulation. This is combined with the P_N+1DG‐P_NDG element pair, which has discontinuous (order N+1) representation for velocity and discontinuous (order N) representation for pressure. A number of test cases are used to evaluate the method's ability to model fracture flow. The first is used to verify the performance of the element pair on structured and unstructured meshes of different resolution. Multiphase flow is then modelled in a range of idealised and simple fracture patterns. Solutions with sharp saturation fronts and computational economy in terms of mesh size are illustrated.     

The Enlightenment,Pyropolitics, and the Problem of Evil

Abstract

In this article I argue that the ideal of rationality espoused by the European Enlightenment hinges on the separation of the light of reason from heat, with which it had been conjoined in the classical element of fire at least since the Greek antiquity. As a result, evil, from the standpoint of the Enlightenment is tantamount to heat without light, while evil, critically viewed outside this tradition, inheres in the absolute separation between the two aspects of the life-giving, albeit ineluctably dangerous, fire.     

Hybrid finite element–finite volume discretization of complex geologic structures and a new simulation workflow demonstrated on fractured rocks

The generation of computational meshes of complex geological objects is a challenge: their shape needs to be retained, resolution has to adapt to local detail, and variations of material properties in the objects have to be represented. Also mesh refinement and adaptation must be sufficient to resolve variations in the computed variable(s). Here, we present an unstructured hybrid finite element, node‐centred finite‐volume discretization suitable for solving fluid flow, reactive transport, and mechanical partial differential equations on a complex geometry with inhomogeneous material domains. We show that resulting meshes accurately capture free‐form material interfaces as defined by non‐uniform rational B‐spline curves and surfaces. The mesh discretization error is analysed for the elliptic pressure equation and an error metric is introduced to guide mesh refinement. Finite elements and finite volumes are represented in parametric space and integrations are conducted numerically. Subsequently, integral properties are mapped to physical space using Jacobian transformations. This method even retains its validity when the mesh is deformed. The resulting generic formulation is demonstrated for a transport calculation performed on a complex discrete fracture model.     

Fluid flow in shear zones: insights from the geometry and evolution of ore bodies at Renco gold mine,Zimbabwe

The geometry of mineral deposits can give insights into fluid flow in shear zones. Lode gold ore bodies at Renco Mine, in the Limpopo Belt, Zimbabwe, occur as siliceous breccias and mylonites within amphibolite facies shear zones that dip either gently or steeply. The two sets of ore bodies formed synchronously from hydrothermal fluids. The ore bodies are oblate, but have well‐defined long axes. Larger ore bodies are more oblate. High‐grade gold ore shoots have long axes that plunge down dip; this direction is perpendicular to the long axes of the low‐grade ore bodies. The centres of the high‐grade ore bodies align within the low‐grade ore bodies along strike in both gently and steeply dipping groups. The range of sizes and shapes of the ore bodies are interpreted as a growth sequence. Geometrical models are proposed for the gently and steeply dipping ore bodies, in which individual ore bodies grow with long axes plunging down dip, and merge to form larger, more oblate ore bodies. The models show that when three or more ore bodies coalesce, the long axis of the merged ore body is perpendicular to the component ore bodies, and that ore bodies in the deposit may have a range of shapes due to both growth of individual ore bodies, and their coalescence. The long axes of the high‐grade ore bodies are parallel to the shear directions of both the gently and steeply dipping dip slip shear zones, which were the directions of greatest permeability and fluid flow. The larger, lower grade bodies, which may have formed by coalescence, are elongate perpendicular to these directions.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

Heat transport at boiling,near‐critical conditions

Two‐phase flow and near‐critical phenomena are likely to enhance energy transport in high‐temperature hydrothermal systems. We present a series of two‐dimensional simulations of two‐phase flow of pure water at near‐critical conditions. The results show that at near‐critical conditions, two‐phase convection can be more efficient in transporting energy than single‐phase convection. The highest heat fluxes are attained when two‐phase heat‐pipes form near the bottom boundary, recharging the root of the upflow zone and thereby enabling the formation of broad upflow regions. When the system becomes more vapor‐dominated, it loses this ability, upflow zones become narrower and the energy efficiency drops to more moderate values.     

The significance of pockmarks to understanding fluid flow processes and geohazards

Underwater gas and liquid escape from the seafloor has long been treated as a mere curiosity. It was only after the advent of the side‐scan sonar and the subsequent discovery of pockmarks that the scale of fluid escape and the moon‐like terrain on parts of the ocean floor became generally known. Today, pockmarks ranging in size from the ‘unit pockmark’ (1–10 m wide, < 0.6 m deep) to the normal pockmark (10–700 m wide, up to 45 m deep) are known to occur in most seas, oceans, lakes and in many diverse geological settings. In addition to indicating areas of the seabed that are ‘hydraulically active’, pockmarks are known to occur on continental slopes with gas hydrates and in association with slides and slumps. However, possibly their potentially greatest significance is as an indicator of deep fluid pressure build‐up prior to earthquakes. Whereas only a few locations containing active (bubbling) pockmarks are known, those that become active a few days prior to major earthquakes may be important precursors that have been overlooked. Pockmark fields and individual pockmarks need to be instrumented with temperature and pressure sensors, and monitoring should continue over years. The scale of such research calls for a multinational project in several pockmark fields in various geological settings.