Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact of Boundary-Layer Schemes,Boundary Conditions and Spin-Up

We evaluated the performance of the three-dimensional Weather Research and Forecasting (WRF) mesoscale model, specifically the performance of the planetary boundary-layer (PBL) parametrizations. For this purpose, Cabauw tower observations were used, with the study extending beyond the third GEWEX Atmospheric Boundary-Layer Study (GABLS3) one-dimensional model intercomparison. The WRF model (version 3.4.1) contains 12 different PBL parametrizations, most of which have been only partially evaluated. The GABLS3 case offers a clear opportunity to evaluate model performance, focusing on time series of near-surface weather variables, radiation and surface flux budgets, vertical structure and the nighttime inertial oscillation. The model results revealed substantial differences between the PBL schemes. Generally, non-local schemes tend to produce higher temperatures and higher wind speeds than local schemes, in particular, for nighttime. The WRF model underestimates the 2-m temperature during daytime (about \(2\) K) and substantially underestimates it at night (about \(4\) K), in contrast to the previous studies where modelled 2-m temperature was overestimated. Considering the 10-m wind speed, during the night turbulent kinetic energy based schemes tend to produce lower wind speeds than other schemes. In all simulations the sensible and latent heat fluxes were well reproduced. For the net radiation and the soil heat flux we found good agreement with daytime observations but underestimations at night. Concerning the vertical profiles, the selected non-local PBL schemes underestimate the PBL depth and the low-level jet altitude at night by about 50 m, although with the correct wind speed. The latter contradicts most previous studies and can be attributed to the revised stability function in the Yonsei University PBL scheme. The local, turbulent kinetic energy based PBL schemes estimated the low-level jet altitude and strength more accurately. Compared to the observations, all model simulations show a similar structure for the potential temperature, with a consistent cold bias ( \(\approx \) 2 K) in the upper PBL. In addition to the sensitivity to the PBL schemes, we studied the sensitivity to technical features such as horizontal resolution and domain size. We found a substantial difference in the model performance for a range of 12, 18 and 24 h spin-up times, longer spin-up time decreased the modelled wind speed bias, but it strengthened the negative temperature bias. The sensitivity of the model to the vertical resolution of the input and boundary conditions on the model performance is confirmed, and its influence appeared most significant for the non-local PBL parametrizations.     

Evaluation of the Reduction in Uncertainty Obtained by Conditioning a 3D Stochastic Channel to Multiwell Pressure Data

A stochastic channel embedded in a background facies is conditioned to data observed at wells. The background facies is a fixed rectangular box. The model parameters consist of geometric parameters that describe the shape, size, and location of the channel, and permeability and porosity in the channel and nonchannel facies. We extend methodology previously developed to condition a stochastic channel to well-test pressure data, and well observations of the channel thickness and the depth of the top of the channel. The main objective of this work is to characterize the reduction in uncertainty in channel model parameters and predicted reservoir performance that can be achieved by conditioning to well-test pressure data at one or more wells. Multiple conditional realizations of the geometric parameters and rock properties are generated to evaluate the uncertainty in model parameters. The ensemble of predictions of reservoir performance generated from the suite of realizations provides a Monte Carlo estimate of the uncertainty in future performance predictions. In addition, we provide some insight on how prior variances, data measurement errors, and sensitivity coefficients interact to determine the reduction in model parameters obtained by conditioning to pressure data and examine the value of active and observation well data in resolving model parameters.     

Controls of δS and δO in dissolved sulphate: Learning from a detailed survey in the Llobregat River (Spain)

The S and O isotopic composition of dissolved SO₄, used as a tracer for SO₄ sources, was applied to the water of the Llobregat River system (NE Spain). The survey was carried out at 30 sites where surface water was sampled on a monthly basis over a period of 2a. The concentration of dissolved SO₄ varied from 20 to 1575 mg L⁻¹. Sulphur isotopic compositions clustered in two populations: one – 93% of the samples – had positive values with a mode of +9‰; the other had negative values and a mode of −5‰. Data for δ¹⁸O_SO4 showed a mean value of +11‰, with no bi-modal distribution, though lower values of δ¹⁸O corresponded to samples with negative δ³⁴S. These values can not be explained solely by the contribution of bedrock SO₄ sources: that is, sulphide oxidation and the weathering of outcrops of sulphates, though numerous chemical sediments exist in the basin. Even in a river with a high concentration of natural sources of dissolved SO₄, such as the Llobregat River, the δ³⁴S values suggest that dissolved SO₄ is controlled by a complex mix of both natural and anthropogenic sources. The main anthropogenic sources in this basin are fertilizers, sewage, potash mine effluent and power plant emissions. Detailed river water sampling, together with the chemical and isotopic characterisation of the main anthropogenic inputs, allowed determination of the influence of redox processes, as well as identification of the contribution of natural and anthropogenic SO₄ sources and detection of spatial variations and seasonal changes among these sources. For instance, in the Llobregat River the input of fertilisers is well marked seasonally. Minimum values of δ³⁴S are reported during fertilization periods – from January to March – indicating a higher contribution of this source. The dual isotope approach, δ³⁴S and δ¹⁸O, is useful to better constrain the sources of SO₄. Moreover, in small-scale studies, where the inputs are well known and limited, the mixing models can be enhanced and the contribution of the different sources can be quantified to some extent.     

Cratering and modification of wet‐target craters: Projectile impact experiments and field observations of the Lockne marine‐target crater (Sweden)

