EPANET Input Files of New York tunnels and Pacific City used in a metamodel-based optimization study

Metamodels have proven be very useful when it comes to reducing the computational requirements of Evolutionary Algorithm-based optimization by acting as quick-solving surrogates for slow-solving fitness functions. The relationship between metamodel scope and objective function varies between applications, that is, in some cases the metamodel acts as a surrogate for the whole fitness function, whereas in other cases it replaces only a component of the fitness function. This paper presents a formalized qualitative process to evaluate a fitness function to determine the most suitable metamodel scope so as to increase the likelihood of calibrating a high-fidelity metamodel and hence obtain good optimization results in a reasonable amount of time. The process is applied to the risk-based optimization of water distribution systems; a very computationally-intensive problem for real-world systems. The process is validated with a simple case study (modified New York Tunnels) and the power of metamodelling is demonstrated on a real-world case study (Pacific City) with a computational speed-up of several orders of magnitude.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Mid-Pleistocene (800-400 ka) records of benthic and planktic foraminfer d18O and d13C, planktic foraminifer based SST, and IRD of IODP Site 306-U1314

Stable carbon and oxygen isotopes from benthic and planktic foraminifers, planktic foraminifer assemblages and ice rafted debris from the North Atlantic Site U1314 (Integrated Ocean Drilling Program Expedition 306) were examined to investigate orbital and millennial-scale climate variability in the North Atlantic and its impact on global circulation focusing on the development of glacial periods during the mid-Pleistocene (ca 800-400 ka). Glacial initiations were characterized by a rapid cooling (6-10 °C in less than 7 kyr) in the mean annual sea surface temperature (SST), increasing benthic d18O values and high benthic d13C values. The continuous increase in benthic d18O suggests a continuous ice sheet growth whereas the positive benthic d13C values indicate that the flow of the Iceland Scotland Overflow water (ISOW) was vigorous. Strong deep water formation in the Norwegian Greenland Sea promoted a high transfer of freshwater from the ocean to the continents. However, low SSTs at Site U1314 suggest a subpolar gyre cooling and freshening that may have reduced deep water formation in the Labrador Sea during glacial initiations. Once the 3.5 per mil threshold in the benthic d18O record was exceeded, ice rafting started and ice sheet growth was punctuated by millennial-scale waning events which returned to the ocean part of the freshwater accumulated on the continents. Ice-rafting events were associated with a rapid reduction in the ISOW (benthic d13C values dropped 0.5-1 per mil) and followed by millennial-scale warmings. The first two millennial-scale warm intervals of each glacial period reached interglacial temperatures and were particularly abrupt (6-10 °C in ~3 kyr). Subsequent millennial-scale warm events were cooler probably because the AMOC was rather reduced as suggested by the low benthic d13C values. These two abrupt warming events that occurred at early glacial periods were also observed in the Antarctic temperature and CO2 records, suggesting a close correlation between both Hemispheres. The comparison of the sea surface proxies with the benthic d18O record (as the Southern sign) indicates the presence of a millennial-scale seesaw pattern similar to that seen during the Last Glacial period.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.