Search

Karst geology is characterized by the presence of sinkholes and voids, which may pose significant risk to existing infrastructure. Sinkhole formation is often observed near active quarries, where dewatering operations can alter regional groundwater flow patterns leading to subsidence and increased void formation. In these areas, identifying locations which may be susceptible to sinkhole formation requires an ability to map dissolution features within the rock. Traditional geotechnical explorations alone are not well-suited to this effort as they only provide subsurface information at discrete points and therefore may miss voids within the rock. Geophysical methods offer a means to produce continuous profiles of the rock surface and possible locations for voids but interpreting the results of these tests in karstic geology can be challenging. This study uses 2D electrical resistivity and seismic surveys at a site with previous sinkhole activity along an interstate in central Alabama. The site is adjacent to a limestone quarry. Resistivity data is collected using 2D dipole-dipole and strong-gradient arrays. The seismic data is processed using a full waveform inversion (FWI) technique. Subsurface profiles interpreted from the geophysical surveys are then compared to borehole data from previous site investigations. Results from the geophysical surveys are found to be consistent with borehole data regarding variation of bedrock depth and identification of possible sinkhole features. Potential limitations and sources of error pertaining to each survey type are considered.

Fox Glacier/Te Moeka o Tuawe is a fast-responding maritime glacier that has undergone multiple advance and retreat phases during recent decades. Here we use a combination of repeat photography, Str...

Glaciers are unique indicators of climate change. While recent global-scale glacier decline has been attributed to anthropogenic forcing, direct links between human-induced climate warming and extreme glacier mass-loss years have not been documented. Here we apply event attribution methods to document this at the regional scale, targeting the highest mass-loss years (2011 and 2018) across New Zealand’s Southern Alps. Glacier mass balance is simulated using temperature and precipitation from multiple climate model ensembles. We estimate extreme mass loss was at least six times (2011) and ten times (2018) (>90% confidence) more likely to occur with anthropogenic forcing than without. This increased likelihood is driven by present-day temperatures ~1.0 °C above the pre-industrial average, confirming a connection between anthropogenic emissions and high annual ice loss. These results suggest that as warming and extreme heat events continue and intensify, there will be an increasingly visible human fingerprint on extreme glacier mass-loss years in the coming decades. Detecting a human role in a given year of extreme glacier mass loss is difficult at regional scales. Event attribution methods estimate that two extreme mass-loss years in the New Zealand Southern Alps, 2011 and 2018, were at least six and ten times more likely with anthropogenic climate warming.

The land ice contribution to global mean sea level rise has not yet been predicted 1 using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models 2–8, but primarily used previous-generation scenarios 9 and climate models 10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios 11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

Mountain glacier moraine sequences and their chronologies allow us to evaluate the timing and climate conditions that underpin changes in the equilibrium line altitudes (ELAs), which can provide valuable information on the paleoclimatology of understudied regions such as tropical East Africa. However, moraine sequences are inherently discontinuous, and the precise climate conditions that they represent can be ambiguous due to the sensitivity of mountain glaciers to temperature, precipitation, and other environmental variables. Here, we used a two-dimensional (2-D) ice-flow and mass-balance model to simulate glacier extents and ELAs in the Rwenzori Mountains in East Africa over the past 31,000 yr (31 k.y.), including the Last Glacial Maximum (LGM), late glacial period, and the Holocene Epoch. We drove the glacier model with two independent, continuous temperature reconstructions to simulate possible glacier length changes through time. Model input paleoclimate values came from branched glycerol dialkyl glycerol tetraether (brGDGT) temperature reconstructions from alpine lakes on Mount Kenya for the last ~31 k.y., and precipitation reconstructions for the LGM came from various East African locations. We then compared the simulated fluctuations with the positions and ages (where known) of the Rwenzori moraines. The simulated glacier extents reached within 1.1 km of the dated LGM moraines in one valley (93% of the full LGM extent) when forced by the brGDGT temperature reconstructions (maximum cooling of 6.1 °C) and a decrease in precipitation (−10% than modern amounts). These simulations suggest that the Rwenzori glaciers required a cooling of at least 6.1 °C to reach the dated LGM moraines. Based on the model output, we predict an age of 12–11 ka for moraines located halfway between the LGM and modern glacier extents. We also predict ice-free conditions in the Rwenzori Mountains for most of the early to middle Holocene, followed by a late Holocene glacier readvance within the last 2000 yr.

