ALBERT

Description: The controls that affect the degree of spatio-temporal change to foredunes following scarping are reviewed herein. As sea-levels rises and climate changes, dune scarping will become more common. Thus, it is critical to understand what factors contribute to the magnitude of scarping, and what effect this has on dune systems to better manage coastal erosion into the future. Scarping occurs when foredunes are partially eroded by waves, generally during periods of high water level. The controls on the degree and magnitude of scarping examined include water level, foredune vegetation cover and species present, plant root mass, height and volume of the foredune, the original foredune morphology, surfzone–beach type, and compaction of sediment. Water-level height and duration of high water is the most significant control as it determines the elevation at which wave action can erode the dune and, therefore, the extent of scarping and dune volumetric loss. Higher plant density, greater rooting depth, high root mass, and greater compaction aid in reducing the degree of scaping. The presence of large woody debris and wrack may also influence the degree of scarping. The effects scarping has on the morphology of a foredune after the initial erosion event can range from small changes (e.g. minor, small scarps and slight slumping), to moderate changes such as the foredune translating landwards, to large change such as the transition of an entire dune system into a new transgressive dunefield phase. A new model summarising the key controls and their relationship/significance to the magnitude and extent of scarping is presented.

Description: Humans kill and care for animals in a multitude of contexts. These themes – killing and caring – form the focus of this second report on animal geographies research. Most notably, killing and caring take place through conservation and the production and consumption of food. Other realms of recent research include killing through climate change, formal arrangements of care, how animals are made killable, and the significance of the individual and collective. Further to these two major themes, the review identifies other approaches to animal geographies research in recent years; namely, political and relational. Finally, the propensity for humans to kill and care for animals is highlighted by crisis. I open and close the review by considering connections between animals and two contemporary emergencies: the fires that burned south-eastern Australia over summer 2019–20 and the COVID-19 pandemic.

Description: This paper argues for a sustained study of knowledge-intensive or liberal professions in geography. I review existing work in political geography and related fields to identify a gap in the study of knowledge-intensive professions, especially those that are popularly associated with elites. I draw from sociology, anthropology, and international relations to explain why we need to better understand such professions. By the geographical study of professions and their expertise I mean the examination of the places, spatial networks, and travels of ideas that shape these professions and the expertise created therein.

Description: Exploring ancient socio-economic adaptation is a basic issue of human-environment interaction. Xinjiang in northwest China is a region of high geographic diversity. Past human adaptations to this arid marginal area is a current focus of research interest but still lacks in-depth study. This article presents data from the Wupu Cemetery, located in the extremely arid Hami Basin in the eastern Tianshan Mountains. Archaeobotanical analysis is used to reconstruct the local environment niche and the subsistence economy of inhabitants. Radiocarbon dating results indicate the cemetery was occupied between 3000 and 2400 cal BP, during the late Bronze Age and early Iron Age. In total 16 species of the plant remains are identified, including four cereal crops, foxtail millet ( Setaria italica), broomcorn millet ( Panicum miliaceum), naked barley ( Hordeum vulgare var. coeleste), wheat ( Triticum aestivum) and 12 wild types. The riparian plant Populus euphratica and aquatic plant Typha sp. indicate inhabitants lived in an oasis near the cemetery. Environmental interpretation of this data compares well with other seven sites in arid southern Xinjiang. In addition to faunal remains from the site, it is assumed that a flexible system of multi-crop farming and herding was the subsistence pattern around Wupu. This system was widespread across Inner Asia and appears to have played a central role in adapting to different marginal environments during the Bronze Age and Iron Age.

Description: Transcontinental exchange emerged and intensified in northern China since the late fifth millennium BP (Before present), especially in the arc, which was the core area of the eastern part of the trans-Eurasian exchange during the Late Neolithic and the Bronze Age. In the arc, the exchange profoundly affected the human subsistence strategy and human-environment relationship. Relative to the crop patterns and human diets during the Bronze Age in northern China, systematic investigations of zooarcheological data based on broad spatial and temporal framework to understand the influence of introduced livestock and indigenous livestock on human subsistence are lacking. To show the spatial-temporal variation in animal utilization patterns and its relation to prehistoric trans-Eurasian exchange, the zooarcheological data from 40 sites in northern China dated between 5000 and 2500 BP were analyzed. The strategy of animal utilization in northern China changed substantially from 5000 to 2500 BP, with notable spatial features in different chronological phases. From 5000 to 4300 BP, wild mammals and indigenous livestock (pig, dog) use dominated in the arc and the North China Plain (NCP). During 4300–3500 BP, the importance of introduced livestock (cattle, sheep/goat, horse) exceeded that of indigenous livestock in the arc, whereas indigenous livestock continued to dominate in the NCP. Indigenous livestock acted as the most important animal subsistence in northern China, although the exploitation of introduced livestock increased during 3500–2000 BP. These spatio-temporal differences in animal utilization appear to be closely associated with the prehistoric trans-Eurasian exchange, but were also affected by local environment, agriculture development, and climate change.

Description: We test a recent prediction that stable carbon isotope ratios from UK oaks will display age-trends of more than 4‰ per century by measuring 〉5400 carbon isotope ratios from the late-wood alpha-cellulose of individual rings from 18 modern oak trees and 50 building timbers spanning the 9th–21st centuries. After a very short (c.5 years) juvenile phase with slightly elevated values, the number of series that show rising and falling trends is almost equal (33:35) and the average trend is almost zero. These results are based upon measuring and averaging the trends in individual time-series; the ‘mean of the slopes’ approach. We demonstrate that the more conventional ‘slope of the mean’ approach can produce strong but spurious ‘age-trends’ even when the constituent series are flat, with zero slope and zero variance. We conclude that it is safe to compile stable carbon isotope chronologies from UK oaks without de-trending. The isotope chronologies produced in this way are not subject to the ‘segment length curse’, which applies to growth measurements, such as ring width or density, and have the potential to retain very long-term climate signals.

Description: In Australia, the drivers of precolonial fire regimes remain contentious, with some advocating an anthropogenic-dominated regime, and others highlighting the importance of climate, climatic variability or alternatively some nexus between climate and human activity. Here, we explore the inter-relationships between fire, humans and vegetation using macroscopic charcoal, archaeology and palynology over the last ~5430 cal. year BP from Broughton Island, a small, near-shore island located in eastern Australia. We find a clear link between fire and the reduction of arboreal pollen and rainforest indicators on the island, especially at ~4.0 ka and in the last ~1000 years. Similarities with comparable palaeoenvironmental records of fire in the region and a record of strong El Niño (dry, fire-prone) events supports the contention that climate was a significant influence on the fire regimes of Broughton Island. However, two periods of enhanced fire activity, at ~4000 years BP and ~

Description: Wildfire is a ubiquitous disturbance agent in subalpine forests in western North America. Lodgepole pine ( Pinus contorta var. latifolia), a dominant tree species in these forests, is largely resilient to high-severity fires, but this resilience may be compromised under future scenarios of altered climate and fire activity. We investigated fire occurrence and post-fire vegetation change in a lodgepole pine forest over the past 2500 years to understand ecosystem responses to variability in wildfire and climate. We reconstructed vegetation composition from pollen preserved in a sediment core from Chickaree Lake, Colorado, USA (1.5-ha lake), in Rocky Mountain National Park, and compared vegetation change to an existing fire history record. Pollen samples ( n = 52) were analyzed to characterize millennial-scale and short-term (decadal-scale) changes in vegetation associated with multiple high-severity fire events. Pollen assemblages were dominated by Pinus throughout the record, reflecting the persistence of lodgepole pine. Wildfires resulted in significant declines in Pinus pollen percentages, but pollen assemblages returned to pre-fire conditions after 18 fire events, within c.75 years. The primary broad-scale change was an increase in Picea, Artemisia, Rosaceae, and Arceuthobium pollen types, around 1155 calibrated years before present. The timing of this change is coincident with changes in regional pollen records, and a shift toward wetter winter conditions identified from regional paleoclimate records. Our results indicate the overall stability of vegetation in Rocky Mountain lodgepole pine forests during climate changes and repeated high-severity fires. Contemporary deviations from this pattern of resilience could indicate future recovery challenges in these ecosystems.

Description: The aim of this paper is to elucidate the dynamic interplay of visuality, space and power through an analysis of what we call scopic relational spaces (SRS). Our primary claim is that scopic relations are intrinsically spatial relations and scopic practices are spatial practices. We contend that such analysis facilitates the discernment of significant socio-spatial processes and events that are otherwise unrecognized. We suggest that attention to the socio-spatial phenomena under investigation can contribute to debates about the alleged incommensurability of territorial and relational spatial imaginaries.

Description: The Guizhou karst area is one of the largest continuous areas of karst in the humid climate zone and is representative of karst landforms in China. Large portions of the karst system are characterized by extremely shallow soils underlain by weathered bedrock and water deficits are common. Although the distribution of ecosystem productivity is largely related to variations in the temperature and precipitation, the influence of the substrate in karst areas requires further exploration. We explored the relative importance of the bedrock geochemistry (characterized by the concentrations of Ca, Mg and Si) and climatic factors (temperature and precipitation) to explain the spatial variability in gross primary productivity (GPP) with various degrees of water deficit during the time period 2001–2015. Our results show that the impact of bedrock geochemistry is an important parameter in changing the original relationship between climate and the GPP. The bedrock geochemistry functioned as a “regulator” of the relation between climate and the GPP, which strengthened with decreasing climate favourability. The variations in GPP and surface water storage were significantly different when different elements (Ca, Mg or Si) were dominant. The Mg-rich regions showed the greatest annual variations in the GPP, whereas the Si-rich regions had the strongest surface water storage potential to support vegetation growth. The results of our study are important for systematically evaluating the effects of climate on vegetation productivity and provide a benchmark for global vegetation modelling predictions.

Description: In this paper, we identify the ways in which the existing literature has examined financial technology (FinTech). Using the frame of the ‘FinTech Cube’, we examine how FinTech unfolds through the intersections of key actors, technologies and institutions. We demonstrate the relevance of FinTech for two areas of geographical enquiry: i) the reshaping of global production and financial networks, and ii) financial inclusion and poverty reduction in poorer countries. In doing so, we accord particular attention to the significance of FinTech for theoretical and empirical research in economic geography.

Description: The Late Antique Little Ice Age, spanning the period from 536 CE to roughly 560 CE, saw temperatures in the Northern Hemisphere drop by a degree C in less than a decade. This rapid cooling is thought to have caused widespread famine, epidemic disease, and social disruption. The relationship between cooling and social disruption is examined here using a set of high-resolution climate and historical data. A significant link between cooling and social disruption is demonstrated, but it is also demonstrated that the link is highly variable, with some societies experiencing dramatic cooling changing very little, and others experiencing only slight cooling changing dramatically. This points to variation in vulnerability, and serves to establish the Late Antique Little Ice Age as a context within which naturalistic quasi-experiments on vulnerability to climate change might be conducted.

