Broadband imaging via direct inversion of blended dispersed source array data

Although seismic sources typically consist of identical broadband units alone, no physical constraint dictates the use of only one kind of device. We propose an acquisition method that involves the simultaneous exploitation of multiple types of sources during seismic surveys. It is suggested to replace (or support) traditional broadband sources with several devices individually transmitting diverse and reduced frequency bands and covering together the entire temporal and spatial bandwidth of interest. Together, these devices represent a so‐called dispersed source array. As a consequence, the use of simpler sources becomes a practical proposition for seismic acquisition. In fact, the devices dedicated to the generation of the higher frequencies may be smaller and less powerful than the conventional sources, providing the acquisition system with increased operational flexibility and decreasing its environmental impact. Offshore, we can think of more manageable boats carrying air guns of different volumes or marine vibrators generating sweeps with different frequency ranges. On land, vibrator trucks of different sizes, specifically designed for the emission of particular frequency bands, are preferred. From a manufacturing point of view, such source units guarantee a more efficient acoustic energy transmission than today's complex broadband alternatives, relaxing the low‐ versus high‐frequency compromise. Furthermore, specific attention can be addressed to choose shot densities that are optimum for different devices according to their emitted bandwidth. In fact, since the sampling requirements depend on the maximum transmitted frequencies, the appropriate number of sources dedicated to the lower frequencies is relatively small, provided the signal‐to‐noise ratio requirements are met. Additionally, the method allows to rethink the way to address the ghost problem in marine seismic acquisition, permitting to tow different sources at different depths based on the devices' individual central frequencies. As a consequence, the destructive interference of the ghost notches, including the one at 0 Hz, is largely mitigated. Furthermore, blended acquisition (also known as simultaneous source acquisition) is part of the dispersed source array concept, improving the operational flexibility, cost efficiency, and signal‐to‐noise ratio. Based on theoretical considerations and numerical data examples, the advantages of this approach and its feasibility are demonstrated.     

Blended-acquisition design of irregular geometries towards faster,cheaper, safer and better seismic surveying

The application of blended acquisition has drawn considerable attention owing to its ability to improve the operational efficiency as well as the data quality and health, safety and environment performance. Furthermore, the acquisition of less data contributes to the business aspect, while the desired data density is still realizable via subsequent data reconstruction. The use of fewer detectors and sources also minimizes operational risks in the field. Therefore, a combined implementation of these technologies potentially enhances the value of a seismic survey further. One way to encourage this is to minimize any imperfection in deblending and data reconstruction during processing. In addition, one may derive survey parameters that enable a further improvement in these processes as introduced in this study. The proposed survey design workflow iteratively performs the following steps to derive the survey parameters responsible for source blending as well as the spatial sampling of detectors and sources. The first step is the application of blending and sampling operators to unblended and well-sampled data. We then apply closed-loop deblending and data reconstruction. The residue for a given design from this step is evaluated and subsequently used by genetic algorithms to simultaneously update the survey parameters related to both blending and spatial sampling. The updated parameters are fed into the next iteration until they satisfy the given termination criteria. We also propose a repeated encoding sequence to form a parameter sequence in genetic algorithms, making the size of problem space manageable. The results of the proposed workflow are outlined using blended dispersed source array data incorporating different scenarios that represent acquisition in marine, transition zone and land environments. Clear differences attributed solely to the parameter design are easily recognizable. Additionally, a comparison among different optimization schemes illustrates the ability of genetic algorithms along with a repeated encoding sequence to find better solutions within a computationally affordable time. The optimized parameters yield a notable enhancement in the deblending and data reconstruction quality and consequently provide optimal acquisition scenarios.     

A GCM study of future climate response to aerosol pollution reductions

We use the global atmospheric GCM aerosol model ECHAM5-HAM to asses possible impacts of future air pollution mitigation strategies on climate. Air quality control strategies focus on the reduction of aerosol emissions. Here we investigate the extreme case of a maximum feasible end-of-pipe abatement of aerosols in the near term future (2030) in combination with increasing greenhouse gas (GHG) concentrations. The temperature response of increasing GHG concentrations and reduced aerosol emissions leads to a global annual mean equilibrium temperature response of 2.18 K. When aerosols are maximally abated only in the Industry and Powerplant sector, while other sectors stay with currently enforced regulations, the temperature response is 1.89 K. A maximum feasible abatement applied in the Domestic and Transport sector, while other sectors remain with the current legislation, leads to a temperature response of 1.39 K. Increasing GHG concentrations alone lead to a temperature response of 1.20 K. We also simulate 2–5% increases in global mean precipitation among all scenarios considered, and the hydrological sensitivity is found to be significantly higher for aerosols than for GHGs. Our study, thus highlights the huge potential impact of future air pollution mitigation strategies on climate and supports the need for urgent GHG emission reductions. GHG and aerosol forcings are not independent as both affect and are influenced by changes in the hydrological cycle. However, within the given range of changes in aerosol emissions and GHG concentrations considered in this study, the climate response towards increasing GHG concentrations and decreasing aerosols emissions is additive.     

