Numerical method of moments for solute transport in a nonstationary flow field conditioned on hydraulic conductivity and head measurements

The unconditional stochastic studies on groundwater flow and solute transport in a nonstationary conductivity field show that the standard deviations of the hydraulic head and solute flux are very large in comparison with their mean values (Zhang et al. in Water Resour Res 36:2107–2120, 2000; Wu et al. in J Hydrol 275:208–228, 2003; Hu et al. in Adv Water Resour 26:513–531, 2003). In this study, we develop a numerical method of moments conditioning on measurements of hydraulic conductivity and head to reduce the variances of the head and the solute flux. A Lagrangian perturbation method is applied to develop the framework for solute transport in a nonstationary flow field. Since analytically derived moments equations are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. Instead of using an unconditional conductivity field as an input to calculate groundwater velocity, we combine a geostatistical method and a method of moment for flow to conditionally simulate the distributions of head and velocity based on the measurements of hydraulic conductivity and head at some points. The developed theory is applied in several case studies to investigate the influences of the measurements of hydraulic conductivity and/or the hydraulic head on the variances of the predictive head and the solute flux in nonstationary flow fields. The study results show that the conditional calculation will significantly reduce the head variance. Since the hydraulic head measurement points are treated as the interior boundary (Dirichlet boundary) conditions, conditioning on both the hydraulic conductivity and the head measurements is much better than conditioning only on conductivity measurements for reduction of head variance. However, for solute flux, variance reduction by the conditional study is not so significant.     

The Impact of the Herbicide Diuron on Photosynthesis in Three Species of Tropical Seagrass

The impact and recovery from exposure to the herbicide diuron [DCMU; 3-(3′,4′-dichlorophenyl)-1,1-dimethylurea] was assessed for three tropical seagrasses, maintained in outdoor aquaria over a 10-day period. Photosynthetic stress was detected using chlorophyll a fluorescence, measured with a Diving-PAM (pulse amplitude modulated fluorometer). Exposure to 10 and 100 μg l⁻¹ diuron resulted in a decline in effective quantum yield (ΔF/F_m^′) within 2 h of herbicide exposure in Cymodocea serrulata, Halophila ovalis and Zostera capricorni. Effective quantum yield also declined over the first 24 h of exposure in H. ovalis at even lower diuron concentrations (0.1 and 1.0 μg l⁻¹). Effective quantum yield in H. ovalis and Z. capricorni was significantly depressed at all diuron concentrations (0.1–100 μg l⁻¹) after 5 days exposure, whereas effective quantum yield in C. serrulata was only significantly lower in plants exposed to highest diuron concentrations (10 and 100 μg l⁻¹). Effective quantum yield depression was present 5 days after plants exposed to 10 and 100 μg l⁻¹ diuron were returned to fresh seawater. These results indicate that exposure to herbicide concentrations present in nearshore Queensland sediments present a potential risk to seagrass functioning.     

Numerical study on tide-driven submarine groundwater discharge and seawater recirculation in heterogeneous aquifers

There are many factors affecting submarine groundwater discharge (SGD). However, systematic study of the influences of these factors is still limited. In this study, numerical modeling is performed to quantitatively explore the influences of various factors on SGD in a coastal aquifer. In such locations, tidal and terrestrial hydraulic gradients are the primary forces driving fresh and salt water movement. Unlike steady-state flow, dynamic fresh and salt water mixing at the near-shore seafloor may form an intertidal mixing zone (IMZ) near the surface. By constructing a general SGD model, the effects of various model components such as boundary conditions, model geometry and hydraulic parameters are systematically studied. Several important findings are obtained from the study results: (1) Previous studies have indicated there will be a freshwater discharge tube between the classic transition zone and the IMZ. However, this phenomenon may become unclear with the increase of heterogeneity and anisotropy of the medium’s conductivity field. (2) SGD and IMZ are both more sensitive to the vertical anisotropy ratio of hydraulic conductivity (K_x/K_z) than to the horizontal ratio (K_x/K_y). (3) Heterogeneity of effective porosity significantly affects SGD and IMZ. (4) Increase of the storage coefficient decreases fresh water discharge but increases mixing salt water discharge and total SGD. The increase will also change the shape of the IMZ. (5) Variation of dispersivities does not affect SGD, but significantly changes the distributions of the IMZ and the whole mixing zone. These findings will be helpful to the sampling design of field studies of SGD and to the application of dynamic SGD models to field sites for model development and calibration.     

