ALBERT

Beschreibung: Seismic facies analysis can generate a map to describe the spatial distribution characteristics of reservoirs, and therefore plays a critical role in seismic interpretation. To analyse the characteristics of the horizon of interest, it is usually necessary to extract seismic waveforms along the target horizon using a selected time window. The inaccuracy of horizon interpretation often produces some inconsistent phases and leads to inaccurate classification. Therefore, the developed adaptive phase K-means algorithm proposed a sliding time window to extract seismic waveforms. However, setting the maximum offset of the sliding window is difficult in a real data application. A value that is too large may cause the cross-layer problem, whereas a value that is too small reduces the flexibility of the algorithm. To address this disadvantage, this paper proposes a robust K-means (R-K-means) algorithm with a Gaussian-weighted sliding window for seismic waveform classification. The used weights punish those windows distant from the interpretation horizon in the objective function, consequently producing a smaller range of horizon adjustments even when using relatively large maximum offsets and benefitting the generation of stable and reliable seismic facies maps. The application of real seismic data from the F3 block proves the effectiveness of the proposed algorithm.

Beschreibung: This paper successfully applied the geoengineering workflow for integrated well-test analysis to characterise fluid flow in a newly discovered fractured reservoir in the Barents Sea. A reservoir model containing fractures and matrix was built and calibrated using this workflow to match complex pressure transients measured in the field. We outline different geological scenarios that could potentially reproduce the pressure response observed in the field, highlighting the challenge of non-uniqueness when analysing well-test data. However, integrating other field data into the analysis allowed us to narrow the range of uncertainty, enabling the most plausible geological scenario to be taken forward for more detailed reservoir characterisation and history matching. The results provide new insights into the reservoir geology and the key flow processes that generate the pressure response observed in the field. This paper demonstrates that the geoengineering workflow used here can be applied to better characterise naturally fractured reservoirs. We also provide reference solutions for interpreting well-tests in fractured reservoirs where troughs in the pressure derivative are recognisable in the data.

Beschreibung: As the robustness for the wave equation-based inversion methods, wave equation migration velocity analysis (WEMVA) is stable for overcoming the multipathing problem and has become popular in recent years. As a rapidly developed method, differential semblance optimisation (DSO) is convenient to implement and can automatically detect the moveout existing in common image gathers (CIGs). However, by implementing in the image domain with the target of minimising moveouts and improving coherence of the CIGs, the DSO method often suffers from imaging artefacts caused by uneven illumination and irregular observation geometry, which may produce poor velocity updates with artefact contamination. To deal with this issue, in this paper, by introducing Wiener-like filters, we modify the conventional image matching-based objective function to a new one by introducing the quadratic Wasserstein metric technique. The new misfit function measures the distance of two distributions obtained by the convolutional filters and target functions. With the new misfit function, the adjoint sources and the corresponding gradients are improved. We apply the new method to two numerical examples and one field dataset. The corresponding results indicate that the new method is robust to compensate low frequency components of velocity models.

Beschreibung: Atlantic salmon (Salmo salar) populations throughout the North Atlantic have declined in recent decades largely due to reduced marine survival, yet our understanding of marine distribution patterns and migratory routes remains limited. Here, we assigned archived individual samples (n = 3891) collected over a half century (1968–2018) throughout the North Atlantic to region of origin using range-wide genetic assignment. In the Northwest Atlantic, the distribution of assignments reinforced the importance of the Labrador Sea as an aggregation area, with 73% of all reporting groups detected. Moreover, individuals from six European reporting groups were identified in the Northwest Atlantic, and detections decreased with decreasing latitude spanning an area from Greenland to southern Newfoundland. In the Northeast Atlantic, six North American reporting groups were detected in samples from around the Faroe Islands. Based on the distribution of samples, estimates of trans-Atlantic migration distance averaged 3861 and 2889 km for North American and European salmon respectively. Our analysis highlights the widespread importance of the Labrador Sea and Faroe Islands to the species marine distribution patterns, and the prevalence of long-distance trans-Atlantic migration. Ultimately, the results suggest that environmental conditions experienced by many Atlantic salmon populations span much of the North Atlantic Ocean.

Beschreibung: Heterotrophic protists are essential components of the marine ecosystem, yet they are often excluded from monitoring programmes. With limited resources, monitoring strategies need to be optimised considering both scientific knowledge and available resources. In doing so, it is crucial to understand how sampling frequency affects the value of the data. We analysed 11 years of weekly heterotrophic protist time-series data from Station L4 in the Western English Channel to explore how different sampling intervals impact data quality. In the L4 dataset, comprising 55 protist taxa, the reduction of sampling frequency from weekly to four times a year at specific seasons decreased the number of taxa encountered by 38% for ciliates and 29% for heterotrophic dinoflagellates while the mean annual biomass or its mean variation were not affected. Furthermore, when samples were taken only four times a year, biomass peaks of the ten most important taxa were often missed. The primary motivator for this study was furthering the development of the heterotrophic protist monitoring in temperate and subarctic marine areas, e.g. the Baltic Sea. Based on our findings, we give recommendations on sampling frequency to optimise the value of heterotrophic protist monitoring.

Beschreibung: The Lafaye orbicular body was emplaced in the Villatange tonalite-granodiorite unit of the Guéret magmatic complex (Massif Central, France). It consists of plagioclasic orbicules (4–35 cm diameter) embedded in homogeneous cordierite granodiorite. Orbicule cores consist mostly of residual metasedimentary xenoliths or autolithic plagioclasic cumulates. Rims (0.7–8 cm thickness) are single- or multi-layered; layers, mostly comb-textured, comprise alternating sheets dominated by cordierite (XFe = 0.32–0.37) or plagioclase (mostly An25–30). Additional mineral phases are minor biotite (XFe = 0.52; AlVI = 0.58–0.92 atoms per formula unit) and interstitial quartz. Plagioclase and cordierite morphologies (needle-like, skeletal, branching or fan-shaped) indicate growth under high initial supersaturation. However, the final polyhedral shapes and primary zoning of many individual plagioclase crystals, as well as evidence of partial recrystallization, imply significant textural maturation. Whole-rock major and trace element data (A/CNK = 1.12–1.46) and Sr and Nd isotopic compositions (εNd(355 Ma)  = −8.6 to −7.4; 87Sr/86Sr(355 Ma) = 0.7110–0.7147) suggest that the parental magma of the orbicules resulted from bulk assimilation of aluminous metasediments by a Villatange-type granodioritic magma. Heterogeneous nucleation and growth of plagioclase and cordierite around xenoliths/autoliths are interpreted in terms of (1) adiabatic decompression of magma pulses ascending in dykes leading to superheating and resorption of early solids, and (2) volatile exsolution, inducing undercooling, supersaturation, and rim crystallization. The variability of layers (number, thickness, mineral distribution, and texture) is considered to result from oscillatory crystallization combined with variable plagioclase growth rates linked to changes in the degree of supersaturation as a function of the extent of melt degassing, itself linked to magma transfer dynamics.

Beschreibung: The high fluxes of iron minerals associated with aeolian dry deposition may result in an anomalously high reactive iron content and the rapid reoxidation of hydrogen sulfide in sediments. This will prevent the formation of pyrite and result in a ‘cryptic’ sulfur cycle. We studied the cycling of iron and sulfur in deep water (〉800 m water depth) sediments of the Red Sea and its northern extension, the Gulf of Aqaba. We found that the reactive iron content in the surface sediments of the Gulf of Aqaba and the Red Sea is high, whereas the amount of sulfur-bound iron is very low and decreases with the water depth. The presence of trace amounts of pyrite and zero-valent sulfur, as well as the isotopic compositions of sulfate and pyrite, which are consistent with sulfate reduction under substrate-limiting conditions, suggest that cryptic sulfur cycling is likely to be a result of the rapid reoxidation of hydrogen sulfide rather than the suppression of microbial sulfate reduction. The low amount of reactive iron and high organic carbon content in the sediments of the Shaban Deep, which are overlain by hypersaline hydrothermal brines, result in a non-cryptic sulfur cycle characterized by the preservation of pyrite in the sediments.Supplementary material: Data for sulfur speciation, iron speciation and sulfur isotope composition are available at https://doi.org/10.6084/m9.figshare.c.5508155Thematic collection: This article is part of the Sulfur in the Earth system collection available at: https://www.lyellcollection.org/cc/sulfur-in-the-earth-system

Beschreibung: La Réunion Island includes two major volcanic systems. About 0·5 Myr ago, Piton des Neiges volcano declined, while Piton de la Fournaise volcano grew on its flank. Since then the Piton de la Fournaise shield volcano has produced homogeneous lavas with chemical compositions transitional between alkali and tholeiitic basalts. In April 2007, the volcano emitted a very small volume of trachytic pumice during its largest historical eruption. We conducted a comprehensive petrological and geochemical study of the pumice to understand the occurrence of such silicic melt in the feeding system of this highly active basaltic volcano. Isotopes of Sr, Nd, Pb and O, together with trace elements, indicate that the trachyte is genetically related to the La Réunion mantle plume and derives from crystallization of a typical basalt. The trachyte chemistry records a long and complex history of differentiation and outgassing. The extensive depletion of moderately volatile elements (F, Cl, B, Cs, Cu, Li) and less volatile uranium is consistent with exsolution of dense fluids at depths of several kilometres. Lithium isotopes point to closed-system degassing during the very late stages of crystallization. U-series isotopes and radiogenic 208Pb*/206Pb* constrain the age of U loss to between 0·4 and 2·1 Ma. This age is as old as or older than the Piton de la Fournaise shield edifice. The 2007 trachyte could thus be a liquid remnant of an extinct volcano, such as Piton des Neiges or Les Alizés (Piton de la Fournaise proto-volcano). It could also result from partial melting of an old syenite intrusion or remobilization of interstitial melts not fully solidified. Thermal modelling indicates that the sustained heat flux from hot basaltic magmas rising from the mantle can maintain temperatures above 800 °C in the central feeding system, and prevent total solidification of magmas trapped in this hot core.

Beschreibung: Fisheries management is usually supported by technical and financial measurements (i.e. logbooks and market data), which are helpful for ecological or economic assessments. Yet this information is not able to address social heterogeneity and fisher motivations, which are key to understanding fisher behaviour. This case study of the demersal segment in the Netherlands shows that combining quantitative analysis of logbooks with qualitative data collected by engaging with fishers can capture both fishing activity and its motivations, generating a more social understanding of fisher behaviour. A métier analysis of logbook data describes five dominant fishing practices among the selected segment. Twenty-five in-depth interviews with fishers along with focus groups including other experts identify three social factors that influence fisher behaviour in the Dutch demersal fleet: business structure, working rhythm, and polyvalence. The results show that motivations for fisher behaviour are more complex than complying with regulations or seeking profit: social factors also influence fishing activity. Furthermore, these social factors have real implications for the impacts of management measures on both the fishing communities and the environment, especially in times of change. These results are useful for management strategy development or evaluation because they are feasibly observable through existing data collection protocols.

Beschreibung: Harbour seals are surveyed aerially when they haul-out to moult in August. The proportion of the population hauled out throughout the year is related to temporal, environmental, and meteorological variables. Thus, monitoring is conducted under predefined ranges of conditions. Effects of variation within these ranges are rarely reviewed. We used linear models to assess effects of time, date and weather on the difference between counts predicted by a population growth model and observed counts, based on a 30-year time-series. Our top-ranked model explained 34.4% of the variance. Survey date and its interaction with survey year were the most important variables, with higher counts earlier in August, particularly early in the time series, where surveys may not have been timed optimally to capture the peak in the moult. Cloud cover, wind speed, temperature, and interactions between these were of lesser importance; there were fewer seals on land during cloudy, windy days and on clear, warm days. These effects of weather are likely related to temperature regulation. Power analyses suggested that correction for survey conditions would allow detection of a one percentage point annual change in population growth rate with 80% power 4 years sooner than without taking survey conditions into account.

