Big Data Analytics for Earth Sciences: the EarthServer approach

Big Data Analytics is an emerging field since massive storage and computing capabilities have been made available by advanced e-infrastructures. Earth and Environmental sciences are likely to benefit from Big Data Analytics techniques supporting the processing of the large number of Earth Observation datasets currently acquired and generated through observations and simulations. However, Earth Science data and applications present specificities in terms of relevance of the geospatial information, wide heterogeneity of data models and formats, and complexity of processing. Therefore, Big Earth Data Analytics requires specifically tailored techniques and tools. The EarthServer Big Earth Data Analytics engine offers a solution for coverage-type datasets, built around a high performance array database technology, and the adoption and enhancement of standards for service interaction (OGC WCS and WCPS). The EarthServer solution, led by the collection of requirements from scientific communities and international initiatives, provides a holistic approach that ranges from query languages and scalability up to mobile access and visualization. The result is demonstrated and validated through the development of lighthouse applications in the Marine, Geology, Atmospheric, Planetary and Cryospheric science domains.

Big Data Analytics for Earth Sciences: the EarthServer approach

Big Data Analytics is an emerging field since massive storage and computing capabilities have been made available by advanced e-infrastructures. Earth and Environmental sciences are likely to benefit from Big Data Analytics techniques supporting the processing of the large number of Earth Observation datasets currently acquired and generated through observations and simulations. However, Earth Science data and applications present specificities in terms of relevance of the geospatial information, wide heterogeneity of data models and formats, and complexity of processing. Therefore, Big Earth Data Analytics requires specifically tailored techniques and tools. The EarthServer Big Earth Data Analytics engine offers a solution for coverage-type datasets, built around a high performance array database technology, and the adoption and enhancement of standards for service interaction (OGC WCS and WCPS). The EarthServer solution, led by the collection of requirements from scientific communities and international initiatives, provides a holistic approach that ranges from query languages and scalability up to mobile access and visualization. The result is demonstrated and validated through the development of lighthouse applications in the Marine, Geology, Atmospheric, Planetary and Cryospheric science domains.

Statistical Techniques In Stereology

Stereology typically concerns estimation of properties of a geometric structure from plane section information. This paperprovides a brief review of some statistical aspects of this rapidly developing field, with some reference to applications in the earth sciences. After an introductory discussion of the scope of stereology, section 2 briefly mentions results applicable when no assumptions can be made about the stochastic nature of the sampled matrix, statistical considerations then arising solelyfrom the ‘randomness’ of the plane section. The next two sections postulate embedded particles of specific shapes, the particular case of spheres being discussed in some detail. References are made to results for ‘thin slices’ and other prob-ing mechanisms. Randomly located convex particles, of otherwise arbitrary shape, are discussed in section 5 and the review concludes with a specific application of stereological ideas to some data on neolithic mining.

The Impact of Offshore Wind Farms on Marine Ecosystems: A Review Taking an Ecosystem Services Perspective

The rapid increase in renewable energy generation from wind has increased concerns about the impacts that wind arrays have on the marine environment and what these impacts mean for society. One method for identifying the impacts of offshore wind farms (OWFs) on human welfare is through the assessment and valuation of ecosystem services. Using an ecosystem services approach, this paper reviews the impacts of OWFs on the ecosystem services delivered by marine environments. During the construction phase, supporting services such as reduced energy capture and nutrient cycling are changed due to the introduction of hard substrate and the reduction in soft sediment habitat at turbine bases. This may lead to changes in all other ecosystem services, both negative and positive. Quantifying these changes, however, is a challenge partly due to data limitations and a lack of clear understanding of the impacts of OWFs on the marine ecosystems. Scientific effort needs to quantitatively explore the impacts of OWFs on ecosystem functionality and the gathering of data that enables the assessment of changes to ecosystem services. Data needed to better quantify and value the impacts of OWFs on ecosystem services are suggested. The development of methods which integrate socioeconomic valuation of ecosystem services into the evaluation of renewable energy devices compliments efforts in assessing the environmental impacts and should enable a holistic assessment of the impact of renewable energy production and greenhouse gas mitigation technologies on the U. K. carbon footprint.