Description: The first aim of this corrigendum is to point out and correct calculation errors on solar depression angle and azimuth angle in Bertolin and Domínguez-Castro (2020a, 2020b). The second aim is to recognize that these calculations are correct in Dalin (2020). The third aim is to analyze the chances of Antonio Colla to observe the noctilucent cloud (NLC) taking into account the correct calculations of the twilight sky arc determining the illuminated area of an NLC and the uncertainties in the Colla’s observation report.

Description: In this paper, we propose that there is a politics of encounters centered on the body at play in seeking asylum and refuge, and that it is critical to study how it unfolds from the point of view of both governing and agency. Building on existing work that looks at the role of embodiment in the political struggles of refugees, and leaning on Helmuth Plessner’s original thinking about social embodiment, we develop a theoretical understanding of this political dynamic, illustrating how it can help us make sense of power relations and forms of governance and (latent) resistance involved in it.

Description: In this first report, I overview research on storytelling of the histories of GIS. I suggest that these efforts reveal important pathways for thinking the current moment and conditioning technoscience futures. I argue that these stories around technology illustrate various disciplinary crucibles around technical practices and knowledge work. To report on progress in GIScience is also to narrate the moments of change in the discipline, whether marked as quantitative, qualitative, radical, critical, black, feminist, Marxist, indigenous, environmental, postcolonial, queer, or minor, among others. I suggest ‘minor GISciences’ to encourage attentiveness to the historical conditions and social implications of geographic technologies.

Description: Territorial stigmatisation emerged over the past decade as the prevailing concept for understanding the phenomenon of disreputable places and the processes by which they are produced. Following the work of Loïc Wacquant, several studies have articulated its centrality to the neoliberal restructuring of capital and the state. Departing from Wacquant, several have also examined the various forms of resistance to territorial stigmatisation. In reviewing this literature, this paper argues and outlines how territory – paradoxically under-theorised in the literature to date – can clarify the production of territorial stigmatisation, the obfuscations and legitimations it performs, and resistance and contestation.

Description: This article identifies new geographies of commemoration. These draw on non-representational perspectives that foreground the experiential aspects of commemorative sites and events, drawing these together with methodologies alive to the subtle but excessive feelings that arise in such settings. It argues that one aim of commemoration – to reinforce the contours of national identity – is disrupted by a focus on the experiential world because of the unpredictable and excessive nature of sensory and affective feelings. New research in geographies of commemoration also draws together different temporalities, holding the potential to unsettle and complicate national narratives.

Description: The discovery that orbital variations are the driving force behind Quaternary climate change provides an impetus to set local and regional records of environmental change into the global context, a principle that has been strongly embraced by Quaternary scientists working in New Zealand. Their major achievements and significant current initiatives are reviewed here. The importance of the New Zealand Quaternary stems from its geographical context: a climatically sensitive, remote oceanic, southern location spanning 17 degrees of the mid-latitudes; an obliquely convergent plate boundary setting resulting in a high mountain range athwart the prevailing westerlies, active volcanism, a youthful and dynamic landscape, and mountains high enough to maintain glaciers today; and a remarkably short prehistory. The resultant records show marked environmental changes due not only to climatic oscillations but also to vigorous, active tectonism and volcanism. The Taupo Volcanic Zone, containing the world's strongest concentration of youthful rhyolitic volcanoes, has produced at least 10 000 km3of magma in the last 2 Ma. Climatic interpretations of records from marine sediments in the New Zealand region, together with several long sequences of alternating marine and terrestrial sediments, indicate broad synchrony with Northern Hemisphere events (within limitations of dating), although there are differences in detail for shorter-term climatic events. It is not yet certain that glacial advances coincided precisely with those in the Northern Hemisphere or were of similar duration. Late Cainozoic glaciation commenced c. 2.6-2.4 Ma but the record of glacial deposits is fragmentary and poorly dated except for the most recent events. The Last (Otira) Glaciation, from c. 100-10 ka, was characterized by at least five glacial advances including during the Last Glacial Maximum from 25 to 15 ka, when snowlines fell by 600-800 m. New Zealand evidence for cooling during the Younger Dryas stade is equivocal whilst isotopic records from speleothems, and other data, indicate warmer and wetter conditions from 10-7 ka, broadly conforming with records from mid-latitude Northern Hemisphere locations. Future advances will require sampling at shorter timescales, improvements in the accuracy and precision of existing dating methods and the development of new ones, extension of palaeoecological techniques to cover the full potential of new Zealand's diverse biota, and a stronger emphasis on quantification of palaeoclimatic parameters.

Description: Affecting an area of ca. 800 000 km2and killing up to 100 000 people, the Lisbon earthquake of 1755 is probably the greatest seismic disaster to have struck western Europe. The shock waves of the earthquake placed a temporary brake on the emerging rationalism of the European Enlightenment and attempts to explain the disaster in terms of human sinfulness coloured many contemporary accounts. Notwithstanding these difficulties, through careful archival research it has proved possible to obtain relatively value-free accounts of most aspects of the earthquake and to use these not only to model the physical characteristics of and damage caused by the earthquake, but also to consider the implications for present day hazard assessment and urban planning. This paper reviews the progress that has been made in: identifying source and faulting mechanisms; the processes involved in the generation and impact of tsunamis; damage caused to different types of building and the use being made of historical earthquakes of different sizes - of which the 1755 event is the largest - in defining future hazard scenarios for Lisbon and other areas of Iberia.

Description: This paper reports the first detailed palaeogeographical analysis of the environmental context of late Mesolithic shell midden sites in the lower Tagus area and focuses on the lower Muge valley, which contains an internationally significant Mesolithic record. The lower Muge valley fill comprises buried estuarine and fluvial environments contemporary with Mesolithic settlement. Holocene environmental and palaeogeographic changes influenced Mesolithic settlement-subsistence and midden accumulation. The sudden appearance of large late Mesolithic shell middens throughout Portugal represents a process of increased visibility and preferential preservation of the archaeological record. Prior to ~6100 cal. BC, aggrading valley floor environments did not occupy the entire width of the present lower Tagus floodplain and any sites located in the early Holocene valley are currently deeply buried. Shell midden occupation on terrace levels followed the establishment of aggrading estuarine environments, containing productive shell beds, near the mouth of the lower Muge valley at ~6100 cal. BC. The critical factors in site choice appear to have been the nearby presence of (i) rich shell resources and (ii) freshwater environments. Long-term site occupation and (semi-)sedentary behaviour was favoured by the local presence, for over 2000 years, of rich resources from estuarine, freshwater and open woodland environments. Site abandonment (~5300—4800 cal. BC) coincided with the regional establishment of an open landscape (~5000 cal. BC) and the contraction of local estuarine environments (~5555—3800 cal. BC). The associated gradual decrease in resources and cultural interaction with the expanding early Neolithic communities may have influenced Mesolithic site abandonment.

Description: We present the first step to quantitatively reconstruct vegetation cover in central Europe. Modern vegetation and pollen deposition were compared for 20 small to medium sized lakes and their catchments on the Swiss Plateau, a relatively flat region between the Jura Mountains and the Alps. To correct for the pollen dispersal bias in pollen assemblages, vegetation abundance was distance-weighted using three different approaches. The Relevant Source Area of Pollen (RSAP) and pollen productivity of 13 plant taxa were estimated using three different submodels of the Extended R-Value model (ERV-model). RSAP was 800 m regardless of the applied distance-weighting or ERV submodel. Pollen Productivity Estimates (PPE) varied from 10 to 〈 0.1 among pollen taxa and differed slightly between the models. Relative to grasses most trees were higher pollen producers and some were equal producers, whereas the herb taxa showed lower PPE. Generally, PPE from lowland Switzerland differ from those found in other European regions. Sampling strategies of vegetation and pollen samples are a likely cause for this variation. However, pollen productivity is also influenced by regionally different factors, such as climate, vegetation structure, geology and soil types. In addition, differences at genus or species level may occur between areas. Our comparison between the different regions in Europe shows that PPE of one region may not be directly applicable to other regions.

Description: Wetlands are poorly documented features of many landscapes, and there is often little understanding of the geomorphological controls on their origin, development and characteristics. This paper addresses the apparent paradox of wetlands in drylands, focusing particularly on the geomorphology and sedimentology of wetlands in southern Africa. Drylands are characterized by high (but variable) levels of aridity, reflecting low ratios between precipitation and potential evapotranspiration, so wetlands can only exist where there are locally positive surface water balances for all or part of the year. Most moderate to large wetlands in drylands are thus maintained by river inflows that combine with other factors that serve to impede drainage or reduce infiltration, including faulting, rock outcrops, swelling soils, and ponding by tributary or aeolian sediments. Together with variations in sediment supply, vegetation communities, and levels of animal activity, this promotes a diverse range of wetlands that span a continuum from permanently inundated, to seasonally inundated, to ephemerally inundated. In detail, every wetland has a unique range of geomorphological and sedimentological characteristics but, at a general level, the dryland setting can be shown to impart some distinctive features. By comparison with humid region (tropical and temperate) wetlands, we propose that many wetlands in drylands are characterized by: 1) more frequent and/or longer periods of desiccation; 2) channels that commonly decrease in size and even disappear downstream; 3) higher levels of chemical sedimentation; 4) more frequent fires that reduce the potential for thick organic accumulations and promote aeolian activity; and 5) longer timescales of development that may extend far back into the Pleistocene. Additional studies of wetlands in different drylands may reveal other distinctive characteristics. Correct identification of the factors giving rise to wetlands, and improved understanding of the geomorphological and sedimentological processes governing their development, is vital for the design of sustainable management guidelines for these diverse yet fragile habitats.

Description: The glaciers and snowfields of the Southern Alps of New Zealand are the most significant in the Southern Hemisphere outside Antarctica and South America. The most substantial data on Southern Hemisphere glacier fluctuations come from New Zealand. The nature and behaviour of New Zealand's glaciers are also of wider scientific interest, because they are highly sensitive, high input-output systems that represent the temperate, maritime end of the glacier process-behaviour continuum. The areal extent and volume of glaciers and snow are outlined and an assessment is made of their scientific relevance and of their importance as resources and hazards. The main themes and progress of research on glaciers and snow, including snow avalanches, are reviewed. Glacier research has concentrated on only a few key glaciers and has focused on understanding glacier change. Main topics covered in this review relate to this focus and include fluctuations in termini, other mass balance signals and response to climate variability. Research on mass balance processes, glacier dynamics and glacier hydrology is also outlined. Seasonal snow has received less attention until recently. The main emphasis has been on quantification and past variability and its contribution to river flow, particularly in the most important hydroelectric power catchments of the South Island. Some field measurements have been made of the energy balance over snow. Research on snow avalanches has grown as the demands of winter recreation and alpine tourism have increased the hazard. Research first concentrated on production of avalanche atlases for the most hazardous areas and on quantifying the nature of the hazard. Subsequently, there has been a shift towards more process studies that are related to avalanche formation and runout distance. The main gaps in research on glaciers and snow are identified and key areas for future work proposed. There is an urgent need, in particular, for glacier mass-balance measurements. Extensive data on snow structure need to be synthesized. Satellite imagery should be used for monitoring of seasonal snow. Snow melt during northwest storms needs to be better defined. A more developed engineering approach is required for the study of snow avalanches. New Zealand offers exciting possibilities for the study of cryospheric processes, including response to future climate change.

