Search

We combine high-resolution in situ data (acoustic Doppler current profiler (ADCP), Scanfish, and surface drifters) and remote sensing to investigate the physical characteristics of a major filament observed in the Benguela upwelling region. The 30–50 km wide and about 400 km long filament persisted for at least 40 d. Mixed-layer depths were less than 40 m in the filament and over 60 m outside of it. Observations of the Rossby number Ro from the various platforms provide the spatial distribution of Ro for different resolutions. Remote sensing focuses on geostrophic motions of the region related to the mesoscale eddies that drive the filament formation and thereby reveals |Ro|<0.1. Ship-based measurements in the surface mixed layer reveal 0.5<|Ro|<1, indicating the presence of unbalanced, ageostrophic motions. Time series of Ro from triplets of surface drifters trapped within the filament confirm these relatively large Ro values and show a high variability along the filament. A scale-dependent analysis of Ro, which relies on the second-order velocity structure function, was applied to the latter drifter group and to another drifter group released in the upwelling zone. The two releases explored the area nearly distinctly and simultaneously and reveal that at small scales ( km) Ro values are twice as large in the filament in comparison to its environment with Ro>1 for scales smaller than ∼500 m. This suggests that filaments are hotspots of ageostrophic dynamics, pointing to the presence of a forward energy cascade. The different dynamics indicated by our Ro analysis are confirmed by horizontal kinetic energy wavenumber spectra, which exhibit a power law k−α with α∼5/3 for wavelengths 2π/k smaller than a transition scale of 15 km, supporting significant submesoscale energy at scales smaller than the first baroclinic Rossby radius (Ro1∼30 km). The detected transition scale is smaller than those found in regions with less mesoscale eddy energy, consistent with previous studies. We found evidence for the processes which drive the energy transfer to turbulent scales. Positive Rossby numbers 𝒪(1) associated with cyclonic motion inhibit the occurrence of positive Ertel potential vorticity (EPV) and stabilize the water column. However, where the baroclinic component of EPV dominates, submesoscale instability analysis suggests that mostly gravitational instabilities occur and that symmetric instabilities may be important at the filament edges.

Abstract To evaluate regional axonal-related parameters as a function of disease stage in primary open angle glaucoma (POAG) and visual field (VF) sensitivity. Spectral domain optical coherence tomography was used to acquire 20° scans of POAG (n = 117) or healthy control (n = 52) human optic nerve heads (ONHs). Region specific and mean nerve fibre layer (NFL) thicknesses, border NFL and peripapillary NFL, minimum rim width (MRW)/ area (MRA) and prelamina thickness; and volume were compared across POAG disease stages and with visual field sensitivity. Differences identified between early glaucoma (EG), preperimetric glaucoma (PG) and control (C) ONHs included thinner PG prelamina regions than in controls (p

Particulate matter emissions from wildfires affect climate, weather and air quality. However, existing global and regional aerosol emission estimates differ by a factor of up to 4 between different methods. Using a novel approach, we estimate daily total particulate matter (TPM) emissions from large wildfires in North American boreal and temperate regions. Moderate Resolution Imaging Spectroradiometer (MODIS) fire location and aerosol optical thickness (AOT) data sets are coupled with HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) atmospheric dispersion simulations, attributing identified smoke plumes to sources. Unlike previous approaches, the method (i) combines information from both satellite and AERONET (AErosol RObotic NETwork) observations to take into account aerosol water uptake and plume specific mass extinction efficiency when converting smoke AOT to TPM, and (ii) does not depend on instantaneous emission rates observed during individual satellite overpasses, which do not sample night-time emissions. The method also allows multiple independent estimates for the same emission period from imagery taken on consecutive days. Repeated fire-emitted AOT estimates for the same emission period over 2 to 3 days of plume evolution show increases in plume optical thickness by approximately 10 % for boreal events and by 40 % for temperate emissions. Inferred median water volume fractions for aged boreal and temperate smoke observations are 0.15 and 0.47 respectively, indicating that the increased AOT is partly explained by aerosol water uptake. TPM emission estimates for boreal events, which predominantly burn during daytime, agree closely with bottom-up Global Fire Emission Database (GFEDv4) and Global Fire Assimilation System (GFASv1.0) inventories, but are lower by approximately 30 % compared to Quick Fire Emission Dataset (QFEDv2) PM2. 5, and are higher by approximately a factor of 2 compared to Fire Energetics and Emissions Research (FEERv1) TPM estimates. The discrepancies are larger for temperate fires, which are characterized by lower median fire radiative power values and more significant night-time combustion. The TPM estimates for this study for the biome are lower than QFED PM2. 5 by 35 %, and are larger by factors of 2.4, 3.2 and 4 compared with FEER, GFED and GFAS inventories respectively. A large underestimation of TPM emission by bottom-up GFED and GFAS indicates low biases in emission factors or consumed biomass estimates for temperate fires.

