Prediction of river water temperature: a comparison between a new family of hybrid models and statistical approaches

River water temperature is a key physical variable controlling several chemical, biological and ecological processes. Its reliable prediction is a main issue in many environmental applications, which however is hampered by data scarcity, when using data‐demanding deterministic models, and modelling limitations, when using simpler statistical models. In this work we test a suite of models belonging to air2stream family, which are characterized by a hybrid formulation that combines a physical derivation of the key equation with a stochastic calibration of parameters. The air2stream models rely solely on air temperature and streamflow, and are of similar complexity as standard statistical models. The performances of the different versions of air2stream in predicting river water temperature are compared with those of the most common statistical models typically used in the literature. To this aim, a dataset of 38 Swiss rivers is used, which includes rivers classified into four different categories according to their hydrological characteristics: low‐land natural rivers, lake outlets, snow‐fed rivers and regulated rivers. The results of the analysis provide practical indications regarding the type of model that is most suitable to simulate river water temperature across different time scales (from daily to seasonal) and for different hydrological regimes. A model intercomparison exercise suggests that the family of air2stream hybrid models generally outperforms statistical models, while cross‐validation conducted over a 30‐year period indicates that they can be suitably adopted for long‐term analyses. Copyright © 2016 John Wiley & Sons, Ltd.     

Long‐term monitoring of a mercury contaminated estuary (Ria de Aveiro,Portugal): the effect of weather events and management in mercury transport

The main aim of this research was to assess the mercury transport from an estuarine basin with a background of anthropogenic contamination during a spring tidal cycle (year 2009) and compare it with two previous tidal cycles (years 1994 and 1999), as part of a long‐term monitoring program. Results showed that effective mercury transport occurs both in the dissolved and particulate fractions (0.18 and 0.20 kg per tidal cycle, respectively), and despite an overall decrease in environmental contamination, results more than double previous findings on particulate transport in the system. These findings result essentially from changes in the tidal prism (net export of 2 million m³ of water), given that both dissolved and particulate concentrations did not increase over time. Hydrodynamic simulations were performed to evaluate the effect of physical disturbance (dredging) and weather events (increased freshwater flow) in these processes, and results suggest the increased freshwater flow into the system as the main forcing function for the mercury transport increment. These results highlight the importance of long‐term monitoring programs, since despite an overall improvement in local contamination levels, the enhancement of transport processes through hydrological changes increases environmental pressure away from the contamination source. Copyright © 2012 John Wiley & Sons, Ltd.     

Modeling deepwater seabed steady-state thermal fields around buried pipeline including trenching and backfill effects

Deepwater pipelines are designed to transport mixtures of oil and gas, and their associated impurities at wellhead temperatures that can be in excess of 149 °C (∼300 °F or 422 K) while the external temperature maybe in the range of 5 °C (∼41 °F or 278 K). Depending on the circumstances these pipelines may be buried for physical protection or for additional thermal insulation using robotic trenching equipment. This results in a complex cut and backfill geometry in the seafloor in addition to altering the thermal properties of the backfill. A two-dimensional boundary element model was developed specifically to address to investigate the local steady-state thermal field in the near field of the pipeline. The model allows one to account for the complex geometries in the near field associated with this burial technique, site-specific multi-layered soil conditions and the seawater adjacent to the seafloor. A parametric study was preformed to evaluate effects of the thermal power loss, burial depth, pipe diameter and soil thermal conductivity on the thermal field in the near field of a buried pipeline. The numerical examples illustrate the influence of the backfill thermal property on the temperature at the pipe wall, that the pipe diameter controls the required output thermal power needed to maintain the desired pipe wall temperature, and the importance of pipeline burial depth on seabed temperature distribution above the pipeline.     

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.     

Simulation of site-scale water fluxes in desert and natural oasis ecosystems of the arid region in Northwest China