Beschreibung: Sargassum forests play an important role in coastal waters as habitats for marine organisms, including commercial species. However, human activities have negatively affected their distribution causing a worldwide decline of Sargassum forests. Mapping and monitoring the distribution and biomass of these habitats using acoustic remote sensing techniques is key for their conservation. Nonetheless, most researches based on acoustic remote sensing methods focus on estimations of macrophyte area and its canopy height, and less researches reporting 3D visualization of these habitats. This study demonstrates the use of high-resolution multibeam echosounder (MBES) bathymetric data to visualize the 3D structure of Sargassum forests. Comparing acoustic data and underwater camera photos collected in field surveys, we identified Sargassum individuals as vertical clusters of contiguous sounding points with a base close to the sea bottom in the sounding data of the MBES. Using this criterion, we could distinguish Sargassum echoes, visualize the 3D structure of Sargassum forests and estimate the number of Sargassum individuals in the survey area. Using the relation between thallus length and dry weight of sampled Sargassum plants, standing stock and biomass could be estimated assuming the thallus length was the height of Sargassum plants identified with the MBES.

Beschreibung: This paper reviews application of the Daily Egg Production Method (DEPM) to sardine (Sardinops sagax) off southern Australia between 1995 and 2019. Coefficients of variation (CVs) of estimates of spawning biomass (SB) were reduced from 23–59% to 8–12% by: (i) estimating mean daily egg production (P0), spawning fraction (S), and sex ratio (R) from all historical data rather than annually; and (ii) combining batch fecundity (F) and female weight (W) into a single parameter, relative fecundity (F′ = F^/W). Total daily egg production was estimated most precisely from annual estimates of spawning area (A) and estimates of P0 obtained from historical data. Both S and R were estimated most precisely from historical data. Estimating W and F from historical data did not increase precision. F′ had lower CVs than both W and F, and was stable across years and a wide range of W. Findings demonstrate that A can be converted into a precise estimates of SB using estimates of P0, S, R, and F′ obtained from historical data. However, the possibility that DEPM parameters may change in the future cannot be discounted. Future monitoring should include annual estimation of P0 and periodic (e.g. 3–5 years) re-estimation of adult parameters.

Beschreibung: In his 1953 paper, Harald Sverdrup argued that the development of a spring bloom in the ocean depends on the juxtaposition of two depth horizons: the mixed-layer depth and the critical depth. Mixed-layer depth shallower than the critical depth favours phytoplankton growth in the layer and vice versa. However, mathematically, Sverdrup left the problem unsolved in the form of a transcendental equation. In spite of the high number of citations that this paper has garnered, the solution to this equation has not been found, until now. In this work, we present an analytical solution for the critical depth, as originally defined by Sverdrup. The paper opens with the definition of the critical depth and the description of the Lambert W function. The analytical solution for critical depth follows. Sverdrup’s original model is extended to include the effect of light attenuation by phytoplankton and the analytical solution for steady-state biomass in the mixed layer is derived. The expression for mixed-layer production at steady state is also presented. Two novel variants of the critical depth are defined: the optically uncoupled critical depth and the optically coupled critical depth. It is demonstrated that at steady state the optically coupled critical depth equals the mixed-layer depth and that the irradiance at the base of the mixed layer equals the irradiance at the optically uncoupled critical depth. Competitive exclusion is demonstrated to hold and the optically uncoupled critical depth is linked to the critical light intensity in multi-species competition. Finally, a conservation principle for the critical depth is found.

Beschreibung: Fishery bycatch poses a serious threat to seabird populations globally. Traditional haul-only post-capture observations are inadequate and inefficient to document seabird bycatch due to the substantial bycatch loss known to occur. Pre-capture observations offer an alternative by documenting seabird interactions leading up to captures. Based on the long-term large-scale dedicated field observations, this study revealed significant risk factors for the pre-capture stages of the seabird bycatch process in pelagic longline fisheries using Bayesian methods. Rough sea conditions were found to correlate with more seabirds following fishing vessels. Species identity, density effect, inter-species effect, and sea condition were found to significantly affect how frequently seabirds aggregated around a fishing vessel engage in bait-taking interactions. Intra-species competition was found to be the dominant type of density effect. Moreover, a web of inter-species interactions was identified to facilitate the bait-taking of superior competitors at the expense of inferior ones. The findings of this study are relevant to fishery managers in updating current data collection protocols to alleviate data issues caused by bycatch loss, to conservation biologists in quantifying bycatch risks for susceptible seabird populations, and in aiding the design and evaluation of bycatch mitigation measures.

Beschreibung: Fishing and translocation of marine species for use in aquaculture is widespread. Corkwing, goldsinny, and ballan wrasse (Symphodus melops, Ctenolabrus rupestris, and Labrus bergylta) are fished on the Swedish west coast for use as cleaner-fish in Norwegian salmon farms. Here, we aim to provide knowledge and recommendations to support ecosystem-based management for wrasse fisheries in Sweden. We compared fished and non-fished areas to test if current fishery levels have led to stock depletion. To gain insight on the role of wrasse in the algal belt trophic chain, we analysed the gut contents of goldsinny and corkwing using metabarcoding. Finally, we analysed the trophic interactions of wrasse and potential prey in a mesocosm study. We could not detect any signs of stock depletion or altered size structure in fished areas compared to the protected control area. Gut analyses confirmed both goldsinny and corkwing as non-specialized, omnivorous opportunists and revealed, with 189 prey taxa detected, a broader spectrum of prey than previously known. Common prey items included mesoherbivores such as small gastropods and crustaceans, but also insects and algae. We conclude that there are no visible signs of stock depletion at the current removal level of wrasses by the fishery. However, this emerging fishery should be closely monitored for potential cascading effects on the algal belt ecosystem, and our study could provide a baseline for future monitoring.

