80 years after aerial photography revealed thousands of aligned oval depressions on the USA’s Atlantic Coastal Plain, the geomorphology of the “Carolina bays” remains enigmatic. Geologists and astronomers alike hold that invoking a cosmic impact for their genesis is indefensible. Rather, the bays are commonly attributed to gradualistic fluvial, marine and/or aeolian processes operating during the Pleistocene era. The major axis orientations of Carolina bays are noted for varying statistically by latitude, suggesting that, should there be any merit to a cosmic hypothesis, a highly accurate triangulation network and suborbital analysis would yield a locus and allow for identification of a putative impact site. Digital elevation maps using LiDAR technology offer the precision necessary to measure their exquisitely-carved circumferential rims and orientations reliably. To support a comprehensive geospatial survey of Carolina bay landforms (Survey) we generated about a million km2 of false-color hsv-shaded bare-earth topographic maps as KML-JPEG tile sets for visualization on virtual globes. Considering the evidence contained in the Survey, we maintain that interdisciplinary research into a possible cosmic origin should be encouraged. Consensus opinion does hold a cosmic impact accountable for an enigmatic Pleistocene event - the Australasian tektite strewn field - despite the failure of a 60-year search to locate the causal astroblem. Ironically, a cosmic link to the Carolina bays is considered soundly falsified by the identical lack of a causal impact structure. Our conjecture suggests both these events are coeval with a cosmic impact into the Great Lakes area during the Mid-Pleistocene Transition, at 786 ka ± 5 k. All Survey data and imagery produced for the Survey are available on the Internet to support independent research. A table of metrics for 50,000 bays examined for the Survey is available from an on-line Google Fusion Table: https://goo.gl/XTHKC4 . Each bay is also geospatially referenceable through a map containing clickable placemarks that provide information windows displaying that bay’s measurements as well as further links which allows visualization of the associated LiDAR imagery and the bay’s planform measurement overlay within the Google Earth virtual globe: https://goo.gl/EHR4Lf .

Earth and environmental models are relied upon to investigate system responses that cannot otherwise be examined. In simulating physical processes, models have adjustable parameters which may, or may not, have a physical meaning. Determining the values to assign to these model parameters is an enduring challenge for earth and environmental modellers. Selecting different error metrics by which the models results are compared to observations will lead to different sets of calibrated model parameters, and thus different model results. Furthermore, models may exhibit ‘equifinal’ behaviour, where multiple combinations of model parameters lead to equally acceptable model performance against observations. These decisions in model calibration introduce uncertainty that must be considered when model results are used to inform environmental decision-making. This presentation focusses on the uncertainties that derive from the calibration of a four parameter lumped catchment hydrological model (GR4J). The GR models contain an inbuilt automatic calibration algorithm that can satisfactorily calibrate against four error metrics in only a few seconds. However, a single, deterministic model result does not provide information on parameter uncertainty. Furthermore, a modeller interested in extreme events, such as droughts, may wish to calibrate against more low flows specific error metrics. In a comprehensive assessment, the GR4J model has been run with 500,000 Latin Hypercube Sampled parameter sets across 303 catchments in the United Kingdom. These parameter sets have been assessed against six error metrics, including two drought specific metrics. This presentation compares the two approaches, and demonstrates that the inbuilt automatic calibration can outperform the Latin Hypercube experiment approach in single metric assessed performance. However, it is also shown that there are many merits of the more comprehensive assessment, which allows for probabilistic model results, multi-objective optimisation, and better tailoring to calibrate the model for specific applications such as drought event characterisation. Modellers and decision-makers may be constrained in their choice of calibration method, so it is important that they recognise the strengths and limitations of their chosen approach.