The study of water fluxes is important to better understand hydrological cycles in arid regions. Data-driven machine learning models have been recently applied to water flux simulation. Previous studies have built site-scale simulation models of water fluxes for individual sites separately, requiring a large amount of data from each site and significant computation time. For arid areas, there is no consensus as to the optimal model and variable selection method to simulate water fluxes. Using data from seven flux observation sites in the arid region of Northwest China, this study compared the performance of random forest (RF), support vector machine (SVM), back propagation neural network (BPNN), and multiple linear regression (MLR) models in simulating water fluxes. Additionally, the study investigated inter-annual and seasonal variation in water fluxes and the dominant drivers of this variation at different sites. A universal simulation model for water flux was built using the RF approach and key variables as determined by MLR, incorporating data from all sites. Model performance of the SVM algorithm (R² = 0.25–0.90) was slightly worse than that of the RF algorithm (R² = 0.41–0.91); the BPNN algorithm performed poorly in most cases (R² = 0.15–0.88). Similarly, the MLR results were limited and unreliable (R² = 0.00–0.66). Using the universal RF model, annual water fluxes were found to be much higher than the precipitation received at each site, and natural oases showed higher fluxes than desert ecosystems. Water fluxes were highest during the growing season (May–September) and lowest during the non-growing season (October–April). Furthermore, the dominant drivers of water flux variation were various among different sites, but the normalized difference vegetation index (NDVI), soil moisture and soil temperature were important at most sites. This study provides useful insights for simulating water fluxes in desert and oasis ecosystems, understanding patterns of variation and the underlying mechanisms. Besides, these results can make a contribution as the decision-making basis to the water management in desert and oasis ecosystems.     

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.     

Role of horizontal density advection in seasonal deepening of the mixed layer in the subtropical Southeast Pacific

The mechanisms behind the seasonal deepening of the mixed layer(ML) in the subtropical Southeast Pacific were investigated using the monthly Argo data from 2004 to 2012. The region with a deep ML(more than 175 m) was found in the region of(22?–30?S, 105?–90?W), reaching its maximum depth(～200 m) near(27?–28?S, 100?W) in September. The relative importance of horizontal density advection in determining the maximum ML location is discussed qualitatively. Downward Ekman pumping is key to determining the eastern boundary of the deep ML region. In addition, zonal density advection by the subtropical countercurrent(STCC) in the subtropical Southwest Pacific determines its western boundary, by carrying lighter water to strengthen the stratification and form a "shallow tongue" of ML depth to block the westward extension of the deep ML in the STCC region. The temperature advection by the STCC is the main source for large heat loss from the subtropical Southwest Pacific. Finally, the combined effect of net surface heat flux and meridional density advection by the subtropical gyre determines the northern and southern boundaries of the deep ML region: the ocean heat loss at the surface gradually increases from 22?S to 35?S, while the meridional density advection by the subtropical gyre strengthens the stratification south of the maximum ML depth and weakens the stratification to the north. The freshwater flux contribution to deepening the ML during austral winter is limited. The results are useful for understanding the role of ocean dynamics in the ML formation in the subtropical Southeast Pacific.     

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

利用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月青藏高原地面加热场强度增强(减弱),有利于(不利于)来自北方的冷空气和南方的暖湿气流在东北上空交汇,且上升气流增强(减弱),水汽输送充沛(减少),导致该地区春季首场透雨出现的时间偏早(晚)。     

1971—2010年京津冀大城市热岛效应多时间尺度分析

利用1971—2010年均一化的京津冀区域逐日气温数据与质量控制后的2011年自动站逐时气温数据,分析了北京、天津和石家庄热岛效应的多尺度时间变化特征。结果表明,三个城市平均、最高和最低气温的热岛效应呈非对称性特征,最强为最低气温的热岛效应,其次为平均气温的热岛效应,最弱为最高气温的热岛效应。北京平均气温的热岛效应最强,其次为天津,石家庄相对较弱,石家庄平均气温的热岛效应近40年呈显著上升趋势,每10年达0.13℃。石家庄最高气温的热岛效应最强,其次为北京,最小为天津,近40年北京最高气温的热岛效应呈缓慢上升趋势,每10年增加0.06℃,石家庄变化不明显,天津呈微弱下降趋势。最低气温的热岛在北京最强,其次为天津,最小为石家庄,近40年最低气温热岛效应天津呈明显上升趋势,每10年增加0.18℃,其次为石家庄,北京呈微弱下降趋势。三个城市的平均气温、最低气温的热岛效应季节变化通常表现为夏季较弱,冬季最强。三个城市最高气温的热岛效应季节变化差异较大,北京10月热岛效应最弱,其他月份变化不大;天津热岛效应6月最弱,在1或12月最强;石家庄4和5月热岛效应最强,10月热岛效应最弱。由2011年自动站数据得到的平均气温热岛效应与1971—2010年的40年平均得到的平均气温的热岛效应季节变化具有类似的规律。2011年自动站热岛效应在一天中表现为白天热岛强度较低,而夜间热岛强度较高。