Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial patterns”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Characterization of Aerosol Episodes in the Greater Mediterranean Sea Area from Satellite Observations (2000-2007)

An algorithm able to identify and characterize episodes of different aerosol types above sea surfaces of the greater Mediterranean basin (GMB), including the Black Sea and the Atlantic Ocean off the coasts of Iberia and northwest Africa, is presented in this study. Based on this algorithm, five types of intense (strong and extreme) aerosol episodes in the GMB are identified and characterized using daily aerosol optical properties from satellite measurements, namely MODIS-Terra, Earth Probe (EP)-TOMS and OMIAura. These aerosol episodes are: (i) biomass-burning/urban-industrial (BU), (ii) desert dust (DD), (iii) dust/sea-salt (DSS), (iv) mixed (MX) and (v) undetermined (UN). The identification and characterization is made with our algorithm using a variety of aerosol properties, namely aerosol optical depth (AOD), Angstrom exponent (a), fine fraction (FF), effective radius (reff) and Aerosol Index (AI). During the study period (2000e2007), the most frequent aerosol episodes are DD, observed primarily in the western and central Mediterranean Sea, and off the northern African coasts, 7 times/year for strong episodes and 4 times/year for extreme ones, on average. The DD episodes yield 40% of all types of strong aerosol episodes in the study region, while they account for 71.5% of all extreme episodes. The frequency of occurrence of strong episodes exhibits specific geographical patterns, for example the BU are mostly observed along the coasts of southern Europe and off the Atlantic coasts of Portugal, the MX episodes off the Spanish Mediterranean coast and over the Adriatic and northern Aegean Sea, while the DSS ones over the western and central Mediterranean Sea. On the other hand, the extreme episodes for all but DD aerosol display more patchy spatial patterns. The strong episodes exhibit AOD at 550 nm as high as 1.6 in the southernmost parts of central and eastern Mediterranean Sea, which rise up to 5 for the extreme, mainly DD and DSS, episodes. Although more than 90% of all aerosol episodes last 1 day, there are few cases, mainly extreme DD episodes, which last up to 4 days. Independently of their type, the Mediterranean aerosol episodes occur more frequently in spring (strong and extreme episodes) and summer (strong episodes) and most rarely during winter. A significant year by year variability of Mediterranean aerosol episodes has been identified, more in terms of their frequency than intensity. An analysis of 5-day back trajectories for the most extreme episodes provides confidence on the obtained results of the algorithm, based on the revealed origin and track of air masses causing the episodes. The 25 and 6% of all strong and extreme episodes, respectively, are MX, thus highlighting the co-existence of different aerosol types in the greater Mediterranean. The intensity of both MX and DSS episodes exhibits similar patterns to those of DD strong ones, indicating that desert dust is a determinant factor for the intensity of aerosol episodes in the Mediterranean, including DSS and MX episodes.

Atmospheric Chemistry↗

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗

Alternative Vegetation States in Tropical Forests and Savannas: The Search for Consistent Signals in Diverse Remote Sensing Data

Globally, the spatial distribution of vegetation is governed primarily by climatological factors (rainfall and temperature, seasonality, and inter-annual variability). The local distribution of vegetation, however, depends on local edaphic conditions (soils and topography) and disturbances (fire, herbivory, and anthropogenic activities). Abrupt spatial or temporal changes in vegetation distribution can occur if there are positive (i.e., amplifying) feedbacks favoring certain vegetation states under otherwise similar climatic and edaphic conditions. Previous studies in the tropical savannas of Africa and other continents using the MODerate Resolution Imaging Spectroradiometer (MODIS) vegetation continuous fields (VCF) satellite data product have focused on discontinuities in the distribution of tree cover at different rainfall levels, with bimodal distributions (e.g., concentrations of high and low tree cover) interpreted as alternative vegetation states. Such observed bimodalities over large spatial extents may not be evidence for alternate states, as they may include regions that have different edaphic conditions and disturbance histories. In this study, we conduct a systematic multi-scale analysis of diverse MODIS data streams to quantify the presence and spatial consistency of alternative vegetation states in Sub-Saharan Africa. The analysis is based on the premise that major discontinuities in vegetation structure should also manifest as consistent spatial patterns in a range of remote sensing data streams, including, for example, albedo and land surface temperature (LST). Our results confirm previous observations of bimodal and multimodal distributions of estimated tree cover in the MODIS VCF. However, strong disagreements in the location of multimodality between VCF and other data streams were observed at 1 km scale. Results suggest that the observed distribution of VCF over vast spatial extents are multimodal, not because of local-scale feedbacks and emergent bifurcations (the definition of alternative states), but likely because of other factors including regional scale differences in woody dynamics associated with edaphic, disturbance, and/or anthropogenic processes. These results suggest the need for more in-depth consideration of bifurcation mechanisms and thus the likely spatial and temporal scales at which alternative states driven by different positive feedback processes should manifest.