Abstract— Marine impacts are one category of crater formation in volatile targets. At target water depths exceeding the diameter of the impactor, the zones of vaporization, melting, and excavation of the standard land‐target cratering model develop partially or entirely in the water column. The part of the crater that has a potential of being preserved (seafloor crater) may to a great extent be formed by material emplacement and excavation processes that are very different from land‐target craters. These processes include a high‐energy, water‐jet‐driven excavation flow. At greater water depths, the difference in strength of the target layers causes a concentric crater to evolve. The crater consists of a wide water cavity with a shallow excavation flow along the seabed surrounding a nested, deeper crater in the basement. The modification of the crater is likewise influenced by the water through its forceful resurge to fill the cavity in the water mass and the seafloor. The resurge flow is strongly erosive and incorporates both ejecta and rip‐up material from the seabed surrounding the excavated crater. A combination of field observations and impact experiments has helped us analyze the processes affecting the zone between the basement crater and the maximum extent of the water cavity. The resurge erosion is facilitated by fragmentation of the upper parts of the solid target caused by a) spallation and b) vibrations from the shallow excavation flow and, subsequently, c) the vertical collapse of the water cavity rim wall. In addition, poorly consolidated and saturated sediments may collapse extensively, possibly aided by a violent expansion of the pore water volume when it turns into a spray during passage of the rarefaction wave. This process may also occur at impacts into water‐saturated targets without an upper layer of seawater present. Our results have implications for impacts on both Earth and Mars, and possibly anywhere in the solar system where volatiles exist/have existed in the upper part of the target.     

Subsidence and thermal history of an inverted Late Jurassic‐Early Cretaceous extensional basin (Cameros,North‐central Spain) affected by very low‐ to low‐grade metamorphism

The Cameros Basin (North Spain) is a Late Jurassic‐Early Cretaceous extensional basin, which was inverted during the Cenozoic. It underwent a remarkable thermal evolution, as indicated by the record of anomalous high temperatures in its deposits. In this study, the subsidence and thermal history of the basin is reconstructed, using subsidence analysis and 2D thermal modelling. Tectonic subsidence curves provide evidence of the occurrence of two rapid subsidence phases during the syn‐extensional stage. In the first phase (Tithonian‐Early Berriasian), the largest accommodation space was formed in the central sector of the basin, whereas in the second (Early Barremian‐Early Albian), it was formed in the northern sector. These rapid subsidence phases could correspond to relevant tectonic events affecting the Iberian Plate at that time. By distinguishing between the initial and thermal subsidence and defining their relative magnitudes, Royden's (1986) method was used to estimate the heat flow at the end of the extensional stage. A maximum heat flow of 60–65 mW/m² is estimated, implying only a minor thermal disturbance associated with extension. In contrast with these data, very high vitrinite reflectance, anomalously distributed in some case with respect to the typical depth‐vitrinite reflectance relation, was measured in the central‐northern sector of the basin. Burial and thermal data are used to construct a 2D thermal basin model, to elucidate the role of the processes involved in sediment heating. Calibration of the thermal model with the vitrinite reflectance (%Ro) and fluid inclusion (FI) data indicates that in the central and northern sectors of the basin, an extra heat source, other than a typical rift, is required to explain the observed thermal anomalies. The distribution of the %Ro and FI values in these sectors suggests that the high temperatures and their distribution are related to the circulation of hot fluids. Hot fluids were attributed to the hydrothermal metamorphic events affecting the area during the early post‐extensional and inversion stages of the basin.     

Carbon isotopic fractionation during a mesocosm bloom experiment dominated by Emiliania huxleyi: Effects of CO2 concentration and primary production

We investigated the effect of CO₂ and primary production on the carbon isotopic fractionation of alkenones and particulate organic matter (POC) during a natural phytoplankton bloom dominated by the coccolithophore Emiliania huxleyi. In nine semi-closed mesocosms (∼11 m³ each), three different CO₂ partial pressures (pCO₂) in triplicate represented glacial (∼180 ppmv CO₂), present (∼380 ppmv CO₂), and year 2100 (∼710 ppmv CO₂) CO₂ conditions. The largest shift in alkenone isotopic composition (4-5‰) occurred during the exponential growth phase, regardless of the CO₂ concentration in the respective treatment. Despite the difference of ∼500 ppmv, the influence of pCO₂ on isotopic fractionation was marginal (1-2‰). During the stationary phase, E. huxleyi continued to produce alkenones, accumulating cellular concentrations almost four times higher than those of exponentially dividing cells. Our isotope data indicate that, while alkenone production was maintained, the interaction of carbon source and cellular uptake dynamics by E. huxleyi reached a steady state. During stationary phase, we further observed a remarkable increase in the difference between δ¹³C of bulk organic matter and of alkenones spanning 7-12‰. We suggest that this phenomenon is caused mainly by a combination of extracellular release of ¹³C-enriched polysaccharides and subsequent particle aggregation induced by the production of transparent exopolymer particles (TEP).     

Abundance, diversity and succession of aquatic Coleoptera and Heteroptera in a cluster of artificial ponds in the North German Lowlands

Our study describes and evaluates environmental influences on assemblages of aquatic Coleoptera and Heteroptera in artificial ponds situated near Lake Steinhude in Lower Saxony (Germany). We determined temporal dynamics and colonization patterns for 14 ponds of different age. In total, we recorded 4941 individuals that represented 87 species of aquatic beetles and bugs. Between 30 and 40 species were found in most of the ponds. Heteropteran species of the families Corixidae and Notonectidae acted as pioneer species in new ponds, while aquatic coleopterans predominated in older ponds. The results of Canonical Correspondence Analyses (CCA) showed that among the key factors affecting community structure were land use, vegetation cover, water chemistry and the age of the ponds. We found that the distribution of adjacent ponds on areas with different land use has a positive influence on the diversity and abundance of the aquatic insect fauna.