Multiphase lattice Boltzmann simulations for porous media applications

Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this article, we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.     

Iron isotope fractionation during planetary differentiation

The Fe isotope composition of samples from the Moon, Mars (SNC meteorites), HED parent body (eucrites), pallasites (metal and silicate) and the Earth's mantle were measured using high mass resolution MC-ICP-MS. These high precision measurements (δ⁵⁶Fe ≈ ± 0.04‰, 2 S.D.) place tight constraints on Fe isotope fractionation during planetary differentiation.Fractionation during planetary core formation is confined to < 0.1‰ for δ⁵⁶Fe by the indistinguishable Fe isotope composition of pallasite bulk metal (including sulfides and phosphides) and olivine separates. However, large isotopic variations (≈ 0.5‰) were observed among pallasite metal separates, varying systematically with the amounts of troilite, schreibersite, kamacite and taenite. Troilite generally has the lightest (δ⁵⁶Fe ≈ − 0.25‰) and schreibersite the heaviest (δ⁵⁶Fe ≈ + 0.2‰) Fe isotope composition. Taenite is heavier then kamacite. Therefore, these variations probably reflect Fe isotope fractionation during the late stage evolution and differentiation of the S- and P-rich metal melts, and during low-temperature kamacite exsolution, rather than fractionation during silicate-metal separation.Differentiation of the silicate portion of planets also seems to fractionate Fe isotopes. Notably, magmatic rocks (partial melts) are systematically isotopically heavier than their mantle protoliths. This is indicated by the mean of 11 terrestrial peridotite samples from different tectonic settings (δ⁵⁶Fe = + 0.015 ± 0.018‰), which is significantly lighter than the mean of terrestrial basalts (δ⁵⁶Fe = + 0.076 ± 0.029‰). We consider the peridotite mean to be the best estimate for the Fe isotope composition of the bulk silicate Earth, and probably also of bulk Earth. The terrestrial basaltic mean is in good agreement with the mean of the lunar samples (δ⁵⁶Fe = + 0.073 ± 0.019‰), excluding the high-Ti basalts. The high-Ti basalts display the heaviest Fe isotope composition of all rocks measured here (δ⁵⁶Fe ≈ + 0.2‰). This is interpreted as a fingerprint of the lunar magma ocean, which produced a very heterogeneous mantle, including the ilmenite-rich source regions of these basalts.Within uncertainties, samples from Mars (SNC meteorites), HED (eucrites) and the pallasites (average olivine + metal) have the same Fe isotope compositions as the Earth's mantle. This indicates that the solar system is very homogeneous in Fe isotopes. Its average δ⁵⁶Fe is very close to that of the IRMM-014 standard.     

Insights from a Qualitative Study of Rural Communes: Physical and Social Dimensions of Place

Existing research on people–place dynamics often focuses on either physical environment or social dimensions, but many theoretical models posit that both factors are important. I argue that rural communes provide unique settings to explore both social and environmental aspects of place. Based on qualitative interview data with members of two communes, I address emotional connection to the land, changing social relationships over time, and the interplay of physical and social aspects of place. Results demonstrate how the connections to land and social ties are both prevalent. However, social interaction has decreased while commitment to sharing and protecting the physical environment has endured.     

Double-difference Relocation of the Aftershocks of the Tecom��n, Colima, Mexico Earthquake of 22 January 2003