The only complete inventory of New Zealand glaciers was based on aerial photography starting in 1978. While there have been partial updates using 2002 and 2009 satellite data, most glaciers are still represented by the 1978 outlines in contemporary global glacier databases. The objective of this project is to establish an updated glacier inventory for New Zealand. We have used Landsat 8 OLI satellite imagery from February and March 2016 for delineating clean glaciers using a semi-automatic band ratio method and debris-covered glaciers using a maximum likelihood classification. The outlines have been checked against Sentinel-2 MSI data, which have a higher resolution. Manual post processing was necessary due to misclassifications (e.g. lakes, clouds), mapping in shadowed areas, and combining the clean and debris-covered parts into single glaciers. New Zealand glaciers cover an area of 794 ± 34 km2 in 2016 with a debris-covered area of 10%. Of the 2918 glaciers, seven glaciers are >10 km2 while 71% is 2. The debris cover on those largest glaciers is >40%. Only 15 glaciers are located on the North Island. For a selection of glaciers, we were able to calculate the area reduction between the 1978 and 2016 inventories.

Clinical trials provide evidence essential for progress in health care, and as the complexity of medical care has increased, the demand for such data has dramatically expanded. Conducting clinical trials has also become more complicated, evolving to meet increasing challenges in delivering clinical care and meeting regulatory requirements. Despite this, the general approach to data collection remains the same, requiring that researchers submit clinical data in response to study treatment protocols, using precisely defined data structures made available in study-specific case report forms. Currently, research data management is not integrated within the patient’s clinical care record, creating added burden for clinical staff and opportunities for error. During the past decade, the electronic health record has become standard across the US healthcare system and is increasingly used to collect and analyze data reporting quality metrics for clinical care delivery. Recently, electronic health record data have also been used to address clinical research questions; however, this approach has significant drawbacks due to the unstructured and incomplete nature of current electronic health record data. This report describes steps necessary to use the electronic health record as a tool for conducting high-quality clinical research.

Quantifying historic changes in glacier size and mass balance is important for understanding how the cryosphere responds to climate variability and change. Airborne photogrammetry enables glacier extent and equilibrium line altitudes (ELAs) to be monitored for more glaciers at lower cost than traditional mass-balance programs and other remote-sensing techniques. Since 1977, end-of-summer-snowlines, which are a proxy for annual ELAs, have been recorded for 50 glaciers in the Southern Alps of New Zealand using oblique aerial photographs. In this study, we use structure from motion photogrammetry to estimate the camera parameters, including position, for historic photographs, which we then use to measure glacier change. We apply this method to a small maritime New Zealand glacier (Brewster Glacier, 1670–2400 m a.s.l.) to derive annual ELA and length records between 1981 and 2017, and quantify the uncertainties associated with the method. Our length reconstruction shows largely continuous terminus retreat of 365 ± 12 m for Brewster Glacier since 1981. The ELA record, which compares well with glaciological mass-balance data measured between 2005 and 2015, shows pronounced interannual variability. Mean ELAs range from 1707 ± 6 to 2303 ± 5 m a.s.l., with the highest ELAs occurring in the last decade.

Glaciers offer the potential to reconstruct past climate over timescales from decades to millennia. They are found on nearly every continent, and at the Last Glacial Maximum, glaciers were larger in all regions on Earth. The physics of glacier-climate interaction are relatively well understood, and glacier models can be used to reconstruct past climate from geological evidence of past glacier extent. This can lead to significant insights regarding past, present, and future climate. For example, glacier modeling has demonstrated that the near-ubiquitous global pattern of glacier retreat during the last few centuries resulted from a global-scale climate warming of ∼1°C, consistent with instrumental data and climate proxy records. Climate reconstructions from glaciers have also demonstrated that the tropics were colder at the Last Glacial Maximum than was originally inferred from sea surface temperature reconstructions. Future efforts to reconstruct climate from glaciers may provide new constraints on climate sensitivity to CO2 forcing, polar amplification of climate change, and more.

Recent model studies suggest that interannual climatic variability could be confounding the interpretation of glacier fluctuations as climate signals. Paleoclimate interpretations of moraine positions and associated cosmogenic exposure ages may have large uncertainties if the glacier in question was sensitive to interannual variability. Here we address the potential for interannual temperature and precipitation variability to cause large shifts in glacier length during the Holocene. Using a coupled ice-flow and mass-balance model, we simulate the response of Cameron Glacier, a small mountain glacier in New Zealand's Southern Alps, to two types of climate forcing: equilibrium climate and variable climate. Our equilibrium results suggest a net warming trend from the Early Holocene ( 10.69 ± 0.41 ka; 2.7 °C cooler than present) to the Late Holocene (CE 1864; 1.3 °C cooler than present). Interannual climatic variability cannot account for the Holocene glacier fluctuations in this valley. Future studies should consider local environmental characteristics, such as a glacier's climatic setting and topography, to determine the magnitude of glacier length changes caused by interannual variability.