Atmospheric bridge on orbital time scales

Large-scale atmospheric patterns are examined on orbital timescales using a climate model which explicitly resolves the atmosphere–ocean–sea ice dynamics. It is shown that, in contrast to boreal summer where the climate mainly follows the local radiative forcing, the boreal winter climate is strongly determined by modulation of circulation modes linked to the Arctic Oscillation/North Atlantic Oscillation (AO/NAO) and the Northern/Southern Annular Modes. We find that during a positive phase of the AO/NAO the convection in the tropical Pacific is below normal. The related atmospheric circulation provides an atmospheric bridge for the precessional forcing inducing a non-uniform temperature anomalies with large amplitudes over the continents. We argue that this is important for mechanisms responsible for multi-millennial climate variability and glacial inception.     

Evaluation of a Finite-Element Sea-Ice Ocean Model (FESOM) set-up to study the interannual to decadal variability in the deep-water formation rates

The characteristics of a global set-up of the Finite-Element Sea-Ice Ocean Model under forcing of the period 1958–2004 are presented. The model set-up is designed to study the variability in the deep-water mass formation areas and was therefore regionally better resolved in the deep-water formation areas in the Labrador Sea, Greenland Sea, Weddell Sea and Ross Sea. The sea-ice model reproduces realistic sea-ice distributions and variabilities in the sea-ice extent of both hemispheres as well as sea-ice transport that compares well with observational data. Based on a comparison between model and ocean weather ship data in the North Atlantic, we observe that the vertical structure is well captured in areas with a high resolution. In our model set-up, we are able to simulate decadal ocean variability including several salinity anomaly events and corresponding fingerprint in the vertical hydrography. The ocean state of the model set-up features pronounced variability in the Atlantic Meridional Overturning Circulation as well as the associated mixed layer depth pattern in the North Atlantic deep-water formation areas.     

Estimating groundwater recharge for an arid karst system using a combined approach of time‐lapse camera monitoring and water balance modelling

Groundwater is the principal water resource in semi‐arid and arid environments. Therefore, quantitative estimates of its replenishment rate are important for managing groundwater systems. In dry regions, karst outcrops often show enhanced recharge rates compared with other surface and sub‐surface conditions. Areas with exposed karst features like sinkholes or open shafts allow point recharge, even from single rainfall events. Using the example of the As Sulb plateau in Saudi Arabia, this study introduces a cost‐effective and robust method for recharge monitoring and modelling in karst outcrops. The measurement of discharge of a representative small catchment (4.0 · 10⁴ m²) into a sinkhole, and hence the direct recharge into the aquifer, was carried out with a time‐lapse camera. During the monitoring period of two rainy seasons (autumn 2012 to spring 2014), four recharge events were recorded. Afterwards, recharge data as well as proxy data about the drying of the sediment cover are used to set up a conceptual water balance model. The model was run for 17 years (1971 to 1986 and 2012 to 2014). Simulation results show highly variable seasonal recharge–precipitation ratios between 0 and 0.27. In addition to the amount of seasonal precipitation, this ratio is influenced by the interannual distribution of rainfall events. Overall, an average annual groundwater recharge for the doline (sinkhole) catchment on As Sulb plateau of 5.1 mm has estimated for the simulation period. Copyright © 2015 John Wiley & Sons, Ltd.     

Carbon sequestration and environmental effects of afforestation with Pinus radiata D. Don in the Western Cape, South Africa

