Studies of grand minima in sunspot cycles by using a flux transport solar dynamo model

We propose that grand minima in solar activity are caused by simultaneous fluctuations in the meridional circulation and the Babcock-Leighton mechanism for the poloidal field generation in the flux transport dynamo model. We present the following results: (a) fluctuations in the meridional circulation are more effective in producing grand minima; (b) both sudden and gradual initiations of grand minima are possible; (c) distributions of durations and waiting times between grand minima seem to be exponential; (d) the coherence time of the meridional circulation has an effect on the number and the average duration of grand minima, with a coherence time of about 30 yr being consistent with observational data. We also study the occurrence of grand maxima and find that the distributions of durations and waiting times between grand maxima are also exponential, like the grand minima. Finally we address the question of whether the Babcock-Leighton mechanism can be operative during grand minima when there are no sunspots. We show that an alpha-effect restricted to the upper portions of the convection zone can pull the dynamo out of the grand minima and can match various observational requirements if the amplitude of this alpha-effect is suitably fine-tuned.

Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Seasonal Deuterium Excess in a Tien Shan Ice Core: Influence of Moisture Transport and Recycling in Central Asia

Stable water isotope (delta(18)O, deltaD) data from a high elevation (5100 masl) ice core recovered from the Tien Shan Mountains, Kyrgyzstan, display a seasonal cycle in deuterium excess (d = deltaD - 8* delta(18)O) related to changes in the regional hydrologic cycle during 1994 - 2000. While there is a strong correlation (r(2) = 0.98) between delta(18)O and dD in the ice core samples, the regression slope (6.9) and mean d value (23.0) are significantly different than the global meteoric water line values. The resulting time-series ice core d profile contains distinct winter maxima and summer minima, with a yearly d amplitude of similar to 15 - 20parts per thousand. Local-scale processes that may affect d values preserved in the ice core are not consistent with the observed seasonal variability. Data from Central Asian monitoring sites in the Global Network of Isotopes in Precipitation (GNIP) have similar seasonal d changes. We suggest that regional-scale hydrological conditions, including seasonal changes in moisture source, transport, and recycling in the Caspian/ Aral Sea region, are responsible for the observed spatial and temporal d variability.

Spatial Patterns in Seasonal and Interannual Variability of Chlorophyll and Sea Surface Temperature in the California Current

Physical forcing and biological response within the California Current System (CCS) are highly variable over a wide range of scales. Satellite remote sensing offers the only feasible means of quantifying this variability over the full extent of the CCS. Using six years (1997-2003) of daily SST and chlorophyll imagery, we map the spatial dependence of dominant temporal variability at resolutions sufficient to identify recurrent mesoscale circulation and local pattern associated with coastal topography. Here we describe mean seasonal cycles and interannual variation; intraseasonal variability is left to a companion paper ( K. R. Legaard and A. C. Thomas, manuscript in preparation, 2006). Coastal upwelling dictates seasonality along north-central California, where weak cycles of SST fluctuate between spring minima and late summer maxima and chlorophyll peaks in early summer. Off northern California, chlorophyll maxima are bounded offshore by the seasonally recurrent upwelling jet. Seasonal cycles differ across higher latitudes and in the midlatitude Southern California Bight, where upwelling winds are less vigorous and/or persistent. Seasonality along south-central Baja is strongly affected by processes other than upwelling, despite year-round upwelling-favorable winds. Interannual variation is generally dominated by El Nino and La Nina conditions. Interannual SST variance is greatest along south-central Baja, although interannual variability constitutes a greater fraction of total variance inshore along southern Oregon and much of California. Patterns of interannual chlorophyll variance are consistent with dominant forcing through the widespread depression and elevation of the nutricline during El Nino and La Nina, respectively. Interannual variability constitutes a greater fraction of total chlorophyll variance offshore.

Isotopic signature of specific biomarkers derived from anaerobic methanotrophic archaea (ANME groups) and sulphate-reducing bacteria (SRB) at different cold seep provinces