Glaciers experienced worldwide retreat during the twentieth and early twenty first centuries, and the negative trend in global glacier mass balance since the early 1990s is predominantly a response to anthropogenic climate warming. The exceptional terminus advance of some glaciers during recent global warming is thought to relate to locally specific climate conditions, such as increased precipitation. In New Zealand, at least 58 glaciers advanced between 1983 and 2008, and Franz Josef and Fox glaciers advanced nearly continuously during this time. Here we show that the glacier advance phase resulted predominantly from discrete periods of reduced air temperature, rather than increased precipitation. The lower temperatures were associated with anomalous southerly winds and low sea surface temperature in the Tasman Sea region. These conditions result from variability in the structure of the extratropical atmospheric circulation over the South Pacific. While this sequence of climate variability and its effect on New Zealand glaciers is unusual on a global scale, it remains consistent with a climate system that is being modified by humans., Many New Zealand glaciers advanced during recent global warming, bucking a worldwide trend of glacier retreat. Here, the authors show that these glacier advances were forced by a sequence of unusually cool years in the New Zealand region, rather than a period of increased precipitation.

Recognising the scarcity of glacier mass-balance data in the Southern Hemisphere, a mass-balance measurement programme was started at Brewster Glacier in the Southern Alps of New Zealand in 2004. Evolution of the measurement regime over the 11 years of data recorded means there are differences in the spatial density of data obtained. To ensure the temporal integrity of the dataset a new geostatistical approach is developed to calculate mass balance. Spatial co-variance between elevation and snow depth allows a digital elevation model to be used in a co-kriging approach to develop a snow depth index (SDI). By capturing the observed spatial variability in snow depth, the SDI is a more reliable predictor than elevation and is used to adjust each year of measurements consistently despite variability in sampling spatial density. The SDI also resolves the spatial structure of summer balance better than elevation. Co-kriging is used again to spatially interpolate a derived mean summer balance index using SDI as a co-variate, which yields a spatial predictor for summer balance. The average glacier-wide surface winter, summer and annual balances over the period 2005–15 are 2484, −2586 and −102 mm w.e., respectively, with changes in summer balance explaining most of the variability in annual balance.

Geological records of mountain glacier fluctuations provide useful evidence for tracing the magnitude and rate of past temperature change. In this study, we present air temperature reconstructions for the last glacial termination in New Zealand derived using snowline reconstructions and numerical glacier modelling. We target the Arthur's Pass moraines in the Otira River catchment, which have previously been dated to the Lateglacial using cosmogenic 10Be. Recalculation of these exposure ages using a locally calibrated 10Be production rate indicates that these moraines formed ca. 16–14 ka. Our glacier modelling experiments and snowline reconstructions exhibit good agreement and show that the Arthur's Pass moraines formed in a climate that was 2.2–3.5 °C colder than present. Combining our results with other, proximal glacier records shows that ice in this catchment retreated ca. 50 km from the coastal plain to the main divide during the interval 17–15 ka, in response to a temperature increase of at least ca. 3 °C. Over half of this retreat occurred after the glacier had withdrawn from an overdeepened basin. Thus, we conclude that temperature increase was the primary driver of widespread and rapid glacier retreat in New Zealand at the onset of the last glacial termination.

The occurrence of terrestrial palynomorphs in Quaternary marine sedimentary sequences allows for direct land–sea correlations and provides a means for transferring Marine Isotope Stage chronologies to terrestrial records that extend beyond the range of radiocarbon dating. Both of these important applications require an implicit assumption that the lag between pollen release and final deposition on the seafloor – here referred to as source-to-sink residence time – is negligible in relation to the chronological resolution of the sedimentary sequence. Most studies implicitly assume zero lag, and where studies do take palynomorph residence time into account, its magnitude is rarely quantified. In Westland, New Zealand, fluvial transport is the main source of delivery of terrestrial pollen offshore to the adjacent East Tasman Sea. We radiocarbon-dated organic matter carried and deposited by contemporary Westland rivers that drain catchments with varying degrees of disturbance. The ages obtained ranged widely from essentially modern (i.e., − 57 ± 22 cal yr BP) to 3583 ± 188 cal yr BP, suggesting that precisely constraining the residence time in this region is unlikely to be achieved. We also compared the timing of four palynomorph events characterising Westland's late Pleistocene, along with the well-dated Kawakawa/Oruanui Tephra (KOT), between marine core MD06-2991 and four terrestrial records from Westland. Critically, all palynomorph events and the KOT are chronologically indistinguishable with respect to the independently dated marine and terrestrial records, supporting the general principle of transferring the marine chronology onto the terrestrial records in this setting. In other regions, particularly those lacking the high soil production and erosion rates that characterise Westland, we suggest that similar tests of marine residence time should be conducted before assumptions of zero or negligible lag are invoked.