A three-step methodology to assess the carbon sequestration and the environmental impact of afforestation projects in the framework of the Flexible Mechanisms of the Kyoto Protocol (Joint Implementation and Clean Development Mechanism) was developed and tested using a dataset collected from the Jonkershoek forest plantation, Western Cape, South Africa, which was established with Pinus radiata in former native fynbos vegetation and indigenous forest. The impact of a change in land use was evaluated for a multifunctional, a production and a non-conversion scenario. First, the carbon balance was modelled with GORCAM and was expressed as (1) C sequestration in tC ha⁻¹ year⁻¹ in soil, litter, and living biomass according to the rules of the first commitment period of the Kyoto Protocol, and (2) CO₂ emission reductions in tC ha⁻¹ year⁻¹, which includes carbon sequestered in the above-mentioned pools and additionally in wood products, as well as emission reductions due to fossil fuel substitution. To estimate forest growth, three data sources were used: (1) inventory data, (2) growth simulation with a process-based model, and (3) yield tables. Second, the effects of land use change were assessed for different project scenarios using a method related to Life Cycle Assessment (LCA). The method uses 17 quantitative indicators to describe the impact of project activities on water, soil, vegetation cover and biodiversity. Indicator scores were calculated by comparing indicator values with reference values, estimated for the climax vegetation. The climax vegetation is the site-specific ecosystem phase with the highest exergy content and the highest exergy flow dissipation capacity. Third, the land use impact per functional unit of 1 tC sequestered was calculated by combining the results of step 1 and step 2. The average baselines to obtain carbon additionality are 476 tC ha⁻¹ for indigenous forest and 32 tC ha⁻¹ for fynbos. Results show that the influence of the growth assessment method on the magnitude of C sequestration and hence on the environmental impact per functional unit is large. When growth rate is assessed with the mechanistic model and with the yield table, it is overestimated in the early years and underestimated in the long term. The main conclusion of the scenario analysis is that the production forest scenario causes higher impacts per functional unit than the multifunctional scenario, but with the latter being less efficient in avoiding CO₂ emissions. The proposed method to assess impacts on diverse components of the ecosystem is able to estimate the general tendency of the adverse and positive effects of each scenario. However, some indicators, more specifically about biodiversity and water balance, could be improved or reinterpreted in light of specific local data about threat to biodiversity and water status.     

Interpretation of single-well tracer tests using fractional-flow dimensions. Part 2: A case study

Generalized equations using fractional-flow dimensions were derived to estimate the Darcy and groundwater-flow velocities obtained from the point-dilution and the single-well injection-withdrawal field tests. Flow velocities can only be estimated from single-well tests if the kinematic porosity or the hydraulic conductivity and hydraulic gradient are known a priori. A pumping test performed on the boreholes will yield an estimate of the fractional-flow dimension and the extent of the flow region by applying the generalized radial flow (GRF) model of Barker [Barker JA (1988) A generalized radial flow model for hydraulic tests in fractured rock. Water Resour Res 24(10):1796–1804]. These parameters are used in the generalised equations for the single-well tracer tests to estimate the flow-through area and, therewith, the Darcy or flow velocity. The generalized procedure, described in detail in Part 1 of this paper, is applied to two boreholes on the Campus Test Site located at the University of the Free State, Bloemfontein, South Africa, and it is shown that the two independent tests (i.e. the point-dilution and the single-well injection-withdrawal tests) yield similar estimates of the natural-seepage velocity in the aquifer. The estimated natural-flow velocity obtained by using fractional dimensions is about two times higher than the velocity estimated by using the standard method (i.e. flow dimension n=2, flow thickness equal to length of the sealed-off section of the borehole). Electronic Publication     

Temporal and spatial differences in salinity and water chemistry in SW Florida estuaries: Effects of human-impacted watersheds

Three estuaries near Naples, Florida with variably modified watersheds have been investigated to understand the chemical consequences of altering drainage patterns. Blackwater River (near natural drainage, control site), Henderson Creek (moderately modified watershed), and Faka-Union Canal (severe channelization) were sampled for temperature, salinity, δ¹⁸O, δ¹³C of dissolved inorganic carbon (DIC), molality of CO₂ (ΣCO₂), and Mg:Ca and Sr:Ca ratios between freshwater and marine water end members over a 17-mo period. Carbon isotope composition followed similar seasonal patterns as salinity. Freshwater and seawater end members are more negative than the global average, likely reflecting equilibration with local carbon sources derived from mangrove leaf litter and groundwater. δ¹³C responds to differences in primary productivity between estuaries. Henderson Creek has higher primary productivity than Blackwater River (probable due to higher sewage input and agricultural runoff) and has more positive δ¹³C and lower ΣCO₂. δ¹⁸O is affected by seasonal input of freshwater from atmospheric precipitation, evaporation, and groundwater. Late summer and fall rains lower the δ¹⁸O of estuarine water, whereas evaporative conditions in the dry season elevate δ¹⁸O to values that can be more positive upstream than those from the Gulf of Mexico (estuarine inversion). Evaporation produces water in the Gulf of Mexico that is >1‰ more positive than the global sea surface average most of the year. The very negative δ¹⁸O values in Blackwater River and Henderson Creek likely reflect atmospheric and groundwater contribution. Mg:Ca and Sr:Ca ratios of Gulf water from all three estuaries are similar to global averages at low latitudes. Freshwater end members among estuaries are different in that Blackwater River has higher ratios, suggesting a groundwater contribution. Dolomitic rocks in the subsurface likely provide a source of Mg ions.