Reactive Transport Process of Earthquake-induced Hydrochemical Changes in Guanding Thermal Spring,Western Sichuan,China

Earthquake-related hydrochemical changes in thermal springs have been widely observed;however,quantitative modeling of the reactive transport process is absent.In the present study,we apply reactive transport simulation to capture the hydrochemical responses in a thermal spring following the Wenchuan Ms 8.0 and Lushan Ms 7.0 earthquakes.We first constrain deep reservoir geothermal fluid compositions and temperature by multicomponent geothermometry,and then a reactive geochemical transport model is constructed to reproduce the hydrochemical evolution process.The results show that the recharge from the shallow aquifer increases gradually until it reaches a peak because of the permeability enhancement caused by the Lushan earthquake,which may be the mechanism to explain the earthquake-related hydrochemical responses.In contrast to the postseismic effect of the Wenchuan earthquake,the chemical evolution can be considered as hydrochemical anomalies related to the Lushan earthquake.This study proves that the efficient simulation of reactive transport processes is useful for investigating earthquake-related signals in hydrochemical time series.

Investigating the elliptic anisotropy of identified particles in p-Pb collisions with a multi-phase transport model

The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.

Can We Estimate the Impacts of Fixed Links Using Current Transport Models?Experiences from Concept Evaluations of Fixed Links on E39

Transport analysis and impact evaluations are important input for decisions about infrastructure projects.The impacts on transport from fjord crossing tunnels or bridges are the foundation for the cost benefit analysis,and also the basis for estimating the income from toll collection.Based on experiences from concept evaluations of several fixed link projects on E39,and an ongoing overall analysis,we question the results from transport analysis made by the official tools for such analysis:the RTM(regional transport model)which estimates the demand for trips below 10 km,the NTM(national transport model),for trips of 10 km or more,and the freight transport model.Both the NTM and the freight transport model are integrated in the RTM in the net assignment stage.We will demonstrate strengths and weaknesses in the transport models by showing contra intuitive or questionable results using the model as it is.The following questions arose as the initial results from the transport model were presented:Are the transport models able to capture immediate as well as long-term impacts?How would different assumptions about the monetary costs on these projects affect the forecasted demand and the cost benefit analysis?Are there other and wider ranges of impacts,if the analysis covers the total coastal highway as a whole,compared to evaluating impacts of each fixed link project individually?Do we have enough data to include transport effects of wider impacts of the fixed link projects?We had to deal with these questions in the concept evaluations carried out for the various fixed links project and in the current overall evaluation.We would like to suggest improvements in the analysis tools and emphasize requirements for knowledge about impacts of fixed links projects.

1-D coupled non-equilibrium sediment transport modeling for unsteady flows in the discontinuous Galerkin framework

A high-resolution, 1-D numerical model has been developed in the discontinuous Galerkin framework to simulate 1-D flow behavior, sediment transport, and morphological evaluation under unsteady flow conditions. The flow and sediment concentration variables are computed based on the one-dimensional shallow water flow equations, while empirical equations are used for entrainment and deposition processes. The sediment transport model includes the bed load and suspended load components. New formulations for Harten-Lax-van Leer （HLL） and Harten-Lax-van Contact （HLLC） are presented for shallow water flow equations that include the bed load and suspended load fluxes. The computational results for the flow and morphological changes after two dam break events are compared with the physical model tests. Results show that the modified HLL and HLLC formulations are robust and can accurately predict morphological changes in highly unsteady flows.

Reactive Transport Modeling of Long-Term CO2 Sequestration Mechanisms at the Shenhua CCS Demonstration Project,China

Carbon dioxide injection into deep saline aquifers results in a variety of strongly coupled physical and chemical processes. In this study, reactive transport simulations using a 2-D radial model were performed to investigate the fate of the injected CO2, the effect of CO2-water-rock interactions on mineral alteration, and the long-term CO2 sequestration mechanisms of the Liujiagou Formation sandstone at the Shenhua CCS（carbon capture and storage） pilot site of China. Carbon dioxide was injected at a constant rate of 0.1 Mt/year for 30 years, and the fluid flow and geochemical transport simulation was run for a period of 10 000 years by the TOUGHREACT code according to the underground conditions of the Liujiagou Formation. The results show that different trapping phases of CO2 vary with time. Sensitivity analyses indicate that plagioclase composition and chlorite presence are the most significant determinants of stable carbonate minerals and CO2 mineral trapping capacity. For arkosic arenite in the Liujiagou Formation, CO2 can be immobilized by precipitation of ankerite, magnesite, siderite, dawsonite, and calcite for different mineral compositions, with Ca（2＋）, Mg（2＋）, Fe（2＋） and Na＋ provided by dissolution of calcite, albite（or oligoclase） and chlorite. This study can provide useful insights into the geochemistry of CO2 storage in other arkosic arenite（feldspar rich sandstone） formations at other pilots or target sites.

