Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

A long‐lived metamorphic history in the contact aureole of the Mooselookmeguntic pluton revealed by in situ dating of monazite grains preserved as inclusions in staurolite porphyroblasts

Emplacement of the Mooselookmeguntic pluton, located in the western Maine region of the northern Appalachians, was thought to have occurred towards the end of the Acadian deformation at around 370 Ma. Crystallization ages from different parts of the pluton suggest a more sequential emplacement history over a period of c. 20 Myr. Foliation inflection/intersection axes (FIAs) within porphyroblasts from its aureole reveal at least five periods of garnet and staurolite growth. The orientation of FIAs in both garnet and staurolite porphyroblasts trend successively from ESE–WNW, NNW–SSE, E–W, ENE–WSW to NE–SW. Electron probe microanalysis dating of monazite grains included in staurolite porphyroblasts containing one of these five periods of FIA development reveals a succession of apparent ages from 410 Ma to 345 Ma. A similar spread of crystallization ages can be observed for plutons from Maine and adjacent regions. This succession indicates that deformation and metamorphism began well before and continued long after what is classically regarded as the Acadian orogeny. The thermal structure of the orogen progressively evolved to enable pluton emplacement, and it continued to develop afterwards with magmatic fluids still forming at depth.     

Structural geometry and controls on basement‐involved deformation in the northern Flinders Ranges,Adelaide Fold Belt,South Australia

In the northern Flinders Ranges, Neoproterozoic and Cambrian sedimentary rocks were deformed and variably metamorphosed during the ca 500 Ma Cambro‐Ordovician Delamerian Orogeny. Balanced and restored structural sections across the northern Flinders Ranges show shortening of about 10–20%. Despite the presence of suitable evaporitic detachment horizons at the basement‐cover interface, the structural style is best interpreted to be thick‐skinned involving basement with only a minor proportion of the overall shortening accommodated along stratigraphically controlled detachments. Much of the contractional deformation was localised by the inversion of former extensional faults such as the Norwest and Paralana Faults, which both controlled the deposition of Neoproterozoic cover successions. As such, both faults represent major, long‐lived structures which effectively define the present boundaries of the northern Flinders Ranges with the Gawler Craton to the west and the Curnamona Craton to the east. The most intense deformation, which resulted in exhumation of the basement along the Paralana Fault to form the Mt Painter and Babbage Inliers, coincides with extremely high heat flows related to extraordinarily high heat‐production rates in the basement rocks. High heat flow in the northern Flinders Ranges suggests that the structural style not only reflects the pre‐Delamerian basin architecture but is also a consequence of the reactivation of thermally perturbed, weakened basement.     

Variability in bedform morphology and kinematics with transport stage

Bedform geometry is widely recognized to be a function of transport stage. Bedform aspect ratio (height/length) increases with transport stage, reaches a maximum, then decreases as bedforms washout to a plane bed. Bedform migration rates are also linked to bedform geometry, in so far as smaller bedforms in coarser sediment tend to migrate faster than larger bedforms in finer sediment. However, how bedform morphology (height, length and shape) and kinematics (translation and deformation) change with transport stage and suspension have not been examined. A series of experiments is presented where initial flow depth and grain size were held constant and the transport stage was varied to produce bedload dominated, mixed‐load dominated and suspended‐load dominated conditions. The results show that the commonly observed pattern in bedform aspect ratio occurs because bedform height increases then decreases with transport stage, against a continuously increasing bedform length. Bedform size variability increased with transport stage, leading to less uniform bedform fields at higher transport stage. Total translation‐related and deformation‐related sediment fluxes all increased with transport stage. However, the relative contribution to the total flux changed. At the bedload dominated stage, translation‐related and deformation‐related flux contributed equally to the total flux. As the transport stage increased, the fraction of the total load contributed by translation increased and the fraction contributed by deformation declined because the bedforms got bigger and moved faster. At the suspended‐load dominated transport stage, the deformation flux increased and the translation flux decreased as a fraction of the total load, approaching one and zero, respectively, as bedforms washed out to a plane bed.     

Upper-Tropospheric Environment–Tropical Cyclone Interactions over the Western North Pacific: A Statistical Study

Based on 25-year(1987–2011) tropical cyclone(TC) best track data, a statistical study was carried out to investigate the basic features of upper-tropospheric TC–environment interactions over the western North Pacific. Interaction was defined as the absolute value of eddy momentum flux convergence(EFC) exceeding 10 m s~(-1)d~(-1). Based on this definition, it was found that 18% of all six-hourly TC samples experienced interaction. Extreme interaction cases showed that EFC can reach~120 m s~(-1)d~(-1) during the extratropical-cyclone(EC) stage, an order of magnitude larger than reported in previous studies.Composite analysis showed that positive interactions are characterized by a double-jet flow pattern, rather than the traditional trough pattern, because it is the jets that bring in large EFC from the upper-level environment to the TC center. The role of the outflow jet is also enhanced by relatively low inertial stability, as compared to the inflow jet. Among several environmental factors, it was found that extremely large EFC is usually accompanied by high inertial stability, low SST and strong vertical wind shear(VWS). Thus, the positive effect of EFC is cancelled by their negative effects. Only those samples during the EC stage, whose intensities were less dependent on VWS and the underlying SST, could survive in extremely large EFC environments, or even re-intensify. For classical TCs(not in the EC stage), it was found that environments with a moderate EFC value generally below ~25 m s~(-1)d~(-1) are more favorable for a TC's intensification than those with extremely large EFC.     