Anaerobic methane-oxidizing microbial communities in sediments at cold methane seeps are important factors in controlling methane emission to the ocean and atmosphere. Here, we investigated the distribution and carbon isotopic signature of specific biomarkers derived from anaerobic methanotrophic archaea (ANME groups) and sulphate-reducing bacteria (SRB) responsible for the anaerobic oxidation of methane (AOM) at different cold seep provinces of Hydrate Ridge, Cascadia margin. The special focus was on their relation to in situ cell abundances and methane turnover. In general, maxima in biomarker abundances and minima in carbon isotope signatures correlated with maxima in AOM and sulphate reduction as well as with consortium biomass. We found ANME-2a/DSS aggregates associated with high abundances of sn-2,3-di-O-isoprenoidal glycerol ethers (archaeol, sn-2-hydroxyarchaeol) and specific bacterial fatty acids (C16:1omega5c, cyC17:0omega5,6) as well as with high methane fluxes (Beggiatoa site). The low to medium flux site (Calyptogena field) was dominated by ANME-2c/DSS aggregates and contained less of both compound classes but more of AOM-related glycerol dialkyl glycerol tetraethers (GDGTs). ANME-1 archaea dominated deeper sediment horizons at the Calyptogena field where sn-1,2-di-O-alkyl glycerol ethers (DAGEs), archaeol, methyl-branched fatty acids (ai-C15:0, i-C16:0, ai-C17:0), and diagnostic GDGTs were prevailing. AOM-specific bacterial and archaeal biomarkers in these sediment strata generally revealed very similar d13C-values of around -100 per mill. In ANME-2-dominated sediment sections, archaeal biomarkers were even more 13C-depleted (down to -120 per mill), whereas bacterial biomarkers were found to be likewise 13C-depleted as in ANME-1-dominated sediment layers (d13C: -100 per mill). The zero flux site (Acharax field), containing only a few numbers of ANME-2/DSS aggregates, however, provided no specific biomarker pattern. Deeper sediment sections (below 20 cm sediment depth) from Beggiatoa covered areas which included solid layers of methane gas hydrates contained ANME-2/DSS typical biomarkers showing subsurface peaks combined with negative shifts in carbon isotopic compositions. The maxima were detected just above the hydrate layers, indicating that methane stored in the hydrates may be available for the microbial community. The observed variations in biomarker abundances and 13C-depletions are indicative of multiple environmental and physiological factors selecting for different AOM consortia (ANME-2a/DSS, ANME-2c/DSS, ANME-1) along horizontal and vertical gradients of cold seep settings.

Sediment, isotopic and geochemical record during Marine Isotope Stages 29 to 34 of IODP Site 339-U1387

Centennial-to-millennial scale records from IODP Site U1387, drilled during IODP Expedition 339 into the Faro Drift at 558 m water depth, now allow evaluating the climatic history of the upper core of the Mediterranean Outflow (MOW) and of the surface waters in the northern Gulf of Cadiz during the early Pleistocene. This study focuses on the period from Marine Isotope Stage (MIS) 29 to 34, i.e. the interval surrounding extreme interglacial MIS 31. Conditions in the upper MOW reflect obliquity, precession and millennial-scale variations. The benthic d18O signal follows obliquity with the exception of an additional, smaller d18O peak that marks the MIS 32/31 transition. Insolation maxima (precession minima) led to poor ventilation and a sluggish upper MOW core, whereas insolation minima were associated with enhanced ventilation and often also increased bottom current velocity. Millennial-scale periods of colder sea-surface temperatures (SST) were associated with short-term maxima in flow velocity and better ventilation, reminiscent of conditions known from MIS 3. A prominent contourite layer, coinciding with insolation cycle 100, was formed during MIS 31 and represents one of the few contourites developing within an interglacial period. MIS 31 surface water conditions were characterized by an extended period (1065-1091 ka) of warm SST, but SST were not much warmer than during MIS 33. Interglacial to glacial transitions experienced 2 to 3 stadial/interstadial cycles, just like their mid-to-late Pleistocene counterparts. Glacial MIS 30 and 32 recorded periods of extremely cold (< 12°C) SST that in their climatic impact were comparable to the Heinrich events of the mid and late Pleistocene. Glacial MIS 34, on the other hand, was a relative warm glacial period off southern Portugal. Overall, surface water and MOW conditions at Site U1387 show strong congruence with Mediterranean climate, whereas millennial-scale variations are closely linked to North Atlantic circulation changes.

Middle to late Miocene oxygen isotope stratigraphy of ODP site 1085 (SE Atlantic): new constrains on Miocene climate variability and sea-level fluctuations