Description: Mediterrranean regions are characterized by high spatiotemporal heterogeneity of vegetation patterns. Understanding the dynamic nature of these environments requires detailed data for wide regions regarding changes in their phyto-ecology, biomass and productivity. This article assesses the current status of satellite remote sensing in this field of application. Mapping the five main life-forms (physiognomic classes) in Mediterranean regions (forests, woodlands, scrub, dwarf shrubs and herbaceous growth) has attracted major attention in recent years. Methodologies developed for this purpose are based on the spectral, temporal and spatial (textural) information domains provided by satellite data. Wide regional vegetation mapping was achieved using phenological classification of vegetation indices derived mainly from NOAA AVHRR images. More detailed mapping was conducted with multispectral techniques in local areas using mainly Landsat TM images. Assessments of multispectral and multi-temporal categories have shown limitations in their applicability over wide regions due mainly to the heterogeneity of Mediterranean regions. This heterogeneity cannot be regarded as a simple mixing of life-forms over large areas but, rather, the formation of transitional zones of varying mixtures resulting from disturbance and recovery cycles. Productivity and biomass monitoring has been found to be an active methodological development due to the introduction of new off-nadir viewing sensors in the visible and infrared spectral bands, and because of the development of methodologies for the retrieval of biophysical information from Synthetic Aperature Radar (SAR) data. Studies of ecosystem evolution using satellite data were conducted mainly in the fields of fire disturbance and desertification. Further progress in the remote sensing of Mediterranean vegetation ecology requires a better synergy of sensors, methods and ancillary data.

Description: Pollen analysis of a 48 m AMS radiocarbon-dated sediment sequence from the Guadiana estuary provides the first record of Lateglacial and Holocene vegetation history in the Algarve province of Portugal. This paper focuses on the record of terrestrial pollen taxa, which document a series of forest expansions and declines during the period 13 000 cal. BP to 1600 cal. BP and provide insights into climate evolution in southwestern Iberia. The main vegetation phases identified in the Guadiana valley are (1) Lateglacial interstadial (Allerød) forest with Quercus and Pinus under a temperate, moist, continental climate; (2) a Younger Dryas forest decline ( Quercus) and expansion of pinewoods, xeric scrub and open ground habitats (with Juniperus , Artemisia, Ephedra distachya type, Centaurea scabiosa type) under arid and cold conditions; (3) an early Holocene forest/scrub/open-ground vegetation mosaic developing under a warm, dry and continental climate; (4) a maximum of Quercus forest and thermomediterranean evergreen taxa ( Olea, Phillyrea, Pistacia) reflecting a warm, moist oceanic climate between c. 9000 cal. BP and c. 5000 cal. BP; and (5) the expansion of shrublands with Cistaceae and Ericaceae under a drier climatic regime and increasing anthropogenic activity since c. 5000 cal. BP. Holocene episodes of maximum climatic aridity are identified in the record of xerophytic taxa ( Juniperus, Artemisia, Ephedra distachya type) centred around 10 200 cal. BP, 7800 cal. BP, 4800 cal. BP, 3100 cal. BP and 1700 cal. BP. Regional comparisons suggest a correlation of arid phases across southern Iberia and northwest Africa, which can be related to abrupt North Atlantic coolings (Bond events).