Savanna↗

An Example of Aerosol Pattern Variability Over Bright Surface Using High Resolution MODIS MAIAC: The Eastern and Western Areas of the Dead Sea and Environs

The extreme rate of evaporation of the Dead Sea (DS) has serious implications for the surrounding area, including atmospheric conditions. This study analyzes the aerosol properties over the western and eastern parts of the DS during the year 2013, using MAIAC (Multi-Angle Implementation of Atmospheric Correction) for MODIS, which retrieves aerosol optical depth (AOD) data at a resolution of 1km. The main goal of the study is to evaluate MAIAC over the study area and determine, for the first time, theprevailing aerosol spatial patterns. First, the MAIAC-derived AOD data was compared with data from three nearby AERONET sites (Nes Ziona - an urban site, and Sede Boker and Masada - two arid sites), and with the conventional Dark Target (DT) and Deep Blue (DB) retrievals for the same days and locations, on a monthly basis throughout 2013. For the urban site, the correlation coefficient (r) for DTDB products showed better performance than MAIAC (r = 0.80, 0.75, and 0.64 respectively) year-round. However, in the arid zones, MAIAC showed better correspondence to AERONET sites than the conventional retrievals (r = 0.58 - 0.60 and 0.48 - 0.50 respectively). We investigated the difference in AOD levels, and its variability, between the Dead Sea coasts on a seasonal basis and calculated monthlyseasonal AOD averages for presenting AOD patterns over arid zones. Thus, we demonstrated that aerosol concentrations show a strong preference for the western coast, particularly during the summer season. This preference, is most likely a result of local anthropogenic emissions combined with the typical seasonal synoptic conditions, the Mediterranean Sea breeze, and the region complex topography. Our results also indicate that a large industrial zone showed higher AOD levels compared to an adjacent reference-site, i.e., 13 percent during the winter season.

AOD↗

Intraseasonal oscillations in the global atmosphere. I - Northern Hemisphere and tropics

Oscillatory modes in the Northern Hemisphere and in the tropics were examined systematically. The 700 mb heights were used to analyze extratropical oscillations, and the outgoing longwave radiation to study tropical oscillations in convection. All datasets were band-pass filtered to focus on the intraseasonal (IS) band of 10-120 days. Leading spatial patterns of variability were obtained by applying empirical orthogonal function analysis to these IS data. The leading principal components were subjected to singular spectrum analysis. In the Northern Hemisphere, there are two important modes of oscillation with periods near 48 and 23 days, respectively. The 48-day mode is the most important of the two. It has both traveling and standing components, and is dominated by a zonal wavenumber two. The 23-day mode has the spatial structure and propagation properties described by Branstator and by Kushnir (1987). In the tropics, the 40-50 day oscillation documented by Madden and Julian (1972), Weickmann (1983), Lau, and their colleagues, dominates the Indian and Pacific oceans from 60 deg E to the date line. From 170 deg W to 90 deg W, however, a 24-28 day oscillation is equally strong. The extratropical modes are often independent of, and sometimes lead, the tropical modes.

Ghil, Michael↗

Comparing Tropopause-Penetrating Convection Identifications Derived from NEXRAD and GOES over the Contiguous United States