Currently available soil volumetric water content (VWC) sensors have several drawbacks that pose certain challenges for implementation on large scale for farms. Such issues include cost, scalability, maintenance, wires running through fields, and single-spot resolution. The development of a passive soil moisture sensing system utilizing Radio Frequency Identification (RFID) would allay many of these issues. The type of passive RFID tags discussed in this paper currently cost between 8 to 15 cents retail per tag when purchased in bulk. An incredibly cheap, scalable, low-maintenance, wireless, high-resolution system for sensing soil moisture would be possible if such tags were introduced into the agricultural world. This paper discusses both the use cases as well as examines one implementation of the tags. In 2015, RFID tag manufacturer SmarTrac started selling RFID moisture sensing tags for use in the automotive industry to detect leaks during quality assurance. We place those tags in soil at a depth of 4 inches and compared the moisture levels sensed by the RFID tags with the relative permittivity (εr) of the soil as measured by an industry-standard probe. Using an equation derived by Topp et al, we converted to VWC. We tested this over a wide range of moisture conditions and found a statistically significant, correlational relationship between the sensor values from the RFID tags and the probe’s measurement of εr. We also identified a possible function for mapping vales from the RFID tag to the probe bounded by a reasonable margin of error.

The outer heliosphere is an interesting region characterized by the interaction between the solar wind and the interstellar neutral atoms. Having accomplished the mission to Pluto in 2015 and currently on the way to the Kuiper Belt, the New Horizons spacecraft is following the footsteps of the two Voyager spacecraft that previously explored this region lying roughly beyond 30 AU from the Sun. We model the three-dimensional, time-dependent solar wind plasma flow to the outer heliosphere using our own software Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS), which, in addition to the thermal solar wind plasma, takes into account charge exchange of the solar wind protons with interstellar neutral atoms and treats nonthermal ions (i.e., pickup ions) born during this process as a separate fluid. Additionally, MS-FLUKSS allows us to model turbulence generated by pickup ions. We use MS-FLUKSS to investigate the evolution of plasma and turbulent fluctuations along the trajectory of the New Horizons spacecraft using plasma and turbulence parameters from OMNI data as time-dependent boundary conditions at 1 AU for the Reynolds-averaged MHD equations. We compare the model with in situ plasma observations by New Horizons, Voyager 2, and Ulysses. We also compare the model pickup proton parameters with those derived from the Ulysses-SWICS data.

Examination of the radiation budget at the surface of the Earth shows that there are three factors affecting the surface temperature; the amount of solar radiation absorbed by the atmosphere and by the surface respectively, and the amount of leakage of infrared radiation emitted from the surface directly into space. If there were no leakage, the upwelling infrared radiation from the Earth’s surface would be equal to the incoming solar radiation absorbed by the atmosphere plus twice the solar radiation absorbed by the surface. This results from the summation of a sequence of equal upward and downward re-emissions of infrared radiation absorbed by the atmosphere following the initial absorption of solar radiation. At current levels of solar absorption, this would result in total upwelling radiation of approximately 398.6 W/m2, or a maximum surface temperature of 16.4°C. Allowing for leakage of infrared radiation through the atmospheric window, the resulting emission from the Earth’s surface is reduced to around 396 W/m2, corresponding to the current average global surface temperature of around 15.9°C. Absorption of solar and infrared radiation by greenhouse gases is determined by the absorption bands for the respective gases and their concentrations. Absorption of incoming solar radiation is largely by water vapor and ozone, and an increase in absorption would reduce not increase the surface temperature. Moreover, it is probable that all emitted infrared radiation that can be absorbed by greenhouse gases, primarily water vapor, with a small contribution from carbon dioxide and ozone, is already fully absorbed, and the leakage of around 5.5 % corresponds to the part of the infrared red spectrum that is not absorbed by greenhouse gases. The carbon dioxide absorption bands, which represent a very small percentage of the infrared spectrum, are most likely fully saturated. In these circumstances, increased concentrations of greenhouse gases, and carbon dioxide in particular, will have no effect on the emitted radiation. The surface temperature is probably at the thermodynamic limit for the current luminosity of the sun. Satellite based measurements since 1979 suggest that any global warming over the past 150 years may be due to an increase in total solar irradiance, which we are still a decade or two from being able to confirm.

