Analytical Solutions to the Linearized Boussinesq Equation for Assessing the Effects of Recharge on Aquifer Discharge1

Boggs, Kevin G., Robert W. Van Kirk, Gary S. Johnson, Jerry P. Fairley, and P. Steve Porter, 2010. Analytical Solutions to the Linearized Boussinesq Equation for Assessing the Effects of Recharge on Aquifer Discharge. Journal of the American Water Resources Association (JAWRA) 46(6):1116–1132. DOI: 10.1111/j.1752-1688.2010.00479.x Abstract: There is a need to develop a general understanding of how variations in aquifer recharge are reflected in discharge. Analytical solutions to the linearized Boussinesq equation governing flow in an unconfined aquifer provide a unified mathematical framework to quantify relationships among lag time, attenuation and distance between aquifer recharge and discharge and the effect of an up-gradient no-flow boundary. We applied this framework to three types of recharge: (1) instantaneous, (2) periodic, and (3) constant rate for a finite duration. When the temporal scale of recharge exceeds the diffusive aquifer time scale, recharge will be reflected in discharge quickly and with little attenuation. When aquifer time scale is large, most recharge events are shorter in scale than that of the aquifer, resulting in large attenuation. Attenuation is more sensitive to boundary effects than lag time, and boundary effects increase as recharge time scale increases. Boundary effects can often be ignored when the recharge source is farther than 1/3 of the domain length away from the no-flow boundary. We illustrate analytical results with application to the economically critical Eastern Snake River Plain Aquifer in Idaho. In this aquifer, detectable annual and decadal cycles in discharge can result from recharge no farther than 20 and 60 km away from the discharge point, respectively. The effects of more distant, long-term recharge can be detected only after a time lag of several decades.     

Uncertainty in health risks due to anthropogenic primary fine particulate matter from different source types in Finland

The emission-exposure and exposure-response (toxicity) relationships are different for different emission source categories of anthropogenic primary fine particulate matter (PM_2.5). These variations have a potentially crucial importance in the integrated assessment, when determining cost-effective abatement strategies. We studied the importance of these variations by conducting a sensitivity analysis for an integrated assessment model. The model was developed to estimate the adverse health effects to the Finnish population attributable to primary PM_2.5 emissions from the whole of Europe. The primary PM_2.5 emissions in the whole of Europe and in more detail in Finland were evaluated using the inventory of the European Monitoring and Evaluation Programme (EMEP) and the Finnish Regional Emission Scenario model (FRES), respectively. The emission-exposure relationships for different primary PM_2.5 emission source categories in Finland have been previously evaluated and these values incorporated as intake fractions into the integrated assessment model. The primary PM_2.5 exposure-response functions and toxicity differences for the pollution originating from different source categories were estimated in an expert elicitation study performed by six European experts on air pollution health effects. The primary PM_2.5 emissions from Finnish and other European sources were estimated for the population of Finland in 2000 to be responsible for 209 (mean, 95% confidence interval 6–739) and 357 (mean, 95% CI 8–1482) premature deaths, respectively. The inclusion of emission-exposure and toxicity variation into the model increased the predicted relative importance of traffic related primary PM_2.5 emissions and correspondingly, decreased the predicted relative importance of other emission source categories. We conclude that the variations of emission-exposure relationship and toxicity between various source categories had significant impacts for the assessment on premature deaths caused by primary PM_2.5.     

Spatial variability of carbonaceous aerosols and associated source tracers in two cites in the Midwestern United States

Semi-continuous and 24-h averaged measurements of fine carbonaceous aerosols were made concurrently at three sites within each of two U.S. Midwestern Cities; Detroit, Michigan and Cleveland, Ohio; during two, one-month intensive campaigns conducted in July of 2007 and January & February of 2008. A comparison of 24-h measurements revealed substantial intra-urban variability in carbonaceous aerosols consistent with the influence of local sources, and excesses in both PM_2.5 organic carbon (OC) and elemental carbon (EC) were identified at individual sites within each city. High time-resolved black carbon (BC) measurements indicated that elemental carbon concentrations were higher at sites adjacent to freeways and busy surface streets, and temporal patterns suggested that excess EC at sites adjacent to freeways was dominated by mobile source emissions while excesses in EC away from traffic corridors was dominated by point/area source emissions. The site-to-site variability in OC concentrations was approximately 7% within the neighborhood scale (0.5–4 km) and between 4 and 27% at the urban scale (4–100 km). In contrast, measurements of organic source tracers, in conjunction with a Chemical Mass Balance (CMB) source-apportionment model, indicated that the spatial variation in the contribution of both mobile and stationary sources to PM_2.5 OC often exceeded the variation in OC mass concentration by a factor of 3 or more. Markers for mobile sources, biomass smoke, natural gas, and coal combustion differed by as much as 60% within the neighborhood scale and by greater than 200% within the urban scale. The observations made during this study suggest that the urban excess of carbonaceous aerosols is much more complex than has been previously reported and that a more rigorous, source-oriented approach should be taken in order to assess the risk associated with exposure to carbonaceous aerosols within the industrialized environments of the Midwestern United States.     