Overshooting cloud tops (OTs) are a result of powerful updrafts that rapidly transport air from the troposphere to the lower stratosphere. Aside from the well-known relationships between strong updrafts and severe weather, the boundaries of updrafts, including the OT, are turbulent, which leads to irreversible mixing of air between the stratosphere and troposphere, altering the composition of both layers. In order to better understand the impact of OTs, knowing precisely when and where they occur is vital. Within the last decade remote-sensing observations from satellites and radars have been used to identify OTs, but the results thus far have not been entirely consistent due to challenges with accurately identifying OTs with these sensor datasets. This study compares OTs identified using NEXRAD reflectivity data from the GridRad dataset with an approach that uses GOES infrared (IR) data. The GridRad method estimates the altitude of echo tops directly from the radar pointing geometry, merging a combination of radars that enables products every 5 mins, with 2 km spatial and 1 km vertical grid spacing. The GOES satellite method uses spatial patterns in tropopause-relative IR brightness temperature to derive an OT likelihood. The study region covers a large part of the contiguous United States during selected active convection dates in 2017 and includes data from GOES-13/16. A large sample of OTs observed by NEXRAD and GOES-13/16 satellites are analyzed to better understand NEXRAD and IR observations of OTs, quantify agreement between the OT detection methods, and demonstrate how the 4x increase in spatial resolution from GOES-13 to GOES-16 impacts OT algorithm detection performance. The results show that for nearly time-matched scenes, with many severe storms, GOES-13 and -16 detection accuracy is quite comparable; however, some missed detections are attributed to poorer GOES-13 image resolution. It was found that GOES-13 OT regions are ~2 K warmer than GOES-16 on average, but such differences exceed 5 K for the most intense OT regions. It is found that false GOES IR-based detections are unavoidable, given that cold cloud pixels significantly colder than the tropopause from nearby or recently decayed OTs can persist within anvils for 15 minutes or more.

John W. Cooney↗

Global Assessment of the Capability of Satellite Precipitation Products to Retrieve Landslide-Triggering Extreme Rainfall Events

Rainfall-induced landsliding is a global and systemic hazard that is likely to increase with the projections of increased frequency of extreme precipitation with current climate change. However, our ability to understand and mitigate landslide risk is strongly limited by the availability of relevant rainfall measurements in many landslide prone areas. In the last decade, global satellite multisensor precipitation products (SMPP) have been proposed as a solution, but very few studies have assessed their ability to adequately characterize rainfall events triggering landsliding. Here, we address this issue by testing the rainfall pattern retrieved by two SMPPs (IMERG and GSMaP) and one hybrid product [Multi-Source Weighted-Ensemble Precipitation (MSWEP)] against a large, global database of 20 comprehensive landslide inventories associated with well-identified storm events. We found that, after converting total rainfall amounts to an anomaly relative to the 10-yr return rainfall R*, the three products do retrieve the largest anomaly (of the last 20 years) during the major landslide event for many cases. However, the degree of spatial collocation of R* and landsliding varies from case to case and across products, and we often retrieved R* > 1 in years without reported landsliding. In addition, the few (four) landslide events caused by short and localized storms are most often undetected. We also show that, in at least five cases, the SMPP’s spatial pattern of rainfall anomaly matches landsliding less well than does ground-based radar rainfall pattern or lightning maps, underlining the limited accuracy of the SMPPs. We conclude on some potential avenues to improve SMPPs’ retrieval and their relation to landsliding.

Extreme events↗

Northern Late Winter Planetary Waves: MRO/MARCI Observations and Mars Climate Model Simulations

As does Earth, Mars presents pronounced global atmospheric circulation patterns. Solar differential heating drives mean meridional overturning (Hadley) circulations which are deep and intense, are hemispherically asymmetric, and where a cross-equatorial single cell dominates. Within middle and high latitudes, thermally indirect eddy-driven (Ferrel) circulation cells have been indicated. Differently, however, large-amplitude orography on planetary and continental scales on Mars can force very non-Earth-like hemispheric circulation patterns. Recent observations from the Mars Reconnaissance Orbiter, "Mars Color Imager" (MARCI) instrument are utilized that emphasize water ice clouds in ultra-violet (UV) wavelengths, and these measurements have been binned into "daily global maps" (DGMs) of water-ice cloud optical depth. The presence of large-scale, extratropical quasi-stationary atmospheric wave disturbances in middle and late winter of the northern hemisphere have been found to be present in such DGMs. In combination with such observations, a full-physics Mars global climate model (NASA ARC marsgcm 2.1) is applied to place the observations into context. During late northern winter, it is found that strong, forced Rossby modes (i.e., planetary waves) exist, and with direct correlation to columnintegrated cloud opacity undulating spatial patterns. At this season, zonal wavenumber s = 2 dominates (in contrast to wavenumber s = 1), consistent with MGS/TES analyses at this particular season (Banfield et al., 2003). Large-scale, planetary waves dictate the "coherence" of the northern polar vortex. Fundamentally, such forced planetary waves influence the polar vortex's impermeability (wave-induced) to tracer transport (e.g., dust and water-ice aerosol) and temporal mean water vapor spatial variations. The large-scale dynamical features of such planetary waves will be highlighted and discussed.