The net surface water flux (evaporation minus precipitation minus runoff, E-P-R) of the Atlantic Ocean is approximately 0.4 – 0.6 Sv (1 Sv = 10^9 kg s-1) larger than that of the Pacific Ocean, as shown in atmospheric and oceanic reanalyses and by oceanographic estimates. This asymmetry is linked to the asymmetry in sea surface salinity and the existence of the Atlantic Meridional Overturning Circulation. It is shown that the reason for the asymmetry in E-P-R is greater precipitation per unit area over the Pacific south of 30N, while evaporation rates are similar over both basins. It is further argued that the Pacific Ocean is anomalous compared to the Atlantic and Indian Oceans in terms of atmospheric moisture flux convergence and precipitation across the tropics and subtropics. To clarify the mechanism by which water vapour is exported out of the Atlantic basin and imported into the Pacific, we use an air mass trajectory model driven by ERA-Interim reanalysis. Using 12-hourly releases of 14-day back trajectories on the boundaries of ocean drainage basins over the period 2010-2014, we are able to partition the atmospheric moisture fluxes between basins according to their origins (i.e. last contact with the boundary layer). We show that at most a quarter of the E-P-R asymmetry is explained by higher moisture export to the Arctic and Southern basins from the Atlantic than from the Pacific. The main contributions come from differences in the longitudinal atmospheric transport of moisture between the Atlantic, Indian and Pacific basins. In particular, during the Asian summer monsoon the recurvature of the low level flow in the Somali Jet results in a much weaker westward moisture transport from the Indian into the Atlantic basin than across Central America (where it is similar to the zonal average) while there is stronger eastward transport from the Indian to Pacific basins. The net effect is stronger moisture convergence into the Pacific, but weaker into the Atlantic. In contrast to previous thinking, the role of the moisture flux across Central America in the asymmetry, albeit significant, is not dominant.

By solving Laplace’s tidal equations along the equatorial Pacific thermocline, assuming a delayed-differential effective gravity forcing due to a combined lunar+solar (lunisolar) stimulus, we are able to precisely match ENSO periodic variations over wide intervals. The underlying pattern is difficult to decode by conventional means such as spectral analysis, which is why it has remained hidden for so long, despite the excellent agreement in the time-domain. What occurs is that a non-linear seasonal modulation with monthly and fortnightly lunar impulses along with a biennially-aligned “see-saw” is enough to cause a physical aliasing and thus multiple folding in the frequency spectrum. So, instead of a conventional spectral tidal decomposition, we opted for a time-domain cross-validating approach to calibrate the amplitude and phasing of the lunisolar cycles. As the lunar forcing consists of three fundamental periods (draconic, anomalistic, synodic), we used the measured Earth’s length-of-day (LOD) decomposed and resolved at a monthly time-scale [1] to align the amplitude and phase precisely. Even slight variations from the known values of the long-period tides will degrade the fit, so a high-resolution calibration is possible. Moreover, a narrow training segment from 1880-1920 using NINO34/SOI data is adequate to extrapolate the cycles of the past 100 years (see attached figure). To further understand the biennial impact of a yearly differential-delay, we were able to also decompose using difference equations the historical sea-level-height readings at Sydney harbor to clearly expose the ENSO behavior. Finally, the ENSO lunisolar model was validated by back-extrapolating to Unified ENSO coral proxy (UEP) records dating to 1650. The quasi-biennial oscillation (QBO) behavior of equatorial stratospheric winds derives following a similar pattern to ENSO via the tidal equations, but with an emphasis on draconic forcing. This improvement in ENSO and QBO understanding has implications for vastly simplifying global climate models due to the straightforward application of a well-known and well-calibrated forcing. [1] Na, Sung-Ho, et al. “Characteristics of Perturbations in Recent Length of Day and Polar Motion.” Journal of Astronomy and Space Sciences 30 (2013): 33-41.