Beschreibung: Target strength model inputs including morphometry, material properties, lipid composition, and in situ orientations were measured for sub-Arctic krill (Euphausia pacifica, Thysanoessa spinifera, T. inermis, and T. raschii) in the eastern Bering Sea (EBS, 2016) and Gulf of Alaska (GOA, 2017). Inter-species and -regional animal lengths were significantly different (F1,680 = 114.10, p 

Beschreibung: Atlantic cod (Gadus morhua) is one of the most commercially important fish species in the North Atlantic. Environmental factors, such as water temperatures, influence growth of individuals over time, thus forming population-specific growth patterns across climatic regions. Here we develop an integrative approach to investigate the role of temperature in shaping geographic differences of cod growth in the Celtic Sea, North Sea, Iceland, and Barents Sea. We combine a physiology-based growth model and 50-years observational temperature data of 0.5 × 0.5° spatial resolution to simulate continuous growth of cod. The model generated weight-at-age data for the period 1959–2007 which we compared to observational data from fishery-independent scientific surveys. In the Celtic and the northern North Sea, simulated growth matches well observational data. We also show that relatively warm temperatures in the Celtic Sea facilitate maximum growth rates; future warming is likely to have a negative impact on growth of these cod stocks. Growth simulations in Icelandic waters and the Barents Sea are less consistent with local observational data. More complex growth patterns in these regions are probably shaped by ontogenetic shifts in temperature regimes, feeding conditions and physiological adaptations. These findings should stimulate further research on critical processes to be considered in population-specific projections of growth of cod and productivity.

Beschreibung: In designing and performing surveys of animal abundance, monitoring programs often struggle to determine the sampling intensity and design required to achieve their objectives, and this problem greatly increases in complexity for multispecies surveys with inherent trade-offs among species. To address these issues, we conducted a multispecies stratified random survey design optimization using a spatiotemporal operating model and a genetic algorithm that optimizes both the stratification (defined by depth and longitude) and the minimum optimal allocation of samples across strata subject to prespecified precision limits. Surveys were then simulated under those optimized designs and performance was evaluated by calculating the precision and accuracy of a resulting design-based abundance index. We applied this framework to a multispecies fishery-independent bottom trawl survey in the Gulf of Alaska, USA. Incorporating only spatial variation in the optimization failed to produce population estimates within the prespecified precision constraints, whereas including additional spatiotemporal variation ensured that estimates were both unbiased and within prespecified precision constraints. In general, results were not sensitive to the number of strata in the optimized solutions. This optimization approach provides an objective quantitative framework for designing new, or improving existing, survey designs for many different ecosystems.

Beschreibung: Changes in animal movement and behaviour at fine scales (tens of metres) in immediate proximity to tidal stream turbine structures are largely unknown and have implications for risks of animal collision with turbine blades. This study used upward-facing multibeam echosounder data to detect and track animal movement comprising fish, diving seabirds, and marine mammals. Measurements over spring-neap tidal cycles at a turbine structure (no blades present) are compared to a neighbouring reference area with no structure and comparable conditions, with measurements consecutive in time to maximize comparability. The majority of tracked animals (93.4% around turbine structure and 99.1% without turbine structure) were observed swimming against the flow, with 87.5% and 97.8%, respectively, making ground and showing capability of manoeuvring in tidal stream flow speeds. Track tortuosity increased around the turbine structure compared to the reference site, particularly in the wake and at low flow speeds, indicating animal station-holding or milling behaviour. These data also evidence the benefits of multibeam echosounders to measure animal movement through larger measurement volumes rather than relying on single-beam echosounders to measure animal presence alone, including to avoid large biases overestimating the size of schools swimming against the flow measured by time-in-beam.

Beschreibung: The Antarctic Peninsula is one of the most rapidly warming regions on earth, and it is likely that the abundance and distribution of marine predators will change as a result.Procellariiform seabirds are highly mobile predators, which target specific habitat characteristics associated with underlying distributions of prey and areas of increased prey availability. We use ship surveys and hurdle models, to estimate the summer distribution and relative density of 11 seabird species within the northern Antarctic Peninsula marine ecosystem. Models differed among species; however, sea surface temperature and depth were frequently associated with seabird occurrence and had the greatest explanatory power across many species. Null models based on observation data were better at predicting seabird density than models that included environmental covariates. This suggests that the main driver of distribution patterns is the broad-scale habitat features, and fine-scale aggregations within these ranges are harder to predict. Our seabird distribution models reflect known habitat associations, species hotspots, and community organization relative to oceanic and coastal marine processes. Application of species distribution models will benefit the assessments of critical habitat and potential responses to climate change and anthropogenic disturbance, which will provide insight into how species may change in polar ecosystems.

Beschreibung: Many data-limited fish stocks worldwide require management advice. Simple empirical management procedures have been used to manage data-limited fisheries but do not necessarily ensure compliance with maximum sustainable yield objectives and precautionary principles. Genetic algorithms are efficient optimization procedures for which the objectives are formalized as a fitness function. This optimization can be included when testing management procedures in a management strategy evaluation. This study explored the application of a genetic algorithm to an empirical catch rule and found that this approach could substantially improve the performance of the catch rule. The optimized parameterization and the magnitude of the improvement were dependent on the specific stock, stock status, and definition of the fitness function. The genetic algorithm proved to be an efficient and automated method for tuning the catch rule and removed the need for manual intervention during the optimization process. Therefore, we conclude that the approach could also be applied to other management procedures, case-specific tuning, and even data-rich stocks. Finally, we recommend the phasing out of the current generic ICES “2 over 3” advice rule in favour of case-specific catch rules of the form tested here, although we caution that neither works well for fast-growing stocks.

Beschreibung: Ocean grabbing occurs when traditional users, such as small-scale fishers, are pushed aside by new development activities. This grabbing must be prevented to avoid sea uses that maintain or increase social inequity. In this paper, we show that in tropical Atlantic countries, such as Brazil and Senegal, examples of ocean grabbing already occur. In this context, we analyse if Maritime Spatial Planning (MSP) may be an opportunity to limit ocean grabbing or, to the contrary, poses a risk to increase it. MSP calls for an ecosystem approach that requires integrated coastal and marine management and involves stakeholders in developing a shared vision of the future, where society and environment are preserved. However, recent studies have shown that MSP is a process to be used cautiously to ensure equitable decisions. Meanwhile, the concept is spreading worldwide including in tropical Atlantic countries. We highlight that context matters and the specificities of the tropical Atlantic must be taken into account when deploying MSP processes. In the tropical Atlantic context, there is increased imbalances of stakeholder power, traps from decision support tools, and a need for adaptive management. These specific features must be addressed when deploying MSP in a way to avoid ocean grabbing.

Beschreibung: The influence of oceanographic variables on assemblages of meso- and bathypelagic fish was investigated along a Northeastern Atlantic Ocean transect (Cape Verde to the Bay of Biscay) during May 2019. Fish were collected using a mrozooplankton trawl during daylight hours at ten stations. Along the transect, 17 hydrographic stations were also performed with a CTD (Conductivity, Temperature, and Depth). A total of 130 fish taxa were identified. The dominant family was Gonostomatidae, with four species (Cyclothone braueri, Cyclothone microdon, Cyclothone pseudopallida, and Cyclothone pallida) being responsible of more than 78% of the total density. The most frequent species that appeared to be ubiquitous were C. braueri and C. pseudopallida, while Myctophidae was the most diverse family. Multivariate analyses revealed two clusters related with the latitudinal gradient. The fish community in the southern stations (25–37°N) was more diverse than in the northern stations (42–48°N). Temperature from 300 to 700 m depth explained 65% of variation in terms of density and 58% in terms of biomass, both statistically significant. The investigated variation in the deep-pelagic ecosystems on a large spatial scale gives essential information to ecosystem management approhes and marine spatial planning.

Beschreibung: In the European Union (EU), the ecosystem approach to fisheries is implemented through several directives and polices, which are overarched by the Marine Strategy Framework Directive (MSFD). The MSFD requires the assessment of the environmental status of exploited fish and shellfish stocks, among others, to advise the EU Common Fisheries Policy (CFP) on sustainable catch options. The stock assessments for the CFP and the MSFD are supposed to be congruent, yet they differ in several substantial subjects. While the annual CFP assessment is based on two indicators, fishing mortality (F) and spawning stock biomass (SSB), in its Descriptor 3 (D3) the MSFD requires the assessment of three criteria (F, SSB, and age or size structure) within a six year period. Assessing exploited fish and shellfish stocks under the MSFD by using existing CFP assessments has therefore been a half-completed task, which had to be fulfilled by the member states of the EU. This paper suggests six easy steps, referred to as “quick wins” (QW), which are based on existing information from CFP stock assessments. The implementation of these six QW would allow for assessments of exploited fish and shellfish stocks that are compliant to D3. These QW are to (i) assess length/age structure within a stock, (ii) analyse and assess selectivity by fisheries, (iii) use all available information to assess a stock, (iv) use response indicators to assess environmental targets, (v) provide integrated stock-specific advice, and (vi) provide assessments with a mid-term perspective over a-period of six years. International Council for the Exploration of the Sea (ICES) combines the infrastructure and expertise to produce stock-specific assessment products tailored to the requirements of MSFD D3. Thereby, ICES could provide a valuable service to EU member states in the north-east Atlantic region by providing scientifically validated, quality-assured, and MSFD-compliant single-stock assessment products. This would be a big advancement towards implementing the ecosystem approach to fisheries management within Europe.

Beschreibung: Managers, stakeholders, and scientists recognize the need for collaborative, transparent, integrated approaches to complex resource management issues, and frameworks to address these complex issues are developing. Through the course of 2019, the Mid-Atlantic Fishery Management Council developed a conceptual model of ecosystem linkages and risks for summer flounder, a species of recreational and commercial fisheries importance. The proximal aim of the model was to develop a list of integrated management questions that could be refined and addressed through a future quantitative management strategy evaluation. As such, this conceptual model served as a scoping tool. However, the true value of the conceptual model lays elsewhere: familiarizing resource managers historically focused on single-species management with the potential utility of an ecosystem approach to management. This paper details the goals and development of the conceptual model and situates this process in the broader context of best practices for collaborative open science and scientific reproducibility. Further, it highlights a successful path by which the shift towards ecosystem-based management can be actuated.

Beschreibung: Understanding the competitive interactions of ecological similar species is essential to determine their role and niche in the ecosystem. Using both conventional dietary methods and stable isotope analysis, we examined the feeding ecology, trophic position and possible niche partitioning of Northeast Atlantic mackerel (Scomber scombrus), Icelandic summer spawning (ISS) herring (Clupea harengus), and Norwegian spring spawning (NSS) herring from Icelandic waters during the spring and summer 2012 and 2014. The stomach analysis showed differences in prey preferences among the species during summer, where mackerel diets were almost exclusively calanoid copepods, while herring ate larger zooplankton, i.e. euphausiids and amphipods. Analysis of isotopic diet contribution of mackerel and herring using Bayesian mixing models, representing mainly the spring consumption, revealed that euphausiids were the main dietary contributor of both mackerel and herring. Fish prey were also an important dietary source for ISS herring, and amphipods were important to the diet of NSS herring. In addition, mackerel encompassed a broad isotopic niche, whereas the two herring stocks had narrower and separated niches from each other. The results from this analysis reveal new and holistic information into the diet of these species, which have various ecological implications.

Beschreibung: Theoretical size at the onset of maturity (TSOM) for female Norway lobster was estimated by a new methodology based on the probability distributions of mature individuals built on physiological maturity measures. Onset of maturity using TSOM varied from 18.4- to 33.7-mm carapace length for the Irish functional management units (FUs). These estimates showed a significant negative linear relationship (R2 = 0.60) with population density and a significant positive linear relationship with average size in females (R2 = 0.84). The size class at which 50% of the females are sexually mature (L50) was linked to the new TSOM metric by a significant positive linear relationship (R2 = 0.40). This set of linear relationships ultimately allowed TSOM and L50 to be estimated without a requirement for maturity stages to be distinguished. As well as contributing to the stock assessment and management of Nephrops (e.g. in data-limited FUs) and its potential for the calibration of more routinely used estimates, TSOM might be applied in new species and meta-analyses where size of maturity data are scarce. This new metric also better-defines the maturity process since, taken together, TSOM, L50, and smallest berried female represent sequential maturity events: (i) onset of maturity, (ii) 50% mature (from gonad staging), and (iii) berried females.

Beschreibung: Although gelatinous zooplankton are an important component of marine ecosystems, gelatinous mesozooplankton that are

Beschreibung: The carapace length (CL) at which American lobster (Homarus americanus) females reach maturity can be used to evaluate egg production, growth patterns, and the overall health of lobster stocks. The female maturity datasets used to represent Gulf of Maine (GOM) lobsters in the 2015 Atlantic States Marine Fisheries Commission American Lobster Stock Assessment were collected in the 1990s by the Maine Department of Marine Resources at two coastal sites. Many studies have demonstrated an inverse relationship between temperature and the size at maturity in female lobsters, and GOM waters have warmed significantly over this period. To update these GOM maturity datasets, we used ovarian staging to determine the maturity status of over 1200 females from fives sites over 3 years. Broad application of this methodology in tandem with key growth measurements on females 50–120 mm CL allowed us to characterize reproductive development and generate maturity ogives (proportion mature at a given CL). We observed a latitudinal gradient in the size at maturity across this coastal region of the GOM and quantified a decrease in this size over 25 years. These findings have implications for future stock assessment approaches and management measures implemented to sustain this valuable fishery.

Beschreibung: The carbonate chemistry in coastal waters is more variable compared with that of open oceans, both in magnitude and time scale of its fluctuations. However, knowledge of the responses of coastal phytoplankton to dynamic changes in pH/pCO2 has been scarcely documented. Hence, we investigated the physiological performance of a coastal isolate of the coccolithophore Emiliania huxleyi (PML B92/11) under fluctuating and stable pCO2 regimes (steady ambient pCO2, 400 μatm; steady elevated pCO2, 1200 μatm; diurnally fluctuating elevated pCO2, 600–1800 μatm). Elevated pCO2 inhibited the calcification rate in both the steady and fluctuating regimes. However, higher specific growth rates and lower ratios of calcification to photosynthesis were detected in the cells grown under diurnally fluctuating elevated pCO2 conditions. The fluctuating pCO2 regime alleviated the negative effects of elevated pCO2 on effective photochemical quantum yield and relative photosynthetic electron transport rate compared with the steady elevated pCO2 treatment. Our results suggest that growth of E. huxleyi could benefit from diel fluctuations of pH/pCO2 under future-projected ocean acidification, but its calcification was reduced by the fluctuation and the increased concentration of CO2, reflecting a necessity to consider the influences of dynamic pH fluctuations on coastal carbon cycles associated with ocean global changes.

Beschreibung: Millions of people’s livelihoods rely on artisanal fisheries. However, in many regions fishers are increasingly facing ciguatera poisoning, a seafood-borne illness. The toxin, produced by benthic dinoflagellates, can spread through marine food webs and to humans by direct consumption. Ciguatera risk can play a major role in fisher’s activities but has never been considered in any marine spatial plans thus far. To fill this gap, we examined if integrating ciguatera in systematic conservation plans could affect these decisions. We developed through map-based interviews, a novel seven-step framework to collect and map local knowledge on ciguatera risk and fisheries activities with two innovations: (i) better mapping of fishing grounds by combining geomorphological habitat and fishing gear information, and (ii) integrating ciguatera risk directly into systematic spatial planning designs and scenarios conceived to maximize benthic habitat conservation while minimizing impacts to fishers. The approach is illustrated for Raivavae Island, in French Polynesia, Pacific Ocean. We found that integrating ciguatera significantly improved prioritization solutions with a 24–38% decrease of costs to fishers compared with scenarios based solely on fishery data. This framework was designed for scientists and managers to optimize the implementation of conservation plans and could be generalized to ciguatera-prone areas.

Beschreibung: Deep-sea marine fishes support important fisheries but estimates of their distributions are often incomplete as the data behind them may reflect fishing practices, access rights, or political boundaries, rather than actual geographic distributions. We use a simple suitable habitat model based on bottom depth, temperature, and salinity to estimate the potential distribution of Greenland halibut (Reinhardtius hippoglossoides). A large presence-only dataset is examined using multivariate kernel densities to define environmental envelopes, which we link to spatial distribution using a pan-Arctic oceanographic model. Occurrences generally fit the model well, although there were gaps in the predicted circum-Arctic distribution likely due to limited survey activity in many of the ice-covered seas around the Arctic Ocean. Bottom temperature and depth were major factors defining model fit to observations, but other factors, such as ecosystem interactions and larval drift could also influence distribution. Model predictions can be tested by increasing sampling effort in poorly explored regions and by studying the connectivity of putative populations. While abundances of Greenland halibut in the High Arctic are currently low, some areas are predicted to be suitable habitat for this species, suggesting that on-going sea-ice melt may lead to fisheries expansion into new areas.

Beschreibung: On the continental shelf, New Zealand bryozoans dominate ecologically-important three-dimensional benthic habitat providing structural complexity which hosts a wide variety of fauna, including economically valuable species (e.g. oysters and blue cod). The association between these species and bryozoan-dominated biogenic habitat commonly results in trawling damage to the benthos; eliminating pressure from destructive fishing practices could support bryozoan regrowth. In 2002, a voluntary fishing restriction was designated over part of one such bryozoan-dominated biogenic habitat, but the efficacy of this ban has not been assessed. Statutory marine protected areas have been proposed in the same area; they aim to achieve bryozoan protection and recovery. A comparison of the bryozoan fauna using benthic images from 2003 and 2019 from the Otago shelf, South Island, New Zealand/Te Waipounamu, Aotearoa shows that a 17-year voluntary fishing ban has been ineffective at increasing the cover and richness of bryozoan-dominated biogenic habitat on the Otago shelf. This study provides a baseline for future monitoring and develops useful metrics for recovery, while calling for further research on the biology and ecology of habitat-forming bryozoans.

Beschreibung: Ecosystem-based management is generally viewed as one of the most promising avenues for addressing the various anthropogenic pressures facing the world’s marine ecosystems. These approaches have been developed to varying degrees by individual countries or international organisations, but there remain a large proportion of marine ecosystems, particularly in developing regions, that have not yet been the subject of such research. In these areas, lack of effective regulation and the often high importance of the marine environment in providing food and economic opportunities, together create conditions where marine resources and habitats come under unsustainable levels of pressure. Here, we present a data-limited assessment approach to discern marine ecological patterns, in this case for the exclusive economic zone of Vietnam. By combining data from environmental and biological surveys from the Vietnamese national survey dataset and local oceanographic models, we have identified a series of 12 candidate ecological production units, delineated by their environmental characteristics, and the key commercial species that exist within them. These units are suggested as a possible foundation for a spatial management structure in the Vietnamese exclusive economic zone including considerations such as placement of marine protected areas, or ecological boundaries of key areas of socio-economic importance.

Beschreibung: This paper explores the possibility of using the ensemble modelling paradigm to fully capture assessment uncertainty and improve the robustness of advice provision. We identify and discuss advantages and challenges of ensemble modelling approaches in the context of scientific advice. There are uncertainties associated with every phase in the stock assessment process: data collection, assessment model choice, model assumptions, interpretation of risk, up to the implementation of management advice. Additionally, the dynamics of fish populations are complex, and our incomplete understanding of those dynamics and limited observations of important mechanisms, necessitate that models are simpler than nature. The aim is for the model to capture enough of the dynamics to accurately estimate trends and abundance, and provide the basis for robust advice about sustainable harvests. The status quo approach to assessment modelling has been to identify the “best” model and generate advice from that model, mostly ignoring advice from other model configurations regardless of how closely they performed relative to the chosen model. We discuss and make suggestions about the utility of ensemble models, including revisions to the formal process of providing advice to management bodies, and recommend further research to evaluate potential gains in modelling and advice performance.

Beschreibung: Climate effects on marine fish depend on life stage, particularly when life stages differ in habitat utilization. In the present study, we investigated life stage-dependent responses of lesser sandeel (Ammodytes marinus) to temperature at contrasting geographical scales. We related population density and individual growth to temperature and found different temperature responses between the first and the second years of life. During the first year of life, fish size was the single most important factor influencing sandeel abundances, indicating a positive relationship between growth and survival. In contrast, during the second year of life, autumn bottom temperature was negatively correlated with sandeel abundance, suggesting elevated mortality in warm years. Southerly areas, experiencing higher temperatures in general, were also the areas showing the strongest response to temperature. This study sheds light on how warming impacts population dynamics of one of the most important forage fishes in the North Sea and supports the discussion of underlying mechanisms.

Beschreibung: Recent ichnological analysis conducted in two sections (Rodiles and Lastres) of the Asturian Basin revealed the presence of Halimedides Lorenz von Liburnau 1902, which occurs just above the black shales related to the end of the T-OAE. Halimedides is associated with the recovery of the tracemaker community after the reestablishment of favourable, oxic, conditions. The appearance of Halimedides after the T-OAE event, previously not registered, supports the close relationship of the tracemaker with oxygen conditions, as occurs in other anoxic events including the Cretaceous OAE-1a and OAE-2. Also a relation between morphometric and paleoenvironmental parameters is observed, occurring larger and densely chambered specimens in darker, weakly oxygenated facies, while smaller and sparsely chambered forms are registered in lighter, better oxygenated sediments.

Beschreibung: The geometry, distribution, and rock properties (i.e. porosity and permeability) of turbidite reservoirs, and the processes associated with turbidity current deposition, are relatively well known. However, less attention has been given to the equivalent properties resulting from laminar sediment gravity-flow deposition, with most research limited to cogenetic turbidite-debrites (i.e. transitional flow deposits) or subsurface studies that focus predominantly on seismic-scale mass-transport deposits (MTDs). Thus, we have a limited understanding of the ability of sub-seismic MTDs to act as hydraulic seals and their effect on hydrocarbon production, and/or carbon storage. We investigate the gap between seismically resolvable and sub-seismic MTDs, and transitional flow deposits on long-term reservoir performance in this analysis of a small (

Beschreibung: Micritic sediments containing dark, discrete, organic-rich burrows, situated in a light grey background carbonate mud, were deposited over a broad geographic area in deep-shelf, bathyal, and basinal environments in the western margin of the Tethys Ocean during the Early and Middle Jurassic. These hemipelagic deposits represent a distinct depositional regime marked by low-energy, soft-bottom and only locally dysoxic environments. Still, it is unclear whether the trace fossil assemblages occurring in these deposits pertain to a network of several community types -the ichnotaxa differing from basin to basin- or else a single community of environmentally broad-ranging, burrow-producing species. Lower Jurassic trace-fossil assemblages are found in the Western Carpathians and in the Subbetic, Betic Cordillera, i.e. in basins separated by more than 2,000 km in their original palaeogeographic areas. The stereotypical Chondrites and Zoophycos trace-fossil assemblages that occur in the analysed deposits share two ichnogenera of distinctive morphology (Lamellaeichnus, Teichichnus). Agglutinated foraminifera Bathysiphon occurs together with the described trace fossil assemblage and determines epibenthic palaeoenvironmental conditions. In the Western Carpathians, a Lamellaeichnus-dominated assemblage alternates with a Zoophycos-dominated assemblage in small, metre-scale cycles in the upper Pliensbachian, and the proportion of the Zoophycos assemblage increases stratigraphically upwards, probably owing to reduced basin ventilation during the early Toarcian. Within the southern Iberian Palaeomargin, represented by the Betic Cordillera, Zoophycos is scarce in the facies.

Beschreibung: Recently, time-lapse seismic (4D seismic) has been steadily used to demonstrate the relation between field depletion and 4D seismic response, subsequently to have more efficient field management. A key component of the reservoir monitoring is the knowledge of fluid movement and pressure variations. This information is vital to assist infill drillings and a trustworthy source to update reservoir models, consequently improving model-based reservoir management and decision-making process. However, in practice the 4D seismic interpretation of reservoirs with multipart production regime possesses ambiguities through different levels of uncertainty. Complex nature of some 4D seismic signals emphasizes the roles of competing effects, geology, rock and fluid interactions. Hence, a reliable 4D interpretation requires an interdisciplinary approach entailing data analysis and insights from geophysics, engineering and geology. In this research, a step-wise workflow was introduced to reduce uncertainties in the 4D seismic interpretation and provide diagnoses to perform better reservoir surveillance. In parallel, the workflow expresses the use of engineering data analysis to conduct a consistent interpretation and encompasses the 3D and 4D seismic attributes with engineering data analysis. This study is implemented in a Brazilian heavy-oil offshore field where production started in 2013. The field experienced intense production activity up to 2016, making the deep-water field an ideal candidate to explore the challenges in interpreting complex 4D signals. Beyond these challenges, significant understanding of reservoir behavior is obtained and suggestions are made to improve the reservoir simulation model, which could support reservoir engineers with data assimilation applications.

Beschreibung: We have applied multiscale deformable-layer tomography (DLT) to build a laterally varying velocity model, using a single-offset vertical seismic profile (VSP) data set acquired for a salt proximity survey in southern Texas. The purpose of the VSP survey is to delineate the 2D salt flank using the P-wave reflections. Previous study has identified an anhydrate layer as the cap rock of the salt dome. The large impedance contrasts of this anhydrite layer generate strong downgoing P (sediment)-S (anhydrite)-P (salt) waves recorded by downhole geophones. Incidentally, the P-S-P-waves have similar traveltimes as those of the P-wave salt flank reflections, thus contaminating the imaging of the salt flank. Identifying shear-mode contamination requires an accurate velocity model of anhydrite. However, the extremely poor coverage of the single-offset VSP greatly challenges tomographic techniques to determine the lateral velocity variation. We tackle this problem using multiscale DLT, which characterizes the velocity field by a set of deformable layers. We constrain the layer velocities using the check-shot data and invert for the geometric variation. The inverted model indicates that the anhydrite layer has a “thick-thin-thick” lateral variation with offset, and the S-wave in the anhydrite layer helps in imaging the P-S-P-waves along the well track. The estimated anhydrite layer geometry is validated by the kinematic accuracies of P-waves in the data domain and P-S-P-waves in the image domain. Some in-salt dipping structures are determined by multiscale DLT as well. This field data example indicates that multiscale DLT is feasible for estimating velocities using VSP data of the single-offset situation. An accurate velocity model is the key for modeling and adaptive subtraction of the shear-mode contamination related to the salt geometry.

Beschreibung: Amplitude variation with incidence angle (AVA) analysis is an essential tool for discriminating lithology in the hydrocarbon reservoirs. Compared with the traditional AVA inversion using only P-wave information, joint AVA inversion using PP and PS seismic data provides better estimation of rock properties (e.g., density, P- and S-wave velocities). At present, the most used AVA inversions depend on the approximations of Zoeppritz equations (e.g., Shuey and Aki-Richards approximations), which are not suitable for formations with strong contrast interfaces and seismic data with large incidence angles. Based on the previous derivation of accurate Jacobian matrix, we find that the sign of each partial derivative of reflection coefficient with respect to P-, S-wave velocities and density changes across the interface, represents good indicator for the reflection interfaces. Accordingly, we propose an adaptive stratified joint PP and PS AVA inversion using the accurate Jacobian matrix that can automatically obtain the layer information and can be further used as a constraint in the inversion of in-layer rock properties (density, P- and S-wave velocities). Due to the use of the exact Zoeppritz equations and accurate Jacobian matrix, this proposed inversion method is more accurate than traditional AVA inversion methods, has higher computational efficiency and can be applied to seismic wide-angle reflection data or seismic data acquired for formations with strong contrast interfaces. The model study shows that this proposed inversion method works better than the classical Shuey and Aki-Richards approximations at estimating reflection interfaces and in-layer rock properties. It also works well in handling a part of the complex Marmousi 2 model and real seismic data.

Beschreibung: Summary Landslides can cause devastating damage. In particular, heavy rainfall-triggered landslides pose a chain of natural hazards. However, such events are often difficult to detect, leaving the physical processes poorly understood. Here we apply a novel surface-wave detector to detect and locate landslides during the transit of Typhoon Talas 2011. We identify multiple landslides triggered by Typhoon Talas, including a landslide in the Tenryu Ward, Shizuoka prefecture, Japan, ∼400 km east from the typhoon track. The Tenryu landslide displaced a total volume of 1.2 − −1.5 × 106 m. The landslide is much smaller than those detected by using globally recorded surface waves, yet the event generated coherent seismic signals propagating up to 3000 km away. Our observations show that attributes of small and large landslides may follow the same empirical scaling relationships, indicating possible invariant failure mechanisms. Our results also suggest an alerting technology to detect and locate landslides with a sparse seismic network.

Beschreibung: Summary We studied the broad-band spectra of the 8 largest earthquakes that have occurred in Chile in the last 25 years using strong-motion records and 1-Hz high-rate GNSS (cGNSS) data. To avoid the numerical instability problem with the double integration of the accelerograms, we computed velocity spectra integrating the acceleration time series in the spectral domain and compared them to time-differentiated the cGNSS displacement records. To compute the velocity spectrum, we used a multitaper algorithm so as to provide stability over the entire spectral band. We found that the velocity spectra of records obtained close to the main rupture of the earthquakes are different from classical Aki and Brune spectra. The velocity spectrum of large events in Chile presents a flat trend at low frequencies produced by the near-field waves. This trend converges at low frequencies to the static displacement as determined from GNSS data. For different magnitude earthquakes, we observe a transition in the ground-velocity spectrum from a decay of ${f^{ - 1}}$ at high frequencies and a flat trend at low frequencies to a more classical model with a peak at the corner frequency. The source-station distance influences the shape of the velocity spectrum at low frequencies, but there is no simple rule for the records available at present. At intermediate frequencies, the spectra are controlled by surface waves and S waves. We found a transition in the velocity spectrum for the 2014 Iquique earthquake, which indicates a change in the decay of the spectrum for stations at distances greater than ∼200 km. Finally, we show that the flat low-frequency trend of the velocity spectra determined from accelerograms, and the peak ground-displacement (PGD) determined from GNSS data scales with the moment to the power 2/3.

Beschreibung: Early Paleozoic trondhjemites, gneissic granites and alkali granites in southern Yunnan preserve important records of the tectonic evolution of the Prototethyan Ocean and regional correlations. Zircon ages suggest that these granitoids were emplaced from 476 to 436 Ma. The trondhjemites are characterized by high Na2O and low K2O contents, with εNd(t) values of −1.9 to −3.5 and εHf(t) values of −2.8 to +3.9. The trondhjemites were derived from an amphibolite source with a juvenile mafic component. The gneissic granites belong to the metaluminous low-K calc-alkaline series with an εNd(t) value of −6.2 and εHf(t) values of −5.0 to −0.4. The alkali granites belong to the high-K calc-alkaline series and yield εNd(t) values of −10.1 to −10.7 and εHf(t) values of −7.9 to −2.3. The gneissic granites were derived from an ‘ancient' lower mafic crust, whereas the alkali granites were derived from a meta-sedimentary source. These granitoids were formed during the subduction of the Prototethyan Ocean beneath the Simao Block and can be compared with similar igneous rocks from the Truong Son and Tam Ky-Phuoc Son zones in southern Laos. Our study, along with Early Paleozoic igneous suites from southern Laos, central Vietnam and the Malay Peninsula, suggests an arc–back-arc system along the northern margin of Gondwana.Supplementary material: Tables of zircon U–Pb and in-situ Hf and geochemical data are available at https://doi.org/10.6084/m9.figshare.c.5322386

Beschreibung: Observation of basin-scale networks of sandstone intrusions are described from subsurface studies and outcrop locations. Regional scale studies are prevalent in the volume and two new regionally significant subsurface sand injection complexes are described. Higher resolution studies, both outcrop and subsurface, show the small-scale complexity but high level of connectedness of sandstone intrusions. Discordance with bedding at all scales is diagnostic of sandstone intrusions. The propensity of hydraulic fractures to develop and fill with fluidised sand in a broad range of host rocks is demonstrated by examples from metamorphic and magmatic basement, and lignite. Terminology used to describe sandstone intrusions and other elements of sand injection complexes is diverse.

Beschreibung: Full waveform inversion, a high-resolution seismic imaging method, is known to require sufficiently accurate initial models to converge toward meaningful estimations of the subsurface mechanical properties. This limitation is due to the non-convexity of the least-squares distance with respect to kinematic mismatch. We propose a comparison of five misfit functions promoted recently to mitigate this issue: adaptive waveform inversion, instantaneous envelope, normalized integration, and two methods based on optimal transport. We explain which principles these methods are based on and illustrate how they are designed to better handle kinematic mismatch than a least-squares misfit function. By doing so, we can exhibit specific limitations of these methods in canonical cases. We further assess the interest of these five approaches for application to field data based on a synthetic Marmousi case study. We illustrate how adaptive waveform inversion and the two methods based on optimal transport possess interesting properties, making them appealing strategies applicable to field data. Another outcome is the definition of generic tools to compare misfit functions for full-waveform inversion.

Beschreibung: Gravity and magnetic data have unique advantages for studying the lateral extents of geological bodies. There is a class of methods for edge recognition called the maximum-edge-recognition methods that use their extreme values to locate the edges of geological bodies. These methods include the total horizontal derivative, the analytic signal amplitude, the theta map, and the normalized standard deviation. These are all first-order derivative-based techniques. There are also higher-order derivative-based methods that are derived from the first-order filters, for example, the total horizontal derivative of the tilt angle. We present an edge recognition filter that is based on the idea of the normalized vertical derivatives of existing methods. For each maximum-edge-recognition method, we first calculate its nth-order vertical derivative and then use thresholding to locate its peaks. The peak values are subsequently normalized by the values of the original maximum-edge-recognition method. Testing on synthetic and real data shows that the normalized vertical derivatives of the maximum-edge-recognition methods have higher accuracy, better lateral resolution and are more interpretable than existing techniques, and thus are a worthwhile addition to the set of edge-detection tools for potential-field data.

Beschreibung: Conventional finite-difference frequency-domain (FDFD) methods can describe wave attenuation and velocity dispersion more easily than time-domain methods. However, there are significant challenges associated with computational costs for solving the linear system when frequency-domain methods are applied in models with large dimensions or fine-scale property variations. Direct-iterative solvers and parallel strategies attempt a tradeoff between memory and time costs. We follow the general framework of heterogeneous multiscale method and develop a multiscale FDFD approach to solve the Helmholtz equation with lower memory and time costs. To achieve this, the discrete linear system approximating the Helmholtz equation is constructed on a coarse mesh, making its dimension much smaller than that of conventional methods. The coefficient matrix in the linear system of dimension-reduction captures fine-scale heterogeneity in the media by coupling fine- and coarse-scale meshes. Several test models are used to verify the accuracy of our multiscale method and investigate potential sources of error. Numerical results demonstrate that our method accurately approximates the wavefields of fine-scale solutions at low frequencies of the source, and could produce solutions with small errors by reducing the size of the coarse mesh cells at high frequencies as well. Comparisons of computational costs with conventional FDFD methods show that the proposed multiscale method significantly reduces computation time and memory consumption.

Beschreibung: Summary Seismic imaging techniques such as elastic full waveform inversion (FWI) have their spatial resolution limited by the maximum frequency present in the observed waveforms. Scales smaller than a fraction of the minimum wavelength cannot be resolved, and only a smoothed, effective version of the true underlying medium can be recovered. These finite-frequency effects are revealed by the upscaling or homogenization theory of wave propagation. Homogenization aims at computing larger scale effective properties of a medium containing small-scale heterogeneities. We study how this theory can be used in the context of FWI. The seismic imaging problem is broken down in a two-stage multiscale approach. In the first step, called homogenized full waveform inversion (HFWI), observed waveforms are inverted for a smooth, fully anisotropic effective medium, that does not contain scales smaller than the shortest wavelength present in the wavefield. The solution being an effective medium, it is difficult to directly interpret it. It requires a second step, called downscaling or inverse homogenization, where the smooth image is used as data, and the goal is to recover small-scale parameters. All the information contained in the observed waveforms is extracted in the HFWI step. The solution of the downscaling step is highly non-unique as many small-scale models may share the same long wavelength effective properties. We therefore rely on the introduction of external a priori information, and cast the problem in a Bayesian formulation. The ensemble of potential fine-scale models sharing the same long wavelength effective properties is explored with a Markov chain Monte Carlo algorithm. We illustrate the method with a synthetic cavity detection problem: we search for the position, size and shape of void inclusions in a homogeneous elastic medium, where the size of cavities is smaller than the resolving length of the seismic data. We illustrate the advantages of introducing the homogenization theory at both stages. In HFWI, homogenization acts as a natural regularization helping convergence toward meaningful solution models. Working with fully anisotropic effective media prevents the leakage of anisotropy induced by the fine scales into isotropic macro-parameters estimates. In the downscaling step, the forward theory is the homogenization itself. It is computationally cheap, allowing us to consider geological models with more complexity (e.g. including discontinuities) and use stochastic inversion techniques.

Beschreibung: Summary Particle Image Velocimetry (PIV), a method based on image cross-correlation, is widely used for obtaining velocity fields from time series of images of deforming objects. Rather than instantaneous velocities, we are interested in reconstructing cumulative deformation, and use PIV-derived incremental displacements for this purpose. Our focus is on analogue models of tectonic processes, which can accumulate large deformation. Importantly, PIV provides incremental displacements during analogue model evolution in a spatial reference (Eulerian) frame, without the need for explicit markers in a model. We integrate the displacements in a material reference (Lagrangian) frame, such that displacements can be integrated to track the spatial accumulative deformation field as a function of time. To describe cumulative, finite deformation, various strain tensors have been developed, and we discuss what strain measure best describes large shape changes, as standard infinitesimal strain tensors no longer apply for large deformation. PIV or comparable techniques have become a common method to determine strain in analogue models. However, the qualitative interpretation of observed strain has remained problematic for complex settings. Hence, PIV-derived displacements have not been fully exploited before, as methods to qualitatively characterize cumulative, large strain have been lacking. Notably, in tectonic settings, different types of deformation - extension, shortening, strike-slip - can be superimposed. We demonstrate that when shape changes are described in terms of Hencky strains, a logarithmic strain measure, finite deformation can be qualitatively described based on the relative magnitude of the two principal Hencky strains. Thereby, our method introduces a physically meaningful classification of large 2D strains. We show that our strain type classification method allows for accurate mapping of tectonic structures in analogue models of lithospheric deformation, and complements visual inspection of fault geometries. Our method can easily discern complex strike-slip shear zones, thrust faults and extensional structures and its evolution in time. Our newly developed software to compute deformation is freely available and can be used to post-process incremental displacements from PIV or similar autocorrelation methods.

Beschreibung: Sue Bowler remembers the Guernsey-born astronomer

Beschreibung: Jacqueline Mitton examines one strand in the life and research of the pioneering astronomer

Beschreibung: Robert Massey ponders the potential post-pandemic world for those working in astronomy and geophysics, following an online meeting in January

Beschreibung: A digest of some of the big stories from the past two months

Beschreibung: This study discusses the petrological and geochemical features of two monogenetic Miocene volcanoes, Mount Early and Sheridan Bluff, which are the above-ice expressions of Earth's southernmost volcanic field located at c. 87° S on the East Antarctic Craton. Their geochemistry is compared to basalts from the West Antarctic Rift System to test affiliation and resolve mantle sources and cause of melting beneath East Antarctica. Basaltic lavas and dykes are olivine-phyric and comprise alkaline (hawaiite and mugearite) and subalkaline (tholeiite) types. Trace element abundances and ratios (e.g. La/Yb, Nb/Y, Zr/Y) of alkaline compositions resemble basalts from the West Antarctic rift and ocean islands (OIB), while tholeiites are relatively depleted and approach the concentrations levels of enriched mid-ocean ridge basalt (E-MORB). The magmas evolved by fractional crystallization with contamination by crust; however, neither process can adequately explain the contemporaneous eruption of hawaiite and tholeiite at Sheridan Bluff. Our preferred scenario is that primary magmas of each type were produced by different degrees of partial melting from a compositionally similar mantle source. The nearly simultaneous generation of lower degrees of melting to produce alkaline types and higher degrees of melting forming tholeiite was most likely to have been facilitated by the detachment and dehydration of metasomatized mantle lithosphere.

Beschreibung: Two small monogenetic volcanoes are exposed at Mount Early and Sheridan Bluff, in the upper reaches of Scott Glacier. In addition, the presence of abundant fresh volcanic detritus in moraines at two other localities suggests further associated volcanism, now obscured by the modern Antarctic ice sheet. One of those occurrences has been attributed to a small subglacial volcano only c. 200 km from South Pole, making it the southernmost volcano in the world. All of the volcanic outcrops in the Scott Glacier region are grouped in a newly defined Upper Scott Glacier Volcanic Field, which is part of the McMurdo Volcanic Group (Western Ross Supergroup). The volcanism is early Miocene in age (c. 25–16 Ma), and the combination of tholeiitic and alkaline mafic compositions differs from the more voluminous alkaline volcanism in the West Antarctic Rift System. The Mount Early volcano was erupted subglacially, when the contemporary ice was considerably thicker than present. By contrast, lithologies associated with the southernmost volcano, currently covered by 1.5 km of modern ice, indicate that it was erupted when any associated ice was either much thinner or absent. The eruptive setting for Sheridan Bluff is uncertain and is still being investigated.

Beschreibung: The Erebus Volcanic Province is the largest Neogene volcanic province in Antarctica, extending c. 450 km north–south and 170 km wide east–west. It is dominated by large central volcanoes, principally Mount Erebus, Mount Bird, Mount Terror, Mount Discovery and Mount Morning, which have sunk more than 2 km into underlying sedimentary strata. Small submarine volcanoes are also common, as islands and seamounts in the Ross Sea (Terror Rift), and there are many mafic scoria cones (Southern Local Suite) in the Royal Society Range foothills and Dry Valleys. The age of the volcanism ranges between c. 19 Ma and present but most of the volcanism is

Beschreibung: The Hudson Strait Platform and basins Tectono-Sedimentary Element (HSPB TSE) is part of a major topographic feature that connects Hudson Bay and Foxe Basin with the Labrador Sea in the Canadian Arctic. The Paleozoic succession (Ordovician to Silurian) unconformably overlies the Precambrian basement and reaches a maximum preserved thickness of less than 600 m on the islands. High-resolution marine seismic data indicate that the offshore part of the Hudson Strait is underlain by several fault-controlled sub-basins with a half-graben geometry. The sedimentary succession in the sub-basins is thicker than the one preserved in nearby islands and includes an upper sedimentary package for which the nature and age remain poorly constrained. Upper Ordovician source rocks have been mapped onshore. Known potential reservoir rocks consist of Ordovician clastics and Ordovician-Silurian reefs and dolostones.

Beschreibung: Summary It is generally accepted that melt extraction from the mantle at mid-ocean ridges is concentrated in narrow regions of elevated melt fraction called channels. Two feedback mechanisms have been proposed to explain why these channels grow by linear instability: shear flow of partially molten mantle and reactive flow of the ascending magma. These two mechanisms have been studied extensively, in isolation from each other, through theory and laboratory experiments as well as field and geophysical observations. Here, we develop a consistent theory that accounts for both proposed mechanisms and allows us to weigh their relative contributions. We show that interaction of the two feedback mechanisms is insignificant and that the total linear growth rate of channels is well-approximated by summing their independent growth rates. Furthermore, we explain how their competition is governed by the orientation of channels with respect to gravity and mantle shear. By itself, analysis of the reaction-infiltration instability predicts the formation of tube-shaped channels. We show that with the addition of even a small amount of extension in the horizontal, the combined instability favours tabular channels, consistent with the observed morphology of dunite bodies in ophiolites. We apply the new theory to mid-ocean ridges by calculating the accumulated growth and rotation of channels along streamlines of the solid flow. We show that reactive flow is the dominant instability mechanism deep beneath the ridge axis, where the most unstable orientation of high-porosity channels is sub-vertical. Channels are then rotated by the solid flow away from the vertical. The contribution of the shear-driven instability is confined to the margins of the melting region. Within the limitations of our study, the shear-driven feedback does not appear to be responsible for significant melt focusing or for the shallowly dipping seismic anisotropy that has been obtained by seismic inversions.

Beschreibung: The Devonian Period was a crucial interval in the evolution of plants. During its 60 myr duration, it witnessed the successive evolution of roots, wood, trees and forests, and many of the biogeomorphic phenomena that operate in modern terrestrial environments came online for the first time. The Old Red Sandstone (ORS) of Svalbard consists of a near-continuous Silurian to Late Devonian record of land plant-colonized sedimentary environments and provides a perfect natural laboratory to aid understanding of the facies signatures and evolution of these phenomena. Here we describe and illustrate a catalogue of ORS features that provide evidence for the stepwise appearance of novel plant-sediment interactions, including: preserved plant material and rooting structures, early large woody debris accumulations, cannel coal deposits, and the oldest known vegetation-induced sedimentary structures, in addition to vegetation-influenced motifs of elevated mudrock content and complex alluvial sand bodies. These characteristics are combined to reconstruct changes to non-marine environments in this Devonian ‘landscape factory’. In addition to tectonic and climate influences, plant evolution first served as a control on the construction of the sedimentary record during this period and has persisted as a fundamental influence on Earth surface processes and landforms ever since.

Beschreibung: Gaussberg is a nunatak composed of lamproite pillow lava situated on the coast of East Antarctica. It is the most isolated Quaternary volcanic centre in Antarctica but it is important palaeoenvironmentally and petrologically out of all proportion to its small size. The edifice has a likely low, shield-like, morphology c. 1200 m high and possibly up to 10 km wide, which is unusually large for a lamproite construct. Gaussberg was erupted subglacially at 56 ± 5 ka, which places it late in the last glacial, close to the peak of marine isotope stage 3. The coeval ice sheet was c. 1300 m thick, and c. 420 m has been removed from the ice surface since Gaussberg erupted. Lamproite is a rare ultrapotassic mantle-derived magma, and Gaussberg is one of two type examples worldwide. Although traditionally considered as related in some way to the Kerguelen plume, it is more likely that the Gaussberg magma is a product of a separate magmatic event. It is ascribed to the storage and long-term (Gy) isolation of sediment emplaced by subduction in the Transition Zone of the deep mantle, followed by entrainment and subsequent melting in a plume.

Beschreibung: Summary To understand earth processes, geoscientists infer subsurface earth properties such as electromagnetic resistivity or seismic velocity from surface observations of electromagnetic or seismic data. These properties are used to populate an earth model vector, and the spatial variation of properties across this vector sheds light on the underlying earth structure or physical phenomenon of interest, from groundwater aquifers to plate tectonics. However, to infer these properties the spatial characteristics of these properties need to be known in advance. Typically, assumptions are made about the length scales of earth properties, which are encoded a priori in a Bayesian probabilistic setting. In an optimisation setting, appeals are made to promote model simplicity together with constraints which keep models close to a preferred model. All of these approaches are valid, though they can lead to unintended features in the resulting inferred geophysical models owing to inappropriate prior assumptions, constraints or even the nature of the solution basis functions. In this work it will be shown that in order to make accurate inferences about earth properties, inferences can first be made about the underlying length scales of these properties in a very general solution basis. From a mathematical point of view, these spatial characteristics of earth properties can be conveniently thought of as “properties” of the earth properties. Thus, the same machinery used to infer earth properties can be used to infer their length scales. This can be thought of as an “infer to infer” paradigm analogous to the “learning to learn” paradigm which is now commonplace in the machine learning literature. However, it must be noted that (geophysical) inference is not the same as (machine) learning, though there are many common elements which allow for cross-pollination of useful ideas from one field to the other, as is shown here. A non-stationary trans-dimensional Gaussian Process (TDGP) is used to parameterise earth properties, and a multi-channel stationary TDGP is used to parameterise the length scales associated with the earth property in question. Using non-stationary kernels, i.e., kernels with spatially variable length scales, models with sharp discontinuities can be represented within this framework. As GPs are multi-dimensional interpolators, the same theory and computer code can be used to solve geophysical problems in 1D, 2D and 3D. This is demonstrated through a combination of 1D and 2D non-linear regression examples and a controlled source electromagnetic (CSEM) field example. The key difference between this and previous work using TDGP is generalised nested inference and the marginalisation of prior length scales for better posterior subsurface property characterisation.

Beschreibung: Summary The Z-transform of a complex time signal (or the analytic signal of a real signal) is equal to the Z-transform of a prediction error divided by the Z-transform of the prediction error operator. This inverse is decomposed into a sum of partial fractions, which are used to obtain impulse response operators formed by non-causal filters which complex-conjugate symmetric coefficients. The time-components are obtained by convolving the filters with the original signal, and the peak frequencies, corresponding to the poles of the prediction error operator, are used for mapping the time-components into frequency components. For non-stationary signals, this decomposition is done in sliding time windows, and the signal component values, in the middle of each window, are attributed to the peak value of its frequency response which corresponds to the pole of this partial fraction component. The result is an exact, but non-unique, time-frequency representation of the input signal. A sparse signal decomposition can be obtained by summing along the frequency axis in patches with similar characteristics in the time-frequency domain. The peak amplitude frequency of each new time component is obtained by computing a scalar prediction error operator in sliding time windows, resulting in a sparse time-frequency representation. In both cases, the result is a time-frequency matrix where an estimate of the frequency content of the input signal can be obtained by summation over the time variable. The performance of the new method is demonstrated with excellent results on a synthetic time signal, the LIGO gravitational wave signal, and on seismic field data.

Beschreibung: Summary Physical properties of near-surface soil and rock layers play a fundamental role in the seismic site effects analysis, being an essential element of seismic hazard assessment. Site-specific mechanical properties (i.e. shear- and compressional-wave velocities and mass density) can be inferred from surface wave dispersion and horizontal-to-vertical or ellipticity data by non-linear inversion techniques. Nevertheless, results typically exhibit significant inherent non-uniqueness as different models may fit the data equally well. Standard optimization inversion techniques minimize data misfit, resulting in a single representative model, rejecting other models providing similar misfit values. An alternative inversion technique can be formulated in the Bayesian framework, where the posterior probability density on the model space is inferred. This paper introduces an inversion approach of surface wave dispersion and ellipticity data based on a novel multizonal transdimensional Bayesian formulation. In particular, we parameterize one-dimensional layered velocity models by the varying number of Voronoi nuclei, allowing us to treat the number of layers as an unknown parameter of the inverse problem. The chosen parameterization leads to the transdimensional formulation of the model space, sampled by a reversible jump Markov chain Monte Carlo algorithm to provide an ensemble of random samples following the posterior probability density of model parameters. The used type of the sampling algorithm controls a model complexity (i.e. the number of layers) self-adaptively based on the measured data's information content. The method novelty lies in the parsimonious selection of sampling models and in the multizonal formulation of prior assumptions on model parameters, the latter allows including additional site-specific constraints in the inversion. These assumptions may be based on, e.g. stratigraphic logs, standard penetration tests, known water table, and bedrock depth. The multizonal formulation fully preserves the validity of the transdimensional one, as demonstrated analytically. The resultant ensemble of model samples is a discrete approximation of the posterior probability density function of model parameters and associated properties (e.g. VS30, quarter-wavelength average velocity profile, and theoretical SH-wave amplification function). Although the ultimate result is the posterior probability density function, some representative models are selected according to data fit and maximum of the posterior probability density function. We first validate our inversion approach based on synthetic tests and then apply it to field data acquired from the active seismic survey and single-station measurements of ambient vibrations at the SENGL seismic station site in central Switzerland.

Beschreibung: Volcán Quizapu, Chile, is an under-monitored volcano that was the site of two historical eruptions: an effusive eruption in 1846-1847 and a Plinian eruption in 1932, both of which discharged ∼5 km3 (DRE) of lava and/or tephra. The majority of material erupted in both cases is trachydacite, nearly identical for each event. We present H2O-saturated, phase equilibrium experiments on this end-member dacite magma, using a pumice sample from the 1932 eruption as the main starting material. At an oxygen fugacity (fO2) of ∼NNO+0.2, the phase assemblage of An25-30 plagioclase + amphibole + orthopyroxene, without biotite, is stable at 865±10 °C and 110±20 MPa H2O pressure (PH2O), corresponding to ∼4 km depth. At these conditions, experiments also reproduce the quenched glass composition of the starting pumice. At slightly higher PH2O and below 860 °C, biotite joins the equilibrium assemblage. Because biotite is not part of the observed Quizapu phase assemblage, its presence places an upper limit on PH2O. At the determined storage PH2O of ∼110 MPa, H2O undersaturation of the magma with XH20fluid==0.87 would align Ptotal to mineral-based geobarometry estimates of ∼130 MPa. However, XH20fluid=1〈 1 is not required to reproduce the Quizapu dacite phase assemblage and compositions. A second suite of experiments at lower fO2 shows that the stability fields of the hydrous silicates (amphibole and biotite) are significantly restricted at NNO-2 relative to NNO+0.2. Additional observations of Quizapu lava and pumice samples support the existing hypothesis that rapid pre-eruptive heating drove the effusive 1846-1847 eruption, with important refinements. We demonstrate that microlites in the end-member dacite lavas are consistent with in situ crystallization (during ascent), rather than transfer from an andesite. In one end-member dacite lava, newly identified reverse zoning in orthopyroxene and incipient destabilization of amphibole are consistent with small degrees of heating. Our work articulates a clear direction for future Quizapu studies, which are warranted given the active nature of the Cerro Azul-Descabezado Grande volcanic axis.

Beschreibung: Summary Until now, the polar motion resonance (PMR) complex frequency has been determined in the seasonal and retrograde diurnal band of the polar motion. In this study this resonance is studied in the prograde diurnal band, where polar motion is mainly composed of periodic terms caused by the diurnal oceanic tide. The resonance parameters (period and quality factor) are encompassed in the frequency transfer function between generating tidal potential and polar motion, and can be estimated accordingly. To this aim, we gather three published sets of prograde diurnal terms determined from GNSS and VLBI, to which we append our own estimates based upon a processing of the VLBI delays over the period 1990-2020. Then, by fitting the PMR parameters so that the prograde diurnal terms match the corresponding components of the tide generating potential, we obtained a resonance period of about 401 days and an equivalent quality factor of −22, differing from the ones reigning in the seasonal band (PPMR ≈ 431 days; QPMR ≈ 56 − 255) and in the retrograde diurnal band (PPMR ≈ 380 days; QPMR ≈ −10). Our estimates confirm strikingly the theoretical prediction derived from the tidal ocean angular momentum derived from the FES 2014 ocean tide model.

Beschreibung: Calc-alkaline lamprophyres from the western Carpathians occur as dykes and sills in the crystalline complexes (predominantly granites and gneisses) of the Tatric Unit. Some of the lamprophyre dykes have been strongly overprinted by tectonism and hydrothermal fluid ingress. They have a similar mineralogical composition, and they are comprised of clinopyroxene, amphibole, biotite and plagioclase. Based on their modal composition, they can be classified as spessartites and kersantites, and based on their geochemical composition, most of them are of calc-alkali type. Lamprophyres from individual core complexes (e.g. the Malá Fatra Mountains and Nízke Tatry Mountains) exhibit variable Nb, Ta and Sr-Nd isotope signatures. These differences are probably due to compositional variations in the mantle source and/or the lower crust at the site of lamprophyre melt generation, or variable incorporation of crustal material. The age of the lamprophyres is Permian (ca. 265 Ma) based on U-Pb LA-ICP-MS dating of apatite micro-phenocrysts.

Beschreibung: Aspiculate demosponges are rarely described in geological history on account of the absence of spicules that are stable and resistant to degradation. One exception is the exquisite preservation of sponges without any mineralised skeletons discovered in Lagerstätten (e.g. the Burgess Shale). The Chengjiang Biota, an early example of a Burgess Shale-type Biota in South China (Cambrian Series 2, Stage 3), is one of the only examples of convincing aspiculate sponges until now. Here, we describe Vauxia pregracilenta sp. nov. and V. paraleioia sp. nov., as well as two poorly preserved vauxiid specimens (Vauxia sp.) in open nomenclature, from the Chengjiang Biota. V. pregracilenta has a fan-like holdfast and branches in various sizes, as well as a typical two-layered net-like skeleton, without spicules. The endosomal layer is hexagonal, while the dermal layer is sub-rectangular. V. paraleioia is characterised by a two-layered subconical skeleton, with the dermal layer ornamented with vertical surface grooves. The openings of the dermal and endosomal layers of V. paraleioia are both hexagonal but of different sizes. These newly discovered Vauxia species indicate that the aspiculate sponges were diversified in the early Cambrian period. Partial silicification of the fibres of aspiculate Vauxia are confirmed from the Chengjiang Biota.Thematic collection: This article is part of the Advances in the Cambrian Explosion collection available at: https://www.lyellcollection.org/cc/advances-cambrian-explosion

Beschreibung: Summary Bayesian inversion of electromagnetic data produces crucial uncertainty information on inferred subsurface resistivity. Due to their high computational cost, however, Bayesian inverse methods have largely been restricted to computationally expedient 1D resistivity models. In this study, we successfully demonstrate, for the first time, a fully 2D, trans-dimensional Bayesian inversion of magnetotelluric data. We render this problem tractable from a computational standpoint by using a stochastic interpolation algorithm known as a Gaussian process to achieve a parsimonious parametrization of the model vis-a-vis the dense parameter grids used in numerical forward modeling codes. The Gaussian process links a trans-dimensional, parallel tempered Markov chain Monte Carlo sampler, which explores the parsimonious model space, to MARE2DEM, an adaptive finite element forward solver. MARE2DEM computes the model response using a dense parameter mesh with resistivity assigned via the Gaussian process model. We demonstrate the new trans-dimensional Gaussian process sampler by inverting both synthetic and field magnetotelluric data for 2D models of electrical resistivity, with the field data example converging within 10 days on 148 cores, a non-negligible but tractable computational cost. For a field data inversion, our algorithm achieves a parameter reduction of over 32x compared to the fixed parameter grid used for the MARE2DEM regularized inversion. Resistivity probability distributions computed from the ensemble of models produced by the inversion yield credible intervals and interquartile plots that quantitatively show the non-linear 2D uncertainty in model structure. This uncertainty could then be propagated to other physical properties that impact resistivity including bulk composition, porosity and pore-fluid content.

Beschreibung: Summary Ambient wavefield data acquired on existing (so-called “dark fiber”) optical fiber networks using distributed acoustic sensing (DAS) interrogators allow users to conduct a wide range of subsurface imaging and inversion experiments. In particular, recorded low-frequency (

Beschreibung: The Jurassic and Cretaceous sedimentary infill of the Irish and Celtic Sea basins is intimately associated with the breakup of the supercontinent Pangea, and the opening of the Atlantic margin. Previous basin studies have constrained tectonism, basin uplift and sediment composition, but sediment provenance and routing have not received detailed consideration. Current hypotheses for basin infill suggest localised sediment sourcing throughout the Jurassic and Cretaceous, despite a dynamic tectonic and paleoenvironmental history spanning more than 100 million years. We present detrital zircon, white mica and apatite geochronology alongside heavy mineral data from five basins. Findings reveal that basin infill derived predominantly from distal sources with lesser periods of local sourcing. We deduce that tectonically induced marine transgression and regression events had a first-order control on distal versus proximal sedimentary sourcing. Additionally, tectonism which uplifted the Fastnet Basin region during the Middle–Late Jurassic recycled basin sediments into the connected Celtic and Irish Sea Basins. Detrital geochronology and heavy mineral evidence support three distinct provenance switches throughout the Jurassic and Cretaceous in these basins. Overall an integrated multi-proxy provenance approach provides novel insights to tectonic and environmental controls on basin infill as demonstrated in the Irish and Celtic Sea Basins.

Beschreibung: This contribution is an overview of the Early Jurassic dinoflagellate cysts of the Lusitanian Basin in Portugal, with particular emphasis on the effects of the Jenkyns Event (Toarcian Oceanic Anoxic Event - T-OAE) on the evolution of this planktonic group. We review and discuss data from 214 samples from six Lower Jurassic successions (upper Sinemurian-upper Toarcian) in the Lusitanian Basin. The late Pliensbachian radiation of dinoflagellate cysts was well recognised in this basin. The pre-Jenkyns Event interval is highly productive, with maximum abundance and species richness values. However, this palaeoenvironmental perturbation severely affected the evolution of this group for the remainder of the Early Jurassic. The prolonged recovery of the dinoflagellates in the Toarcian following the Jenkyns Event is not typical of the northern regions (Arctic and Boreal realms), where new species began to evolve earlier compared with southern European basins.

Beschreibung: Summary The stability of partial differential equations determines the properties of their solutions. This study focuses on the stability analysis of the equations describing wave propagation in fluids-saturated porous media. We briefly introduce the stability analysis method for the wave propagation equations and discuss the adverse effects on the solutions. In this way, the first part of this paper is mainly devoted to the analysis of the Tuncay and Corapcioglu's (TC) model, which describes the dynamic behaviour of porous media saturated with two immiscible fluids. It is pointed out that the TC model allows spatially bounded but time-exponentially exploding solutions and may yield unstable numerical results. Based on the deduced unstable factors, we construct a stable equivalent fluid (SEF) model. We rigorously analyze the stability of the SEF model using the energy method. For predicting the influence of saturation on wave velocity, the robustness of this model is preserved due to its consistency with the original TC model. Furthermore, the numerical simulations of the wave fields show that the results of the TC model exponentially increase with time after the initial effective wave signal, which does not occur in the SEF model curves. This indicates the necessity of considering the stability from a mathematical point of view during the construction of physical model. It could be useful to merge the mathematical stability theory with the geophysical wave propagation modelling theory.

Beschreibung: In the Southern Ribeira Belt (Brazil), the Itatins Complex, composed of ortho- and para-derived migmatites record granulite-facies metamorphism. Our work focuses on the ortho-derived migmatite, which displays features indicative of partial melting ranging from in-situ melting (cm-thick lenses and veins) to accumulations of melt as leucosomes and patches of diatexite. The unit comprises a biotite felsic granulite (metagranodiorite) with a residual neosome and a leucocratic neosome, generally referred to as leucosome. Phase equilibrium modelling using a modified bulk composition from the geochemical standard JG-1 (granodiorite) was carried out to evaluate whether a protolith of granodiorite bulk composition could produce the mineral assemblages observed in the residual neosome and the compositions of melt and residue produced by the partial melting. Isopleth thermobarometry of plagioclase, biotite and orthopyroxene points to P-T conditions of 0.7-1.0 GPa and 845-872 ºC. The modified granodiorite bulk composition generated just 10 mol % melt at these conditions, far lower than the amount estimated from outcrops (around 32 %). A T-XH2O pseudosection shows that increasing H2O content enhances the fertility of the bulk composition and it can become capable of producing melt volumes similar to that observed in outcrop. The melt compositions determined for partial melting of the modified granodiorite bulk composition in open- and closed-systems using THERMOCALC are granitic, and peraluminous. The melts have SiO2 contents between 52.02 mol % and 77.27 mol. % and show XFe values of 0.68 and 0.86. In K2O vs. Na2O+CaO, K2O vs. XFe and K2O-(FeO+MgO)-(Na2O+CaO) space natural samples (5 from the residual neosome and 2 from leucosome) from the Itatins Complex plot in the direction of the trends for melt and residue compositions obtained from modelling of the modified granodiorite bulk composition. The differences between modeled compositions and real rocks are interpreted to result from the interaction of melt with its residual counterpart. Chondrite-normalized REE patterns from the residual neosome show considerably greater depletion of HREE than geochemical standard JG-1 (granodiorite), and the two analyzed leucosomes show lower HREE contents. Zircon rims have REE patterns interpreted to indicate metamorphic recrystallization. U-Pb zircon Concordia ages indicate that the protolith formed during the Paleoproterozoic (2137 Ma) but that metamorphic recrystallization is Neoproterozoic (between 632 and 601 Ma). The metamorphic age pre-dates the intense granitic magmatism that affected the Southern Ribeira Belt between 600 and 570 Ma.

Beschreibung: Summary Uncovering the distribution of magnitudes and arrival times of aftershocks is a key to comprehending the characteristics of earthquake sequences, which enables us to predict seismic activities and conduct hazard assessments. However, identifying the number of aftershocks immediately after the main shock is practically difficult due to contaminations of arriving seismic waves. To overcome this difficulty, we construct a likelihood based on the detected data, incorporating a detection function to which Gaussian process regression (GPR) is applied. The GPR is capable of estimating not only the parameters of the distribution of aftershocks together with the detection function, but also credible intervals for both the parameters and the detection function. The property that the distributions of both the Gaussian process and aftershocks are exponential functions leads to an efficient Bayesian computational algorithm to estimate hyperparameters. After its validation through numerical tests, the proposed method is retrospectively applied to the catalog data related to the 2004 Chuetsu earthquake for the early forecasting of the aftershocks. The results show that the proposed method stably and simultaneously estimates distribution parameters and credible intervals, even within t ≤ 3h after the main shock.

Beschreibung: Late Permian to Early Jurassic strata in northern Peru allows us to carry out a seismo-stratigraphic, litho-tectonic and chemostratigraphic analysis connecting the Andean-Amazonian foreland basins of Huallaga, Ucayali, southern Marañón, and the Eastern Cordillera. This analysis and data integration from Ecuador to western Brazil and southern Peru and Bolivia, allow us to redefine the timing of the major documented tectonic phases and corresponding palaeogeographies of western Gondwana from the late Permian to Triassic. Three litho-tectonic sequences and four associated deformation stages are recognized: 1) A sequence, tectonic relaxation, during late Permian; 2) A-B intra-sequence, folding-and-thrusting attributed to a continuation in time of the Gondwanide Orogeny, during the Early to Middle Triassic; 3) B sequence, rifting, attributed to Gondwana breakup during the Middle and Late Triassic; and 4) C Sequence, thermal sag, during the Late Triassic. Evaporites and carbonates (A sequence) dominated a low subsidence basin with southern restricted marine inflow at the Permian-Triassic boundary. A novel palaeogeographic model for these evaporites suggests that this saline basin extended up to 50,000 km2 in a restricted environment area with a potential bullseye pattern. The last pulse of the Gondwanide Orogeny and associated fold and thrust belt (A-B intra-sequence) exhumed previous the sequence generating emerged areas with little to no sedimentation. Red beds (B sequence) characterize the rifting stage, representing the syn-depositional infill of continental grabens, likely extending to the Acre Basin in Brazil. Finally, during the thermal sag, a marine inflow likely from the northwestern part of Peru generated sedimentation of carbonates and evaporites (C Sequence) to the west and east of the Peruvian margin. This sediment differentiation was, in part, controlled by the existence of pre-existing grabens associated to the previous rifting stage. This interpretation, together with other evaporitic occurrences attributed here to a Late Triassic epoch in south and north Peru and west Brazil, suggest the existence of an evaporitic basin filling an undeformed area of probably ca. 170,000 km2. It is therefore suggestive of the existence of a Late Triassic (Norian to Rhaetian; 217 to 204 Ma) salt giant controlled by thermal sag in western Gondwana. Our results are of great relevance for any future interpretation related to mass extinctions, paleoclimatic analysis and ocean dynamics during the Permian and Triassic as well as natural resources distribution between Ecuador and Bolivia.