Application of a data assimilation method via an ensemble Kalman filter to reactive urea hydrolysis transport modeling

With growing importance of water resources in the world, remediations of anthropogenic contaminations due to reactive solute transport become even more important. A good understanding of reactive rate parameters such as kinetic parameters is the key to accurately predicting reactive solute transport processes and designing corresponding remediation schemes. For modeling reactive solute transport, it is very difficult to estimate chemical reaction rate parameters due to complex processes of chemical reactions and limited available data. To find a method to get the reactive rate parameters for the reactive urea hydrolysis transport modeling and obtain more accurate prediction for the chemical concentrations, we developed a data assimilation method based on an ensemble Kalman filter (EnKF) method to calibrate reactive rate parameters for modeling urea hydrolysis transport in a synthetic one-dimensional column at laboratory scale and to update modeling prediction. We applied a constrained EnKF method to pose constraints to the updated reactive rate parameters and the predicted solute concentrations based on their physical meanings after the data assimilation calibration. From the study results we concluded that we could efficiently improve the chemical reactive rate parameters with the data assimilation method via the EnKF, and at the same time we could improve solute concentration prediction. The more data we assimilated, the more accurate the reactive rate parameters and concentration prediction. The filter divergence problem was also solved in this study.     

Mid-Mountain Clouds at Whistler During the Vancouver 2010 Winter Olympics and Paralympics

A comprehensive study of mid-mountain clouds and their impacts on the Vancouver 2010 Winter Olympics and Paralympics is presented. Mid-mountain clouds were frequently present on the Whistler alpine venue, as identified in an extensive archive of webcam images over a 45-day period from February 5 to March 21, 2010. These clouds posed serious forecast challenges and had significant impacts on some Olympic and Paralympic alpine skiing competitions. Under fair weather conditions, a diurnal upslope (anabatic) flow can work in concert with a diurnal temperature inversion aloft to produce a localized phenomenon known as “Harvey’s Cloud” at Whistler. Two detailed case studies in this paper suggest that mid-mountain clouds can also develop in the area as a result of a moist valley flow interacting with a downslope flow descending from the mountaintop. A southerly inflow through the Sea-to-Sky corridor can be channeled by the local topography into a westerly upslope flow toward Whistler Mountain, resulting in orographic clouds on the alpine venue. Under favorable circumstances, these clouds are trapped to the mid-mountain zone by the leeward subsidence of an elevated southerly flow. The presence of the downslope subsidence was manifested by a distinguished dry layer observed on the top of the mid-mountain clouds in both cases. It is the subsidence-induced adiabatic warming that imposes a strong buoyant suppression to trap the mid-mountain cloud. On the other hand, the subsidence-induced dry layer has the potential to trigger evaporative instability to periodically breakup the mid-mountain cloud.     

Mineralogy and geothermometry of high-temperature rhyolites from the central and western Snake River Plain

Voluminous mid-Miocene rhyolitic ash-flow tuffs and lava flows are exposed along the northern and southern margins of the central and western Snake River Plain. These rhyolites are essentially anhydrous with the general mineral assemblage of plagioclase ±sanidine ± quartz + augite + pigeonite ± hypersthene ± fayalitic olivine + Fe-Ti oxides + apatite + zircon which provides an opportunity to compare feldspar, pyroxene, and Fe-Ti oxide equilibration temperatures for the same rocks. Estimated pyroxene equilibration temperatures (based on the geothermometers of Lindsley and coworkers) range from 850 to 1000°C, and these are well correlated with whole-rock compositions. With the exception of one sample, agreement between the two-pyroxene thermometers tested is well within 50°C. Fe-Ti oxide geothermometers applied to fresh magnetite and ilmenite generally yield temperatures about 50 to 100°C lower than the pyroxene temperatures, and erratic results are obtained if these minerals exhibit effects of subsolidus oxidation and exsolution. Results of feldspar thermometry are more complicated, and reflect uncertainties in the thermometer calibrations as well as in the degree of attainment of equilibrium between plagioclase and sanidine. In general, temperatures obtained using the Ghiorso (1984) and Green and Usdansky (1986) feldspar thermometers agree with the pyroxene temperatures within the respective uncertainties. However, uncertainties in the feldspar temperatures are the larger of the two (and exceed ±60°C for many samples). The feldspar thermometer of Fuhrman and Lindsley (1988) produces systematically lower temperatures for many of the samples studied. The estimated pyroxene temperatures are considered most representative of actual magmatic temperatures for these rhyolites. This range of temperatures is significantly higher than those for rhyolites from many other suites, and is consistent with the hypothesis that the Snake River Plain rhyolitic magmas formed by partial fusion of relatively dry (e.g. granulitic) crustal lithologies.     