Quantitative palaeoclimate reconstructions provide data for evaluating the mechanisms of past, natural climate variability. Geometries of former mountain glaciers, constrained by moraine mapping, afford the opportunity to reconstruct palaeoclimate, due to the close relationship between ice extent and local climate. In this study, we present results from a series of experiments using a 2-D coupled energy balance–ice flow model that investigate the palaeoclimate significance of Last Glacial Maximum moraines within nine catchments in the central North Island, New Zealand. We find that the former ice limits can be simulated when present-day temperatures are reduced by between 4 and 7 °C, if precipitation remains unchanged from present. The spread in the results between the nine catchments is likely to represent the combination of chronological and model uncertainties. The majority of catchments targeted require temperature decreases of 5.1 to 6.3 °C to simulate the former glaciers, which represents our best estimate of the temperature anomaly in the central North Island, New Zealand, during the Last Glacial Maximum. A decrease in precipitation of up to 25 % from present, as suggested by proxy evidence and climate models, increases the magnitude of the required temperature changes by up to 0.8 °C. Glacier model experiments using reconstructed topographies that exclude the volume of post-glacial (

Daily time lapse imagery and pixel tracking was used to monitor and track spatial and temporal changes in sediment-mantled hillslopes, during and immediately after glacier retreat in the Fox Glacier/Te Moeka o Tuawe valley, New Zealand. Observations from 2014 to 2018 of the Fox Glacier and surrounding hillslopes show hillslope failure is primarily coincident with, and triggered by, glacier retreat with rainfall accelerating movement of the hillslope. During glacier retreat, failure of the hillslope primarily occurred through sediment sliding, internal deformation of the sediment and occasional surficial debris falls, flows and avalanches delivering approximately 9.2 M m3 of sediment directly to the underside of the glacier over the study period with a maximum daily averaged movement of 0.4 m per day of the main sediment mass. Following debuttressing, hillslope failure became dominated by localised rainfall-induced debris flows which initiated gullying of the sediment-mantled slope. Ongoing instability of the slope and associated movement is maintained by toe erosion from the Fox River and melting dead ice, while continued rapid failure is facilitated through localised debris flows. The tracking of temporal and spatial changes of sediment-mantled hillslopes during glacier retreat has shown broad-scale hillslope response to occur quickly within days of rainfall or accelerated glacier retreat, particularly during summer. Debris flows, commonly thought to be a dominant erosion process within paraglacial environments, only occur following complete debuttressing and are supply-limited, only occurring after sediment sliding has occurred. Unlike many other case studies, sediment connectivity immediately following glacier retreat is high due to a lack of storage space and high rainfall inducing mass movements, efficiently delivering hillslope sediments to the proglacial stream channel. Attempts to quantify displaced volumes of sediment from paraglacial systems are likely underestimated due to a) a lack of focus on the early and latter stages of debuttressing b) a reliance on debris flows being a primary transport mechanism, and c) sediment being delivered sub-glacially rather than supra-glacially, further enhancing connectivity.

Polychromatic laser light can reduce speckle noise in wavefront-sensing and imaging applications that use direct-detection schemes. To help quantify the achievable reduction in speckle, this paper investigates the accuracy and numerical efficiency of three separate wave-optics methods. Each method simulates the active illumination of extended objects with polychromatic laser light. In turn, this paper uses the Monte Carlo method, the depth-slicing method, and the spectral-slicing method, respectively, to simulate the laser-object interaction. The limitations and sampling requirements of all three methods are discussed. Further, the numerical efficiencies of the methods are compared over a range of conditions. The Monte Carlo method is found to be the most efficient, while spectral slicing is more efficient than depth slicing for well-resolved objects. Finally, Hu's theory is used to quantify method accuracy when possible (i.e., for well-resolved objects). In general, the theory compares favorably to the simulation methods.

Glacier mass balance is more sensitive to warming than cooling, but feedbacks related to the exposure of previously buried firn and ice in very warm years is not generally considered in sensitivity...

Short-term glacier velocity variations typically occur when a water input is accommodated by an increase in the subglacial water pressure. Although these velocity variations have been well documented on many glaciers, few studies have considered them on glaciers where heavy rain and glacier melt occur year-round. This study investigates the relationship between water inputs and glacier velocity on Franz Josef Glacier, New Zealand. We installed six GNSS stations across the lower glacier during austral summer 2010/11 and one station during summer 2012/13. Glacier velocity remained elevated at all stations for ∼7 days following large rain events. During diurnal melt events, we find velocity variations in the early afternoon (12:00–16:00) at 600 m a.s.l. and in the late evening (20:00–01:00) at 400 m a.s.l. We hypothesize that the late-evening velocity variations occurred as an upstream region of high subglacial water pressures and accelerated ice motion propagated downstream. This mechanism may also explain the increased longitudinal compression and transverse extension across the lower glacier during speed-up events. Our results indicate that the subglacial drainage system likely decreases in efficiency upstream and that the water input variability can still cause short-term velocity variations despite the large year-round water inputs.