Coccolithophore surface distributions in the North Atlantic and their modulation of the air-sea flux of CO2 from 10 years of satellite Earth observation data

Coccolithophores are the primary oceanic phytoplankton responsible for the production of calcium carbonate (CaCO3). These climatically important plankton play a key role in the oceanic carbon cycle as a major contributor of carbon to the open ocean carbonate pump (similar to 50 %) and their calcification can affect the atmosphere-to-ocean (air-sea) uptake of carbon dioxide (CO2) through increasing the seawater partial pressure of CO2 (pCO(2)). Here we document variations in the areal extent of surface blooms of the globally important coccolithophore, Emiliania huxleyi, in the North Atlantic over a 10-year period (1998-2007), using Earth observation data from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). We calculate the annual mean sea surface areal coverage of E. huxleyi in the North Atlantic to be 474 000 +/- 104 000 km(2), which results in a net CaCO3 carbon (CaCO3-C) production of 0.14-1.71 Tg CaCO3-C per year. However, this surface coverage (and, thus, net production) can fluctuate inter-annually by -54/+81% about the mean value and is strongly correlated with the El Nino/Southern Oscillation (ENSO) climate oscillation index (r = 0.75, p < 0.02). Our analysis evaluates the spatial extent over which the E. huxleyi blooms in the North Atlantic can increase the pCO(2) and, thus, decrease the localised air-sea flux of atmospheric CO2. In regions where the blooms are prevalent, the average reduction in the monthly air-sea CO2 flux can reach 55%. The maximum reduction of the monthly air-sea CO2 flux in the time series is 155 %. This work suggests that the high variability, frequency and distribution of these calcifying plankton and their impact on pCO(2) should be considered if we are to fully understand the variability of the North Atlantic air-to-sea flux of CO2. We estimate that these blooms can reduce the annual N. Atlantic net sink atmospheric CO2 by between 3-28 %.

Transfer Across the Air-Sea Interface

The efficiency of transfer of gases and particles across the air-sea interface is controlled by several physical, biological and chemical processes in the atmosphere and water which are described here (including waves, large- and small-scale turbulence, bubbles, sea spray, rain and surface films). For a deeper understanding of relevant transport mechanisms, several models have been developed, ranging from conceptual models to numerical models. Most frequently the transfer is described by various functional dependencies of the wind speed, but more detailed descriptions need additional information. The study of gas transfer mechanisms uses a variety of experimental methods ranging from laboratory studies to carbon budgets, mass balance methods, micrometeorological techniques and thermographic techniques. Different methods resolve the transfer at different scales of time and space; this is important to take into account when comparing different results. Air-sea transfer is relevant in a wide range of applications, for example, local and regional fluxes, global models, remote sensing and computations of global inventories. The sensitivity of global models to the description of transfer velocity is limited; it is however likely that the formulations are more important when the resolution increases and other processes in models are improved. For global flux estimates using inventories or remote sensing products the accuracy of the transfer formulation as well as the accuracy of the wind field is crucial.

Developments since 2005 in understanding potential environmental impacts of CO2 leakage from geological storage

This paper reviews research into the potential environmental impacts of leakage from geological storage of CO2 since the publication of the IPCC Special Report on Carbon Dioxide Capture and Storage in 2005. Possible impacts are considered on onshore (including drinking water aquifers) and offshore ecosystems. The review does not consider direct impacts on man or other land animals from elevated atmospheric CO2 levels. Improvements in our understanding of the potential impacts have come directly from CO2 storage research but have also benefitted from studies of ocean acidification and other impacts on aquifers and onshore near surface ecosystems. Research has included observations at natural CO2 sites, laboratory and field experiments and modelling. Studies to date suggest that the impacts from many lower level fault- or well-related leakage scenarios are likely to be limited spatially and temporarily and recovery may be rapid. The effects are often ameliorated by mixing and dispersion of the leakage and by buffering and other reactions; potentially harmful elements have rarely breached drinking water guidelines. Larger releases, with potentially higher impact, would be possible from open wells or major pipeline leaks but these are of lower probability and should be easier and quicker to detect and remediate.