Northern polar vortex↗

Structure and Evolution of the Lunar Procellarum Region as Revealed by GRAIL Gravity Data

The Procellarum region is a broad area on the nearside of the Moon that is characterized by low elevations, thin crust, and high surface concentrations of the heat-producing elements uranium, thorium, and potassium. The Procellarum region has been interpreted as an ancient impact basin approximately 3200 km in diameter, though supporting evidence at the surface would have been largely obscured as a result of the great antiquity and poor preservation of any diagnostic features. Here we use data from the Gravity Recovery and Interior Laboratory (GRAIL) mission to examine the subsurface structure of Procellarum. The Bouguer gravity anomalies and gravity gradients reveal a pattern of narrow linear anomalies that border the Procellarum region and are interpreted to be the frozen remnants of lava-filled rifts and the underlying feeder dikes that served as the magma plumbing system for much of the nearside mare volcanism. The discontinuous surface structures that were earlier interpreted as remnants of an impact basin rim are shown in GRAIL data to be a part of this continuous set of quasi-rectangular border structures with angular intersections, contrary to the expected circular or elliptical shape of an impact basin. The spatial pattern of magmatic-tectonic structures bounding Procellarum is consistent with their formation in response to thermal stresses produced by the differential cooling of the province relative to its surroundings, coupled with magmatic activity driven by the elevated heat flux in the region.

Procellarum↗

Snow Distribution Patterns Revisited: A Physics-Based and Machine Learning Hybrid Approach to Snow Distribution Mapping in the Sub-Arctic

Snowpack distribution in Arctic and alpine landscapes often occurs in repeating, year-to-year patterns due to local topographic, weather, and vegetation characteristics. Previous studies have suggested that with years of observational data, these snow distribution patterns can be statistically integrated into a snow process modeling workflow. Recent advances in snow hydrology and machine learning (ML) have increased our ability to predict snowpack distribution using in-situ observations, remote sensing data sets, and simple landscape characteristics that can be easily obtained for most environments. Here, we propose a hybrid approach to couple a ML snow distribution pattern (MLSDP) map with a physics-based, snow process model. We trained a random forest ML algorithm on tens of thousands of snow survey observations from a subarctic study area on the Seward Peninsula, Alaska, collected during peak snow water equivalent (SWE). We validated hybrid model outputs using in-situ snow depth and SWE observations, as well as a light detection and ranging data set and a distributed temperature profiling sensor data set. When the hybrid results were compared with the physics-based method, the hybrid method more accurately depicted the spatial patterns of the snowpack, areas of drifting snow, and years when no in-situ observations were used in the random forest ML training data set. The hybrid method also showed improvements in root mean squared error at 61% of locations where time-series estimations of snow depth were observed. These results can be applied to any physics-based model to improve the snow distribution patterning to reflect observed conditions in high latitude and high elevation cold region environments.

54 ENVIRONMENTAL SCIENCES↗

Climate change-resilient snowpack estimation in the Western United States

Abstract In the 21st century, warmer temperatures and changing atmospheric circulation will likely produce unprecedented changes in Western United States snowfall 1–3 , with impacts on the timing, amount, and spatial patterns of snowpack 4–7 . The ~900 snow pillow stations are indispensable to water resource management by measuring snow-water equivalent (SWE) 8,9 in strategic but fixed locations 10,11 . However, this network may not be impacted by climate change in the same way as the surrounding area 12 and thus fail to accurately represent unmeasured locations; climate change thereby threatens our ability to measure the effects of climate change on snow. In this work, we show that maintaining the current peak SWE estimation skill is nonetheless possible. We find that explicitly including spatial correlations—either from gridded observations or learned by the model—improves skill at predicting distributed snowpack from sparse observations by 184%. Existing artificial intelligence methods can be useful tools to harness the many available sources of snowpack information to estimate snowpack in a nonstationary climate.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Forced ENSO-Like Hydrological Expressions in Simulations of the Preindustrial and Mid-Holocene