An artificial-neural-network-based surrogate modeling workflow for reactive transport modeling

Process-based reactive transport modeling(RTM)integrates thermodynamic and kinetically controlled fluid-rock interactions with fluid flow through porous media in the subsurface and surface environment.RTM is usually conducted through numerical programs based on the first principle of physical processes.However,the calculation for complex chemical reactions in most available programs is an iterative process,where each iteration is in general computationally intensive.A workflow of neural networkbased surrogate model as a proxy for process-based reactive transport simulation is established in this study.The workflow includes(1)base case RTM design,(2)development of training experiments,(3)surrogate model construction based on machine learning,(4)surrogate model validation,and(5)prediction with the calibrated model.The training experiments for surrogate modeling are generated and run prior to the predictions using RTM.The results show that the predictions from the surrogate model agree well with those from processes-based RTM but with a significantly reduced computational time.The well-trained surrogate model is especially useful when a large number of realizations are required,such as the sensitivity analysis or model calibration,which can significantly reduce the computational time compared to that required by RTM.The benefits are(1)it automatizes the experimental design during the sensitivity analysis to get sufficient numbers and coverage of the training cases;(2)it parallelizes the calculations of RTM training cases during the sensitivity analysis to reduce the simulation time;(3)it uses the neural network algorithm to rank the sensitivity of the parameters and to search the optimal solution for model calibration.

Modeling the circulation and sediment transport in the Beibu Gulf

Water circulation and sediment transport in the Beibu Gulf are important for its environmental protection and resource exploitation.By employing the Regional Ocean Modeling System（ROMS）,we studied the seasonal variation of circulation,sediment transport and long-term morphological evolution in the Beibu Gulf.The simulation results show that the circulation induced by tide and wind is cyclonic both in winter and summer in the gulf and that the wind-driven circulation is stronger in winter than that in summer.The sediment concentration is higher in the Qiongzhou Strait,west of the Hainan Island and the coast of Vietnam and the Leizhou Peninsula.The sediment is transported westwards in winter and eastwards in summer in the Qiongzhou Strait.The west entrance of the Qiongzhou Strait is dominated by westward transport all the year round.The sediment discharged by rivers is deposited near the river mouths.The simulated result demonstrates that the sediment transport is mainly controlled by tidal induced bottom resuspension in the Beibu Gulf.Four characteristics are summarized for the distribution patterns of erosion and deposition.（1） The erosion and deposition are insignificant in most area of the gulf.（2） Sediment deposition is more significant in the mouths of Qiongzhou Strait.（3） The erosion is observed in the seabed of Qiongzhou Strait.（4） Erosion and deposition occur alternatively in the west of Hainan Island.

Simulation of transport mechanism of radium isotopes in aquifer on the southern coast of Laizhou Bay

Naturally occurring radium(^(223)Ra,^(224)Ra,^(226)Ra,and^(228)Ra)isotopes have been widely applied as geochemical tracers in marine environments,especially when estimating the submarine groundwater discharge(SGD).In this sense,the influencing factors and transport mechanism of radium isotope activity in aquifers can be key information for SGD estimation.This work evaluates the adsorption/desorption behavior of^(224)Ra and^(226)Ra in the solid-liquid phase through a leaching experiment and analysis of field data.The results suggested that radium isotope activity was positively correlated with salinity and grain size,in the case of abundant sediments.Through ion analysis,we found that the ions(Na^(+),Ca^(2+),Mg^(2+),and Ba^(2+))exchanged with radium isotopes in the process of transport.A 1-D reactive transport model was established to simulate the transport process of radium isotope in aquifers.The model successfully simulated the variation of radium isotope desorption activity with salinity and was subsequently verified in the field.This study contributes to the understanding of the geochemical behavior of radium isotopes in aquifers and provides guidance for selecting a suitable groundwater endmember in SGD estimation.

Modeling the Synergetic Effect of Various Factors on Chloride Transport in Nonsaturated Concrete