Beschreibung: In December 1903, Williams Speirs Bruce, leader of the Scottish National Antarctic Expedition, arrived in Buenos Aires and contacted Francisco Moreno, director of the Museo de La Plata to request his assistance. Bruce asked Moreno to be an intermediary with the Argentine government and to facilitate Bruce's wish for Argentina to take over the meteorological station that the Scottish expedition had established on Laurie Island (South Orkney Islands). Moreno was please to provide the necessary assistance and was instrumental in Bruce achieving his ambition. As a gesture of appreciation, before leaving Buenos Aires Bruce presented a small collection of Laurie Island rock specimens to Moreno as a donation to the Museo de La Plata. This donation initiated the museum's Antarctic collection.

Beschreibung: Ed Cliver, Lyndsay Fletcher and Hugh Hudson are looking for a photograph of Richard Carrington. Can you help?

Beschreibung: Richard Ellis recalls the life of the eminent cosmic-ray physicist

Beschreibung: Claire Davies recounts a successful outreach project based on Saturn and Jupiter's 2020 conjunction, a celebration of closeness in a year dominated by distancing

Beschreibung: First-arrival traveltime tomography (FATT) is used to delineate shallow velocity structures to identify static effects in oil exploration as well as to characterize the near surface for geotechnical purposes. Because FATT is generally used for land seismic data processing, it becomes necessary to consider irregular topography especially when performing wave-based tomography. However, the standard Cartesian finite-difference method cannot properly handle irregular topography. Hence, the embedded boundary method (EBM) is incorporated into the frequency-domain damped-wave equation in order to correctly describe irregular topography. The developed modeling algorithm is used to calculate first-arrival traveltimes and to perform FATT. The EBM-based modeling algorithm accurately describes the irregular surfaces of numerical velocity models on a regular mesh by exploiting the mirror image principle. The accuracy of the EBM-based traveltime calculation is validated by using two homogeneous velocity models with dipping and complex surfaces. The validation results demonstrate that the proposed algorithm is unaffected by the staircase approximation. The FATT is then applied to synthetic and real data to demonstrate the applicability of the developed algorithm to velocity models with complex topography. For the real data example, the inverted velocity model is used to apply static corrections. The processing results demonstrate an improvement in the continuity of seismic events, thus confirming the accuracy of the developed FATT method.

