Locating regional seismic events with global optimization based on interval arithmetic

More than 1000 seismic events in northern Europe at distances of up to 400  km from the detecting network are located using an optimization method in which the global minimum of the traveltime function residuals is searched for using an Interval Arithmetic (IA) method. Epicentres are determined using P waves detected by the Finnish national seismic network: up to 15 stations were used in the analysis. The IA results coincide with locations provided by the University of Helsinki bulletins with a median location bias of 7.6  km.
  A second data set of 59 explosions in the Siilinjärvi mine in central Finland was examined in detail, because the locations of the explosions were known exactly. In this case, the median difference of IA locations was 3.8  km from the average location of mine explosions, while all 59 events were located within 9  km of the 'true' epicentres. The corresponding median error of the University of Helsinki locations was smaller (3.2  km), but some Helsinki locations were well over 10  km from the mine. The convergence towards the global optimum using interval arithmetic was fast when compared with the conventional least-squares approaches for epicentre determinations.     

Intra-City Variation in Urban Morphology and Turbulence Structure in Helsinki, Finland

Most atmospheric boundary-layer theories are developed over vegetative surfaces and their applicability at urban sites is questionable. Here, we study the intra-city variation of turbulence characteristics and the applicability of boundary-layer theory using building-morphology data across Helsinki, and eddy-covariance data from three sites: two in central Helsinki (400 m apart) and one 4 km away from the city centre. The multi-site measurements enable the analysis of the horizontal scales at which quantities that characterize turbulent transport vary: (i) Roughness characteristics vary at a 10-m scale, and morphometric estimation of surface-roughness characteristics is shown to perform better than the often used rule-of-thumb estimates (average departures from the logarithmic wind profile are 14 and 44 %, respectively). (ii) The drag coefficient varies at a 100-m scale, and we provide an updated parametrization of the drag coefficient as a function of z/z _H (the ratio of the measurement height to the mean building height). (iii) The transport efficiency of heat, water vapour and CO₂ is shown to be weaker the more heterogeneous the site is, in terms of sources and sinks, and strong scalar dissimilarity is observed at all sites. (iv) Atmospheric stability varies markedly even within 4 km across the city: the median difference in nocturnal sensible heat fluxes between the three sites was over 50W m^?2. Furthermore, (v) normalized power spectra and cospectra do not vary between sites, and they follow roughly the canonical theory as developed over vegetated terrain.     

Ejective and Sweeping Motions Above a Peatland and Their Role in Relaxed-Eddy-Accumulation Measurements and Turbulent Transport Modelling

The three turbulent velocity components, water vapour (\(\text {H}_2\text {O}\)), carbon dioxide (\(\text {CO}_{2}\)), and methane (\(\text {CH}_{4}\)) concentration fluctuations are measured above a boreal peatland and analyzed using conditional sampling and quadrant analysis. The overarching question to be addressed is to what degree lower-order cumulant expansion methods describe transport efficiency and the relative importance of ejections and sweeps to momentum, \(\text {CH}_{4}\), \(\text {CO}_{2}\) and \(\text {H}_2\text {O}\) fluxes across a range of atmospheric flow regimes. The patchy peatland surface creates distinctly different source and sink distributions for the three scalars in space and time thereby adding to the uniqueness of the set-up. The measured and modelled fractional contributions to the momentum flux show that sweep events dominate over ejections in agreement with prior studies conducted in the roughness sublayer. For scalar fluxes, ejections dominate the turbulent fluxes over sweeps. While ejective motions persist longer for momentum transport, sweeping events persist longer for all three scalars. Third-order cumulant expansions describe many of the results detailed above, and the results are surprising given the highly non-Gaussian distribution of \(\text {CH}_{4}\) turbulent fluctuations. Connections between the asymmetric contributions of sweeps and ejections and the flux-transport term arising in scalar turbulent-flux-budget closure are derived and shown to agree reasonably well with measurements. The proposed model derived here is much simpler than prior structural models used to describe laboratory experiments. Implications of such asymmetric contributions on, (i) the usage of the now proliferating relaxed-eddy-accumulation method in turbulent flux measurements, (ii) the constant-flux assumption, and (iii) gradient-diffusion closure models are presented.     

On the Factors Controlling the Snow Surface and 2-m Air Temperatures Over the Arctic Sea Ice in Winter

Factors controlling the magnitudes of, and short-term variations in, the potential temperatures of the snow surface and the air at the height of 2 m θ_S and θ_{2 m} over Arctic sea ice in winter are analysed. The study addresses the winters of 1986–1987 and 1987–1988, and is based on the temperature, wind, and cloud observations made by Russian drifting ice stations. It also relies on the ERA40 re-analyses of the European Centre for Medium-Range Weather Forecasts, which were utilised to calculate the lateral heat advection at the sites of the ice stations. The cloud cover and wind speed were more important than the heat advection in controlling the magnitudes of θ_{2 m} and θ_S, while on a time scale of 24 h, during steady forcing conditions, the heat advection was the most important factor affecting the changes in θ_S and θ_{2 m}. During changing conditions, and considering individual factors separately, the monthly mean 24-h temperature changes were less than ± 5 °C: the effect of the cloud cover was the largest, and that of the heat advection was the smallest. When simultaneous changes in the three factors were analysed, the seasonal mean temperature changes were even of the order of ±15 °C, with the strongest warming events exceeding 35 K in a single day. The difference θ_S − θ_{2 m} reached its lowest seasonal mean values during conditions of clear skies (−1.3 °C), light winds (−1.3 °C) and warm-air advection (−0.8 °C). θ_S and θ_{2 m} followed each other closely, even during major synoptic-scale temperature variations.     