The development of high-power fiber lasers is of great interest due to the advantages they offer relative to other laser technologies. Currently, the maximum power from a reportedly single-mode fiber amplifier stands at 10 kW. Though impressive, this power level was achieved at the cost of a large spectral linewidth, making the laser unsuitable for coherent or spectral beam combination techniques required to reach power levels necessary for airborne tactical applications. An effective approach in limiting the SBS effect is to insert an electro-optic phase modulator at the low-power end of a master oscillator power amplifier (MOPA) system. As a result, the optical power is spread among spectral sidebands; thus raising the overall SBS threshold of the amplifier. It is the purpose of this work to present a comprehensive numerical scheme that is based on the extended nonlinear Schrodinger equations that allows for accurate analysis of phase modulated fiber amplifier systems in relation to the group velocity dispersion and Kerr nonlinearities and their effect on the coherent beam combining efficiency. As such, we have simulated a high-power MOPA system modulated via filtered pseudo-random bit sequence format for different clock rates and power levels. We show that at clock rates of ≥30 GHz, the combination of GVD and self-phase modulation may lead to a drastic drop in beam combining efficiency at the multi-kW level. Furthermore, we extend our work to study the effect of cross-phase modulation where an amplifier is seeded with two laser sources.

In this paper we report the generation of flat top optical spectrum using an arbitrary waveform generator to increase the SBS threshold in high power optical fiber amplifiers. The optical spectrum consists of a number of discrete spectral lines, ranging from 16 to 380, within the bandwidth of 2GHz, corresponding to line spacing between 133 MHz and 5 MHz. These discrete spectral lines correspond to a PRBS pattern of n = 4 to n = 8. The SBS threshold and coherence properties of the flat top spectrum are measured and compared to that of the filtered PRBS in a kilowatt class fiber amplifier. It is experimentally demonstrated that for large frequency line spacing, the flat top spectrum significantly outperforms the corresponding filtered PRBS, but as the line spacing is decreased to less than the Brillouin bandwidth, the two modulation waveforms have similar enhancement factors in the SBS threshold due to the enhanced crosstalk between neighboring frequency components.

Hybrid microstructured fibers, utilizing both air holes and high index cladding structures, provide important advantages over conventional fiber including robust fundamental mode operation with large core diameters (>30μm) and spectral filtering (i.e. amplified spontaneous emission and Raman suppression). This work investigates the capabilities of a hybrid fiber designed to suppress stimulated Brillouin scattering (SBS) and modal instability (MI) by characterizing these effects in a counter-pumped amplifier configuration as well as interrogating SBS using a pump-probe Brillouin gain spectrum (BGS) diagnostic suite. The fiber has a 35 μm annularly gain tailored core, the center doped with Yb and the second annulus comprised of un-doped fused silica, designed to optimize gain in the fundamental mode while limiting gain to higher order modes. A narrow-linewidth seed was amplified to an MI-limited 820 W, with near-diffraction-limited beam quality, an effective linewidth ~ 1 GHz, and a pump conversion efficiency of 78%. Via a BGS pump-probe measurement system a high resolution spectra and corresponding gain coefficient were obtained. The primary gain peak, corresponding to the Yb doped region of the core, occurred at 15.9 GHz and had a gain coefficient of 1.92×10-11 m/W. A much weaker BGS response, due to the pure silica annulus, occurred at 16.3 GHz. This result demonstrates the feasibility of power scaling hybrid microstructured fiber amplifiers

Knowledge of the ice thickness distribution of glaciers and ice caps is an important prerequisite for many glaciological and hydrological investigations. A wealth of approaches has recently been presented for inferring ice thickness from characteristics of the surface. With the Ice Thickness Models Intercomparison eXperiment (ITMIX) we performed the first coordinated assessment quantifying individual model performance. A set of 17 different models showed that individual ice thickness estimates can differ considerably – locally by a spread comparable to the observed thickness. Averaging the results of multiple models, however, significantly improved the results: On average over the 21 considered test cases, comparison against direct ice thickness measurements revealed deviations in the order of 10 ± 24% of the mean ice thickness (1-sigma estimate). For models relying on multiple data sets -- such as surface ice velocity fields, surface mass balance, or rates of ice thickness change – the results highlighted the sensitivity to input data consistency. Together with the requirement of being able to handle large regions in an automated fashion, the capacity of better accounting for uncertainties in the input data will be a key for an improved next generation of ice thickness estimation approaches.