Can benthic community structure be used to predict the process of bioturbation in real ecosystems

Disentangling the roles of environmental change and natural environmental variability on biologically mediated ecosystem processes is paramount to predict future marine ecosystem functioning. Bioturbation, the biogenic mixing of sediments, has a regulating role in marine biogeochemical processes. However, our understanding of bioturbation as a community level process and of its environmental drivers is still limited by loose use of terminology, and a lack of consensus about what bioturbation is. To help resolve these challenges, this empirical study investigated the links between four different attributes of bioturbation (bioturbation depth, activity and distance, and biodiffusive transport); the ability of an index of bioturbation (BPc) to predict each of them; and their relation to seasonality, in a shallow coastal system – the Western Channel Observatory, UK. Bioturbation distance depended on changes in benthic community structure, while the other three attributes were more directly influenced by seasonality in food availability. In parallel, BPc successfully predicted bioturbation distance but not the other attributes of bioturbation. This study therefore highlights that community bioturbation results from this combination of processes responding to environmental variability at different time-scales. However, community level measurements of bioturbation across environmental variability are still scarce, and BPc is calculated using commonly available data on benthic community structure and the functional classification of invertebrates. Therefore, BPc could be used to support the growth of landscape scale bioturbation research, but future uses of the index need to consider which bioturbation attributes the index actually predicts. As BPc predicts bioturbation distance, estimated here using a random-walk model applicable to community settings, studies using either of the metrics should be directly comparable and contribute to a more integrated future for bioturbation research.

Evidence of mass failure in the Hess Deep Rift from multi-resolutional bathymetry data

New regional swath and near-bottom bathymetric data provide constraints on shallow structures at the Hess Deep Rift, an oceanic rift that exposes the crust and upper mantle of fast-spreading oceanic lithosphere created at the East Pacific Rise. These data reveal the presence of a lobate structure with a length of ~ 4 km and a width of ~ 6 km south of an Intrarift Ridge, north of Hess Deep. The lobe consists of a series of concentric benches that are widest in the center of the lobe and narrower at the edges, with a dominant bench separating two distinct morphologic regions in the lobe. There are two end-member possible interpretations of this feature: 1) the lobate structure represents a mass failure with little translation that contains coherent blocks that preserve rift-related lineaments; or 2) it represents degraded tectonic structures, and the lobate form is accounted for by, for example, two intersecting faults. We favor the slump interpretation because it more readily accounts for the lobate form of the feature and the curved benches and based on the presence of other similar lobes in this region. In the slump model, secondary structures within the benches may indicate radial spreading during or after failure. The large lobate structure we identify south of the Intrarift Ridge in Hess Deep is one of the first features of its kind identified in an oceanic rift, and illustrates that mass failure may be a significant process in these settings, consistent with the recognition of their importance in mid-ocean ridges, oceanic islands, and continental rifts. Understanding the structure of the Hess Deep Rift is also important for reconstructing the section of fast-spreading oceanic crust exposed here.

Defining scenarios of future vectors of change in marine life and associated economic sectors

Addressing the multitude of challenges in marine policy requires an integrated approach that considers the multitude of drivers, pressures, and interests, from several disciplinary angles. Scenarios are needed to harmonise the analyses of different components of the marine system, and to deal with the uncertainty and complexity of the societal and biogeophysical dynamics in the system. This study considers a set of socio-economic scenarios to (1) explore possible futures in relation to marine invasive species, outbreak forming species, and gradual changes in species distribution and productivity; and (2) harmonise the projection modelling performed within associated studies. The exercise demonstrates that developing interdisciplinary scenarios as developed in this study is particularly complicated due to (1) the wide variety in endogeneity or exogeneity of variables in the different analyses involved; (2) the dual role of policy decisions as variables in a scenario or decisions to be evaluated and compared to other decisions; and (3) the substantial difference in time scale between societal and physical drivers.

iMarNet: an ocean biogeochemistry model intercomparison project within a common physical ocean modelling framework

Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.

Modelling the Stoichiometric Regulation of C-Rich Toxins in Marine Dinoflagellates