On 22 January 2003, the M _w?=?7.6 Tecomán earthquake struck offshore of the state of Colima, Mexico, near the diffuse triple junction between the Cocos, Rivera, and North American plates. Three-hundred and fifty aftershocks of the Tecomán earthquake with magnitudes between 2.6 and 5.8, each recorded by at least 7 stations, are relocated using the double difference method. Initial locations are determined using P and S readings from the Red Sismológica Telemétrica del Estado de Colima (RESCO) and a 1-D velocity model. Because only eight RESCO stations were operating immediately following the Tecomán earthquake, uncertainties in the initial locations and depths are fairly large, with average uncertainties of 8.0?km in depth and 1.4?km in the north?Csouth and east?Cwest directions. Events occurring between 24 January and 31 January were located using not only RESCO phase readings but also additional P and S readings from 11 temporary stations. Average uncertainties decrease to 0.8?km in depth, 0.3?km in the east?Cwest direction, and 0.7?km in the north?Csouth direction for events occurring while the temporary stations were deployed. While some preliminary studies of the early aftershocks suggested that they were dominated by shallow events above the plate interface, our results place the majority of aftershocks along the plate interface, for a slab dipping between approximately 20° and 30°. This is consistent with the slab positions inferred from geodetic studies. We do see some upper plate aftershocks that may correspond to forearc fault zones, and faults inland in the upper plate, particularly among events occurring more than 3?months after the mainshock.     

Galaxy groups in the 2dF Galaxy Redshift Survey: the catalogue

We construct a galaxy groups catalogue from the public 100-K data release of the 2dF Galaxy Redshift Survey. The group identification is carried out using a slightly modified version of the group-finding algorithm developed by Huchra & Geller. Several tests using mock catalogues allow us to find the optimal conditions to increase the reliability of the final group sample. A minimum number of four members, an outer number density enhancement of 80 and a linking radial cut-off of 200 km s⁻¹ are the best obtained values from the analysis. Using these parameters, approximately 90 per cent of groups identified in real space have a redshift space counterpart. On the other hand, the level of contamination in redshift space reaches 30 per cent, including ∼6 per cent of artificial groups and ∼24 per cent of groups associated with binaries or triplets in real space. The final sample comprises 2209 galaxy groups covering the sky region described by Colless et al.spanning over the redshift range of  0.003 ≤ z ≤ 0.25  with a mean redshift of 0.1.     

Overcoming unreported violence using place-based ambulance data: The case for mapping hotspots based on health data for crime prevention initiatives

A key concern in crime analysis is the “hidden crime” problem. Crime events unaccounted for in police records limit the external validity of official statistics and, more importantly, hinder the ability of the police to manage crime and utilize their resources effectively. The problem is exacerbated in proactive initiatives aimed at curbing violence through hotspot policing, where inaccuracies and imprecision, or, worse, no data at all, diminish prevention efforts. Previous studies have sought to overcome the data problem by juxtaposing police records with ambulance data on assault callouts and have found profound disparities. Specifically, researchers matched “crime hotspots” with “ambulance hotspots” (rather than individual events) because patient confidentiality considerations have prevented health professionals from sharing subject-level data with the police. However, health services can safely share spatial data on wider areas that do not disclose personal information. We build on this line of inquiry by analyzing data from the Thames Valley, United Kingdom, and observing spatial hotspots of different sizes. The results demonstrate that while the police and ambulance services attend to the same communities and similar types of facilities, the police are “blinded” to the location of nearly 8 out of 10 assaults. The incongruency is shown even with severe assaults, but to a lesser extent. We then simulate the reduction in injuries if the police had access to health data at different spatial levels and show that even under the most conservative set of assumptions, such an approach can prevent between 113 and 116 violent injuries each year that might otherwise require hospitalization.     

~(40)Ar/~(39)Ar Age and Geological Significance of Lamprophyres in Pishan on the SW Margin of the Tarim Terrane, NW China

The dike swarm of lamprophyres at Pishan on the SW margin of the Tarim terrane is composed of ultra-potassic lamprophyres that have intruded into the Proterozoic hornblende schists. 40Ar/39Ar dating data suggest that the isochron ages for phlogopite and whole rock separates of the lamprophyres are 231.7±0.3 Ma and 228.5±0.3 Ma, respectively. Geochemical data indicate that the dikes were derived from a subduction-related or metasomatized mantle. During the Middle-Late Triassic the western Kunlun area entered the post-orogeny period, and thus the occurrence of the dike swarm in the study area may have been controlled by underplating.     

