Identification of CO2 disposal locations in an ocean general circulation model of the North Pacific

A basinwide ocean general circulation model of the North Pacific is used to identify which location is more efficient for ocean CO2 sequestration in the North Pacific. Four injection depths at each one of fifteen locations are chosen. In terms of effectiveness index (EI) and escape factor (EF), it is clear that the effectiveness increases with increasing latitude at the end of the 50 a injection period. Sitebysite differences in the EI can be over 9% for the 1 000 m injection depth in the western North Pacific at the end of 50 a of continuous injection. The difference is much larger for the 500 m injection. The difference decreases with increasing injection depth. However, the sitebysite difference is small for the injection in the eastern North Pacific. The sequestration is more efficient for the injection in the east than in the west. For the 500 m injection depth, the difference in effectiveness between the west and the east is over 10% at the end of 50 a injection period. The largest concentration of sequestered CO2 increases with increasing injection depth. For the injection in both the western and central North Pacific, the largest exchange flux always appears to be at about 42°N, 150°E, whereas for the injection in the eastern area the large flux appears to be in the equatorial region (120°W).     

On the importance of soil moisture in calibration of rainfall–runoff models: two case studies

Streamflow modelling results from the GR4H and PDM hydrological models were evaluated in two Australian sub-catchments, using (1) calibration to streamflow and (2) joint-calibration to streamflow and soil moisture. Soil moisture storage in the models was evaluated against soil moisture observations from field measurements. The PDM had the best performance in terms of both streamflow and soil moisture estimations during the calibration period, but was outperformed by GR4H during validation. It was also shown that the soil moisture estimation was improved significantly by joint-calibration for the case where streamflow and soil moisture estimations were poor. In other cases, addition of the soil moisture constraint did not degrade the results. Consequently, it is recommended that GR4H be used, in preference to the PDM, in the foothills of the Murrumbidgee catchment or other Australian catchments with semi-arid to sub-humid climate, and that soil moisture data be used in the calibration process.     

Numerical simulation of the detachment dynamics in North China Basin

ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｄｅｔａｃｈｍｅｎｔｄｙｎａｍｉｃｓｉｎＮｏｒｔｈＣｈｉｎａＢａｓｉｎＤｏｎｇ－Ｎｉｎｇ；ＺＨＡＮＧ（张东宁）ａｎｄＲｏｎｇ－ＳｈｅｎｇＺＥＮＧ（曾融生）（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，Ｓ...     

The role of salinity in fluvio-deltaic morphodynamics: A long-term modelling study

Salinity difference between terrestrial river discharge and oceanic tidal water plays a role in modifying the local flow field and, as a consequence, estuarine morphodynamics. Although widely recognized, recent numerical studies exploring the long-term morphological evolution of river-influenced estuaries with two-dimensional, depth-averaged models have mostly neglected salinity. Using a three-dimensional morphodynamic model, we aim to gain more insight into the effect of salinity on the morphodynamics of fluvio-deltaic systems. Model results indicate that the resultant estuarine morphology established after 600 years differs remarkably when a salinity gradient is included. A fan-shaped river-mouth delta exhibits less seaward expansion and is cut through by narrower channels when salinity is included. The inclusion of salinity tends to generate estuarine circulation, which favours landward sediment transport and hence limits the growth of the delta while enhancing the development of intertidal areas. The formation of deltaic channel–shoal patterns resulting from morphodynamic evolution tends to strengthen salinity stratification, which is characterized by an increased gradient Richardson number. The direction of the depth-averaged residual sediment transport over a tide may be opposite to the direction of residual velocity, indicating the significant influence of baroclinic effects on the net sediment transport direction (and hence morphological change). The effect of salinity on morphological evolution becomes less profound when the strength of tidal or fluvial forcing is dominant over the other. The effects of sediment type and flocculation, which are particularly important when salinity gradients are present, are also discussed. Overall, this study highlights that neglecting salinity to simulate long-term estuarine morphodynamics requires more careful justification, particularly when the environment is characterized by fine sediment types (favouring suspended transport), and relatively large river discharge and estuarine depth (favouring baroclinic effects). © 2020 John Wiley & Sons, Ltd.     

A patch hierarchy approach to modeling surface and subsurface hydrology in complex flood‐plain environments