Ambient TSP concentration and dustfall in major cities of China: Spatial distribution and temporal variability

Based on environmental monitoring data in 93 major cities and meteorological records at 398 weather stations in China from 1981 to 2007, total suspended particle (TSP) concentration, the intensity of dustfall, and sand and dust storm frequency (F_d) were analysed. During the past 27 years, the annual average TSP concentration (C_TSP) in 93 cities was 402 μg m^?3. Annual average C_TSP decreased from the north to the south and from inland to the coast areas with a peak value of 628.8 μg m^?3 in Lanzhou. In the 1980s, 1990s and 2000s, annual average C_TSP was 628.7, 319.2, and 250.1 μg m^?3, respectively. Annual average intensity of dustfall (I_d) was 240.5 t km^?2 a^?1, decreased from northern to southern China and from inland to the coast areas with the maximum value of 717.2 t km^?2 a^?1 in Baotou. In the 1980s, 1990s and 2000s, annual average I_d was 334.8, 220.9, 146 t km^?2 a^?1 respectively. Annual average I_d in the Loess Plateau region was commonly higher than 200 t km^?2 a^?1. The annual average F_d decreased from arid regions in northwestern China to humid areas in southeastern China with two sand and sand storm centers existing in Xinjiang Taklamakan Desert and western Inner Mongolia. The annual average F_d in the 1980s, 1990s, 2000s was 16, 8, 6 days respectively, decreased steadily from 18 days in 1981–5 days in 2007. Annual average I_d had a positive linear relation to annual average C_TSP (R² = 0.96). Annual average F_d had a positive relation with annual average C_TSP (R² = 0.97) as well as annual average I_d (R² = 0.94). TSP was the chief pollutant influencing Air Pollution Index (API) in northern China in spring and winter seasons. Sand and dust storm might be a major factor affecting the temporal variability and spatial distribution of TSP and dustfall in China.     

Assimilation of conventional and satellite wind observations in a mesoscale atmospheric model for studying atmospheric dispersion

A mesoscale atmospheric model PSU/NCAR MM5 is used to provide operational weather forecasts for a nuclear emergency response decision support system on the southeast coast of India. In this study the performance of the MM5 model with assimilation of conventional surface and upper-air observations along with satellite derived 2-d surface wind data from QuickSCAT sources is examined. Two numerical experiments with MM5 are conducted: one with static initialization using NCEP FNL data and second with dynamic initialization by assimilation of observations using four dimensional data assimilation (FDDA) analysis nudging for a pre-forecast period of 12 h. Dispersion simulations are conducted for a hypothetical source at Kalpakkam location with the HYSPLIT Lagrangian particle model using simulated wind field from the above experiments. The present paper brings out the differences in the atmospheric model predictions and the differences in dispersion model results from control and assimilation runs. An improvement is noted in the atmospheric fields from the assimilation experiment which has led to significant alteration in the trajectory positions, plume orientation and its distribution pattern. Sensitivity tests using different PBL and surface parameterizations indicated the simple first order closure schemes (Blackadar, MRF) coupled with the simple soil model have given better results for various atmospheric fields. The study illustrates the impact of the assimilation of the scatterometer wind and automated weather stations (AWS) observations on the meteorological model predictions and the dispersion results.     

Verification of diffusive and pumped samplers for volatile organic compounds using a controlled atmosphere test facility

There is a requirement to verify the performance of sorbent-based passive or active samplers and to extend their use, where possible, to monitor volatile organic compounds (VOCs) that are known to be photochemical ozone pre-cursors or are relevant to the activities of the petrochemical industry. We report measurements of the 14-day diffusive uptake rates for the VOCs: i-butane (2-methyl propane), n-butane, i-pentane (2-methyl butane), n-pentane, n-hexane, benzene, toluene, and m-xylene (at environmental level concentrations) for industry standard axial samplers (Perkin–Elmer-type samplers) containing the sorbents Carbopack-X, -Z, -B or Tenax-TA. We also present data on back-diffusion, blank levels, and storage for the above sorbents, and describe the simultaneous use of the sorbent Carbopack-X for pumped sampling of certain VOCs. The results were obtained by dosing samplers in a controlled atmosphere test facility (CATFAC) operating under well-defined conditions of concentration, nominal temperature of 20 °C, wind speed of 0.5 m s^?1, and relative humidities of 0% and 80%. Field measurements were also obtained to provide supplementary data to support the laboratory study. Results are compared to existing published data, where these are available.     