The middle Miocene delta18O increase represents a fundamental change in earth's climate system due to a major expansion and permanent establishment of the East Antarctic Ice Sheet accompanied by some effect of deepwater cooling. The long-term cooling trend in the middle to late Miocene was superimposed by several punctuated periods of glaciations (Mi-Events) characterized by oxygen isotopic shifts that have been related to the waxing and waning of the Antarctic ice-sheet and bottom water cooling. Here, we present a high-resolution benthic stable oxygen isotope record from ODP Site 1085 located at the southwestern African continental margin that provides a detailed chronology for the middle to late Miocene (13.9-7.3 Ma) climate transition in the eastern South Atlantic. A composite Fe intensity record obtained by XRF core scanning ODP Sites 1085 and 1087 was used to construct an astronomically calibrated chronology based on orbital tuning. The oxygen isotope data exhibit four distinct delta18O excursions, which have astronomical ages of 13.8, 13.2, 11.7, and 10.4 Ma and correspond to the Mi3, Mi4, Mi5, and Mi6 events. A global climate record was extracted from the oxygen isotopic composition. Both long- and short-term variabilities in the climate record are discussed in terms of sea-level and deep-water temperature changes. The oxygen isotope data support a causal link between sequence boundaries traced from the shelf and glacioeustatic changes due to ice-sheet growth. Spectral analysis of the benthic delta18O record shows strong power in the 400-kyr and 100-kyr bands documenting a paleoceanographic response to eccentricity-modulated variations in precession. A spectral peak around 180-kyr might be related to the asymmetry of the obliquity cycle indicating that the response of the dominantly unipolar Antarctic ice-sheet to obliquityinduced variations probably controlled the middle to late Miocene climate system. Maxima in the delta18O record, interpreted as glacial periods, correspond to minima in 100-kyr eccentricity cycle and minima in the 174-kyr obliquity modulation. Strong middle to late Miocene glacial events are associated with 400-kyr eccentricity minima and obliquity modulation minima. Thus, fluctuations in the amplitude of obliquity and eccentricity seem to be the driving force for the middle to late Miocene climate variability.

Benthic and planktic foraminiferal stable isotope record and coccolith abundance spanning the last 500,000 years of ODP Hole 121-758C

The Indian Summer Monsoon (ISM) is an inter-hemispheric and highly variable ocean-atmosphere-land interaction that directly affects the densely populated Indian subcontinent. Here, we present new records of palaeoceanographic variability that span the last 500,000 years from the eastern equatorial Indian Ocean, a relatively under-sampled area of ISM influence. We have generated carbon and oxygen stable isotope records from three foraminiferal species from Ocean Drilling Program Site 758 (5°N, 90°E) to investigate the oceanographic history of this region. We interpret our resultant Dd18O (surface-thermocline) record of upper water-column stratification in the context of past ISM variability, and compare orbital phase relationships in our Site 758 data to other climate and monsoon proxies in the region. Results suggest that upper water-column stratification at Site 758, which is dominated by variance at precession and half-precession frequencies (23, 19 and 11 ka), is forced by both local (5°N) insolation and ISM winds. In the precession (23 ka) band, stratification minima at Site 758 lag northern hemisphere summer insolation maxima (precession minima) by 9 ka, which is consistent with Arabian Sea ISM phase estimates and suggests a common wind forcing in both regions. This phase implicates a strong sensitivity to both ice volume and southern hemisphere insolation forcing via latent heat export from the southern subtropical Indian Ocean. Additionally, we find evidence of possible overprinting of millennial-scale events during glacial terminations in our stratification record, which suggests an influence of remote abrupt climate events on ISM dynamics.

Miocene stable isotope record of IODP Site 320-U1335 and revised stable isotope record of IODP Site 321-U1337

The Miocene Climatic Optimum (MCO; ~16.9 to 14.7 Ma) provides an outstanding opportunity to investigate climate-carbon cycle dynamics during a geologically recent interval of global warmth. We present benthic stable oxygen (d18O) and carbon (d13C) isotope records (5-12 kyr time resolution) spanning the late early to middle Miocene interval (18 to 13 Ma) at Integrated Ocean Drilling Program (IODP) Site U1335 (eastern equatorial Pacific Ocean). The U1335 stable isotope series track the onset and development of the MCO as well as the transitional climatic phase culminating with global cooling and expansion of the East Antarctic ice-sheet at ~13.8 Ma. We integrate these new data with published stable isotope, geomagnetic polarity and X-ray fluorescence (XRF) scanner-derived carbonate records from IODP Sites U1335, U1336, U1337 and U1338 on a consistent, astronomically-tuned timescale. Benthic isotope and XRF scanner-derived CaCO3 records depict prominent 100 kyr variability with 400 kyr cyclicity additionally imprinted on d13C and CaCO3 records, pointing to a tight coupling between the marine carbon cycle and climate variations. Our inter-site comparison further indicates that the lysocline behaved in highly dynamic manner throughout the MCO, with >75% carbonate loss occurring at paleo-depths ranging from ~3.4 to ~4 km in the eastern equatorial Pacific Ocean. Carbonate dissolution maxima coincide with warm phases (d18O minima) and d13C decreases, implying that climate-carbon cycle feedbacks fundamentally differed from the late Pleistocene glacial-interglacial pattern, where dissolution maxima correspond to d13C maxima and d18O minima. Carbonate dissolution cycles during the MCO were, thus, more similar to Paleogene hyperthermal patterns.