Description: A multiproxy record including benthic foraminifera, diatoms and XRF data of a marine sediment core from a SW Greenland fjord provides a detailed reconstruction of the oceanographic and climatic variations of the region during the last 4400 cal. years. The lower part of our record represents the final termination of the Holocene Thermal Maximum. After the onset of the `Neoglaciation' at approximately 3.2 ka cal. BP, the fjord system was subject to a number of marked hydrographical changes that were closely linked to the general climatic and oceanographic development of the Labrador Sea and the North Atlantic region. Our data show that increased advection of Atlantic water (Irminger Sea Water) from the West Greenland Current into the Labrador Sea was a typical feature of Northeast Atlantic cooling episodes such as the `Little Ice Age' and the `European Dark Ages', while the advection of Irminger Sea Water decreased significantly during warm episodes such as the `Mediaeval Warm Period' and the `Roman Warm Period'. Whereas the `Mediaeval Warm Period' was characterized by relatively cool climate as suggested by low meltwater production, the preceding `Dark Ages' display higher meltwater runoff and consequently warmer climate. When compared with European climate, these regional climate anomalies indicate persisting patterns of advection of colder, respectively warmer air masses in the study region during these periods and thus a long-term seesaw climate pattern between West Greenland and Europe.

Description: This article presents an overview of palaeofluvial geomorphology research in southern Africa. For the purposes of this article this includes South Africa, Zimbabwe, Namibia, Lesotho, Swaziland and Botswana. Although interest in fluvial systems has a long history in southern Africa, the scientific study of rivers was initiated by the discovery of the first alluvial diamond along the banks of the Orange River in 1867. Since then, significant progress has been made in unravelling the fluvial history of southern Africa from the early Archaean Ventersdorp Contact Reef River to modern channel process studies. The development of an understanding of palaeofluvial systems has occurred along two main lines. The first was alluvial diamond exploration work undertaken by the large mining houses. The second line was of a more ‘academic’ interest and included determining the impact of superimposition, tectonics, base level and climate changes. The review suggests that southern Africa fluvial systems have shown large-scale changes in drainage pattern, discharge and sediment yield and that these can be related to a complex set of causative factors including the geological template, the Jurassic rifting of Gondwana, tectonic episodes and climate change.

Description: Sedimentological studies of coarse-grained alluvial rivers reveal patterns of bed material sorting at a variety of spatial scales ranging from downstream fining over the length of the long profile to the vertical segregation of a coarse surface layer at the scale of individual particles. This article reviews the mechanisms that sort bed material by size during sediment entrainment, transport and deposition and discusses some of the inter-relationships that exist between patterns and processes of sediment sorting at different spatial and temporal scales. At initiation of motion, sorting can arise from the preferential entrainment of the finer fractions from the heterogeneous bed sediments. Bedload grain-size distributions are modified during transport as different size fractions are routed along different transport pathways under the influence of nonuniform bed topography and associated flow patterns, and during deposition as the variable pocket geometry of the rough bed surface and turbulence intensity of the flow control the size of the particles that deposit. The review highlights the poor understanding of the many feedback linkages that exist between patterns and processes of sediment sorting at different scales and the need for a greater awareness of the spatial and temporal bounds of these linkages.

Description: Wherever people gain their livelihood in mountains and steeplands, the productive capacity of the soils they use is likely to be affected by mass movement erosion. The impacts of mass movement erosion on land productivity are significant but under-rated in the scientific literature. Impacts on cropping are here reported from 15 countries in south and southeast Asia, east Africa, the Caribbean and Melanesia, but accounts are generalized or anecdotal, and do not quantify crop loss or damage attributable to mass movement separately from that due to surface or fluvial erosion. Impacts on pastoral grazing have been studied in New Zealand, where production losses of up to 80% at field scale, and up to 20% at farm scale, have been measured. Studies in the Pacific Northwest coastal forests of North America show plantation forest wood volume declines by 35-50% on eroded sites. Mass movement impacts on production from tropical forests or agroforestry appear to be as yet undocumented.The reasons for lack of documentation are, first, that most soil erosion-productivity research has been done on gently sloping cropland, which is subject to surface rather than mass movement erosion. Secondly, geomorphological research in steeplands has dealt with mass movement as a hazard to human life, settlements and infrastructure -with limited identification of its contribution to sediment loads in rivers, and disregarding its impact on land productivity.We suggest there are many other countries where significant impacts are likely to occur, and that erosion-productivity studies should pay more attention to this type of erosion. Studies should not be restricted to cropland, but also extend to grazing land, plantation forestry, agro-forestry and traditional uses of natural forest as mass movement appears to affect all these forms of land-based production, particularly in densely populated steeplands whether tropical or temperate. Topics needing study are the documentation and costing of productivity losses, ways to reduce mass movement impacts on productivity, and ways to enhance recovery of soil on eroded areas (e.g., revegetation with fertility-building shrubs and legumes).

Description: Remote sensing has demonstrated wide applicability in the area of estimating and mapping forest physical and structural features. Focus in recent years has been directed towards measuring the biophysical/physiological character of forest ecosystems in order to estimate and predict forest ecosystem health and sustainability. The following reviews the relationship between forest condition and reflectance; remote-sensing measurements (and derivatives) that provide biophysical/physiological information; and the potential of hyperspectral sensors in the measurement of these parameters.

Description: This article draws attention to the potential of magnetic resonance imaging (MRI) for studies of patterns and mechanisms of water infiltration into soils, and stresses the vitally important need for collaboration between hydrologists, soil physicists and MRI experts. A brief introduction of the principles of MRI is given. This is followed by a review of the literature relating to nonpreferential infiltration, preferential infiltration exhibiting fingering and preferential infiltration involving a wide range of macropore flow. These differing degrees of complexity of infiltration dynamics require the employment of noninvasive and nondestructive techniques for their detailed investigation. Finally, an overview of applications of MRI to the detection of the spatial and temporal distribution of soil moisture and its changes is given. General conclusions are drawn from previous and current research, and the potential of the application of MRI to infiltration studies is summarized.

Description: The analysis and interpretation of remote sensing data facilitates investigation of land surface complexity on large spatial scales. We introduce here a geometrically high-resolution data set provided by the airborne High Resolution Stereo Camera (HRSC-A). The sensor records digital multispectral and panchromatic stereo bands from which a very high-resolution ground elevation model can be produced. After introducing the basic principles of the HRSC technique and data, applications of HRSC data within the multidisciplinary Research Training Group 437 are presented. Applications include geomorphologic mapping, geomorphometric analysis, mapping of surficial grain-size distribution, rock glacier kinematic analysis, vegetation monitoring and three-dimensional landform visualization. A final evaluation of the HRSC data based on three years of multipurpose usage concludes this presentation. A combination of image and elevation data opens up various possibilities for visualization and three-dimensional analysis of the land surface, especially in geomorphology. Additionally, the multispectral imagery of the HRSC data has potential for land cover mapping and vegetation monitoring. We consider HRSC data a valuable source of high-resolution terrain information with high applicability in physical geography and earth system science.

Description: Following the Holocene thermal maximum, the Greenland ice sheet began to expand, and many local glaciers formed. The maximum extent of most outlet glaciers from the Inland Ice was reached around 100 to 200 years ago. However, some exceptions have been recorded. In southern Greenland, the Neoglacial Narssarssuaq stade marks an older re-advance that exceeded the historical advance. The dating of the Narssarssuaq stade has been much debated. Here we present new radiocarbon dates from a sediment core collected from a small lake within the moraines formed during this stade. A minimum age for the Narssarssuaq stade of around 1200 yr BP is suggested from the dates.

Description: Fourteen hillslope soil profiles were sampled under natural vegetation (ie, grassland or forest) and plantations in the Nilgiri highlands, Southern India. Delta13C ratios were measured at different depths and14C ages determined for six profiles. In these highland soils where the turnover rate of organic matter is extremely low, the δ13C ratios of entire soil profiles have recorded signatures of past land cover. By correlating the data with results previously obtained from peat bogs and with knowledge concerning the history of human settlement, we distinguish three contrasting trajectories of palaeoenvironmental history and landscape change since the last glacial maximum. In the Central Nilgiris, between 18 and 10 ka BP, forest expansion occurred through the conjunction of a wetter climate (the maximum of southwest monsoon-related humidity occurring at c. 11 ka BP) and higher temperatures; since 10 ka BP, the reversal towards grassland vegetation is attributed to drier conditions. In the Western Nilgiris, where strong southwest monsoon winds permanently restrict forest patches to sheltered valley sites, steady but limited expansion of forest from 18 ka BP to the present is recorded and attributed to rising temperatures. The Southern and Eastern Nilgiris, where the northeast monsoon contributes 20% of the annual rainfall, are the less sensitive to fluctuations in the southwest monsoon. In these areas, rapid and extensive expansion of forest occurred mainly as a consequence of higher temperatures from 18 ka BP to the present. Massive deforestation by Badaga cultivators and Europeans planters followed after the sixteenth century AD.As a result, and in contrast with the Western Nilgiris where the land cover mosaic has remained remarkably stable in the last 18 ka BP, the current landscape differs sharply from the land cover pattern detected by the soil record.

Description: The paradigm of landscape ecology describes a landscape as a mosaic of landscape elements including the matrix, patches and corridors. Corridors are described as linear disruptions to the matrix, produced by anthropogenic actions or by streams which produce riparian corridors. Snow avalanches and debris flows are other geomorphic processes that should be considered as geomorphic process corridors rather than as disturbance patches. They possess requisite linearity, and they accomplish the five functions of a corridor: habitat, conduit, filter, source and sink. The definition of corridor in landscape ecology should be modified to embrace the concept of geomorphic process corridors.

Description: A characteristic of beaver ecology is their ability to build dams and, thus, to modify the landscape to increase its suitability for their occupation. This ability gives beaver great significance as a geomorphic agent. In order to review the hydrogeomorphological effects of beaver dam-building activity, this article places a context on the likely distribution and magnitude of beaver activity by considering the spatial and temporal variability of distributions of beaver and the habitat characteristics which might favour the establishment of substantial beaver populations. A description is then given of the nature and potential dimensions of instream structures built by beaver and the environmental conditions under which dam building has been observed to occur. The hydrogeomorphological impact of dam building is then appraised both locally and at the landscape scale, illustrating the very significant process modification caused by beaver. While the European beaver, Castor fiber, is the main focus of this review, it necessarily draws extensively on the much larger literature concerning the North American beaver (Castor canadensis).

Description: With global periglacial geomorphology undergoing significant advancements, it is appropriate to review the past and current status of such research in Africa. A brief historical overview of research outputs and approaches is presented for the respective African regions. Potential future quantitative periglacial research needs and approaches identified for Africa include: the examination of active periglacial processes, the identification of landforms and ground-ice forms, the potential for environmental change and the palaeoenvironmental reconstruction, and the application of periglacial studies.It is demonstrated that while periglacial geomorphology has progressed significantly in southern Africa, there has been little or no advancement elsewhere on the continent over the last two decades. None the less, on a more positive note, it is concluded that Africa has considerable potential in future global periglacial research.

Description: Remote sensing has served as an efficient method of gathering data about glaciers since its emergence. The recent advent of Geographic Information Systems (GIS) and Global Positioning Systems (GPS) has created an effective means by which the acquired data are analysed for the effective monitoring and mapping of temporal dynamics of glaciers. A large number of researchers have taken advantage of remote sensing, GIS and GPS in their studies of glaciers. These applications are comprehensively reviewed in this paper. This review shows that glacial features identifiable from aerial photographs and satellite imagery include spatial extent, transient snowline, equilibrium line elevation, accumulation and ablation zones, and differentiation of ice/snow. Digital image processing (e.g., image enhancement, spectral ratioing and automatic classification) improves the ease and accuracy of mapping these parameters. The traditional visible light/infrared remote sensing of two-dimensional glacier distribution has been extended to three-dimensional volume estimation and dynamic monitoring using radar imagery and GPS. Longitudinal variations in glacial extent have been detected from multi-temporal images in GIS. However, the detected variations have neither been explored nor modelled from environmental and topographic variables. GPS has been utilized independent of remote sensing and GIS to determine glacier ice velocity and to obtain information about glacier surfaces. Therefore, the potential afforded by the integration of nonconventional remote sensing (e.g., SAR interferometry) with GIS and GPS still remains to be realized in glaciology. The emergence of new satellite images will make remote sensing of glaciology more predictive, more global and towards longer terms.

Description: A widening variety of applications is diversifying geomorphometry ( digital terrain modelling), the quantitative study of topography. An amalgam of earth science, mathematics, engineering and computer science, the discipline has been revolutionized by the computer manipulation of gridded terrain heights, or digital elevation models (DEMs). Its rapid expansion continues. This article reviews the remarkable diversity of recent morphometric work in 15 selected topics and discusses their significance and prospects. The quantitative analysis of industrial microsurface topography is introduced to the earth science community. The 14 other topics are Internet access to geomorphometry; global DEMs; DEM modelling of channel networks; self-organized criticality; fractal and wavelet analysis; soil resources; landslide hazards; barchan dunes; harvesting wind energy; sea-ice surfaces; sea-floor abyssal hills; Japanese work in morphometry; and the emerging fields of landscape ecology and image understanding. Closing remarks note reasons for the diversity within geomorphometry, speculate on future trends and recommend creating a unified field of surface representation.

Description: Postdisaster search and rescue is an important application of ultrawideband (UWB) radar systems, which mainly detect trapped victims by their respiratory-motion response. The development of a respiratory-motion detection (RMD) algorithm that can eliminate nonstationary clutter and noise is a challenging task for the application. A new algorithm is proposed to deal with the task in this letter. It uses the multichannel singular spectrum analysis (MSSA) technique to reconstruct the respiratory-motion response detected by a UWB radar. During the reconstruction, the periodicity and range interrelation characteristics of the response are exploited to adaptively identify signal subspaces. The performance of the algorithm is verified both by simulated and real data. The results show its improved performance over the reference algorithms, e.g., a singular-value-decomposition-based algorithm. The adaptive-MSSA-based RMD algorithm has great promise not only in practical use but also for future research of UWB-radar-based human being remote sensing.

Description: Land–sea contamination observed in Soil Moisture and Ocean Salinity (SMOS) brightness temperature images is found to have two main contributions: the floor error inherent of image reconstruction and a multiplicative error either in the antenna temperature or in the visibility samples measured by the correlator. The origin of this last one is traced down to SMOS calibration parameters to yield a simple correction scheme, which is validated against several geophysical scenarios. Autoconsistency rules in interferometric synthesis together with redundant and complementary calibration procedures provide a robust SMOS calibration scheme.

Description: This letter proposes a novel algorithm, which is based on the generalized method of moments (GMM), for the estimation and correction of phase errors induced in synthetic aperture radar (SAR) imagery. The GMM algorithm is used to replace the original phase-estimation kernel in the basic structure of the phase-gradient-autofocus algorithm. Since this novel algorithm does not require the observed signal to be a certain distribution model, it is able to estimate arbitrary phase errors. The GMM algorithm has the ability of estimating range-dependent phase errors, which makes it an efficient estimator. As a result, higher accuracy of the estimated phase errors and a better focused image can be achieved. Excellent results have been obtained in autofocusing and imaging experiments on real SAR data.

Description: Strong clutter reflections of terrain and marine surfaces obscure the contrast between the target-of-interest and clutter (terrain and marine surface reflections) in synthetic aperture radar (SAR) images and consequently hinder the efficiency of image interpretation and analysis. To overcome this problem, this letter proposes an efficient clutter suppression method in SAR images, which is named shedding irrelevant patterns (SIP). The essence is to construct a regression function that can suppress clutter and preserve the target patterns concurrently. We assume that the clutter is irrelevant to the target-of-interest and distinguishable in patterns in terms of image-pixel distribution and intensity (spatial information). Experimental results show the efficiency of the proposed method in both clutter suppression and target pattern preservation.

Description: A concern in hyperspectral image classification is the high number of required training samples. When traditional classifiers are applied, feature reduction (FR) techniques are the most common approaches to deal with this problem. Subspace-based classifiers, which are developed based on high-dimensional space characteristics, are another way to handle the high dimension of hyperspectral images. In this letter, a novel subspace-based classification approach is proposed and compared with basic and improved subspace-based classifiers. The proposed classifier is also compared with traditional classifiers that are accompanied by an FR technique and the well-known support vector machine classifier. Experimental results prove the efficiency of the proposed method, especially when a limited number of training samples are available. Furthermore, the proposed method has a very high level of automation and simplicity, as it has no parameters to be set.

Description: To improve the spatial density of measurement points of persistent-scatterer interferometry, distributed scatterer (DS) should be considered and processed. An important procedure in DS interferometry is the phase triangulation (PT). This letter introduces two modified PT algorithms (i.e., equal-weighted PT and coherence-weighted PT) and analyzes the mathematical relations between different published PT methods (i.e., the maximum-likelihood phase estimator, least squares estimator, and eigendecomposition-based phase estimators). The analysis shows that the above five PT methods share very similar mathematical forms with different weight values in the estimation procedure.

Description: Sparse representation-based classifier and its variants have been widely adopted for hyperspectral image (HSI) classification recently. However, sparse representation is unstable so that similar features might obtain significantly different sparse codes. Despite the instability, we find that the sparse codes follow a class-dependent distribution under the structured dictionary consisting of training samples from all classes. Based on this observation, a novel discriminative feature, sparse code histogram (SCH), is developed for HSI classification. By counting the SCH of each sample from the sparse codes of its spatial neighbors, we can statistically obtain the distribution pattern of sparse codes of the class to which the sample belongs, and then treat the SCH as a new feature for classification. To reduce the possible outliers among the neighbors, a shape-adaptive neighborhood extractor is also employed to enhance the stability of the histogram feature. Experimental results demonstrate that SCH enjoys a strong discriminative power, which can achieve notably better performance than several state-of-the-art methods for HSI classification with limited training samples.

Description: Several detection statistics have been proposed for detecting fine ground disturbances between two synthetic aperture radar (SAR) images, such as vehicle tracks. The standard method involves estimating a local correlation coefficient between images. Other methods have been proposed using various statistical hypothesis tests. One of these alternative methods is a generalized likelihood ratio test (GLRT), which compares a full-correlation image model to a no-correlation image model. In this letter, we expand the GLRT to polarimetric SAR data and derive the appropriate GLRT detection statistics. Additionally, we explore relaxing the equal variance/equal polarimetric covariance assumptions used in previous results and find improved performance on macroscopic scene changes.

Description: For downward-looking linear array 3-D synthetic aperture radar, the resolution in cross-track direction is much lower than the ones in range and azimuth. Hence, superresolution reconstruction algorithms are desired. Since the cross-track signal to be reconstructed is sparse in the object domain, compressive sensing algorithm has been used. However, the imaging processing on the 3-D scene brings large computational loads, which renders challenges in both data acquisition and processing. To cover this shortage, truncated singular value decomposition is utilized to reconstruct a reduced-redundancy spatial measurement matrix. The proposed algorithm provides advantages in terms of computational time while maintaining the quality of the scene reconstructions. Moreover, our results on uniform linear array are generally applicable to sparse nonuniform linear array. Superresolution properties and reconstruction accuracies are demonstrated using simulations under the noise and clutter scenarios.

Description: This letter proposes a signal processing method of passive bistatic radar (PBR) exploiting an uncooperative radar as an illuminator. Compared with other opportunity illuminators, the transmitting signal of a radar usually has a better ambiguity function, which leads to a higher range resolution. Two channels are needed in PBR system. The reference channel is used to estimate radar signal parameters and reconstruct directly propagated signal. The surveillance channel is used to receive scattered wave. An array antenna and a simultaneous multibeam algorithm are necessary in the surveillance channel due to the flexible beam scanning of the uncooperative radar. The procedure of the proposed method is explained in detail, which is then followed by a field experiment. Preliminary results from the field experiment show that the proposed method can be applied to target angle and bistatic range measurement successfully.

Description: In order to achieve 3-D imaging with an airborne down-looking linear-array synthetic aperture radar (LASAR), a uniform virtual antenna array may be obtained by aperture synthesis of the cross-track sparse multiple-input–multiple-output array. However, the actual 3-D imaging quality is unavoidably degraded by errors in the virtual element position. In this letter, we investigate the effects of these errors on the forms and the degrees of image quality degradation by decomposing the error-related stochastic processes via an orthogonal transform based on discrete Legendre polynomials. It should be noted that these analyses are helpful for designing a LASAR system and providing a reference for specifying the requisite precision of measurement devices and calibration methods. Finally, we briefly consider the use of calibration methods to eliminate the effects of errors.

Description: In problems where labeled data are scarce, semisupervised learning (SSL) techniques are an attractive framework that can exploit both labeled and unlabeled data. These approaches typically rely on a smoothness assumption such that examples that are similar in input space should also be similar in label space. In many domains, such as remotely sensed hyperspectral image (HSI) classification, the data violate this assumption. In response, we propose a general method by which a neighborhood graph used in SSL is learned using binary classifiers that are trained to predict whether a pair of pixels shares the same label. Working within the framework of semisupervised neural networks (SSNNs), we show that our approach improves on the performance of the SSNN on two HSI data sets.

Description: In this letter, new models for the spatial correlation of sea clutter texture and intensity are proposed as improved versions of current power law models or exponential decay model. The models for texture have three unknown parameters, and thus can be called triparametric models. The structure of the models is a weighted sum of two components, which can describe the decaying process of the correlation coefficient with spatial lags, as well as the periodic behavior due to the existence of transient coherent structures in sea clutter. Unknown parameters are optimized by the nonlinear least square fit method. Models for sea clutter intensity can be obtained through a linear transform for uncorrelated speckle based on the compound-Gaussian representation of sea clutter. The proposed models are validated and compared with current models using S- and C-band measured sea clutter data. Analysis results indicate the effectiveness of the proposed models in that they can describe the behavior of spatial correlation coefficients with higher accuracy.

Description: Terrestrial laser scanning (TLS) has become a popular tool for acquiring source data points which can be used to construct digital elevation models (DEMs) for a wide number of applications. A TLS point cloud often has a very fine spatial resolution, which can represent well the spatial variation of a terrain surface. However, the uncertainty in DEMs created from this relatively new type of source data is not well understood, which forms the focus of this letter. TLS survey data representing four terrain surfaces of different characteristics were used to explore the effects of surface complexity and typical TLS data density (in terms of data point spacing) on DEM accuracy. The spatial variation in TLS data can be decomposed into parts corresponding to the signal of spatial variation (of terrain surfaces) and noise due to measurement error. We found a linear relation between the DEM error and the typical TLS data spacings considered (30–100 mm) which arises as a function of the interpolation error, and a constant contribution from the propagated data noise. This letter quantifies these components for each of the four surfaces considered and shows that, for the interpolation method considered here, higher density sampling would not be beneficial.

Description: A novel way to estimate the live fuel moisture content (LFMC) was explored from the ratio of canopy water content (CWC) and foliage dry biomass (FDB). The CWC was estimated using the PROSAIL (PROSPECT + SAIL) radiative transfer model from the Landsat 8 product. A weak constraint 4-D variational data assimilation method was employed to assimilate the temporally estimated leaf area index into a soil-water-atmosphere-plant (SWAP) model for optimizing the model control variables. Then, the SWAP model was reinitialized with this optimum set of control variables, and better prediction of FDB was obtained. Results showed that a high accuracy level was achieved for the estimated CWC ( $R^{2}=0.91$ , $mbox{RMSE}=84.74 mbox{g/m}^2$ ) and FDB ( $R^2=0.88$ , $mbox{RMSE}=48.54 mbox{g/m} ^2$ ) when compared with in situ measured values. However, the accuracy level of estimated LFMC was poor ( $R^2=0.59$ , $mbox{RMSE} =30.85%$ ) . Further analyses find that the estimated LFMC is reliable for low LFMC but challenged for high LFMC, which indicates that the presented method still makes sense to the assessment of wildfire risk since the wildfire generally occurs when the vegetation is in low LFMC condition.

Description: In this letter, we present the use of experimental human micro-Doppler signature data gathered by a multistatic radar system to discriminate between unarmed and potentially armed personnel walking along different trajectories. Different ways of extracting suitable features from the spectrograms of the micro-Doppler signatures are discussed, particularly empirical features such as Doppler bandwidth, periodicity, and others, and features extracted from singular value decomposition (SVD) vectors. High classification accuracy of armed versus unarmed personnel (between 90% and 97% depending on the walking trajectory of the people) can be achieved with a single SVD-based feature, in comparison with using four empirical features. The impact on classification performance of different aspect angles and the benefit of combining multistatic information is also evaluated in this letter.

Description: Automatic urban area detection in remote sensing images is an important application in the field of earth observation. Most of the existing methods employ feature classifiers and thereby contain a data training process. Moreover, some methods cannot detect urban areas in complex scenes accurately. This letter proposes an automatic urban area detection method that uses multiple features that have different resolutions. First, a downsampled low-resolution image is used to segment the candidate area. After the corner points of the urban area are extracted, a weighted Gaussian voting matrix technique is employed to integrate the corner points into the candidate area. Then, the edge features and homogeneous region are extracted by using the original high-resolution image. Using these results as the input, the processes of guided filtering and contrast enhancement can finally detect accurately the urban areas. This method combines multiple features, such as corner, edge, and regional characteristics, to detect the urban areas. The experimental results show that the proposed method has better detection accuracy for urban areas than the existing algorithms.

Description: In marine sciences, time series are often nonlinear and nonstationary. Adequate and specific methods are needed to analyze such series. In this letter, an application of the empirical mode decomposition method (EMD) associated to the Hilbert spectral analysis (HSA) is presented. Furthermore, EMD-based time-dependent intrinsic correlation (TDIC) analysis is applied to consider the correlation between two nonstationary time series. Four temperature time series obtained from automatic measurements in nearshore waters of the Réunion island are considered, recorded every 10 min from July 2011 to January 2012. The application of the EMD on these series and the estimation of their power spectra using the HSA are illustrated. The authors identify low-frequency tidal waves and display the pattern of correlations at different scales and different locations. By TDIC analysis, it was concluded that the high-frequency modes have small correlation, whereas the trends are perfectly correlated.

Description: Detecting vehicles in aerial images provides important information for traffic management and urban planning. Detecting the cars in the images is challenging due to the relatively small size of the target objects and the complex background in man-made areas. It is particularly challenging if the goal is near-real-time detection, i.e., within few seconds, on large images without any additional information, e.g., road database and accurate target size. We present a method that can detect the vehicles on a 21-MPixel original frame image without accurate scale information within seconds on a laptop single threaded. In addition to the bounding box of the vehicles, we extract also orientation and type (car/truck) information. First, we apply a fast binary detector using integral channel features in a soft-cascade structure. In the next step, we apply a multiclass classifier on the output of the binary detector, which gives the orientation and type of the vehicles. We evaluate our method on a challenging data set of original aerial images over Munich and a data set captured from an unmanned aerial vehicle (UAV).

Description: We compare five slope correction methods developed by Walter et al. , Montes et al. , Schleppi et al. , España et al. , and Gonsamo et al. (referred to as WAL, MON, SCH, ESP, and GON, respectively) using artificial fisheye pictures simulated by graphics software and a lookup table (LUT) retrieval method. The LUT is built by simulating the directional gap fraction as a function of leaf area index (LAI) and average leaf inclination angle (ALIA) using the Poisson law. LAI and ALIA estimates correspond to the case of the LUT that provides the lowest root-mean-square error between the observed gap fractions after slope correction and the simulated ones. Three LAI values (1.5, 3.5, and 5.5), four ALIA values (26.8°, 45°, 57.5°, and 63.2°), and three slope angles (0°, 20°, and 50°) constituted 36 samples of random scenes. ESP is recommended because its results are accurate and independent on the leaf angle distribution (LAD), while GON only performs well for spherical LAD. The three other methods present less good performances with underestimation or overestimation of LAI and/or ALIA depending on the LAD, and the recommended order for them is MON, SCH, and WAL.

Description: In this letter, an improved phase correlation (PC) method based on 2-D plane fitting and the maximum kernel density estimator (MKDE) is proposed, which combines the idea of Stone's method and robust estimator MKDE. The proposed PC method first utilizes a vector filter to minimize the noise errors of the phase angle matrix and then unwraps the filtered phase angle matrix by the use of the minimum cost network flow unwrapping algorithm. Afterward, the unwrapped phase angle matrix is robustly fitted via MKDE, and the slope coefficients of the 2-D plane indicate the subpixel shifts between images. The experiments revealed that the improved method can effectively avoid the impact of outliers on the phase angle matrix during the plane fitting and is robust to aliasing and noise. The matching accuracy can reach 1/50th of a pixel using simulated data. The real image sequence tracking experiment was also undertaken to demonstrate the effectiveness of the proposed PC method with a registration accuracy of root-mean-square error better than 0.1 pixels.

Description: Accurately mapping forest carbon density by combining sample plots and remotely sensed images has become popular because this method provides spatially explicit estimates. However, mixed pixels often impede the improvement of the estimation. In this letter, regression modeling and spectral unmixing analysis were integrated to improve the estimation of forest carbon density for the You County of Hunan, China, using Landsat Thematic Mapper images. Linear spectral unmixing with and without a constraint (LSUWC and LSUWOC) and nonlinear spectral unmixing (NSU) were compared to derive the fractions of five endmembers, particularly forests. Stepwise regression, logistic regression, and polynomial regression (PR) with and without the forest fraction used as an independent variable and the product of the forest fraction image and the map from the best model without the forest fraction were compared. The models were developed using 56 sample plots, and their results were validated using 26 test plots. The decomposition of mixed pixels was assessed using higher spatial resolution SPOT images and a corresponding land cover map. The results showed that 1) LSUWC more accurately estimated the endmember fractions than LSUWOC and NSU, 2) PR had the greatest estimation accuracy of forest carbon, and 3) combining regression modeling and spectral unmixing increased the estimation accuracy by 31%–39%, and introducing the forest fraction into the regressions performed better than the product of forest fraction image and the results from PR without the fraction. This implied that the integrations provided great potential in reducing the impacts of mixed pixels in mapping forest carbon.

Description: This letter proposes a multiresolution technique to address the high computational cost in remote sensing image registration. The scale-invariant feature transform is applied to detect keypoints and descriptors, and then, global information combined with descriptors is utilized to establish keypoint mappings. Keypoints are first classified according to their octaves. Then, in the lowest resolution, the keypoints of the largest octave are mapped with descriptors and the global information, giving an initial affine transformation $T_0$ . In the next octave, the keypoints of the second largest octave are mapped by employing $T_0$ to narrow the space of matching keypoints. By this means, the process of establishing keypoint correspondences is conducted from one resolution (octave) to the next as the obtained transformation gets finer until we get to the highest resolution. Due to the high computational expense of computing global information, the proposed technique is important for aligning large-size remote sensing imagery. Experimental results show that the proposed method can achieve a comparable registration accuracy but with a less computational cost than directly building keypoint mappings on images of large size.

Description: Presents corrections to the paper, "Estimation of forest biomass from two-level model inversion of single-pass InSAR data" (Soja, M.J., et al.,Trans. Geosci. Remote Sens., vol. 53, no. 9, pp. 5083???5099, Sep. 2015).

Description: In this letter, we present an efficient parallel implementation of composite kernels in support vector machines (SVMs) for hyperspectral image (HSI) classification. Our implementation makes effective use of commodity graphics processing units (GPUs). Specifically, we port the calculation of composite kernels to GPUs, perform intensive computations based on NVidia's compute unified device architecture, and execute the rest of the operations related with control and small data calculations in the CPU. Our experimental results, conducted using real hyperspectral data sets and NVidia GPU platforms, indicate significant improvements in terms of computational effectiveness, achieving near-real-time performance of spatial–spectral HSI classification for the first time in the literature.

Description: Compared with airborne laser scanning, terrestrial laser scanning (TLS) offers ground-based point cloud data of trees and provides greater potential to accurately estimate tree and stand parameters. However, there is a lack of effective methods to accurately identify locations of individual trees from TLS point cloud data. It is also unknown whether the estimation accuracy of the parameters, including tree height (H), diameter at breast height (DBH), and so on, using TLS can meet the requirement of forest management and planning. In this letter, a novel method to effectively process point cloud data and further determine the locations of individual trees in a stand based on the central coordinates of point cloud data on a defined grid according to the largest DBH was developed. Moreover, a point-cloud-data-based convex hull algorithm and the cylinder method were, respectively, used to estimate DBH and H of individual trees. This study was conducted in a pure Chinese fir plantation of 45 trees located in Huang-Feng-Qiao forest farm, You County of Hunan, China. The comparison of the estimated and observed values showed that the obtained tree locations had errors of less than 20 cm, and the relative root mean square errors for the estimates of both DBH and H were less than 5%. This implies that TLS is very promising for the retrieval of tree and stand parameters in forest stands. For the applications of these methods to mixed forests with a structure of multilayer canopies, further examination is needed.

Description: In this letter, a novel algorithm for attitude measurement based on a 3-D electromagnetic model (3-D em-model) is proposed. The 3-D em-model is established offline based on the geometric structure of the target, and it can be used to predict the scattering features at different target attitudes. In order to measure the attitude of the air target, we design a bistatic step frequency radar system. The directions of the two radars' lines of sight (LOSs) relative to the target are acquired by matching the high-resolution range profiles (HRRPs) from the target echoes to the HRRPs generated from the 3-D em-model. Since the directions of two radars' LOSs relative to the Earth are already known, the absolute attitude of the target can be acquired. The innovative contributions of this letter are as follows: 1) A comprehensive theoretical analysis of air target attitude measurement based on its own 3-D em-model is proposed; 2) the method can be applied to different kinds of air targets such as aircraft, satellite, missile, etc.; 3) the proposed attitude measurement method does not require target motion model in advance; and 4) the proposed algorithm can be applied to any kind of step frequency waveforms. Experiments using both data predicted by a high-frequency electromagnetic code and data measured in the chamber verify the validity of the method.

Description: The Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) missions provide Level-1 brightness temperature (Tb) observations that are used for global soil moisture estimation. However, the nature of these Tb data differs: the SMOS Tb observations contain atmospheric and select reflected extraterrestrial (“Sky”) radiation, whereas the SMAP Tb data are corrected for these contributions, using auxiliary near-surface information. Furthermore, the SMOS Tb observations are multiangular, whereas the SMAP Tb is measured at 40° incidence angle only. This letter discusses how SMOS Tb, SMAP Tb, and radiative transfer modeling components can be aligned in order to enable a seamless exchange of SMOS and SMAP Tb data in soil moisture retrieval and assimilation systems. The aggregated contribution of the atmospheric and reflected Sky radiation is, on average, about 1 K for horizontally polarized Tb and 0.5 K for vertically polarized Tb at 40° incidence angle, but local and short-term values regularly exceed 5 K.

Description: The temporal variations (diurnal and annual) in arboreal $(varepsilon_mathrm{Tree})$ and bare soil $(varepsilon_mathrm{Soil})$ dielectric constants and their correlation with precipitation were examined for several trees in Japan. A significant (1 $sigma$ (standard deviation) and 2 $sigma$ ) $varepsilon_mathrm{Tree}$ increase is observed after rainfall at 89.8% and 90.5% probability. However, rainfall does not always induce significant $varepsilon_mathrm{Tree}$ increases. Rainfall of more than 5 mm/day can induce 1 $sigma$ $varepsilon_mathrm{Tree}$ Tree increase at a 59.6% probability. In order to examine whether the increase in $varepsilon_mathrm{Tree}$ affects the L-band $sigma^{0}$ variation in a forest, the four-year temporal variation of the L-band backscattering coefficient $(sigma^{0})$ was estimated from observations by the Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar. Observed maximum absolute deviations from the mean over the forest area were 1.0 and 1.2 dB for $sigma_{mathrm{HH}}^{0}$ and $sigma_mathrm{HV}^{0}$ , respectively, and 4.0 and 3.0 dB over open land. $sigma^{0}$ and rainfall correlations show that $varepsilon_mathrm{Tree}$ and $sigma_mathrm{Forest}^{0}$ are proportional to precipitation integrated over seven or eight days; $varepsilon_mathrm{Soil}$ and $sigma_mathrm{Open land}^{0}$ are proportional to precipitation integrated over three days. This finding indicates that $varepsilon_mathrm{Tree}$ variations influence $sigma_{mathrm{Forest areas}}^{0}$ . A stronger correlation between $sigma_mathrm{HV}^{0}$ and precipitation is observed in several sites with low $sigma_mathrm{HV}^{0}$ , where less biomass is expected, and several sites with high $sigma_mathrm{HV}^{0}$ , where more biomass is expected. A weaker correlation between $sigma_mathrm{HV}^{0}$ and precipitation is observed for several sites with high $sigma_mathrm{HV}^{0}$ . These differences may be explained by the different contributions of double bounce scattering and potential transpiration, which is a measure of the ability of the atmosphere to remove water from the surface th

Description: The spectral and radiometric quality of airborne imaging spectrometer data is affected by the anisotropic reflectance behavior of the imaged surface. Illumination and observation angle-dependent patterns of surface reflected radiation propagate into products, hinder quantitative assessment of biophysical/biochemical parameters, and decrease the comparability of data from multiple flight lines. The Ross–Li model, originally developed for multiangular observations, can be inverted to estimate and correct for surface anisotropy effects. This requires land cover be stratified into distinct types of scattering behavior. When the observations subsumed in these classes cover a range of view angles, a pseudo multiangular view on the surface can be employed to invert the Ross–Li model. A discrete land cover classification, however, bears the risk of inappropriate scattering correction resulting in spatial artifacts in the corrected data, predominantly in transition regions of two land cover types (e.g., soil and sparse vegetation with varying fractions). We invert the Ross–Li model on continuous land cover fraction layers. We decompose land cover in dominating structural types using linear spectral unmixing. Ross–Li kernel weights and formulations are estimated for each type independently; the correction is then applied pixel-wise according to the fractional distribution. The corrected Airborne Prism EXperiment imaging spectrometer data show significant reduction of anisotropic reflectance effects of up to 90% (average 60% to 75%, $p=0.05$ ), measured in the overlapping regions of adjacent flight lines. No spatial artifacts or spectral irregularities are observed after correction.

Description: Most spaceborne sensors have a tradeoff between high spatial and high temporal resolutions. This tradeoff limits the use of remote sensing data in various applications that require images in both the high spatial and high temporal resolutions. In this paper, we propose a novel technique to create a fine spatial and high temporal resolution images at a ground-based data processing system. Resourcesat-2 is one of the Indian Space Research Organization missions, and it carries the Linear Imaging Self-Scanning sensors (LISS III and LISS IV) and an Advanced Wide-Field Sensor (AWiFS) onboard. The spatial resolution of LISS III is 23.5 m, and that of AWiFS is 56 m. The temporal resolution of LISS III is 24 days, and that of AWiFS is five days. The proposed method creates a synthetic LISS-III image at 23.5-m spatial and five-day temporal resolutions. It is based on the subpixel relationship between a single AWiFS–LISS-III image pair, which is acquired before or after the prediction date. In temporal data composition, spurious spatial discontinuities are inevitable for land-cover type changes. These discontinuities were identified with temporal edge primitives and were smoothed with a spatial-profile-averaging method. A synthetic LISS-III image for time $t_{k}$ is predicted from an AWiFS image at time $t_{k}$ and a single AWiFS–LISS-III image pair at time $t_{0}$ , where $t_{0}ne t_{k}$ . Experimental results demonstrated that the proposed method is superior in terms of the computational efficiency and prediction accuracy with the other existing methods.

Description: Glacier avalanches are natural hazards that could damage infrastructures and threaten lives in high-altitude mountainous terrains. On April 7, 2012, a massive ice avalanche struck a Pakistani base at Gayari sector, Saltoro Valley, and buried/killed 148 soldiers and civilians. Keeping in view the catastrophe, a study was designed with the objectives to: 1) model and simulate the Gayari sector glacier avalanche incident in terms of height, extent, velocity, pressure, and momentum; 2) generate hazard risk assessment of possible other glacier avalanches in the Saltoro Valley through modeling and simulation; and 3) suitability analysis of current camp sites and recommendation of new safe camps sites locations in the Saltoro valley. To simulate the Gayari sector glacier event and other Glacier possible avalanches, a physical process-based rapid mass movements (RAMMS) was used. The RAMMS has two main components, i.e., Voellmy–Salm (VS) model and random kinetic energy, which deals with variables such as avalanche height and the mean avalanche velocity during the course of simulation. The suitability analysis of current camp sites were achieved using weighted overlay analysis with different constraints in ArcGIS Spatial Analyst. The RAMMS model simulation of the Gayari avalanche event predicted a maximum velocity of 74 ms -1 , generating a pressure of 5074 kPa and attaining a height of 45 m, whereas the predicted debris volume on the ground was 3.8145 million m 3 . A good agreement was found between actual debris height and extent, as compared with the RAMMS model output. The RAMMS model simulated all the potential tributary glacier avalanches of Saltoro valley very well. The weighted overlay analyses in ArcGIS Spatial Analyst revealed that the existing camp sites are safe and were not threatened by the glacier avalanche hazard. However, it was recommended that the Gayari camp should not be constructed at the same location and should be relocated- to the proposed safe camp sites identified in this research study. The proposed methodology developed in the current study could be applied in the Siachen conflict zone for avalanche hazard/risk analysis of all the camp sites located in the valley.

Description: Object-oriented remote sensing image classification is becoming more and more popular because it can integrate spatial information from neighboring regions of different shapes and sizes into the classification procedure to improve the mapping accuracy. However, object identification itself is difficult and challenging. Superpixels, which are groups of spatially connected similar pixels, have the scale between the pixel level and the object level and can be generated from oversegmentation. In this paper, we establish a new classification framework using a superpixel-based graphical model. Superpixels instead of pixels are applied as the basic unit to the graphical model to capture the contextual information and the spatial dependence between the superpixels. The advantage of this treatment is that it makes the classification less sensitive to noise and segmentation scale. The contribution of this paper is the application of a graphical model to remote sensing image semantic segmentation. It is threefold. 1) Gradient fusion is applied to multispectral images before the watershed segmentation algorithm is used for superpixel generation. 2) A probabilistic fusion method is designed to derive node potential in the superpixel-based graphical model to address the problem of insufficient training samples at the superpixel level. 3) A boundary penalty between the superpixels is introduced in the edge potential evaluation. Experiments on three real data sets were conducted. The results show that the proposed method performs better than the related state-of-the-art methods tested.

Description: We introduce a practical and accurate model, referred to as “GO4,” to describe near-nadir microwave scattering from the sea surface, and at the same time, we address the issue of the filtered mean square slope (mss) conventionally used in the geometrical optics model. GO4 is a simple correction of this last model, taking into account the diffraction correction induced by the rough surface through what we call an effective mean square curvature (msc). We evaluate the effective msc as a function of the surface wavenumber spectrum and the radar frequency and show that GO4 reaches the same accuracy as the physical optics model in a wide range of incidence and frequency bands with the sole knowledge of the mss and msc parameters. The key point is that the mss entering in GO4 is not the filtered but the total slope. We provide estimation of the effective msc on the basis of classical sea spectrum models. We also evaluate the effective msc from near-nadir satellite data in various bands and show that it is consistent with model predictions. Non-Gaussian effects are discussed and shown to be incorporated in the effective msc. We give some applications of the method, namely, the estimation of the total sea surface mss and the recalibration of relative radar cross sections.

Description: This paper evaluates the radiometric accuracy of observations from the Advanced Technology Microwave Sounder (ATMS) onboard Suomi National Polar-orbiting Partnership and Sondeur Atmospherique du Profil d' Humidité Intropicale par Radiométrie (SAPHIR) onboard Megha-Tropiques through intercalibration and validation versus in situ radiosonde and Global Positioning System Radio Occultation (GPS-RO) observations. SAPHIR and ATMS water vapor channels operate at slightly different frequencies. We calculated the bias due to radiometric errors as the difference between the observed and simulated differences between the two instruments. This difference, which is often referred to as double difference, ranges between 0.3 and 0.7 K, which shows good consistency between the instruments. We used a radiative transfer model to simulate the satellite brightness temperatures (Tbs) using radiosonde and GPS-RO profiles and then compared simulated and observed Tbs. The difference between radiosonde and ATMS Tbs for the middle and upper tropospheric temperature sounding channels was less than 0.5 K at most stations, but the difference between radiosonde and ATMS/SAPHIR Tbs for water vapor channels was between 0.5 and 2.0 K. The larger bias for the water vapor channels is mainly due to several errors in radiosonde humidity observations. The mean differences between the ATMS observations and the Tbs simulated using GPS-RO profiles were 0.2, 0.3, 0.4, 0.2, and −0.2 K for channels 10–14, respectively; and the uncertainty increases from 0.02 K for channel 10 to 0.07 K for channel 14.

Description: This paper presents a completely automatic processing chain for orthorectification of optical pushbroom sensors. The procedure is robust and works without manual intervention from raw satellite image to orthoimage. It is modularly divided in four main steps: metadata extraction, automatic ground control point (GCP) extraction, geometric modeling, and orthorectification. The GCP extraction step uses georeferenced vector roads as a reference and produces a file with a list of points and their accuracy estimation. The physical geometric model is based on collinearity equations and works with sensor-corrected (level 1) optical satellite images. It models the sensor position and attitude with second-order piecewise polynomials depending on the acquisition time. The exterior orientation parameters are estimated in a least squares adjustment, employing random sample consensus and robust estimation algorithms for the removal of erroneous points and fine-tuning of the results. The images are finally orthorectified using a digital elevation model and positioned in a national coordinate system. The usability of the method is presented by testing three RapidEye images of regions with different terrain configurations. Several tests were carried out to verify the efficiency of the procedure and to make it more robust. Using the geometric model, subpixel accuracy on independent check points was achieved, and positional accuracy of orthoimages was around one pixel. The proposed procedure is general and can be easily adapted to various sensors.

Description: This paper presents a new ground moving target imaging (GMTIm) algorithm for airborne synthetic aperture radar (SAR) based on a novel time-frequency representation (TFR), Lv's distribution (LVD). We first analyze generic moving target signatures for a multichannel SAR and then derive the analytical spectrum of a point target moving at a constant velocity by a polar format algorithm for SAR image formation. SAR motion deviation from a predetermined flight track is considered to facilitate airborne SAR applications. LVD, as a recently developed TFR for the analysis of multicomponent linear-frequency-modulated signal, is adopted to represent the target kinematic spectrum in the Doppler centroid frequency and chirp rate domain. As a result, the proposed SAR-GMTIm algorithm is capable of imaging multiple moving targets even when they are located at the same range resolution cell. Some practical issues such as imaging maneuvering targets and small/weak targets are discussed to enhance the applicability of the proposed algorithm. Simulation results with isotropic point moving targets are presented to validate the effectiveness and superiority of the proposed algorithm. Raw data collected by an airborne multichannel SAR are also used to verify the performance improvement made by the proposed algorithm.

Description: Satellite soil moisture estimates have received increasing attention over the past decade. This paper examines the applicability of estimating soil moisture states and soil hydraulic parameters through two particle filter (PF) methods: The PF with commonly used sampling importance resampling (PF-SIR) and the PF with recently developed Markov chain Monte Carlo sampling (PF-MCMC) methods. In a synthetic experiment, the potential of assimilating remotely sensed near-surface soil moisture measurements into a 1-D mechanistic soil water model (HYDRUS-1D) using both the PF-SIR and PF-MCMC algorithms is analyzed. The effects of satellite temporal resolution and accuracy, soil type, and ensemble size on the assimilation of soil moisture are analyzed. In a real data experiment, we first validate the Advanced Microwave Scanning Radiometer--Earth Observing System (AMSR-E) soil moisture products in the Oklahoma Little Washita Watershed. Aside from rescaling the remotely sensed soil moisture, a bias correction algorithm is implemented to correct the deep soil moisture estimate. Both the ascending and descending AMSR-E soil moisture data are assimilated into the HYDRUS-1D model. The synthetic assimilation results indicated that, whereas both updating schemes showed the ability to correct the soil moisture state and estimate hydraulic parameters, the PF-MCMC scheme is consistently more accurate than PR-SIR. For real data case, the quality of remotely sensed soil moisture impacts the benefits of their assimilation into the model. The PF-MCMC scheme brought marginal gains than the open-loop simulation in RMSE at both surface and root-zone soil layer, whereas the PF-SIR scheme degraded the open-loop simulation.

Description: We propose a new deterministic approach for remote sensing retrieval, called modified total least squares (MTLS), built upon the total least squares (TLS) technique. MTLS implicitly determines the optimal regularization strength to be applied to the normal equation first-order Newtonian retrieval using all of the noise terms embedded in the residual vector. The TLS technique does not include any constraint to prevent noise enhancement in the state space parameters from the existing noise in measurement space for an inversion with an ill-conditioned Jacobian. To stabilize the noise propagation into parameter space, we introduce an additional empirically derived regularization proportional to the logarithm of the condition number of the Jacobian and inversely proportional to the L2-norm of the residual vector. The derivation, operational advantages and use of the MTLS method are demonstrated by retrieving sea surface temperature from GOES-13 satellite measurements. An analytic equation is derived for the total retrieval error, and is shown to agree well with the observed error. This can also serve as a quality indicator for pixel-level retrievals. We also introduce additional tests from the MTLS solutions to identify contaminated pixels due to residual clouds, error in the water vapor profile and aerosols. Comparison of the performances of our new and other methods, namely, optimal estimation and regression-based retrieval, is performed to understand the relative prospects and problems associated with these methods. This was done using operational match-ups for 42 months of data, and demonstrates a relatively superior temporally consistent performance of the MTLS technique.

Description: Radar-based remote sensing for measurement of ocean surface waves presents advantages over conventional point sensors such as wave buoys. As its use becomes more widespread, it is important to understand the sensitivity of the extracted wave parameters to the characteristics of the radar and the scatterers. To examine such issues, experiments were performed offshore of the Scripps Institution of Oceanography pier in July 2010. Radar measurements in low wind speeds were performed with a dual-polarized high-resolution X-band pulse-Doppler radar at low grazing angles along with two independent measurements of the surface waves using conventional sensors, a GPS-based buoy, and an ultrasonic array. Comparison between radar cross section (RCS) and Doppler modulations shows peak values occurring nearly in-phase, in contrast with tilt modulation theory. Spectral comparisons between Doppler-based and RCS-based spectra show that Doppler-based spectra demonstrate greater sensitivity to swell-induced modulations, whereas RCS-based spectra show greater sensitivity to small-scale modulations (or generally have more noise at high frequency), and they equally capture energy at the wind wave peak. Doppler estimates of peak period were consistent with the conventional sensors, whereas the RCS differed in assignment of peak period to wind seas rather than swell in a couple of cases. Higher order period statistics of both RCS and Doppler were consistent with the conventional sensors. Radar-based significant wave heights are lower than buoy-based values and contain nontrivial variability of ∼33%. Comparisons between HH and VV polarization data show that VV data more accurately represent the wave field, particularly as the wind speeds decrease.

Description: Polarimetric synthetic aperture radar satellite and ground-based Ku- and X-band scatterometer measurements are used to explore the scattering mechanism for ice in shallow Arctic lakes, wherein strong radiometric responses are seen for floating ice, and low returns are evident where the ice has grounded. Scatterometer measurements confirm that high backscatter is from the ice/water interface, whereas polarimetric decomposition suggests that the dominant scattering mechanism from that interface is single bounce. Using Fresnel equations, a simple model for surface bounce from the ice/water interface is proposed, and its predictions are supported by experimental parameters such as co-pol phase difference, co-pol ratio, and the results of rigorous numerical modeling. Despite early research suggesting double-bounce scattering from columnar air bubbles and the ice/water interface as the dominant scattering mechanism in shallow lakes, this paper strongly supports a single-bounce model.

Description: The on-orbit radiometric calibration of the reflective solar bands (RSBs) of the Visible Infrared Imaging Radiometer Suite (VIIRS) aboard the Suomi National Polar-orbiting Partnership satellite is carried out primarily through observations of a fully illuminated solar diffuser (SD) panel. Accurate knowledge of the solar spectral radiance scattered from the SD is available. The sensor aperture spectral radiance is assumed to be a quadratic polynomial function of a VIIRS detector's background-subtracted response in digital number. The coefficients of the polynomial were initially determined prelaunch. Once on orbit, we assume that these coefficients change uniformly by a common calibration factor, which is referred to as the $F$ -factor. The known solar spectral radiance scattered from a fully illuminated SD allows for the determination of these $F$ -factors. We describe the methodology and the associated algorithms used in the calculation of the RSB $F$ -factors. Our results show that the $F$ -factors change over time, with the largest change occurring at a wavelength of 862 nm (with a value of about 1.55 on day 950 after the satellite launch, relative to its value at the beginning of the launch) . In addition, we estimate the relative error standard deviations of the computed top-of-the-atmosphere reflectance at the detector pixel level. On day 950 of the mission, the relative error standard deviations are all less or equal to 0.016, except for the M11 band (band central wavelength of 2257 nm) , which has a relative error standard deviation of about 0.049 due to a very low signal-to-noise ratio.

Description: We present simplified expressions for the cross-polarized backscatter of a randomly rough surface predicted by the second-order small-slope approximation (SSA2). The simplification is based on appropriate polynomial approximations of the SSA2 kernel function. We obtain numerically efficient expressions for the cross-polarized backscattering amplitude of a deterministic surface in the form of a single space integral involving only the surface elevation and the second (mixed) derivative of the surface elevation. The ensemble average normalized radar cross section is then derived under a Gaussian random process assumption for the surface. The resulting expression has the form of a Kirchhoff integral involving the roughness correlation function and its second- and fourth-order cross-derivatives. Further simplification is achieved for off-nadir observations using a high-frequency approximation; the result is an analytical formula involving only the resonant curvature and the radar-filtered mean square slope in the out-of-plane direction. A numerical validation of the simplified expressions is provided by comparison with exact SSA2 predictions in representative test cases. The dependence of cross-polarized backscattering on the incidence angle as well as wind speed and direction is then investigated for the case of a directional sea surface model. At near nadir incidence, a clear maximum in azimuth of the cross-polarized backscatter is observed for radar look directions 45 $^{circ}$ from the wind direction.

Description: How to produce the difference data of the two temporal images is a crucial factor in image change detection. In this letter, we propose multicontextual mutual information data (MMID) based on the bivariate Gaussian distribution (BGD) for synthetic aperture radar (SAR) image change detection and illustrate their superiorities over the classical difference data. MMID, which are an improved form of image spatial mutual information, are constructed based on the quadrilateral Markov random field (QMRF) and can be factored into the linear combination of the entropies. Then to adapt MMID to the change detection, we construct the 2-D entropies based on the BGD. In this way, MMID are able to capture the intertemporal statistical dependence of the two temporal images and thus can be taken as the feature-level difference data rather than the pixel-level data. The maximum-likelihood method, the automatic threshold method, and the Markov random field method are performed on the MMID of the real two temporal SAR images for the change detection. Experimental results demonstrate the superiorities of MMID over the traditional difference data.

Description: Spectral unmixing has been a popular technique for analyzing remotely sensed hyperspectral images. The goal of unmixing is to find a collection of pure spectral constituents (called endmembers ) that can explain each (possibly mixed) pixel of the scene as a combination of endmembers, weighted by their coverage fractions in the pixel or abundances . Over the last years, many algorithms have been presented to address the three main parts of the spectral unmixing chain: 1) estimation of the number of endmembers; 2) identification of the endmember signatures; and 3) estimation of the per-pixel fractional abundances. However, to date, there is no standardized tool that integrates these algorithms in a unified framework. In this letter, we present HyperMix, an open-source tool for spectral unmixing that integrates different approaches for spectral unmixing and allows building unmixing chains in graphical fashion, so that the end-user can define one or several spectral unmixing chains in fully configurable mode. HyperMix provides efficient implementations of most of the algorithms used for spectral unmixing, so that the tool automatically recognizes if the computer has a graphics processing unit (GPU) available and optimizes the execution of these algorithms in the GPU. This allows for the execution of spectral unmixing chains on large hyperspectral scenes in computationally efficient fashion. The tool is available online from http://hypercomphypermix.blogspot.com.es and has been validated with real hyperspectral scenes, providing state-of-the-art unmixing results.

Description: Seismic signals are nonlinear, and the seismic state-space model can be described as a nonlinear system. The particle filter (PF) method, as an effective method for estimating the state of a nonlinear system, can be applied to deal with seismic random noise attenuation. However, PF suffers from sample impoverishment caused by resampling, which results in serious loss of valid seismic information and leads to inaccurate representation of the reflected signal. To address the impoverishment issue and to further improve the particle quality, we propose a novel method to suppress seismic random noise—the adaptive fission particle filter (AFPF). In AFPF, all the particles undergo a fission process and produce “offspring” particles to maintain particle diversity. To implement the adaptation and to monitor the degree of fission, we apply a fission factor, which takes into account weights that indicate the quality of the particles. This leads to significant improvements in the particle quality, i.e., the proportion of highly weighted particles is increased. The effective seismic information provided by the resulting particles reproduces the true signal more reliably, reducing the bias of PF. In addition, we establish a dynamic state-space model suitable for seismic signals. Experimental results on synthetic records and field data illustrate the superior performance of AFPF in noise attenuation and reflected signal preservation compared with the PF.

Description: The synoptic determination of ocean circulation using the data acquired from space, with a coherent depiction of its turbulent characteristics, remains a fundamental challenge in oceanography. This determination has the potential of revealing all aspects of the ocean dynamic variability on a wide range of spatiotemporal scales and will enhance our understanding of ocean–atmosphere exchanges at superresolution, as required in the present context of climate change. Here, we show a four-year time series of spatial superresolution (4 km) turbulent ocean dynamics generated from satellite data using emerging ideas in signal processing coming from nonlinear physics, low-resolution dynamics, and superresolution oceanic sea surface temperature data acquired from optical sensors. The method at its core consists in propagating across the scales the low-resolution dynamics in a multiresolution analysis computed on adimensional critical transition information.

Description: This paper presents the theory, algorithm, and results of a new bistatic interferometry synthetic aperture radar (InSAR) method. It employs the data acquired in an innovative bistatic configuration, which uses the orbital sensors as transmitters of opportunity and the stationary receivers on the ground, to generate a digital elevation model (DEM). In the bistatic spaceborne/stationary InSAR configuration, the interferometric phase only depends on the target-receiver range, which could not be obtained directly from the measured bistatic range. Therefore, the conventional transforming relationship between the interferometric phase and the topographic height is no longer practical. In order to solve the problem, we introduce a new conversion relationship between the interferometric phase and the topographic height, which is derived by the model of the ellipsoidal projection in the bistatic configuration. Meanwhile, the error analysis of the new conversion is carried out through a simulation. Both the simulated and measured data are used to test and verify the feasibility of the new bistatic InSAR method. In the spaceborne/stationary InSAR experiment, YaoGan-3 (an L-band spaceborne SAR system launched by China) was selected as the transmitter and two stationary receivers were mounted on the top of a tall building. The generated DEM of high quality shows that the presented method performs very well in the bistatic InSAR data process.

Description: We present a nonparametric Bayesian hierarchical model (HDP_IBPs) to classify very high resolution panchromatic satellite images in an unsupervised way, in which the hierarchical Dirichlet process (HDP) and Indian buffet process (IBP) are combined on multiple scenes. The main contribution of this paper is a novel application framework to solve the problems of traditional probabilistic topic models and achieve the effective unsupervised classification of very high resolution (VHR) panchromatic satellite images. In this framework, a VHR satellite image is first oversegmented into basic processing units and divided into a set of subimages. We use the Chinese restaurant franchise process as a construct method of the HDP to capture the latent semantic structures (i.e., classes) and the class proportion (i.e., co-occurrence of topics) for each subimage. Meanwhile, the subimages are grouped into different scenes based on the class proportion. Finally, the IBP is employed to select the most appropriate classes for each subimage from all of the classes based on different scenes in turn. The hierarchical structure of our model transmits the spatial information from the original image to the scene layer implicitly and provides useful cues of classification. The experimental results show that HDP_IBPs outperforms state-of-the-art models in terms of both qualitative and quantitative evaluations.

Description: This paper describes a linear-image-transform-based algorithm for reducing stripe noise, track line artifacts, and motion-induced errors in remote sensing data. Developed for multibeam bathymetry (MB), the method has also been used for removing scalloping in synthetic aperture radar images. The proposed image transform is the composition of an invertible edge detection operator and a fast discrete Radon transform (DRT) due to Götz, Druckmüller, and Brady. The inverse DRT is computed by using an iterative method and exploiting an approximate inverse algorithm due to Press. The edge operator is implemented by circular convolution with a Laplacian point spread function modified to render the operator invertible. In the transformed image, linear discontinuities appear as high-intensity spots, which may be reset to zero. In MB data, a second noise signature is linked to motion-induced errors. A Chebyshev approximation of the original image is subtracted before applying the transform, and added back to the denoised image; this is necessary to avoid boundary effects. It is possible to process data faster and suppress motion-induced noise further by filtering images in nonoverlapping blocks using a matrix representation for the inverse DRT. Processed test images from several MB data sets had less noise and distortion compared with those obtained with standard low-pass filters. Denoising also improved the accuracy in statistical classification of geomorphological type by 10–28% for two sets of invariant terrain features.