Beschreibung: Sumedh Anathpindika reviews some recent results that shed new light on the dynamical evolution of giant molecular clouds and discusses the impact of ambient environment on their ability to form stars

Beschreibung: Leslie Morrison and Catherine Hohenkerk consider the future direction of the leap second

Beschreibung: An RAS Specialist Discussion Meeting organized by Malcolm Druett and Ben Snow on 13 November 2020 drew a broad international audience to explore a tricky layer of the Sun

Beschreibung: This study combines electrical resistivity tomography (ERT), geological information from boreholes and outcrops, and hydrogeologic data to investigate field-scale fault-zone cementation of the Loma Blanca Fault in the Rio Grande Rift. We collected electrical resistivity data from 16 transects and geological samples from 29 boreholes (completed as groundwater wells to 30 m depth) across and around the fault. 2D ERT profiles, whose interpretations are constrained by geological data, indicate: (1) a high resistivity zone in cemented portions of the fault below the water table, and (2) in the unsaturated zone, a low resistivity feature along the cemented portions of the fault. The high resistivity zone below the water table is consistent with a 10% reduction in porosity due to the fault zone cementation. With the same porosity in the unsaturated zone, the low resistivity feature in the cemented fault zone is consistent with saturation 〉0.7, in contrast to saturation 0.2-0.7 for sediment outside of the cemented fault zone. In addition, subsurface samples and ERT profiles delineate a buttress unconformity (i.e., steeply dipping erosional contact) corresponding to a paleo-valley margin. This unconformity truncates the cemented fault zone and separates Pliocene axial-fluvial sand (deposited by an ancestral Rio Grande) from late Quaternary sand and gravel (deposited by the Rio Salado, a Rio Grande tributary). The cemented fault zone in the southern portion of the study area is a hydrogeologic barrier; north of the buttress unconformity, where the cemented fault zone has been removed by erosion, the fault is not a hydrogeologic barrier. The integration of geological, geophysical, and hydrogeological observations was key to developing our understanding of this complex system, and allowed us to demonstrate the utility of ERT in detecting subsurface fault-zone cementation.