A simplified approach to the indirect evaluation of the chemical composition of atmospheric aerosols from PM mass concentrations

The present work shows a straightforward procedure to indirectly estimate the chemical composition at a given site only from the determination of the PM concentrations, and the classification of the days according to the main meteorological scenarios. A previous study based on the mean chemical composition associated to the main meteorological scenarios is required.This experiment has been carried out with data from two monitoring sites in the North-Western Mediterranean, one regional and one suburban background. At both sites, one-year datasets on chemical PM₁₀ composition were obtained. Based on these datasets, the mean PM₁₀ compositions according to the most relevant meteorological situations were calculated for both locations. After that, the reconstruction of the chemical composition for all the days with available PM₁₀ concentrations was completed. Subsequently, the estimated PM₁₀ composition was compared with that determined experimentally. The comparison between the rebuilt and the experimental results was very satisfactory in the case of the regional background site, and relatively replaced in the other case. Furthermore, the validation of the method at the regional background site has been conducted from the reconstruction of a 4-year data base, and subsequent comparison with the experimental chemical composition.Our results show that it is possible to attain a good approach to the chemical composition at regional background sites, where local emission sources are negligible. Conversely, when the local sources rise in importance, i.e., at a suburban background site, the approach is suitable only for those components with a more regional origin and/or those from long range transport of air masses.     

Distributions of radionuclide sorption coefficients (Kd) in sub-surface sediments and the implications for transport calculations

The effect of the spatial variability of K_d on calculations of contaminant travel time in the vadose zone was determined. Depth discrete measurements of K_d were made for a suite of radionuclides (¹⁰⁹Cd, ⁵⁷Co, ⁶⁰Co, ⁸⁵Sr, ¹³⁷Cs, and ⁸⁸Y) utilizing a sediment core from the E-Area at the Savannah River Site. The K_d’s were ordered as ⁸⁵Sr²⁺ < ¹³⁷Cs⁺ < ¹⁰⁹Cd²⁺ < ⁵⁷Co²⁺ = ⁶⁰Co²⁺ << ⁸⁸Y³⁺ and the values generally fell below or near the lowest quartile of values reported in the literature. Correlations were generally weak between soil properties and K_d values. Most importantly, all of the K_d distributions could be reasonably approximated as log-normal. Deterministic and stochastic calculations of contaminant travel time to the water table were made. The deterministic calculations were based on each of three conceptual models of the vadose zone: complete stratification (17 strata, each with a different K_d), two strata (two sections of the vadose zone, each characterized by a single, average K_d), and unstratified (a single zone with an average K_d). Stochastic calculations were based on log-normal fits to the K_d data. The two strata model generally yielded travel times 2× greater than those in the completely stratified model. The unstratified model yielded travel times that were between 3 and 5 times greater than the completely stratified model. The stochastic mean travel times were comparable to those of the two strata model.     

Variation of radiation level and radionuclide enrichment in high background area

Significantly high radiation level and radionuclide concentration along Quilon beach area of coastal Kerala have been reported by several investigators. Detailed gamma radiation level survey was carried out using a portable scintillometer. Detailed studies on radionuclides concentration in different environmental matrices of high background areas were undertaken in the coastal areas of Karunagapalli, Kayankulam, Chavara, Neendakara and Kollam to study the distribution and enrichment of the radionuclides in the region. The absorbed gamma dose rates in air in high background area are in the range 43-17,400 nGyh⁻¹. Gamma radiation level is found to be maximum at a distance of 20 m from the sea waterline in all beaches. The soil samples collected from different locations were analysed for primordial radionuclides by gamma spectrometry. The activity of primordial radionuclides was determined for the different size fractions of soil to study the enrichment pattern. The highest activity of ²³²Th and ²²⁶Ra was found to be enriched in 125-63 μ size fraction. The preferential accumulation of ⁴⁰K was found in <63 μ fraction. The minimum ²³²Th activity was 30.2 Bq kg⁻¹, found in 1000-500 μ particle size fraction at Kollam and maximum activity of 3250.4 Bq kg⁻¹ was observed in grains of size 125-63 μ at Neendakara. The lowest ²²⁶Ra activity observed was 33.9 Bq kg⁻¹ at Neendakara in grains of size 1000-500 μ and the highest activity observed was 482.6 Bq kg⁻¹ in grains of size 125-63 μ in Neendakara. The highest ⁴⁰K activity found was 1923 Bq kg⁻¹ in grains of size <63 μ for a sample collected from Neendakara. A good correlation was observed between computed dose and measured dose in air. The correlation between ²³²Th and ²²⁶Ra was also moderately high. The results of these investigations are presented and discussed in this paper.