The Darwin–Hatherton Glacial system (DHGS) connects the East Antarctic Ice Sheet (EAIS) with the Ross Ice Shelf and is a key area for understanding past variations in ice thickness of surrounding ice masses. Here we present the first detailed measurements of ice thickness and grounding zone characteristics of the DHGS as well as new measurements of ice velocity. The results illustrate the changes that occur in glacier geometry and ice flux as ice flows from the polar plateau and into the Ross Ice Shelf. The ice discharge and the mean basal ice shelf melt for the first 8.5 km downstream of the grounding line amount to 0.24 ± 0.05 km3 a−1 and 0.3 ± 0.1 m a−1, respectively. As the ice begins to float, ice thickness decreases rapidly and basal terraces develop. Constructed maps of glacier geometry suggest that ice drainage from the EAIS into the Darwin Glacier occurs primarily through a deep subglacial canyon. By contrast, ice thins to

A diode pumped (793nm), long-wavelength thulium-doped all-fiber amplifier was demonstrated. Laser gain competition was used to suppress parasitic lasing and an output power of 80 W was generated at 2130nm with 50% optical efficiency.

Ka Roimata o Hine Hukatere (Franz Josef Glacier) is a fast-flowing maritime glacier and its climatological and hydrological drivers are different from those of many previously studied alpine glaciers. The glacier tongue has recently advanced as well as retreated, remains largely snow free, has significant volumes of melt and rainwater inputs throughout the year, and experiences small radiation and air temperature fluctuations over diurnal to seasonal time scales. We discuss measurements of surface velocity made between 2000 and 2012 at annual, seasonal, weekly, and daily time scales together with measurements of glacier geometry change, and calculations of surface water inputs and subglacial water pressure variations derived from a distributed surface mass balance model and a one-dimensional conduit hydrology model, respectively. Annual velocity variations are linked to changes in glacier geometry and advance/retreat cycles with accelerations during thickening and advance and decelerations during...

Polychromatic laser light can reduce speckle contrast in wavefront-sensing and imaging applications that use direct detection schemes. To help quantify the associated reduction in speckle contrast, this study investigates the accuracy and numerical efficiency of three separate wave-optics models that simulate the active illumination of extended objects with polychromatic laser light. The three separate models use spectral slicing, Monte Carlo averaging, and depth slicing, respectively, to simulate the laser-target interaction. The sampling requirements of all three models are discussed. Comparisons to analytical solutions and experimental data are made when possible. In general, the experiments and theory compare favorably with the models.

We report power scaling results of a highly efficient narrow-linewidth monolithic Yb-doped fiber amplifier seeded with two signals, operating at 1038 and 1064 nm. With the appropriate seed power ratio applied, this technique was shown to suppress stimulated Brillouin scattering in conjunction with phase modulation, while generating the output power in predominantly the longer wavelength signal. Notably, the integration of laser gain competition with pseudo-random bit sequence phase modulation, set at a clock rate of 2.5 GHz and utilizing an optimized pattern to match the shortened effective nonlinear length, yielded 1 kW of output power. The beam quality was measured to be near the diffraction limit with no sign of transverse mode instability. Furthermore, the coherent beam combination performance of the amplifier provided a 90% combining efficiency with no indication of spectral broadening when compared to the single-tone case. Overall, the power scaling results represent a significant reduction in spectral linewidth compared to that of commercially available narrow-linewidth Yb-doped fiber amplifiers.