Beschreibung: Tsunami waves were observed along the Bulgarian Black Sea coastline on May 7, 2007. The maximum rise and fall of the sea-level were 1.2 versus 2.0 m, respectively, with wave oscillations between 4-8 minutes. At first submarine landsliding and then later on atmospheric disturbance were suggested as the cause of the tsunami. Numerical modeling by Gusev et al. (2019), assuming a landslide displacing 30 to 60 million m3 material on the slope with a thickness range of over 20-40 m, could produce similar characteristics of the recorded tsunami. In the recent model the landslide initiated on the shelf at a water depth of 100 m with a runout of about 20 km into 1000 m water depth. Subsequent and recent numerical modeling suggested that the failure may have initiated on the slope, anywhere between 200 to 1500 m seafloor depth. The runout of the transported sediments in these latest model was at 1850 m water depth. Just a few years after the tsunami, OMV and its joint venture partners, TOTAL and Repsol, acquired modern deepwater data sets in the same area where the submarine landsliding was assumed to occur. These data sets included multi-beam swath bathymetry area and acquisition of a 3D reflection seismic data. These data sets offer a possibility to establish the presence of speculative submarine landslide responsible for the tsunami, with its geometry and nature. Our results provide direct evidence for the occurrence of large, non-seismic, catastrophic sediment failures along the Bulgarian coast. In this study we illustrate Quaternary submarine landslides on 3D seismic reflection data immediately below than the one responsible for the 2007 event, besides we also briefly point out the potential interpretation pitfall related to sediment waves and mass transport complexes.