Particulate emissions from wildfires impact human health and have a large but uncertain effect on climate. Modelling schemes depend on information about emission factors, emitted particle microphysical and optical properties and ageing effects, while satellite retrieval algorithms make use of characteristic aerosol models to improve retrieval. Ground-based remote sensing provides detailed aerosol characterisation, but does not contain information on source. Here, a method is presented to estimate plume origin land cover type and age for AERONET aerosol observations, employing trajectory modelling using the HYSPLIT model, and satellite active fire and aerosol optical thickness (AOT) observations from Moderate Resolution Imaging Spectroradiometer (MODIS) and Along Track Scanning Radiometer (AATSR). It is applied to AERONET stations located in or near northern temperate and boreal forests for the period 2002–2013. The results from 629 fire attributions indicate significant differences in size distributions and particle optical properties between different land cover types and plume age. Smallest fine mode median radius (Rfv) are attributed to plumes from cropland and/or natural vegetation mosaic (0.143 μm) and grassland (0.157 μm) fires. North American evergreen needleleaf forest emissions show a significantly smaller Rfv (0.164 μm) than plumes from Eurasian mixed forests (0.193 μm) and plumes attributed to the land cover types with sparse tree cover – open shrubland (0.185 μm) and woody savannas (0.184 μm). The differences in size distributions are related to inferred variability in plume concentrations between the land cover types. Significant differences are observed between day and night emissions, with daytime emissions showing larger particle sizes. Smoke is predominantly scattering for all of the classes with median single scattering albedo at 440 nm (SSA(440)) values close to 0.95 except the cropland emissions which have an SSA(440) value of 0.9. Plumes aged for 4 days or older have median Rfv larger by ~0.02 μm compared to young smoke. Differences in size were consistent with a decrease in the Ångström Exponent and increase in the asymmetry parameter. Only an insignificant increase in SSA(λ) with ageing was found.

We develop a method to derive aerosol properties over land surfaces using combined spectral and angular information, such as available from ESA Sentinel-3 mission, to be launched in 2015. A method of estimating aerosol optical depth (AOD) using only angular retrieval has previously been demonstrated on data from the ENVISAT and PROBA-1 satellite instruments, and is extended here to the synergistic spectral and angular sampling of Sentinel-3. The method aims to improve the estimation of AOD, and to explore the estimation of fine mode fraction (FMF) and single scattering albedo (SSA) over land surfaces by inversion of a coupled surface/atmosphere radiative transfer model. The surface model includes a general physical model of angular and spectral surface reflectance. An iterative process is used to determine the optimum value of the aerosol properties providing the best fit of the corrected reflectance values to the physical model. The method is tested using hyperspectral, multi-angle Compact High Resolution Imaging Spectrometer (CHRIS) images. The values obtained from these CHRIS observations are validated using ground-based sun photometer measurements. Results from 22 image sets using the synergistic retrieval and improved aerosol models show an RMSE of 0.06 in AOD, reduced to 0.03 over vegetated targets.

We present a novel technique to infer ground slope angle from waveform LiDAR, known as the independent slope method (ISM). The technique is applied to large footprint waveforms (\(\sim\) mean diameter) from the Ice, Cloud and Land Elevation Satellite (ICESat) Geoscience Laser Altimeter System (GLAS) to produce a slope dataset of near-global coverage at \(0.5^{\circ} \times 0.5^{\circ}\) resolution. ISM slope estimates are compared against high resolution airborne LiDAR slope measurements for nine sites across three continents. ISM slope estimates compare better with the aircraft data (R\(^{2}=0.87\) and RMSE\(=5.16^{\circ}\)) than the Shuttle Radar Topography Mission Digital Elevation Model (SRTM DEM) inferred slopes (R\(^{2}=0.71\) and RMSE\(=8.69^{\circ}\)) ISM slope estimates are concurrent with GLAS waveforms and can be used to correct biophysical parameters, such as tree height and biomass. They can also be fused with other DEMs, such as SRTM, to improve slope estimates.