We report efficient coherent beam combining of five kilowatt-class fiber amplifiers with a diffractive optical element (DOE). Based on a master oscillator power amplifier (MOPA) configuration, the amplifiers were seeded with pseudo random phase modulated light. Each non-polarization maintaining fiber amplifier was optically path length matched and provides approximately 1.2 kW of near diffraction-limited output power (measured M2

We report results from two ~1.5 kW Yb-doped fiber amplifiers with comparable optical to optical efficiencies and linewidths. One amplifier utilized a fiber with a core diameter of 25 μm while the core diameter of the fiber utilized in the other amplifier was 20 μm. Stimulated Brillouin scattering (SBS) suppression in both cases was achieved through pseudo-random bit sequence (PRBS) phase modulation. While the power generated in the larger core fiber was modal instability (MI) limited, no sign of MI was observed in the smaller core fiber. This may allow us to utilize the higher MI threshold fiber to scale further while maintaining sufficiently narrow linewidth for beam combining. Furthermore, in a demonstration of the utility of applying a thermal gradient in conjunction with phase modulation to suppress SBS further, we report on a 1 kW amplifier driven at a PRBS clock rate of 2 GHz. Finally, we compare the coherent beam combining properties of amplifiers seeded with PRBS phase modulated sources to those seeded with white noise sources.

Laser gain competition was used in conjunction with external phase modulation techniques in order to investigate power scaling of narrow-linewidth monolithic Ytterbium-doped fiber amplifiers. In this study, both pseudo-random bit sequence (PRBS) and filtered white noise source (WNS) modulation techniques were separately utilized to drive the external phase modulator for linewidth broadening. The final-stage amplifier was then seeded with the phase modulated narrow-linewidth 1064 nm signal along with a spectrally broader 1038 nm source. Consequently, integration of laser gain competition in conjunction with PRBS phase modulation yields a factor of ∼15 dB in stimulated Brillouin scattering (SBS) threshold enhancement at a clock rate of 2.5 GHz; leading to 1 kilowatt of output power with 85% optical efficiency at 1064 nm. Notably, the combination of PRBS phase modulation with laser gain competition provided superior enhancement in SBS threshold power when compared to the WNS modulated case. The beam quality at maximum power was near the diffraction limit (M2

A low-pass filter is used to suppress the high-frequency components of a PRBS signal and is used to modulate a fiber amplifier. Notably, the coherence length and SBS threshold of a kilowatt-class amplifier are increased.

We present multi-kW power scaling results of narrow line Yb-doped fiber amplifiers. Nonlinear suppression is attained through pseudo-random modulation. Subsequently, 5 pseudorandom modulated kW amplifiers were coherently combined into a 5kW beam.

Understanding relationships between snow accumulation and synoptic climatology is important for assessing the way in which future climate variability will impact on glacier mass balance. However, few studies have as yet examined these relationships. Variability in snow accumulation on mid-latitude glaciers is strongly influenced by atmospheric circulation, orography, and redistribution of snow by wind. Very little is known about these processes in the New Zealand Southern Alps, where it is assumed that west-facing glaciers receive higher snow totals. However, few measurements are available to test this hypothesis. These processes were investigated over a 21-day period in winter 2008 on glaciers located west (Franz Josef Glacier) and east (Tasman Glacier) of the Main Divide of the Southern Alps. We directly measured snow accumulation and considered how it was affected by synoptic weather regime and location with respect to the Main Divide. Both glaciers received ∼75% of their snowfall during troug...

Coverage of ice velocities in the central part of the Southern Alps, New Zealand, is obtained from feature tracking using repeat optical imagery in 2002 and 2006. Precise orthorectification, co-registration and correlation is carried out using the freely available software COSI-Corr. This analysis, combined with short times between image acquisitions, has enabled velocities to be captured even in the accumulation areas, where velocities are lowest and surface features ephemeral. The results indicate large velocities for mountain glaciers (i.e. up to ∼5 m d−1) as well as dynamic changes in some glaciers that have occurred between 2002 and 2006. For the steep and more responsive Fox and Franz Josef Glaciers the speed increased at the glacier snout during the advance period, while the low-angled and debris-covered Tasman Glacier showed no measurable velocity change. Velocity increases on the steeper glaciers are the result of an observed thickening and steepening of the glacier tongues as they moved from a retreat phase in 2002 to an advance phase in 2006. This contrasting behaviour is consistent with historic terminus position changes. The steeper glaciers have undergone several advance/retreat cycles during the observation period (1894 to present), while the low-angled glacier showed little terminus response until retreat resulting from the accelerating growth of a proglacial lake commenced in 1983.

The sensitivity of glaciers to climatic change is key information in assessing the response and sea-level implications of projected future warming. New Zealand glaciers are important globally as an example of how maritime glaciers will contribute to sea-level rise. A spatially distributed energy-balance model is applied to Brewster Glacier, New Zealand, in order to calculate glacier mass balance, run-off and sensitivity to climate change. The model successfully simulates four annual mass-balance cycles. Close to half (52%) of the energy available for melt on the glacier is supplied by turbulent heat fluxes, with radiation less important, except during the winter. Model sensitivity to temperature change is one of the largest reported on Earth, at −2.0 m w.e. a−1 °C−1. In contrast, a 50% change in precipitation is required to offset the mass-balance change resulting from a 1 °C temperature change. Meltwater runoff sensitivity is also very high, increasing 60% with a 1°C warming. The extreme sensitivity of mass balance to temperature change suggests that significant ice loss will occur with even moderate climate warming.

Observations show that glaciers around the world are in retreat and losing mass. Internationally coordinated for over a century, glacier monitoring activities provide an unprecedented dataset of glacier observations from ground, air and space. Glacier studies generally select specific parts of these datasets to obtain optimal assessments of the mass-balance data relating to the impact that glaciers exercise on global sea-level fluctuations or on regional runoff. In this study we provide an overview and analysis of the main observational datasets compiled by the World Glacier Monitoring Service (WGMS). The dataset on glacier front variations (∼42 000 since 1600) delivers clear evidence that centennial glacier retreat is a global phenomenon. Intermittent readvance periods at regional and decadal scale are normally restricted to a subsample of glaciers and have not come close to achieving the maximum positions of the Little Ice Age (or Holocene). Glaciological and geodetic observations (∼5200 since 1850) show that the rates of early 21st-century mass loss are without precedent on a global scale, at least for the time period observed and probably also for recorded history, as indicated also in reconstructions from written and illustrated documents. This strong imbalance implies that glaciers in many regions will very likely suffer further ice loss, even if climate remains stable.

Despite their relatively small total ice volume, mid-latitude valley glaciers are expected to make a significant contribution to global sea-level rise over the next century due to the sensitivity of their mass-balance systems to small changes in climate. Here we use a degree-day model to reconstruct the past century of mass-balance variation at ‘Ka Roimata o Hine Hukatere’ Franz Josef Glacier, New Zealand, and to predict how mass balance may change over the next century. Analysis of the relationship between temperature, precipitation and mass balance indicates that temperature is a stronger control than precipitation on the mass balance of Franz Josef Glacier. The glacier’s mass balance, relative to its 1986 geometry, has decreased at a mean annual rate of 0.02m a–1w.e. between 1894 and 2005. We compare this reduction to observations of terminus advance and retreat, of which Franz Josef Glacier has the best record in the Southern Hemisphere. For the years 2000–05 the relative mass balance ranged from –0.75 to +1.50m a–1w.e., with 2000/01 the only year showing a negative mass balance. In a regionally downscaled Intergovernmental Panel on Climate Change mean warming scenario, the annual relative mass balance will continue to decrease at 0.02m a–1w.e. through the next century.

Pseudo random phase modulation signals have been shown to provide considerable stimulated Brillouin scattering (SBS) suppression in narrow linewidth Yb-doped all-fiber amplifiers. In terms of coherent beam combining, however, pseudo random signals display a linear drop in visibility; leading to pronounced drops in combining efficiencies for small path length deviations. To that end, we report a novel filtered pseudo random modulation approach for enhanced combining efficiency and coherence length performance. Here a low pass radio frequency (RF) filter is used to mitigate the PRBS high frequency components, thereby suppressing the sidelobes in the optical spectrum. This leads to an approximate Gaussian visibility function and improved coherence lengths of up to 27% in a kW class fiber amplifier (954 W). In addition, the spectral sidelobe suppression leads to concurrent SBS threshold enhancement due to a reduction in the spectral overlap between the Rayleigh reflected light and the Stokes shifted light. This reduction in the SBS seeding phenomena leads to ~10% SBS threshold improvements in a kW class fiber amplifier. Theoretical and experimental data is presented to substantiate the improved coherence length and SBS suppression. More importantly, the simultaneous nonlinear SBS suppression and coherence length benefits of the filtered PRBS approach can have a significant impact for high power, narrow linewidth, all-fiber amplifiers.

Ice-penetrating radar and modelling data are presented suggesting the presence of a zone of temperate ice, water ponding or saturated sediment beneath the tongue of Taylor Glacier, Dry Valleys, Antarctica. The proposed subglacial zone lies 3–6 km up-glacier of the terminus and is 400– 1000m across. The zone coincides with an extensive topographic overdeepening to 80m below sea level. High values of residual bed reflective power across this zone compared to other regions and the margins of the glacier require a high dielectric contrast between the ice and the bed and are strongly indicative of the presence of basal water or saturated sediment. Analysis of the hydraulic equipotential surface also indicates strong convergence into this zone of subglacial water flow paths. However, thermodynamic modelling reveals that basal temperatures in this region could not exceed –7˚C relative to the pressure-melting point. Such a result is at odds with the radar observations unless the subglacial water is a hypersaline brine.

Local wind speed variations influence the energy and mass fluxes over snow through snow accumulation, sublimation of drifting and blowing snow, or variations in turbulent fluxes over static snow and ice surfaces. We use idealized model experiments to analyze the sensitivity of turbulent fluxes over static snow surfaces to variations in wind speed under different climatic conditions. We find that the sensitivity (change in the turbulent flux per change of unit wind speed) increases with increasing air temperature and relative humidity. The sensitivity of turbulent fluxes to wind speed is highest when the stability parameter ζ = 1, which occurs at wind speeds typical for glacierized catchments (3-5 m s -1), and exponentially decreases either side of that range. That peak in sensitivity is caused by atmospheric stability corrections in the model, and occurs independently of the flux-profile relationships we tested. Our results quantify the significant effect of local wind speed variations on turbulent fluxes over snow and ice and can be used to estimate potential model uncertainties in different climates, especially for the typical assumption in distributed hydrological models that the wind speed is spatially constant. © 2012 Elsevier Ltd. All rights reserved.