The structure(s), distribution and dynamics of CDOM have been investigated over the last several decades largely through optical spectroscopy (including both absorption and fluorescence) due to the fairly inexpensive instrumentation and the easy-to-gather data (over thousands published papers from 1990-2016). Yet, the chemical structure(s) of the light absorbing and emitting species or constituents within CDOM has only recently being proposed and tested through chemical manipulation of selected functional groups (such as carbonyl and carboxylic/phenolic containing molecules) naturally occurring within the organic matter pool. Similarly, fitting models (among which the PArallel FACtor analysis, PARAFAC) have been developed to better understand the nature of a subset of DOM, the CDOM fluorescent matter (FDOM). Fluorescence spectroscopy coupled with chemical tests and PARAFAC analyses could potentially provide valuable insights on CDOM sources and chemical nature of the FDOM pool. However, despite that applications (and publications) of PARAFAC model to FDOM have grown exponentially since its first application/publication (2003), a large fraction of such publications has misinterpreted the chemical meaning of the delivered PARAFAC ‘components’ leading to more confusion than clarification on the nature, distribution and dynamics of the FDOM pool. In this context, we employed chemical manipulation of selected functional groups to gain further insights on the chemical structure of the FDOM and we tested to what extent the PARAFAC ‘components’ represent true fluorophores through a controlled chemical approach with the ultimate goal to provide insights on the chemical nature of such ‘components’ (as well as on the chemical nature of the FDOM) along with the advantages and limitations of the PARAFAC application.

In case of nuclear accidents like Fukushima disaster, the influence of 137Cs depositional forms (soluble and/or solid forms) on mineral soil of forest environment on its availability have not reported yet. Soluble (137Cs tagged ultra-pure water) and solid (137Cs contaminated litter-OL and fragmented litter-OF) input forms were mixed with the mineral soils collected under Fukushima coniferous and broadleaf forests. The mixtures then incubated under controlled laboratory condition to evaluate the extractability of 137Cs in soil over time in the presence of decomposition process through two extracting reagents- water and ammonium acetate. Results show that extracted 137Cs fraction with water was less than 1% for soluble input form and below detection limit for solid input form. On the same way with acetate reagent, the extracted 137Cs fraction ranged from 46 to 56% for soluble input and 2 to 15% for solid input, implying the nature of 137Cs contamination strongly influences the extractability and hence the mobility of 137Cs in soil. Although the degradation rate of the organic materials has been calculated in the range of 0.18 ± 0.1 to 0.24 ± 0.1 y-1, its impact on 137Cs extractability appeared very weak at least within the observation period, probably due to shorter observation period. Concerning the treatments of solid 137Cs input forms through acetate extraction, relatively more 137Cs has been extracted from broadleaf organic materials mixes (BL-OL & BL-OF) than the coniferous counterparts. This probably is due to the fact that the lignified coniferous organic materials (CED-OL & CED-OF) components tend to retain more 137Cs than that of the broadleaf. Generally, by extrapolating these observations in to a field context, one can expect more available 137Cs fraction in forest soil from wet depositional pathways such as throughfall and stemflow than those attached with organic materials like litter (OL) and its eco-processed forms (OF).

The kinetic theory of twisted wave instability is developed for the nonthermal dusty plasmas in the presence of helical electric field carrying orbital angular momentum (OAM). The Laguerre Gaussian (LG) mode function is employed to decompose the mode into its longitudinal and azimuthal components. The longitudinal component illustrates the spatial while the azimuthal component describes the phasor variation of the wave. The dielectric functions of quasi-electrostatic dust ion acoustic (DIA) and dust acoustic (DA) twisted modes are obtained and solved analytically to calculate the growth rates. The numerical results are also shown pictorially.