The effects on climate of land-cover change, predominantly from the conversion of forests to crops or grassland, are reasonably well understood for low and high latitudes but are largely unknown for temperate latitudes. The main reason for this gap in our knowledge is that there are compensating effects on the energy and water balance that are related to changes in land-surface albedo, soil evaporation and plant transpiration. We analyse how vegetation height affected the response of vegetation during the 2003 European drought using precipitation data, temperature data, normalized difference vegetation index data and a new vegetation height data set obtained from the Geoscience Laser Altimeter System (GLAS) on the Ice, Cloud and land Elevation Satellite (ICESat). At the height of the 2003 drought we find for tall vegetation a significantly smaller decrease in vegetation index and a smaller diurnal temperature (DTR) range, indicating less water stress and drought impacts on tall vegetation. Over Germany for example, 98% of significant correlations showed a smaller anomaly in vegetation index anomaly with greater height, and 95% of significant correlations showed a smaller DTR with greater vegetation height. Over France the equivalent percentages were 94 and 88%, respectively. Vegetation height is likely associated with greater rooting depth, canopy heat capacity or both. Our results suggest that land-surface models can be improved by better estimates of vegetation height and associated with this a more realistic response to drought.

In this paper we provide an overview of new knowledge on oxygen depletion (hypoxia) and related phenomena in aquatic systems resulting from the EU-FP7 project HYPOX ("In situ monitoring of oxygen depletion in hypoxic ecosystems of coastal and open seas, and landlocked water bodies", http://www.hypox.net). In view of the anticipated oxygen loss in aquatic systems due to eutrophication and climate change, HYPOX was set up to improve capacities to monitor hypoxia as well as to understand its causes and consequences. Temporal dynamics and spatial patterns of hypoxia were analyzed in field studies in various aquatic environments, including the Baltic Sea, the Black Sea, Scottish and Scandinavian fjords, Ionian Sea lagoons and embayments, and Swiss lakes. Examples of episodic and rapid (hours) occurrences of hypoxia, as well as seasonal changes in bottom-water oxygenation in stratified systems, are discussed. Geologically driven hypoxia caused by gas seepage is demonstrated. Using novel technologies, temporal and spatial patterns of water-column oxygenation, from basin-scale seasonal patterns to meter-scale sub-micromolar oxygen distributions, were resolved. Existing multidecadal monitoring data were used to demonstrate the imprint of climate change and eutrophication on long-term oxygen distributions. Organic and inorganic proxies were used to extend investigations on past oxygen conditions to centennial and even longer timescales that cannot be resolved by monitoring. The effects of hypoxia on faunal communities and biogeochemical processes were also addressed in the project. An investigation of benthic fauna is presented as an example of hypoxia-devastated benthic communities that slowly recover upon a reduction in eutrophication in a system where naturally occurring hypoxia overlaps with anthropogenic hypoxia. Biogeochemical investigations reveal that oxygen intrusions have a strong effect on the microbially mediated redox cycling of elements. Observations and modeling studies of the sediments demonstrate the effect of seasonally changing oxygen conditions on benthic mineralization pathways and fluxes. Data quality and access are crucial in hypoxia research. Technical issues are therefore also addressed, including the availability of suitable sensor technology to resolve the gradual changes in bottom-water oxygen in marine systems that can be expected as a result of climate change. Using cabled observatories as examples, we show how the benefit of continuous oxygen monitoring can be maximized by adopting proper quality control. Finally, we discuss strategies for state-of-the-art data archiving and dissemination in compliance with global standards, and how ocean observations can contribute to global earth observation attempts.

This work provides a comparison of satellite retrievals of Saharan desert dust aerosol optical depth (AOD) during a strong dust event through March 2006. In this event, a large dust plume was transported over desert, vegetated, and ocean surfaces. The aim is to identify the differences between current datasets. The satellite instruments considered are AATSR, AIRS, MERIS, MISR, MODIS, OMI, POLDER, and SEVIRI. An interesting aspect is that the different algorithms make use of different instrument characteristics to obtain retrievals over bright surfaces. These include multi-angle approaches (MISR, AATSR), polarisation measurements (POLDER), single-view approaches using solar wavelengths (OMI, MODIS), and the thermal infrared spectral region (SEVIRI, AIRS). Differences between instruments, together with the comparison of different retrieval algorithms applied to measurements from the same instrument, provide a unique insight into the performance and characteristics of the various techniques employed. As well as the intercomparison between different satellite products, the AODs have also been compared to co-located AERONET data. Despite the fact that the agreement between satellite and AERONET AODs is reasonably good for all of the datasets, there are significant differences between them when compared to each other, especially over land. These differences are partially due to differences in the algorithms, such as assumptions about aerosol model and surface properties. However, in this comparison of spatially and temporally averaged data, it is important to note that differences in sampling, related to the actual footprint of each instrument on the heterogeneous aerosol field, cloud identification and the quality control flags of each dataset can be an important issue.

We present new coarse resolution (0.5° × 0.5°) vegetation height and vegetation-cover fraction data sets between 60° S and 60° N for use in climate models and ecological models. The data sets are derived from 2003–2009 measurements collected by the Geoscience Laser Altimeter System (GLAS) on the Ice, Cloud and land Elevation Satellite (ICESat), the only LiDAR instrument that provides close to global coverage. Initial vegetation height is calculated from GLAS data using a development of the model of Rosette et al. (2008) with with further calibration on desert sites. Filters are developed to identify and eliminate spurious observations in the GLAS data, e.g. data that are affected by clouds, atmosphere and terrain and as such result in erroneous estimates of vegetation height or vegetation cover. Filtered GLAS vegetation height estimates are aggregated in histograms from 0 to 70 m in 0.5 m intervals for each 0.5° × 0.5°. The GLAS vegetation height product is evaluated in four ways. Firstly, the Vegetation height data and data filters are evaluated using aircraft LiDAR measurements of the same for ten sites in the Americas, Europe, and Australia. Application of filters to the GLAS vegetation height estimates increases the correlation with aircraft data from r = 0.33 to r = 0.78, decreases the root-mean-square error by a factor 3 to about 6 m (RMSE) or 4.5 m (68% error distribution) and decreases the bias from 5.7 m to −1.3 m. Secondly, the global aggregated GLAS vegetation height product is tested for sensitivity towards the choice of data quality filters; areas with frequent cloud cover and areas with steep terrain are the most sensitive to the choice of thresholds for the filters. The changes in height estimates by applying different filters are, for the main part, smaller than the overall uncertainty of 4.5–6 m established from the site measurements. Thirdly, the GLAS global vegetation height product is compared with a global vegetation height product typically used in a climate model, a recent global tree height product, and a vegetation greenness product and is shown to produce realistic estimates of vegetation height. Finally, the GLAS bare soil cover fraction is compared globally with the MODIS bare soil fraction (r = 0.65) and with bare soil cover fraction estimates derived from AVHRR NDVI data (r = 0.67); the GLAS tree-cover fraction is compared with the MODIS tree-cover fraction (r = 0.79). The evaluation indicates that filters applied to the GLAS data are conservative and eliminate a large proportion of spurious data, while only in a minority of cases at the cost of removing reliable data as well. The new GLAS vegetation height product appears more realistic than previous data sets used in climate models and ecological models and hence should significantly improve simulations that involve the land surface.

Remote sensing of aerosol from space is a challenging and typically underdetermined retrieval task, requiring many assumptions to be made with respect to the aerosol and surface models. Therefore, the quality of a priori information plays a central role in any retrieval process (apart from the cloud screening procedure and the forward radiative transfer model, which to be most accurate should include the treatment of light polarization and molecular-aerosol coupling). In this paper the performance of various algorithms with respect to the of spectral aerosol optical thickness determination from optical spaceborne measurements is studied. The algorithms are based on various types of measurements (spectral, angular, polarization, or some combination of these). It is confirmed that multiangular spectropolarimetric measurements provide more powerful constraints compared to spectral intensity measurements alone, particularly those acquired at a single view angle and which rely on a priori assumptions regarding the particle phase function in the retrieval process.

The Geological Survey of Canada (GSC) is making preparations for Canadian participation in GSETT 3 but will be unable to make a formal commitment until the necessary resources have been secured. As Canada is expected to provide at least four alpha stations, and a significant number of beta stations, the financial resources that will be needed are substantial, even though in many respccts the GSC is, with the recent modernization of the Yellowknife array and the ongoing installation of the Canadian National Seismograph Network (CNSN), well positioned to make a significant contribution to GSETT 3. The CNSN currently (October 1993) consists of 17 broad band stations and will grow to 23 and 33 such stations by December 1993 and December 1994 respectively. Some 40 50 short period stations will complete the network. Data from all sites are continuously telemetered in real time to network acquisition centres in Ottawa and Sidney, British Columbia, archived to optical disk, and kept on line in a 72 h ring buffer. Most of the broadband sites could serve as either alpha or beta stations once the necessary software for continuous data transfer, or on request provision, of data from the selected sites has been completed. This software wili be configured so that changes in station selection are easy to implement, and this will provide considerable flexibility to the GSETT 3 planning and operations working groups in selecting the optimum network. Backup stations can be designated in the case of station failures, and the network centre in British Columbia will serve, at least for beta stations, as a backup NDC to that in Ottawa. Data from. the Yellowknife array are collected in Yellowknife and forwarded in ten minute files to Ottawa, where processing is completed and the results archived. This arrangement would not meet the deadlines for receipt of alpha station data at the IDC and new hardware and software will be needed to forward the data more immediately from Yellowknife to Ottawa. Although the procedures and formats for both alpha and beta station data have not yet been agreed upon, or even discussed, by the GSE, it is apparent that new facilities will be re quired in Ottawa to multiplex and reformat data for transmission to the IDC. We anticipate that a dedicated 56 kbaud link will be needed between Ottawa. and Washington.

The Canadian National Seismograph Network currently consists of 5 very-broadband (VBB) and 15 broadband (BB) stations across Canada, supplemented by 6 short period (SP) stations. When it is completed by the end of 1995, a further 1 VBB, 12 BB and over 40 SP stations will have been added. Data from all sites are telemetered in real time to twin network acquisition, processing and archiving centres in Eastern and Western Canada. All data are continuously archived in SEED format on optical disk and access to the most recent three days of data is provided through a mail-based AutoDRM system. Continuous data from the VBB sites are sent to the FDSN Data Management Centre approximately one month after being recorded.

An algorithm for the automatic detection and association of surface waves has been developed and tested over an 18 month interval on broad band data from the Yellowknife array (YKA). The detection algorithm uses a conventional STA/LTA scheme on data that have been narrow band filtered at 20 s periods and a test is then applied to identify dispersion. An average of 9 surface waves are detected daily using this technique. Beamforming is applied to determine the arrival azimuth; at a nonarray station this could be provided by poIarization analysis. The detected surface waves are associated daily with the events located by the short period array at Yellowknife, and later with the events listed in the USGS NEIC Monthly Summaries. Association requires matching both arrival time and azimuth of the Rayleigh waves. Regional calibration of group velocity and azimuth is required. . Large variations in both group velocity and azimuth corrections were found, as an example, signals from events in Fiji Tonga arrive with apparent group velocities of 2.9 3.5 krn/s and azimuths from 5 to + 40 degrees clockwise from true (great circle) azimuth, whereas signals from Kuriles Kamchatka have velocities of 2.4 2.9 km/s and azimuths off by 35 to 0 degrees. After applying the regional corrections, surface waves are considered associated if the arrival time matches to within 0.25 km/s in apparent group velocity and the azimuth is within 30 degrees of the median expected. Over the 18 month period studied, 32% of the automatically detected surface waves were associated with events located by the Yellowknife short period array, and 34% (1591) with NEIC events; there is about 70% overlap between the two sets of events. Had the automatic detections been reported to the USGS, YKA would have ranked second (after LZH) in terms of numbers of associated surface waves for the study period of April 1991 to September 1992.