Diffusion has been systematically described as the main mechanism of chloride transport in reinforced concrete（RC） structure, especially when the concrete is in a saturated state. However, the single mechanism of diffusion is not able to describe the actual chloride ingress in the nonsaturated concrete. Instead, it is dominated by the interaction of diffusion and convection. With the synergetic effects of various factors taken into account, this study aimed to modify and develop an analytical convection- diffusion coupling model for chloride transport in nonsaturated concrete. The model was verified by simulation of laboratory tests and field measurement. The results of comparison study demonstrate that the analytical model developed in this study is efficient and accurate in predicting the chloride profiles in the nonsaturated concrete.

Accessibility Measures in Modal Choice Models: A Proof of Concept Applied to Freight Transport

Modal choice models applied to interregional or international freight transportation network models are often based on rather coarse origin-destination matrices, containing annual transported tonnages between (sub)regions, for instance. Generally, only basic (sometimes constructed) independent variables (transportation costs or transit times) are used because other variables such as shipment sizes, service frequencies, etc. are not available. Using origin-destination matrices and an assignment model, it is also possible to compute spatial accessibility measures that can further be used as additional explanatory variables. Indeed, several published studies have identified network accessibility as an important element in the mode-choice decision. This paper also shows that the inclusion of an accessibility measure in the utility functions of a logit model substantially improves the performance of a transportation network model, both in the modal choice and the assignment levels of the classical four-step model. Consequently, the assignment of the estimated modal demands results in more accurate estimated traffic on the networks. The model presented in this paper is to be considered as a proof of concept because its workflow should further be streamlined to make it easily useable by modelers.

Simulating Transport Capacity, Delivery Speed, and Routing Efficiency to Predict Economic Growth

This study uses a simulation-based approach to investigate the impact of delivery delays due to constraints on transport capacity, transit speed, and routing efficiencies on an economy with various levels of interdependency among firms. The simulation uses object-oriented programming to create specialized production, consumption, and transportation classes. A set of objects from each class is distributed randomly on a 2D plane. A road network is then established between fixed objects using Prim’s MST (Minimum Spanning Tree) algorithm, followed by construction of an all-pair shortest path matrix using the Floyd Warshall algorithm. A genetic algorithm-based vehicle routing problem solver employs the all-pair shortest path matrix to best plan multiple pickup and delivery orders. Production units utilize economic order quantities (EOQ) and reorder points (ROP) to manage inventory levels. Hicksian and Marshallian demand functions are utilized by consumption units to maximize personal utility. The transport capacity, transit speed, routing efficiency, and level of interdependence serve as 4 factors in the simulation, each assigned 3 distinct levels. Federov’s exchange algorithm is used to generate an orthogonal array to reduce the number of combination replays from 34 to just 9. The simulation results of a 9-run orthogonal array on an economy with 6 mining facilities, 12 industries, 8 market centers, and 8 transport hubs show that the level of firm interdependence, followed by transit speed, has the most significant impact on economic productivity. The principal component analysis (PCA) indicates that interdependence and transit speed can explain 90.27% of the variance in the data. According to the findings of this research, a dependable and efficient regional transportation network among various types of industries is critical for regional economic development.

Analytical evaluation of steady-state solute distribution in through- diffusion and membrane behavior test under non-perfectly flushing boundary conditions

The through-diffusion and membrane behavior testing procedure using a closed-system apparatus has been widely used for concurrent measurement of diffusion and membrane efficiency coefficients of low-permeability clay-based barrier materials.However,the common assumption of perfectly flushing conditions at the specimen boundaries could induce errors in analyses of the diffusion coefficients and membrane efficiencies.In this study,an innovative pseudo three-dimensional(3D)analytical method was proposed to evaluate solute distribution along the boundary surfaces of the soil-porous disks system,considering the non-perfectly flushing conditions.The results were consistent with numerical models under two scenarios considering different inflow/outflow positions.The proposed model has been demonstrated to be an accurate and reliable method to estimate solute distributions along the bound-aries.The calculated membrane efficiency coefficient and diffusion coefficient based on the proposed analytical method are more accurate,resulting in up to 50%less relative error than the traditional approach that adopts the arithmetic mean value of the influent and effluent concentrations.The retar-dation factor of the clay specimen also can be calculated with a revised cumulative mass approach.Finally,the simulated transient solute transport matched with experimental data from a multi-stage through-diffusion and membrane behavior test,validating the accuracy of the proposed method.

Converting TREx-RGB green-channel data to 557.7 nm auroral intensity:Methodology and initial results