Geomorphology interacts with surface‐ and ground‐water hydrology across multiple spatial scales. Nonetheless, hydrologic and hydrogeologic models are most commonly implemented at a single spatial scale. Using an existing hydrogeologic computer model, we implemented a simple hierarchical approach to modeling surface‐ and ground‐water hydrology in a complex geomorphic setting. We parameterized the model to simulate ground‐ and surface‐water ?ow patterns through a hierarchical, three‐dimensional, quantitative representation of an anabranched montane alluvial ?ood plain (the Nyack Flood Plain, Middle Fork Flathead River, Montana, USA). Comparison of model results to ?eld data showed that the model provided reasonable representations of spatial patterns of aquifer recharge and discharge, temporal patterns of ?ood‐water storage on the ?ood plain, and rates of ground‐water movement from the main river channel into a large lateral spring channel on the ?ood plain, and water table elevation in the alluvial aquifer. These results suggest that a hierarchical approach to modeling ground‐ and surface‐water hydrology can reproduce realistic patterns of surface‐ and ground‐water ?ux on alluvial ?ood plains, and therefore should provide an excellent ‘quantitative laboratory’ for studying complex interactions between geomorphology and hydrology at and across multiple spatial scales. Copyright © 2004 John Wiley & Sons, Ltd.     

水文信息处理系统研究

水文信息处理系统是在Windows环境下、用Visual Basic 6．0高级语言研制开发的智能水文资料整编软件，引用数据库存储信息，调用Excel输出整编成果。该软件实现了水位、流量、泥沙等要素关系线的自动生成和人工辅助定线，并完成全部整编计算工作，结束了我国水文资料电算整编工作中依靠手工定线整编的历史，解决了计算机定线、推流和成果输出一次完成的技术难题，真正实现了水文资料整编计算机化。     

Flood risk analyses—how detailed do we need to be?

Applied flood risk analyses, especially in urban areas, very often pose the question how detailed the analysis needs to be in order to give a realistic figure of the expected risk. The methods used in research and practical applications range from very basic approaches with numerous simplifying assumptions up to very sophisticated, data and calculation time demanding applications both on the hazard and on the vulnerability part of the risk. In order to shed some light on the question of required model complexity in flood risk analyses and outputs sufficiently fulfilling the task at hand, a number of combinations of models of different complexity both on the hazard and on the vulnerability side were tested in a case study. The different models can be organized in a model matrix of different complexity levels: On the hazard side, the approaches/models selected were (A) linear interpolation of gauge water levels and intersection with a digital elevation model (DEM), (B) a mixed 1D/2D hydraulic model with simplifying assumptions (LISFLOOD-FP) and (C) a Saint-Venant 2D zero-inertia hyperbolic hydraulic model considering the built environment and infrastructure. On the vulnerability side, the models used for the estimation of direct damage to residential buildings are in order of increasing complexity: (I) meso-scale stage-damage functions applied to CORINE land cover data, (II) the rule-based meso-scale model FLEMOps+ using census data on the municipal building stock and CORINE land cover data and (III) a rule-based micro-scale model applied to a detailed building inventory. Besides the inundation depths, the latter two models consider different building types and qualities as well as the level of private precaution and contamination of the floodwater. The models were applied in a municipality in east Germany, Eilenburg. It suffered extraordinary damage during the flood of August 2002, which was well documented as were the inundation extent and depths. These data provide an almost unique data set for the validation of flood risk analyses. The analysis shows that the combination of the 1D/2D model and the meso-scale damage model FLEMOps+ performed best and provide the best compromise between data requirements, simulation effort, and an acceptable accuracy of the results. The more detailed approaches suffered from complex model set-up, high data requirements, and long computation times.     

Middle Proterozoic–early Palaeozoic evolution of central Baltoscandian intracratonic basins: evidence for asthenospheric diapirs