Integral Inequalities and Applications to Partial Differential Equations with “Maxima”

Кремена В. Стефанова - В тази статия са разрешени някои нелинейни интегрални неравенства, които включват максимума на неизвестната функция на две променливи. Разгледаните неравенства представляват обобщения на класическото неравенство на Гронуол-Белман. Значението на тези интегрални неравенства се определя от широките им приложения в качествените изследвания на частните диференциални уравнения с “максимуми” и е илюстрирано чрез някои директни приложения.

Spatial and seasonal variability in elemental content, δ13C, and δ15N ofThalassia testudinum from South Florida and its implications for ecosystem studies

Elemental and isotopic composition of leaves of the seagrassThalassia testudinum was highly variable across the 10,000 km2 and 8 years of this study. The data reported herein expand the reported range in carbon:nitrogen (C:N) and carbon:phosphorus (C:P) ratios and δ13C and δ15N values reported for this species worldwide; 13.2–38.6 for C:N and 411–2,041 for C:P. The 981 determinations in this study generated a range of −13.5‰ to −5.2‰ for δ13C and −4.3‰ to 9.4‰ for δ15N. The elemental and isotope ratios displayed marked seasonality, and the seasonal patterns could be described with a simple sine wave model. C:N, C:P, δ13C, and δ15N values all had maxima in the summer and minima in the winter. Spatial patterns in the summer maxima of these quantities suggest there are large differences in the relative availability of N and P across the study area and that there are differences in the processing and the isotopic composition of C and N. This work calls into question the interpretation of studies about nutrient cycling and food webs in estuaries based on few samples collected at one time, since we document natural variability greater than the signal often used to imply changes in the structure or function of ecosystems. The data and patterns presented in this paper make it clear that there is no threshold δ15N value for marine plants that can be used as an unambiguous indicator of human sewage pollution without a thorough understanding of local temporal and spatial variability.

Estimating present day extreme water level exceedance probabilities around the coastline of Australia : tides, extra-tropical storm surges and mean sea level

The occurrence of extreme water levels along low-lying, highly populated and/or developed coastlines can lead to considerable loss of life and billions of dollars of damage to coastal infrastructure. Therefore it is vitally important that the exceedance probabilities of extreme water levels are accurately evaluated to inform risk-based flood management, engineering and future land-use planning. This ensures the risk of catastrophic structural failures due to under-design or expensive wastes due to over-design are minimised. This paper estimates for the first time present day extreme water level exceedence probabilities around the whole coastline of Australia. A high-resolution depth averaged hydrodynamic model has been configured for the Australian continental shelf region and has been forced with tidal levels from a global tidal model and meteorological fields from a global reanalysis to generate a 61-year hindcast of water levels. Output from this model has been successfully validated against measurements from 30 tide gauge sites. At each numeric coastal grid point, extreme value distributions have been fitted to the derived time series of annual maxima and the several largest water levels each year to estimate exceedence probabilities. This provides a reliable estimate of water level probabilities around southern Australia; a region mainly impacted by extra-tropical cyclones. However, as the meteorological forcing used only weakly includes the effects of tropical cyclones, extreme water level probabilities are underestimated around the western, northern and north-eastern Australian coastline. In a companion paper we build on the work presented here and more accurately include tropical cyclone-induced surges in the estimation of extreme water level. The multi-decadal hindcast generated here has been used primarily to estimate extreme water level exceedance probabilities but could be used more widely in the future for a variety of other research and practical applications.

The use of penumbrae to define and analyse non-standard fields