Using the water isotope- and vapor source distribution (VSD) tracer-enabled Goddard Institute for Space Studies ModelE-R, we examine changing El Nino-Southern Oscillation (ENSO)-like expressions in the hydrological cycle in a suite of model experiments. We apply strong surface temperature anomalies associated with composite observed El Nino and La Nina events as surface boundary conditions to preindustrial and mid-Holocene model experiments in order to investigate ENSO-like expressions in the hydrological cycle under varying boundary conditions. We find distinct simulated hydrological anomalies associated with El Nino-like ("ENSOWARM") and La Nina-like ("ENSOCOOL") conditions, and the region-specific VSD tracers show hydrological differences across the Pacific basin between El Nino-like and La Nina-like events. The application of ENSOCOOL forcings does not produce climatological anomalies that represent the equal but opposite impacts of the ENSOWARM experiment, as the isotopic anomalies associated with ENSOWARM conditions are generally stronger than with ENSOCOOL and the spatial patterns of change distinct. Also, when the same ENSO-like surface temperature anomalies are imposed on the mid-Holocene, the hydrological response is muted, relative to the preindustrial. Mid-Holocene changes in moisture sources to the analyzed regions across the Pacific reveal potentially complex relationships between ENSO-like conditions and boundary conditions. Given the complex impacts of ENSO-like conditions on various aspects of the hydrological cycle, we suggest that proxy record insights into paleo-ENSO variability are most likely to be robust when synthesized from a network of many spatially diverse archives, which can account for the potential nonstationarity of ENSO teleconnections under different boundary conditions.

boundary conditions↗

Quantifying Spatiotemporal Variability of Glacier Algal Blooms and the Impact on Surface Albedo in Southwestern Greenland

Albedo reduction due to light-absorbing impurities can substantially enhance ice sheet surface melt by increasing surface absorption of solar energy. Glacier algae have been suggested to play a critical role in darkening the ablation zone in southwestern Greenland. It was very recently found that the Sentinel-3 Ocean and Land Colour Instrument (OLCI) band ratio R709 nm∕R673 nm can characterize the spatial patterns of glacier algal blooms. However, Sentinel-3 was launched in 2016, and current data are only available over three melting seasons (2016–2019). Here, we demonstrate the capability of the MEdium Resolution Imaging Spectrometer (MERIS) for mapping glacier algae from space and extend the quantification of glacier algal blooms over southwestern Greenland back to the period 2004–2011. Several band ratio indices (MERIS chlorophyll a indices and the impurity index) were computed and compared with each other. The results indicate that the MERIS two-band ratio index (2BDA) R709 nm∕R665 nm is very effective in capturing the spatial distribution and temporal dynamics of glacier algal growth on bare ice in July and August. We analyzed the interannual (2004–2011) and summer (July–August) trends of algal distribution and found significant seasonal and interannual increases in glacier algae close to the Jakobshavn Isbrae Glacier and along the middle dark zone between the altitudes of 1200 and 1400 m. Using broadband albedo data from the Moderate Resolution Imaging Spectroradiometer (MODIS), we quantified the impact of glacier algal growth on bare ice albedo, finding a significant correlation between algal development and albedo reduction over algae-abundant areas. Our analysis indicates the strong potential for the satellite algal index to be used to reduce bare ice albedo biases in regional climate model simulations.

surface albedo↗

Source Decomposition of Eddy-Covariance CO 2 Flux Measurements for Evaluating a High-Resolution Urban CO 2 Emissions Inventory