Distribution of nitrate over sea-salt and soil derived particles — Implications from a field study

A method for determining the distribution of supermicrometer nitrate between size-segregated sea-salt and soil derived particles is presented. The analysis is based on field data from six measurements at a coastal site in southern Finland, and on a theoretical treatment taking into account the transfer of gaseous species onto particle surfaces and their subsequent reaction. Significant amounts of nitrate were found in both the particle types, with the fraction of nitrate associated with soil particles varying from 20–50% in the 1–2 m size to near 90% in particles larger than 10 m. Overall, the nitrate accumulation followed closely the relative abundances of these two particle types. Two overlapping modes in supermicron nitrate mass size distributions could be identified. The lower mode, associated with sea-salt, was located between the surface-area and volume distribution of sodium peaking at about 2–3 m of EAD. The upper mode peaked at 3–5 m and followed more closely the surface-area distribution of calcium in all samples. At our site, the accumulation of nitrate into both particle types was shown to be limited by an effective surface reaction rate rather than by gas-phase diffusion. This rate was estimated to be considerably larger for sea-salt particles. Strong evidence in support of the saturation of nitrate in sea-salt particles were obtained.     

How to Earn the Status of Honest Broker? Scientists’ Roles Facilitating the Political Water Supply Decision-Making Process

“Honest broker” refers to scientists facilitating complex political decisions from a neutral position, aiming to expand the scope of choice available instead of advocating a certain policy. This article reflects on the experiences of three scientists acting as honest brokers while carrying out a multicriteria evaluation to help politicians choose the water supply option for the city of Oulu, Finland. The process was considered successful and many politicians said that they have never before been that well informed when making a decision. The case primarily shows that besides the common knowledge brokerage roles targeting knowledge exchange, the honest brokerage approach relates more to managing social relationships and balancing political games and rapidly evolving discussions. Trustful relationships increased the perceived credibility of the brokers and ensured two-way communication. On the other hand, the brokers faced ethical dilemmas and difficulties of staying impartial while finding out confidential information related to political strategies.     

Contaminant exposure and effects in Baltic ringed and grey seals as assessed by biomarkers

The Baltic Sea ecosystem has suffered from a heavy pollutant load for more than three decades. Persistent organic pollutants (POPs) and heavy metals have been of most concern due to their persistence and toxic properties. Ringed seals (Phoca hispida baltica) and grey seals (Halichoerus grypus) living in the Baltic Sea have been suffering from pathological impairments, including reproductive disturbances, which have resulted in a depressed reproductive capacity. We investigated several biochemical parameters as potential biomarkers for exposure to and effects of the contaminant load in the Baltic seals. Seals from less polluted areas were used as reference material in terms of the pollution load. In both Baltic seal populations, the levels of some biochemical parameters diverged from those in the reference seals, and some of these showed a clear correlation with the individual contaminant load. Of the potential bioindicators, we propose cytochrome P4501A activity and vitamin E levels, in blubber or plasma, as exposure biomarkers for polychlorinated biphenyls (PCBs) and dichlorodiphenyltrichloroethane (DDT) in both species. The arylhydrocarbon receptor-mediated chemical-activated luciferase gene expression (CALUX) response reflects the whole PCB and DDT burden in ringed seals. Retinyl palmitate (vitamin A) levels showed a negative correlation with the individual POP load, and is proposed as potential effect biomarkers for the depletion of the vitamin A stores. As the nutritional levels of both vitamin A and E have an impact on the vitamin levels in the seals, more information on the dietary vitamin levels is needed before any conclusions can be drawn. As the relationship between biochemical parameters and contaminants varied between the two species, species-specific characteristics has to be considered when monitoring the health status and possible toxic effects of the contaminant load in ringed and grey seals.     

Damage–viscoplastic consistency model with a parabolic cap for rocks with brittle and ductile behavior under low‐velocity impact loading

This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical modelling of bit–rock fracture mechanisms in percussive drilling with a continuum approach

This paper presents a numerical method for continuum modelling of the dynamic bit–rock interaction process in percussive drilling. The method includes a constitutive model based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. The interaction between the drill bit and rock is modelled using contact mechanics by treating the bit as a rigid body. As the bit–rock interaction in percussive drilling is a transient event, the method is implemented in explicit dynamics FEM. The rock strength heterogeneity is characterized at the mesoscopic level statistically using the Weibull distribution. The bit–rock interaction is simulated under axisymmetric conditions using cylindrical and hemispherical buttons. The choice of the quite complex constitutive model accounting, e.g. for plastic compaction, viscoplastic shear and tensile failure along with induced damage and rate dependency is justified by numerical simulations. Moreover, the quasi‐static and dynamic cases are compared in plane strain simulations. Finally, some results clarifying the discrepancy of opinions found in the literature concerning the side (lateral) crack formation are obtained. Copyright © 2010 John Wiley & Sons, Ltd.