Toxin production in marine microalgae was previously shown to be tightly coupled with cellular stoichiometry. The highest values of cellular toxin are in fact mainly associated with a high carbon to nutrient cellular ratio. In particular, the cellular accumulation of C-rich toxins (i.e., with C:N > 6.6) can be stimulated by both N and P deficiency. Dinoflagellates are the main producers of C-rich toxins and may represent a serious threat for human health and the marine ecosystem. As such, the development of a numerical model able to predict how toxin production is stimulated by nutrient supply/deficiency is of primary utility for both scientific and management purposes. In this work we have developed a mechanistic model describing the stoichiometric regulation of C-rich toxins in marine dinoflagellates. To this purpose, a new formulation describing toxin production and fate was embedded in the European Regional Seas Ecosystem Model (ERSEM), here simplified to describe a monospecific batch culture. Toxin production was assumed to be composed by two distinct additive terms; the first is a constant fraction of algal production and is assumed to take place at any physiological conditions. The second term is assumed to be dependent on algal biomass and to be stimulated by internal nutrient deficiency. By using these assumptions, the model reproduced the concentrations and temporal evolution of toxins observed in cultures of Ostreopsis cf. ovata, a benthic/epiphytic dinoflagellate producing C-rich toxins named ovatoxins. The analysis of simulations and their comparison with experimental data provided a conceptual model linking toxin production and nutritional status in this species. The model was also qualitatively validated by using independent literature data, and the results indicate that our formulation can be also used to simulate toxin dynamics in other dinoflagellates. Our model represents an important step towards the simulation and prediction of marine algal toxicity.

Global Carbon Budget 2015

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates as well as consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2, and land-cover change (some including nitrogen–carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2005–2014), EFF was 9.0 ± 0.5 GtC yr−1, ELUC was 0.9 ± 0.5 GtC yr−1, GATM was 4.4 ± 0.1 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 3.0 ± 0.8 GtC yr−1. For the year 2014 alone, EFF grew to 9.8 ± 0.5 GtC yr−1, 0.6 % above 2013, continuing the growth trend in these emissions, albeit at a slower rate compared to the average growth of 2.2 % yr−1 that took place during 2005–2014. Also, for 2014, ELUC was 1.1 ± 0.5 GtC yr−1, GATM was 3.9 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and SLAND was 4.1 ± 0.9 GtC yr−1. GATM was lower in 2014 compared to the past decade (2005–2014), reflecting a larger SLAND for that year. The global atmospheric CO2 concentration reached 397.15 ± 0.10 ppm averaged over 2014. For 2015, preliminary data indicate that the growth in EFF will be near or slightly below zero, with a projection of −0.6 [range of −1.6 to +0.5] %, based on national emissions projections for China and the USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the global economy for the rest of the world. From this projection of EFF and assumed constant ELUC for 2015, cumulative emissions of CO2 will reach about 555 ± 55 GtC (2035 ± 205 GtCO2) for 1870–2015, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2015).

Impacts of Climate Modes on Air-Sea Heat Exchange in the Red Sea

The impacts of various climate modes on the Red Sea surface heat exchange are investigated using the MERRA reanalysis and the OAFlux satellite reanalysis datasets. Seasonality in the atmospheric forcing is also explored. Mode impacts peak during boreal winter [December–February (DJF)] with average anomalies of 12–18 W m−2 to be found in the northern Red Sea. The North Atlantic Oscillation (NAO), the east Atlantic–west Russia (EAWR) pattern, and the Indian monsoon index (IMI) exhibit the strongest influence on the air–sea heat exchange during the winter. In this season, the largest negative anomalies of about −30 W m−2 are associated with the EAWR pattern over the central part of the Red Sea. In other seasons, mode-related anomalies are considerably lower, especially during spring when the mode impacts are negligible. The mode impacts are strongest over the northern half of the Red Sea during winter and autumn. In summer, the southern half of the basin is strongly influenced by the multivariate ENSO index (MEI). The winter mode–related anomalies are determined mostly by the latent heat flux component, while in summer the shortwave flux is also important. The influence of the modes on the Red Sea is found to be generally weaker than on the neighboring Mediterranean basin.