Flow‐Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations

This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams.     

Heat budget of the surface mixed layer south of Africa

ARGO hydrographic profiles, two hydrographic transects and satellite measurements of air–sea exchange parameters were used to characterize the properties and seasonal heat budget variations of the Surface Mixed Layer (SML) south of Africa. The analysis distinguishes the Subtropical domain (STZ) and the Subantarctic Zone (SAZ), Polar Frontal Zone (PFZ) and Antarctic Zone (AZ) of the Antarctic Circumpolar Current. While no Subantarctic Mode Water forms in that region, occurrences of deep SML (up to ∼450 m) are observed in the SAZ in anticyclones detached from the Agulhas Current retroflection or Agulhas Return Current. These are present latitudinally throughout the SAZ, but preferentially at longitudes 10–20° E where, according to previous results, the Subtropical Front is interrupted. Likely owing to this exchange window and to transfers at the Subantarctic Front also enhanced by the anticyclones, the SAZ shows a wide range of properties largely encroaching upon those of the neighbouring domains. Heat budget computations in each zone reveal significant meridional changes of regime. While air–sea heat fluxes dictate the heat budget seasonal variability everywhere, heat is mostly brought through lateral geostrophic advection by the Agulhas Current in the STZ, through lateral diffusion in the SAZ and through air–sea fluxes in the PFZ and AZ. The cooling contributions are by Ekman advection everywhere, lateral diffusion in the STZ (also favoured by the ∼10° breach in the Subtropical Front) and geostrophic advection in the SAZ. The latter likely reflects an eastward draining of water warmed through mixing of the subtropical eddies.     

Mitigating the effects of sand dunes on seismic data from the Rub al Khali basin,Saudi Arabia

The Rub Al Khali region in Saudi Arabia is characterized by the presence of sand dunes separated by salt flats, also called Sabkhas. In general, the elevation of dunes in this region varies between 90 and 250 m above sea level. The presence of these sand dunes, along with the rapidly changing surface topography poses challenges for seismic data acquisition and processing. The high contrast in acoustic impedance between the dune base and the underlying formation often results in amplification of seismic waves that are recorded at stations located on the surface of sand dunes. Attempts to address the issue using conventional surface-consistent amplitude scaling methods without reducing these amplification effects generally fail, thus compromising the suitability of the processed data for amplitude versus offset analysis. In this study, we propose a new reference site technique to reduce the effects of sand dune amplification, enabling the production of data sets that are suitable for amplitude versus offset processing. The proposed technique uses a deterministic approach to derive surface-consistent, frequency-dependent de-amplification functions for shots and stations located on the dunes only. Two-dimensional synthetic and field data examples show that the technique significantly reduces the effects of sand dune amplification.     

Influence of wave groups on SSC patterns over vortex ripples

Estimates of spatial and temporal variations in suspended sand concentrations (SSC) made with a multi-transducer Acoustic Backscatter Sensor (ABS) under a repeated wave group over a mobile rippled bed in the wave research flume at the National Hydraulics Laboratory in Ottawa, Canada, reveal an number of complex and intriguing patterns. Ensemble averages of 8 nearly identical wave groups provided much more robust estimates of SSC and allowed a detailed examination of the wave group effects. The largest SSC near the bed (< 0.10 m) occurs in phase with the largest waves in the group. Above approximately 0.10 m elevation, SSC lags behind the near bed SSC by as much as 2–3 waves; introducing significant curvature (on a semi-log plot) to the SSC profile. The log linear segments of the SSC profile grow and decay systematically on the scale of the wave group. The range in lengths of log-linear profile segments ( 0.03–0.355 m) suggest that the boundary layer thickness also fluctuates throughout the passage of the wave group. Furthermore, there are significant variations in the patterns of SSC, which occur under the largest and smallest waves in the group. Under the largest waves vertical bands of alternating high and low SSC produce an intra-wave modulation in the upper water column ( 0.075–0.30 m). The equivalent horizontal excursion of these bands scales to the ripple length. Under the smaller waves the intra-wave modulation of the SSC disappears and is replaced by temporally homogenous suspension that expands vertically through several individual wave cycles. The former pattern of homogenous suspension appears to be associated with growth of a boundary layer due to the persistent uni-directional horizontal flow during this part of the group together with the persistence of antecedent bed generated turbulence and vorticity which maintains the suspension. The latter pattern of bands of high and low SSC indicates a strong temporal and spatial constraint on the SSC (phase coupling) induced by the presence of the bedforms which may be enhanced by strong reversals in both flow and vorticity under the large waves in the group.     