The recently deployed Transition Region Explorer(TREx)-RGB(red-green-blue)all-sky imager(ASI)is designed to capture“true color”images of the aurora and airglow.Because the 557.7 nm green line is usually the brightest emission line in visible auroras,the green channel of a TREx-RGB camera is usually dominated by the 557.7 nm emission.Under this rationale,the TREx mission does not include a specific 557.7 nm imager and is designed to use the RGB green-channel data as a proxy for the 557.7 nm aurora.In this study,we present an initial effort to establish the conversion ratio or formula linking the RGB green-channel data to the absolute intensity of 557.7 nm auroras,which is crucial for quantitative uses of the RGB data.We illustrate two approaches:(1)through a comparison with the collocated measurement of green-line auroras from the TREx spectrograph,and(2)through a comparison with the modeled green-line intensity according to realistic electron precipitation flux measurements from low-Earth-orbit satellites,with the aid of an auroral transport model.We demonstrate the procedures and provide initial results for the TREx-RGB ASIs at the Rabbit Lake and Lucky Lake stations.The RGB response is found to be nonlinear.Empirical conversion ratios or formulas between RGB green-channel data and the green-line auroral intensity are given and can be applied immediately by TREx-RGB data users.The methodology established in this study will also be applicable to the upcoming SMILE ASI mission,which will adopt a similar RGB camera system in its deployment.

Assessment of the Biospheric Contribution to Surface Atmospheric CO_2 Concentrations over East Asia with a Regional Chemical Transport Model

A regional chemical transport model,RAMS-CMAQ,was employed to assess the impacts of biosphere-atmosphere CO2 exchange on seasonal variations in atmospheric CO2 concentrations over East Asia.Simulated CO2 concentrations were compared with observations at 12 surface stations and the comparison showed they were generally in good agreement.Both observations and simulations suggested that surface CO2 over East Asia features a summertime trough due to biospheric absorption,while in some urban areas surface CO2 has a distinct summer peak,which could be attributed to the strong impact from anthropogenic emissions.Analysis of the model results indicated that biospheric fluxes and fossil-fuel emissions are comparably important in shaping spatial distributions of CO2 near the surface over East Asia.Biospheric flux plays an important role in the prevailing spatial pattern of CO2 enhancement and reduction on the synoptic scale due to the strong seasonality of biospheric CO2 flux.The elevation of CO2 levels by the biosphere during winter was found to be larger than 5 ppm in North China and Southeast China,and during summertime a significant depletion （≥ 7 ppm） occurred in most areas,except for the Indo-China Peninsula where positive bioflux values were found.

Numerical Modelling of Dissolving and Driving Exploitation of Potash Salt in the Qarhan Playa——A Coupled Model of Reactive Solute Transport and Chemical Equilibrium in a Multi-component Underground Brine System

Firstly, the macroscopic chemical equilibrium state of a series of chemical reactions between intercrystal brine and its media salt layer （salt deposit） in Qarhan Salt Lake was studied by using the Pitzer theory. The concept of macroscopic solubility product and its relation with accumulated ore dissolving ratio were presented, which are used in the numerical model of dissolving and driving exploitation of potassium salt in Qarhan Salt Lake. And secondly, with a model forming idea of transport model for reacting solutes in the multi-component fresh groundwater system in porous media being a reference, a two-dimensional transport model coupled with a series of chemical reactions in a multi-component brine porous system （salt deposits） was developed by using the Pitzer theory. Meanwhile, the model was applied to model potassium/magnesium transport in Qarhan Salt Lake in order to study the transfer law of solid and liquid phases in the dissolving and driving process and to design the optimal injection/abstraction strategy for dissolving and capturing maximum Potassium/ Magnesium in the mining of salt deposits in Qarhan Salt Lake.

Further developments of a multi-phase transport model for relativistic nuclear collisions

A multi-phase transport(AMPT)model was constructed as a self-contained kinetic theory-based description of relativistic nuclear collisions as it contains four main components:the fluctuating initial condition,a parton cascade,hadronization,and a hadron cascade.Here,we review the main developments after the first public release of the AMPT source code in 2004 and the corre-sponding publication that described the physics details of the model at that time.We also discuss possible directions for future developments of the AMPT model to better study the properties of the dense matter created in relativistic collisions of small or large systems.

GLOBAL EXISTENCE AND ASYMPTOTIC BEHAVIOR FOR AN 1-D COMPRESSIBLE ENERGY TRANSPORT MODEL

In this article, the global existence and the large time behavior of smooth solutions to the initial boundary value problem for a degenerate compressible energy transport model are established.