Lossless compression can reduce climate data storage by 30-40%. In general, further reductions require lossy compression that also reduces precision. Fortunately, geoscientific models and measurements generate false precision (scientifically meaningless data bits) that can be eliminated without sacrificing scientifically meaningful data. We introduce Bit Grooming, a lossy compression algorithm that removes the bloat due to false- precision, those bits and bytes beyond the meaningful precision of the data. We evaluated Bit Grooming against competitors Linear Packing, Layer Packing, and GRIB2/JPEG2000.

Studies of intraseasonal and annual cycles of meteorological variables, using projections of climate change, are nowadays extremely important to improve regional socio-economic planning for countries. This is particularly true in Costa Rica, as Central America has been identified as a climate change hot spot. Today many of the economic activities in the region, especially those related to agriculture, tourism and hydroelectric power generation are linked to the seasonal cycle of precipitation. Changes in rainfall (mm/day) and in the diurnal temperature range ($^{\circ}$C) for the periods 1970-1999 and 2070-2099 were investigated using the NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) constructed using the CMIP5 (Coupled Model Intercomparison Project version 5) data. Differences between the multi-model ensembles of the two prospective scenarios (RCP 4.5 and RCP 8.5) and the retrospective baseline scenario were computed. This study highlights Costa Rica as an inflexion point of the climate change in the region and also suggests an early onset of the rainy season and future drying conditions.

The broadband ocean bottom seismometer (BBOBS) and its new generation system (BBOBS-NX) have been developed in Japan, and we performed several test and practical observations to create and establish a new category of the ocean floor broadband seismology, since 1999. Now, the data obtained by our BBOBS and BBOBS-NX is proved to be adequate for broadband seismic analyses. Especially, the BBOBS- NX can obtain the horizontal data comparable to land sites in longer periods (10 s –). Moreover, the BBOBST-NX is in practical evaluation for the mobile tilt observation that enables dense geodetic monitoring. The BBOBS-NX system is a powerful tool, although, it has intrinsic limitation of the ROV operation. If this system can be used without the ROV, like as the BBOBS, it should lead us a true breakthrough of ocean bottom seismology. Hereafter, the new autonomous BBOBS-NX is noted as NX-2G in short. The main problem to realize the NX-2G is a tilt of the sensor unit on landing, which exceed the acceptable limit (±8°) in about 50%. As we had no evidence at which moment and how this tilt occurred, we tried to observe it during the BBOBST-NX landing in 2015 by attaching a video camera and an acceleration logger. The result shows that the tilt on landing was determined by the final posture of the system at the penetration into the sediment, and the large oscillating tilt more than ±10° was observed in descending. The function of the NX-2G system is based on 3 stage operations as shown in the image. The glass float is aimed not only to obtain enough buoyancy to extract the sensor unit, but also to suppress the oscillating tilt of the system in descending. In Oct. 2016, we made the first in-situ test of the NX-2G system with a ROV. It was dropped from the sea surface with the video camera and the acceleration logger. The ROV was used to watch the operation of the system at the seafloor. The landing looked well and it was examined from the acceleration data. As the maximum tilt in descending was about ±2.5°, the glass float effectively suppressed the oscillating tilt. The extraction of the sensor unit was also succeeded with the total buoyancy of about 75 kgf within about 2.5 minutes. As the final step experiment, the one-year-long observation of this NX-2G system has been started in this April with the BBOBS, to obtain simultaneous data for the noise level evaluation.