We present the comparison of source-partitioned CO2 flux measurements with a high-resolution urban CO 2 emissions inventory (Hestia). Tower-based measurements of CO and 14 C are used to partition net CO 2 flux measurements into fossil and biogenic components. A flux footprint model is used to quantify spatial variation in flux measurements. We compare the daily cycle and spatial structure of Hestia and eddy-covariance derived fossil fuel CO 2 emissions on a seasonal basis. Hestia inventory emissions exceed the eddy-covariance measured emissions by 0.36 µ mol m −2 s −1 (3.2%) in the cold season and 0.62 µ mol m −2 s −1 (9.1%) in the warm season. The daily cycle of fluxes in both products matches closely, with correlations in the hourly mean fluxes of 0.86 (cold season) and 0.93 (warm season). The spatially averaged fluxes also agree in each season and a persistent spatial pattern in the differences during both seasons that may suggest a bias related to residential heating emissions. In addition, in the cold season, the magnitudes of average daytime biological uptake and nighttime respiration at this flux site are approximately 15% and 27% of the mean fossil fuel CO 2 emissions over the same time period, contradicting common assumptions of no significant biological CO 2 exchange in northern cities during winter. This work demonstrates the effectiveness of using trace gas ratios to adapt eddy-covariance flux measurements in urban environments for disaggregating anthropogenic CO 2 emissions and urban ecosystem fluxes at high spatial and temporal resolution.

eddy-covariance flux measurements↗

Daily, 30 m Resolution NDSI Data for the East River Watershed, CO for 2000-2020

This dataset contains daily Normalized Difference Snow Index (NDSI) values at 30 m spatial resolution for the East River watershed in Colorado, USA. The temporal range of these data includes water years 2001-2020. These data were created using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). This model fuses low spatial and high temporal resolution data from MODIS (500 m, daily) with high spatial and low temporal resolution data from Landsat (30 m, 16 days) to create a 30m synthetic daily snow product. This product allows for the analysis of historical snow covered area trends in the East River Watershed at fine spatiotemporal resolutions where it was not available previously. This research was performed as a part of the Department of Energy’s Subsurface Biogeochemical Research Program with the primary intent of better understanding the timing and spatial patterns of water delivery to the Critical Zone in mountain watersheds. Each .zip file contains one "water year" of data (October 1 - September 30; i.e., water year 2010 starts October 1, 2010 and ends September 30, 2011). Each zip file contains the following: STARFM daily Normalized Difference Snow Index (NDSI) fusion data files in GeoTiff format with one layer for each day between Landsat data acquisition dates (i.e., for dates of Landsat acquisition, the Landsat image is included for that date). The study area is located in an area of Landsat path overlap, so Landsat dates acquisitions are every 7-9 days. Landsat NDSI files containing the high spatial (30m), low temporal (7-9 days due to Landsat path overlap) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Dates for which no Landsat data were obtained are included as NoData layers. MODIS NDSI files containing the high temporal (daily), low spatial (500m) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Please note the MODIS data were resampled to 30m pixels for input into the STARFM model. The data have a scale factor of 10,000 and a no data value of -32767. The projection of all datasets is WGS 84 (EPSG: 4326), which has a latitude/longitude based degree resolution of 0.0002694946 X 0.0002694946, and approximates to the 30 m spatial resolution mentioned above. The Layer Index files in .csv format. They contain information for each layer in the above GeoTiff files regarding the corresponding date for each layer, the fraction of pixels in the image that contain valid data (missing data is due to either cloud cover or poor data quality; these values are not percent snow cover). Dates of Landsat overpass are indicated in these files. If no Landsat data were able to be obtained due to cloud cover or lack of Landsat Tier 1 data available on Google Earth Engine, this is also noted.

EARTH SCIENCE > CRYOSPHERE > SNOW/ICE↗

Interannual Variability and Trends of Extratropical Ozone: Northern Hemisphere - Part 1