Abandoned Underground Storage Tank Location Using Fluxgate Magnetic Surveying: A Case Study

In 1993, during the removal of a diesel and a gasoline underground storage tank at the municipal garage of the Village of Kohler, Sheboygan County, Wisconsin, soil testing revealed environmental contamination at the site. A site investigation revealed the possibility of a second on-site source of petroleum contamination. Limited historical data and the present usage of structures within the suspected source area precluded the use of most invasive sampling methods and most geophysical techniques. A fluxgate magnetometer survey, followed by confirmatory excavation, was conducted at the site. The fluxgate magnetometer survey identified nine possible magnetic anomalies within the 18 × 25 m area. The subsequent excavation near the anomalies revealed the presence of five paired and two individual 2000 L underground storage tanks. The fluxgate magnetometer survey, although affected by the proximity of buildings, was able to detect the buried tanks within 3 m of the brick structures, using a 1.5 × 1.5 m sampling array.     

Comparisons of Snowfall Measurements in Complex Terrain Made During the 2010 Winter Olympics in Vancouver

Solid precipitation (SP) intensity ( $ R_{\text{s}} $ ) using four automatic gauges, Pluvio, PARSIVEL (PArticle, SIze and VELocity), FD12P and POSS, and radar reflectivity factor ( $ Z $ ) using the POSS and PARSIVEL were measured at a naturally sheltered station (VOA) located at high level (1,640 m) on the Whistler Mountain in British Colombia, Canada. The R _s and other standard meteorological parameters were collected from March 2009, and from November 2009, to February 2010. The wind speed (ws) measured during this period ranged from 0 to 4.5 ms^?1, with a mean value of 0.5 ms^?1. The temperature varied from 4 to ?17 °C. The SP amount reported by the PARSIVEL was higher than that reported by the Pluvio by more than a factor of 2, while the FD12P and POSS measured relatively smaller amounts, but much closer to that reported by the Pluvio and manual measurements. The dependence of R _s from the PARSIVEL on wind speed was examined, but no significant dependence was found. The PARSIVEL’s precipitation retrieval algorithm was modified and tested using three different snow density size relationships (ρ _s–D) reported in literature. It was found that after modification of the algorithm, the derived R _s amounts using the raw data agreed reasonably well with the Pluvio. Statistical analysis shows that more than 95 % of $ Z_{{h_{\text{poss}} }} $ data measured by POSS appears to correlates well with the reflectivity factors determined using the three ρ _s–D relationships. The automated Pluvio accumulation and manually determined daily SP amount (SPm) measured during five winter months were compared. The mean ratio (MR) and the mean difference (MD), and the correlation coefficient (r) calculated using the data collected using the two methods, were found to be 0.96, 0.4 and 0.6 respectively, indicating respectable agreement between these two methods, with only the Pluvio underestimating the amount by about 4 %.     

Xixth International Congress of Photogrammetry and Remote Sensing

The XIXth International Congress of Photogrammetry and Remote Sensing was held at the Amsterdam RAI International Exhibition and Congress Centre, the Netherlands, from 16th to 23rd July 2000. Reports are given on the Congress as a whole, including the General Assembly, some Technical Commission activities and the Congress Exhibition. Papers from the Congress are published in the International Archives of Photogrammetry and Remote Sensing, 33.