The Ganga plains represent the abode of more than 400 million people and a region of severe anthropogenic disturbance to natural processes. Changing agricultural practices, inefficient use of water, contamination of groundwater systems, and decrease in soil fertility are some of the issues that have affected the long-term resilience of hydrological processes. The quantification of these processes demands a network of hydro-meteorological instrumentation, low-cost sensors, continuous engagement of stakeholders and real time data transmission at a fine interval. We have therefore set up a Critical Zone Observatory (CZO) in a small watershed (21 square Km ) that forms an intensively managed rural landscape consisting of 92% of agricultural land in the Pandu River Basin (a small tributary of the Ganga River). Apart from setting up a hydrometeorological observatory, the major science questions we want to address relate to development of water balance model, understanding the soil-water interaction and estimation of nutrient fluxes in the watershed. This observatory currently has various types of sensors that are divided into three categories: (a) spatially not dense but temporally fine data, (b) spatially dense but temporally not fine data and(c) spatially dense and temporally fine data. The first category represent high cost sensors namely automatic weather stations that are deployed at two locations and provide data at 15 minute interval. The second category includes portable soil moisture, discharge and groundwater level at weekly/ biweekly interval. The third category comprises low-cost sensors including automatic surface and groundwater level sensors installed on open wells to monitor the continuous fluctuation of water level at every 15 minutes. In addition to involving the local communities in data collection (e.g. manual rainfall measurement, water and soil sampling), this CZO also aims to provide relevant information to them for improving their sustainability. The preliminary results show significant heterogeneity in soil type, cropping system, fertilizer application, water quality, irrigation source etc. within a small catchment.

Non-proteinaceous and proteinaceous biological aerosols are abundant within the atmosphere and have the potential to impact the climate through cloud and precipitation formation. In this study, we present the differences in the laboratory-measured freezing capabilities of the non-proteinaceous and proteinaceous biological materials to determine which has more potential to impact the ice nucleation in the clouds. As non-proteinaceous surrogates, we examined multiple cellulose materials (e.g., microcrystalline and nanocrystalline cellulose) whose sizes range from ~100 nm to >100 μm (according to manufacturer report). For proteinaceous proxies, we looked at different gram-negative bacteria, such as Pseudamonas aeruginosa, Escherichia coli, Serratia marcescens, Citrobacter freundii, and Snomax, (which contains P. syringae) that can be found around the proximity of the Texas Panhandle. By using the Cryogenic Refrigeration Applied Freezing Test (CRAFT) system, we estimated immersion freezing efficiency (i.e., ice nucleation activity scaled to a unit of mass) of each sample at the temperatures greater than -30°C. We have observed that not all gram-negative bacteria has high immersion freezing activity, but the few do have a warmer temperature onset (>-20 °C) than the cellulose used. For those that did not exhibit substantial freezing efficiencies, they had similar freezing properties as the broth, in which the bacteria were incubated, as well as the cellulose materials examined. These observations suggest the presence and potential importance of bacterial cellulose in the atmospheric ice nucleation. From here, we need to conduct more in-depth investigation in the effects of a wider variety of atmospherically relevant biological aerosols to get a better understanding of the effects of said aerosols on overall aerosol-cloud interactions. Acknowledgments: K. Cory would like to acknowledge NSF-EAPSI and JSPS Summer Program for the travel fellowship support. N. Hiranuma acknowledges financial aids by the Higher Education Assistance Fund (HEAF), WTAMU Office of Graduate School and Killgore Research Center.