The authors apply principal component analysis (PCA) to the extratropical total column ozone from the combined merged ozone data product and the European Centre for Medium-Range Weather Forecasts assimilated ozone from January 1979 to August 2002. The interannual variability (IAV) of extratropical O-3 in the Northern Hemisphere (NH) is characterized by four main modes. Attributable to dominant dynamical effects, these four modes account for nearly 60% of the total ozone variance in the NH. The patterns of variability are distinctly different from those derived for total O-3 in the tropics. To relate the derived patterns of O-3 to atmospheric dynamics, similar decompositions are performed for the 30 100-Wa geopotential thickness. The results reveal intimate connections between the IAV of total ozone and the atmospheric circulation. The first two leading modes are nearly zonally symmetric and represent the connections to the annular modes and the quasi-biennial oscillation. The other two modes exhibit in-quadrature, wavenumber-1 structures that, when combined, describe the displacement of the polar vortices in response to planetary waves. In the NH, the extrema of these combined modes have preferred locations that suggest fixed topographical and land-sea thermal forcing of the involved planetary waves. Similar spatial patterns and trends in extratropical column ozone are simulated by the Goddard Earth Observation System chemistryclimate model (GEOS-CCM). The decreasing O-3 trend is captured in the first mode. The largest trend occurs at the North Pole, with values similar to-1 Dobson Unit (DU) yr(-1). There is almost no trend in tropical O-3. The trends derived from PCA are confirmed using a completely independent method, empirical mode decomposition, for zonally averaged O-3 data. The O-3 trend is also captured by mode 1 in the GEOS-CCM, but the decrease is substantially larger than that in the real atmosphere.

Yung, Yuk L.↗

Complex Spatial Structure in a Population of Didymopanax pittieri, A Tree of Wind-Exposed Lower Montane Rain Forest

Didymopanax pittieri is a common shade-intolerant tree colonizing treefall gaps in the elfin forests on windswept ridgecrests in the lower montane rain forests of the Cordillera de Tilarain, Costa Rica. All D. pittieri taller than > 0.5 m in a 5.2-ha elfin forested portion of a gridded study watershed in the Monteverde Cloud Forest Preserve were located, mapped, and measured. This local population of D. pittieri is spatially inhomogeneous, in that density increases with increasing wind exposure; D. pittieri are more abundant near ridge crests than lower on windward slopes. The important and ubiquitous phenomenon of spatial inhomogeneity in population density is addressed and corrected for in spatial analyses by the application of the inhomogeneous version of Ripley's K. The spatial patterns of four size classes of D. pittieri (<5 cm dbh, 5-10 cm dbh, 10-20 cm dbh, and> 20 cm dbh) were investigated. Within the large-scale trend in density driven by wind exposure, D. pittieri saplings are clumped at the scale of treefall gaps and at the scale of patches of aggregated gaps. D. pittieri 5-10 cm dbh are randomly distributed, apparently due to competitive thinning of sapling clumps during the early stages of gap-phase regeneration. D. pittieri larger than 10 cm dbh are overdispersed at a scale larger than that of patches of gaps. Natural disturbance can influence the distribution of shade intolerant tree populations at several different spatial scales, and can have discordant effects at different life history stages.

Lawton, Robert M.↗

Tracking the Rates and Mechanisms of Canopy Damage and Recovery Following Hurricane Maria Using Multitemporal Lidar Data

Hurricane Maria, a Category 4 storm, snapped and uprooted canopy trees, removed large branches, and defoliated vegetation across Puerto Rico. The magnitude of forest damages and the rates and mechanisms of forest recovery following Maria provide important benchmarks for understanding the ecology of extreme events. We used airborne Lidar data acquired before (2017) and after Maria (2018, 2020) to quantify landscape-scale changes in forest structure along a 439-ha elevational gradient (100–800 m) in the Luquillo Experimental Forest. Damages from Maria were widespread, with 73% of the study area losing ≥ 1 m in canopy height (mean = −7.1 m). Taller forests at lower elevations suffered more damage than shorter forests above 600 m. Yet only 13.5% of the study area had canopy heights ≤ 2 m in 2018, a typical threshold for forest gaps, highlighting the importance of damaged trees and advanced regeneration on post-storm forest structure. Heterogeneous patterns of regrowth and recruitment yielded shorter and more open forests by 2020. Nearly 45% of forests experienced initial height loss > 1 m (2017–2018) followed by rapid height gain > 1 m (2018–2020), whereas 21.6% of forests with initial height losses showed little or no height gain, and 17.8% of forests exhibited no height changes larger than ± 1 m in either period. Canopy layers < 10 m tall accounted for most increases in canopy height and fractional cover between 2018 and 2020, with gains split evenly between height growth and lateral crown expansion by surviving individuals. These findings benchmark rates of gap formation, crown expansion, and canopy closure following hurricane damage and highlight the diversity of ecosystem impacts from heterogeneous spatial patterns and vertical stratification of forest regrowth following a major disturbance event.

forest structure↗