The three intracratonic sedimentary basins located in central Baltoscandinavia, namely the Bothnian Gulf basin, the Bothnian Sea basin and the Baltic basin, developed in response to Middle Proterozoic and Late Proterozoic tectonic events, separated in time by about 800 Ma. Only the Baltic basin was subsequently affected by Caledonian orogenesis and Mesozoic rifting. Crustal extension was minor or did not take place during the Proterozoic basin evolution phases. However, according to the Moho topography, crustal thinning did take place. This was probably a result of subcrustal magmatism. On a craton-wide scale, the ages of granitoids, which intruded during the Middle Proterozoic basin formation, generally decrease from east to west. This fact, combined with the evidence provided by mantle-derived flood basalt magmatism, points to a moving asthenospheric diapir as the cause for basin development. Asthenospheric upwelling was probably also responsible for the second, Late Proterozoic, basin evolution phase, as evidenced by the lack of crustal thinning and extension, and the occurrence of tholeiitic intrusions. In addition, a Late Proterozoic thermally induced palaeo-high, located at about the position of the intracratonic basins, is compatible with indications from glaciations. As the ages of Late Proterozoic intracratonic basins also decrease from east to west across the craton, the location of asthenospheric diapirism during this time interval was also moving. For the Fennoscandian lithosphere, the presence of fundamental lithospheric weakness zones (e.g. terrane boundaries) might be an explanation for the formation of two generations of basins originating from asthenospheric upwelling at about the same location in the Fennoscandian Shield. The spacing and size of the Proterozoic intracratonic basins suggest that the asthenospheric diapirism was not deep seated. Therefore, sublithospheric convective processes might be the cause for the asthenospheric upwellings. Such processes are related to Rayleigh–Taylor instabilities in the sublithospheric mantle. Emplacement of an asthenospheric diapir causes a thermal bulge at the surface of the lithosphere. Modelling results demonstrate that erosion of the surficial high, succeeded by cooling of the lithosphere, can explain the accumulation of early Palaeozoic sediments in the Bothnian Sea basin, taking into account post-Ordovician vertical and lateral erosion of the basin fill.     

Constraints on the timing and conditions of high‐grade metamorphism,charnockite formation and fluid–rock interaction in the Trivandrum Block,southern India

Incipient charnockites have been widely used as evidence for the infiltration of CO₂‐rich fluids driving dehydration of the lower crust. Rocks exposed at Kakkod quarry in the Trivandrum Block of southern India allow for a thorough investigation of the metamorphic evolution by preserving not only orthopyroxene‐bearing charnockite patches in a host garnet–biotite felsic gneiss, but also layers of garnet–sillimanite metapelite gneiss. Thermodynamic phase equilibria modelling of all three bulk compositions indicates consistent peak‐metamorphic conditions of 830–925 °C and 6–9 kbar with retrograde evolution involving suprasolidus decompression at high temperature. These models suggest that orthopyroxene was most likely stabilized close to the metamorphic peak as a result of small compositional heterogeneities in the host garnet–biotite gneiss. There is insufficient evidence to determine whether the heterogeneities were inherited from the protolith or introduced during syn‐metamorphic fluid flow. U–Pb geochronology of monazite and zircon from all three rock types constrains the peak of metamorphism and orthopyroxene growth to have occurred between the onset of high‐grade metamorphism at c. 590 Ma and the onset of melt crystallization at c. 540 Ma. The majority of metamorphic zircon growth occurred during protracted melt crystallization between c. 540 and 510 Ma. Melt crystallization was followed by the influx of aqueous, alkali‐rich fluids likely derived from melts crystallizing at depth. This late fluid flow led to retrogression of orthopyroxene, the observed outcrop pattern and to the textural and isotopic modification of monazite grains at c. 525–490 Ma.     

Morphodynamic assessment of side channel systems using a simple one‐dimensional bifurcation model and a comparison with aerial images

Side channel construction is a common intervention applied to increase a river's conveyance capacity and to increase its ecological value. Past modelling efforts suggest two mechanisms affecting the morphodynamic change of a side channel: (1) a difference in channel slope between the side channel and the main channel and (2) bend flow just upstream of the bifurcation. The objective of this paper was to assess the conditions under which side channels generally aggrade or degrade and to assess the characteristic timescales of the associated morphological change. We use a one‐dimensional bifurcation model to predict the development of side channel systems and the characteristic timescale for a wide range of conditions. We then compare these results to multitemporal aerial images of four side channel systems. We consider the following mechanisms at the bifurcation to be important for side channel development: sediment diversion due to the bifurcation angle, sediment diversion due to the transverse bed slope, partitioning of suspended load, mixed sediment processes such as sorting at the bifurcation, bank erosion, deposition due to vegetation, and floodplain sedimentation. There are limitations to using a one‐dimensional numerical model as it can only account for these mechanisms in a parametrized manner, but the model reproduces general behaviour of the natural side channels until floodplain‐forming processes become important. The main result is a set of stability diagrams with key model parameters that can be used to assess the development of a side channel system and the associated timescale, which will aid in the future design and maintenance of side channel systems. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.