Multispectral thermal infrared mapping of the 1 October 1988 Kupaianaha flow field,Kilauea volcano,Hawaii

Multispectral thermal infrared radiance measurements of the Kupaianaha flow field were acquired with the NASA airborne Thermal Infrared Multispectral Scanner (TIMS) on the morning of 1 October 1988. The TIMS data were used to map both the temperature and emissivity of the surface of the flow field. The temperature map depicted the underground storage and transport of lava. The presence of molten lava in a tube or tumulus resulted in surface temperatures that were at least 10° C above ambient. The temperature map also clearly defined the boundaries of hydrothermal plumes which resulted from the entry of lava into the ocean. The emissivity map revealed the boundaries between individual flow units within the Kupaianaha field. In general, the emissivity of the flows varied systematically with age but the relationship between age and emissivity was not unique. Distinct spectral anomalies, indicative of silica-rich surface materials, were mapped near fumaroles and ocean entry sites. This apparent enrichment in silica may have resulted from an acid-induced leaching of cations from the surfaces of glassy flows. Such incipient alteration may have been the cause for virtually all of the emissivity variations observed on the flow field, the spectral anomalies representing areas where the acid attack was most intense.     

平面SV波入射下半圆凸起地形地表运动解析解

采用波函数展开法,并借助辅助函数思想,给出了平面SV波入射下半圆凸起地形地表运动的一个解析解,并对解答的收敛性和截断计算精度进行了检验;最后通过若干算例分析了入射频率和入射角度以及凸起地形宽度对地表运动的影响. 数值结果表明,凸起地形对入射SV波具有显著的放大作用, 该放大作用可达自由场的4倍以上. 当入射频率较低时, 位移峰值多位于半空间地表; 而当入射频率较高时,位移峰值多位于凸起地形表面.      

Equilibrium modeling of the beach profile on a macrotidal embayed low tide terrace beach

Eleven-year long time series of monthly beach profile surveys and hourly incident wave conditions are analyzed for a macrotidal Low Tide Terrace beach. The lower intertidal zone of the beach has a pluriannual cycle, whereas the upper beach profile has a predominantly seasonal cycle. An equilibrium model is applied to study the variation of the contour elevation positions in the intertidal zone as a function of the wave energy, wave power, and water level. When forcing the model with wave energy, the predictive ability of the equilibrium model is around 60% in the upper intertidal zone but decreases to 40% in the lower intertidal zone. Using wave power increases the predictive ability up to 70% in both the upper and lower intertidal zones. However, changes around the inflection point are not well predicted. The equilibrium model is then extended to take into account the effects of the tide level. The initial results do not show an increase in the predictive capacity of the model, but do allow the model free parameters to represent more accurately the values expected in a macrotidal environment. This allows comparing the empirical model calibration in different tidal environment. The interpretation of the model free parameter variation across the intertidal zone highlights the behavior of the different zones along the intertidal beach profile. This contributes to a global interpretation of the four model parameters for beaches with different tidal ranges, and therefore to a global model applicable at a wide variety sites.     

Long term variations in the mesospheric mean flow observed at Grahamstown (South Africa) and Adelaide (Australia)

The seasonal and interannual behaviour of monthly mean winds at a height of 90 km recorded at Grahamstown (33.3°S, 26.5°E) and Adelaide (34.5°S, 138.5°E) between 1987 and 1994 are compared. The zonal wind is found to be consistently stronger at Grahamstown and is always eastward, whereas at Adelaide it sometimes reverses. Maxima tend to occur near the solstices, the primary maximum during summer at Grahamstown, in agreement with satellite results, and during winter at Adelaide. The meridional wind also tends to be stronger at Grahamstown, but at both stations is predominantly northward with a maximum in summer and generally not as strong as the zonal component. This seasonal behaviour is reasonably well understood in terms of the interaction of the mean flow with gravity waves propagating up from below, with coriolis forces also playing an important role in the case of the meridional wind. Satellite observations do not generally support the idea that longitudinal differences between the stations could be attributed to the presence of a tropospheric/stratospheric stationary wave. It is suggested that these differences are more probably associated with local effects. Interannual zonal wind patterns at the two sites are similar over the summer months but are less well correlated during the rest of the year. The underlying causes of this variability are not well understood but are most probably global in nature, at least during the summer.