不同气候变化情景下亚洲北山羊（Capra sibirica）的响应（英文）

目前很少见到关于气候变化影响亚洲北山羊物种栖息地的研究。通过调查气候变化对塔吉克斯坦东部亚洲北山羊(Capra sibirica)分布的影响,并采用生态位建模比较了亚洲北山羊的适宜栖息地的当前与未来分布情况。预计到2070年,现有适宜栖息地的18%(2689 km^2)将变得不适宜亚洲北山羊的生存,损失的区域主要位于研究区域的东南部和西北部地区。新的适宜栖息地可能会扩展到当前亚洲北山羊范围之外:到2070年将扩展30%(4595 km^2)的范围,这些区域与亚洲北山羊现有的分布有很强的相关性。东南部的损失与该地区当前大多数的亚洲北山羊栖息地重叠,主要出现在比研究区域海拔低得多的区域(3500–4000 m)。当同时考虑损失和收益时,亚洲北山羊可能会净扩展到新的适宜栖息地。到2070年,亚洲北山羊的平均栖息地增加量约为30%(1379 km^2),表明适宜栖息地已向北部低温栖息地转移。研究结果有助于规划气候变化情景下塔吉克斯坦东部山区对生物多样性保护的潜在影响。应该特别注意东南地区的高地山羊种群,那里的栖息地可能由于气候对山区生态系统的影响而变得不适合该物种继续生存。     

Proposal for an International Molybdenum Isotope Measurement Standard and Data Representation

Molybdenum isotopes are increasingly widely applied in Earth Sciences. They are primarily used to investigate the oxygenation of Earth's ocean and atmosphere. However, more and more fields of application are being developed, such as magmatic and hydrothermal processes, planetary sciences or the tracking of environmental pollution. Here, we present a proposal for a unifying presentation of Mo isotope ratios in the studies of mass‐dependent isotope fractionation. We suggest that the δ^98/95Mo of the NIST SRM 3134 be defined as +0.25‰. The rationale is that the vast majority of published data are presented relative to reference materials that are similar, but not identical, and that are all slightly lighter than NIST SRM 3134. Our proposed data presentation allows a direct first‐order comparison of almost all old data with future work while referring to an international measurement standard. In particular, canonical δ^98/95Mo values such as +2.3‰ for seawater and ?0.7‰ for marine Fe–Mn precipitates can be kept for discussion. As recent publications show that the ocean molybdenum isotope signature is homogeneous, the IAPSO ocean water standard or any other open ocean water sample is suggested as a secondary measurement standard, with a defined δ^98/95Mo value of +2.34 ± 0.10‰ (2s).     

Modeling Fate and Transport of Volatile Organic Compounds (VOCs) Inside Sewer Systems

Hazardous waste site investigations have shown that volatile organic compounds (VOCs) can be transported via sewer pipes and migrate into indoor spaces. Despite field data confirming the presence of this exposure pathway, there is lack of context-based numerical models that provide guidance to characterize and predict VOCs concentration in sewer gas at vapor intrusion sites. Particularly, this poses a challenge when assessing and mitigating risks associated with these exposure pathways. Therefore, a numerical model has been developed to simulate the concentration of VOCs in sewer gas in different stages throughout the sewer lines. The developed model considers various input parameters, including temperature, sewer liquid depth, groundwater depth, and sewer construction characteristics to incorporate local and operational conditions. The model's output is verified using field data from a sewer system constructed near a Superfund site. Moreover, a sensitivity analysis was conducted to evaluate the model's response to variation of the external input parameters. To the best of our knowledge, this study is the first attempt to model VOCs concentration in sewer gas, particularly to address vapor intrusion. The developed model can be used as a numerical tool to support the development of sewer assessment guidelines, risk assessment studies, and mitigation strategies.