Beschreibung: The STACK (Sooner Trend in the Anadarko [Basin] in Canadian and Kingfisher counties) play primarily produces oil and gas from Mississippian strata. The interval consists of interbedded argillaceous mudstones and calcareous siltstones. Such contrast in rock composition is linked directly to the mechanical stratigraphy of the strata. Brittle (calcareous siltstones) and ductile beds (argillaceous mudstones) are related to the sequence-stratigraphic framework at different scales. We use seismic and well log data to estimate and map geomechanical properties distribution and interpret the mechanical stratigraphy of rocks within the Mississippian strata. First, we defined the parasequences that form the main reservoir zones of the Meramecian-Mississippian strata. Once we established the stratigraphic framework, we estimated and compared rock brittleness index using two independent laboratory-based measurements from the core. The first method MIDBI (MIneralogical Derived Brittleness Index), uses mineralogical composition inverted from FTIR (Fourier Transform InfraRed spectroscopy) analyses, whereas the second method MEDBI (MEchanical Derived Brittleness Index), involves measurements of compressional and shear velocities from core plugs.We use the data from core-plug velocity measurements and well logs and an Artificial Neural Network (ANN) approach to establish relationships between geomechanical properties, well logs, and acoustic impedance values. We then applied these relationships to generate 3-D geomechanical models constrained to seismic volumes. The resulting grid distributions illustrate the stratigraphic variability of the properties at the parasequence scale. Overall, brittle strata decrease in thickness and abundance basinward as the frequency of interbedded brittle and ductile zones increases and gradually transitions into thin calcite-cemented siltstones and clay-rich mudstones. Analysis of the production performance of selected horizontal wells drilled within the Mississippian strata shows that the proportion of brittle and ductile rocks along the well path drilled and the drilled area vertical stacking pattern play a significant role in the hydrocarbon production for these Mississippian units.