Aim A new method of penumbral analysis is implemented which allows an unambiguous determination of field size and penumbra size and quality for small fields and other non-standard fields. Both source occlusion and lateral electronic disequilibrium will affect the size and shape of cross-axis profile penumbrae; each is examined in detail. Method A new method of penumbral analysis is implemented where the square of the derivative of the cross-axis profile is plotted. The resultant graph displays two peaks in the place of the two penumbrae. This allows a strong visualisation of the quality of a field penumbra, as well as a mathematically consistent method of determining field size (distance between the two peak’s maxima), and penumbra (full-widthtenth-maximum of peak). Cross-axis profiles were simulated in a water phantom at a depth of 5 cm using Monte Carlo modelling, for field sizes between 5 and 30 mm. The field size and penumbra size of each field was calculated using the method above, as well as traditional definitions set out in IEC976. The effect of source occlusion and lateral electronic disequilibrium on the penumbrae was isolated by repeating the simulations removing electron transport and using an electron spot size of 0 mm, respectively. Results All field sizes calculated using the traditional and proposed methods agreed within 0.2 mm. The penumbra size measured using the proposed method was systematically 1.8 mm larger than the traditional method at all field sizes. The size of the source had a larger effect on the size of the penumbra than did lateral electronic disequilibrium, particularly at very small field sizes. Conclusion Traditional methods of calculating field size and penumbra are proved to be mathematically adequate for small fields. However, the field size definition proposed in this study would be more robust amongst other nonstandard fields, such as flattening filter free. Source occlusion plays a bigger role than lateral electronic disequilibrium in small field penumbra size.

Running shoes increase Achilles tendon load in walking : an acoustic propagation study

Background: Footwear remains a prime candidate for the prevention and rehabilitation of Achilles tendinopathy as it is thought to decrease tension in the tendon through elevation of the heel. However, evidence for this effect is equivocal. Purpose: This study used an acoustic transmission technique to investigate the effect of running shoes on Achilles tendon loading during barefoot and shod walking. Methods: Acoustic velocity was measured in the Achilles tendon of twelve recreationally–active males (age, 31±9 years; height, 1.78±0.06 m; weight, 81.0±16.9 kg) during barefoot and shod walking at matched self–selected speed (3.4±0.7 km/h). Standard running shoes incorporating a 10– mm heel offset were used. Vertical ground reaction force and spatiotemporal parameters were determined with an instrumented treadmill. Axial acoustic velocity in the Achilles tendon was measured using a custom built ultrasonic device. All data were acquired at a rate of 100 Hz during 10s of steady–state walking. Statistical comparisons between barefoot and shod conditions were made using paired t–tests and repeated measure ANOVAs. Results: Acoustic velocity in the Achilles tendon was highly reproducible and was typified by two maxima (P1, P2) and minima (M1, M2) during walking. Footwear resulted in a significant increase in step length, stance duration and peak vertical ground reaction force compared to barefoot walking. Peak acoustic velocity in the Achilles tendon (P1, P2) was significantly higher with running shoes. Conclusions: Peak acoustic velocity in the Achilles tendon was higher with footwear, suggesting that standard running shoes with a 10–mm heel offset increase tensile load in the Achilles tendon. Although further research is required, these findings question the therapeutic role of standard running shoes in Achilles tendinopathy.

The effect of an in–shoe orthotic heel lift on loading of the Achilles tendon during shod walking

- Study Design Controlled laboratory study - Objective To investigate the effect of a 12–mm in–shoe orthotic heel lift on Achilles tendon loading during shod walking using transmission–mode ultrasonography. - Background Orthotic heel lifts are thought to lower tension in the Achilles tendon but evidence for this effect is equivocal. - Methods The propagation speed of ultrasound, which is governed by the elastic modulus and density of tendon and is proportional to the tensile load to which it is exposed, was measured in the right Achilles tendon of twelve recreationally–active males during shod treadmill walking at matched speeds (3.4±0.7 km/h), with and without addition of a heel lift. Vertical ground reaction force and spatiotemporal gait parameters were simultaneously recorded. Data were acquired at 100Hz during 10s of steady–state walking. Statistical comparisons were made using paired t–tests (α=.05). - Results Ultrasound transmission speed in the Achilles tendon was characterized by two maxima (P1, P2) and minima (M1, M2) during walking. Addition of a heel lift to footwear resulted in a 2% increase and 2% decrease in the first vertical ground reaction force peak and the local minimum, respectively (P<.05). Peak ultrasonic velocity in the Achilles tendon (P1, P2, M2) was significantly lower with addition of an orthotic heel lift (P<.05). - Conclusions Peak ultrasound transmission speed in the Achilles tendon was lower with the addition of a 12–mm orthotic heel lift, indicating the heel lift reduced tensile load in the Achilles tendon, thereby counteracting the effect of footwear. These findings support the addition of orthotic heel lifts to footwear in the rehabilitation of Achilles tendon disorders where management aims to lower tension within the tendon. - Level of Evidence Therapy, level 2a