By using datasets of HadISST monthly SST from 1895 to 2014 and 600-year simulations of two CESM model experiments with/without doubling of CO2 concentration, ENSO characteristics are compared pre- and post- global warming. The main results are as follows. Due to global warming, the maximum climatological SST warming occurs in the tropical western Pacific (La Niña-like background warming) and the tropical eastern Pacific (El Niño-like background warming) for observations and model, respectively, resulting in opposite zonal SST gradient anomalies in the tropical Pacific. The La Niña-like background warming induces intense surface divergence in the tropical central Pacific, which enhances the easterly trade winds in the tropical central-western Pacific and shifts the strongest ocean-atmosphere coupling westward, correspondingly. On the contrary, the El Niño-like background warming causes westerly winds in the whole tropical Pacific and moves the strongest ocean-atmosphere coupling eastward. Under the La Niña-like background warming, ENSO tends to develop and mature in the tropical central Pacific, because the background easterly wind anomaly weakens the ENSO-induced westerly wind anomaly in the tropical western Pacific, leading to the so-called "Central Pacific ENSO (CP ENSO)". However, the so-called "Eastern Pacific ENSO (EP ENSO)" is likely formed due to increased westerly wind anomaly by the El Niño-like background warming. ENSO lifetime is significantly extended under both the El Niño-like and the La Niña-like background warmings, and especially, it can be prolonged by up to 3 months in the situation of El Niño-like background warming. The prolonged El Nino lifetime mainly applies to extreme El Niño events, which is caused by earlier outbreak of the westerly wind bursts, shallower climatological thermocline depth and weaker "discharge" rate of the ENSO warm signal in response to global warming. Results from both observations and the model also show that the frequency of ENSO events greatly increases due to global warming, and many more extreme El Niño and La Niña events appear under the El Niño-like and the La Niña-like background warmings, respectively. This study reconciles the phenomena and mechanisms of different characteristics of ENSO changes in observations and models.

Most climate change impacts are linked to terrestrial vegetation productivity, carbon stocks and land use change. Changes in land use and climate drive the dynamics of terrestrial carbon cycle. These carbon cycle dynamics operate at different spatial and temporal scales. Quantification of the spatial and temporal variability of carbon flux has been challenging because land-atmosphere-carbon exchange is influenced by many factors, including but not limited to, land use change and climate change and variability. The study of terrestrial carbon cycle, mainly gross primary product (GPP), net ecosystem exchange (NEE), soil organic carbon (SOC) and ecosystem respiration (Re) and their interactions with land use and climate change, are critical to understanding the terrestrial ecosystem. The main objective of this study was to examine the interactions among land use, climate change and terrestrial carbon cycling in the state of Texas using satellite measurements. We studied GPP, NEE, Re and SOC distributions for five selected major land covers and all ten climate zones in Texas using Soil Moisture Active Passive (SMAP) carbon products. SMAP Carbon products (Res=9 km) were compared with observed CO2 flux data measured at EC flux site on Prairie View A&M University Research Farm. Results showed the same land cover in different climate zones has significantly different carbon sequestration potentials. For example, cropland of the humid climate zone has higher (-228 g C/m2) carbon sequestration potentials than the semiarid climate zone (-36 g C/m2). Also, shrub land in the humid zone and in the semiarid zone showed high (-120 g C/m2) and low (-36 g C/m2) potentials of carbon sequestration, respectively, in the state. Overall, the analyses indicate CO2 storage and exchange respond differently to various land covers, and environments due to differences in water availability, root distribution and soil properties.

The International Monitoring System (IMS) of the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) includes a hydro-acoustic network as one of the monitoring technologies. The underwater part of this network includes six stations each composed of two sets of three hydrophones or triplets, except for HA01 (Cape Leeuwin, Australia) which is composed of a single triplet. The hydro-acoustic network is now complete with the recent installation of the HA04 station located in the Southern Ocean island of Crozet (France). A large number of calls emanating from marine mammals are recorded by the hydrophones, and we present examples where the animals are sufficiently close to attempt a range estimate. We also present examples of scattered arrivals and related interpretations. One striking example of extremely accurate range estimation is obtained for a whale in the neighborhood of the Cape Leeuwin (Australia) HA01 IMS stations. The proximity to the station and in particular to hydrophone H01W2 was first hypothesized because a running cross-correlation computation showed that the apparent velocity of the source was very high and could not be explained by hydro-acoustic waves travelling within the SOFAR channel. Since the far-field, plane wave assumption does not apply anymore in this case, a grid search was implemented to locate the source of the signal with the added assumption that the source is close to the ocean surface. As a further confirmation of the proximity of the source to the hydrophones, and given the expectation that such a source would generate scattering from the ocean floor and from the free surface, reflections are observed and the travel time of the scattered waves confirm the position calculated from the grid search using the direct arrivals.