Effects of surface flux parameterization on the numerically simulated intensity and structure of Typhoon Morakot (2009)

The effects of surface flux parameterizations on tropical cyclone(TC) intensity and structure are investigated using the Advanced Research Weather Research and Forecasting(WRF-ARW) modeling system with high-resolution simulations of Typhoon Morakot(2009).Numerical experiments are designed to simulate Typhoon Morakot(2009) with different formulations of surface exchange coefficients for enthalpy(C_K) and momentum(C_D) transfers,including those from recent observational studies based on in situ aircraft data collected in Atlantic hurricanes.The results show that the simulated intensity and structure are sensitive to C_K and C_D,but the simulated track is not.Consistent with previous studies,the simulated storm intensity is found to be more sensitive to the ratio of C_K/C_D than to C_K or C_D alone.The pressure-wind relationship is also found to be influenced by the exchange coefficients,consistent with recent numerical studies.This paper emphasizes the importance of C_D and C_K on TC structure simulations.The results suggest that C_D and C_K have a large impact on surface wind and flux distributions,boundary layer heights,the warm core,and precipitation.Compared to available observations,the experiment with observed C_D and C_K generally simulated better intensity and structure than the other experiments,especially over the ocean.The reasons for the structural differences among the experiments with different C_D and C_K setups are discussed in the context of TC dynamics and thermodynamics.     

Long-Term Variability of the Wind Field over the Indian Ocean Based on ERA-Interim Reanalysis

A detailed study of long-term variability of winds using 30 years of data from the European Centre for Medium-range Weather Forecasts global reanalysis (ERA-Interim) over the Indian Ocean has been carried out by partitioning the Indian Ocean into six zones based on local wind extrema. The trend of mean annual wind speed averaged over each zone shows a significant increase in the equatorial region, the Southern Ocean, and the southern part of the trade winds. This indicates that the Southern Ocean winds and the southeast trade winds are becoming stronger. However, the trend for the Bay of Bengal is negative, which might be caused by a weakening of the monsoon winds and northeast trade winds. Maximum interannual variability occurs in the Arabian Sea due to monsoon activity; a minimum is observed in the subtropical region because of the divergence of winds. Wind speed variations in all zones are weakly correlated with the Dipole Mode Index (DMI). However, the equatorial Indian Ocean, the southern part of the trade winds, and subtropical zones show a relatively strong positive correlation with the Southern Oscillation Index (SOI), indicating that the SOI has a zonal influence on wind speed in the Indian Ocean. Monsoon winds have a decreasing trend in the northern Indian Ocean, indicating monsoon weakening, and an increasing trend in the equatorial region because of enhancement of the westerlies. The negative trend observed during the non-monsoon period could be a result of weakening of the northeast trade winds over the past few decades. The mean flux of kinetic energy of wind (FKEW) reaches a minimum of about 100?W?m^?2 in the equatorial region and a maximum of about 1500?W?m^?2 in the Southern Ocean. The seasonal variability of FKEW is large, about 1600?W?m^?2, along the coast of Somalia in the northern Indian Ocean. The maximum monthly variability of the FKEW field averaged over each zone occurs during boreal summer. During the onset and withdrawal of monsoon, FKEW is as low as 50?W?m^?2. The Southern Ocean has a large variation of about 1280?W?m^?2 because of strong westerlies throughout the year.     

东北地区春季首场透雨的变化特征及与青藏高原地面加热场强度的关系

利用1958—2012年4—5月东北地区(39°~55°N、118°~135°E)101个站点逐日降水资料、青藏高原地区(25°~40°N、73.75°~103.75°E)JRA-55的地面感热和潜热通量月平均再分析资料以及NCEP/NCAR-I大气环流场的月平均再分析资料,分析了春播期首场透雨出现日期的时空变化特征及其与透雨量和播种期降水量间的关系,以及对青藏高原地面加热场强度异常的响应及其可能机制。结果表明:透雨日期自1958年以来在东北地区的西北和东南大部分区域呈现略微偏晚的趋势;中部有略微偏早的趋势。春播期首场透雨出现时间偏早(晚)的地方,首场透雨量小(大),春播期总降水量多(少)。同时,4月青藏高原地面加热场强度增强(减弱),有利于(不利于)来自北方的冷空气和南方的暖湿气流在东北上空交汇,且上升气流增强(减弱),水汽输送充沛(减少),导致该地区春季首场透雨出现的时间偏早(晚)。     

海河流域夏季降水与北半球环流异常的相关

使用1981—2013年6—8月的月平均海河流域站点降水资料以及NCEP/NCAR再分析资料,采用SVD方法分析了北半球高度场环流异常与海河流域夏季降水的关系。结果表明夏季海河流域的一致型降水异常主要受热带地区高度场异常的影响,而北半球中高纬的高度异常分布对海河流域东南部降水异常有重要作用。热带及中高纬地区的高度场异常与副热带西风急流和流域上空的水汽输送有一定的对应关系,急流强度和位置的变化以及水汽辐合异常造成海河流域降水异常,尤其中高纬高度异常的这种对应关系更加明显,并且两个地区的高度异常分布对应的海河流域异常水汽输送有不同来源。