237 research outputs found

    On The Development of a Dynamic Contrast-Enhanced Near-Infrared Technique to Measure Cerebral Blood Flow in the Neurocritical Care Unit

    Get PDF
    A dynamic contrast-enhanced (DCE) near-infrared (NIR) method to measure cerebral blood flow (CBF) in the neurocritical care unit (NCU) is described. A primary concern in managing patients with acquired brain injury (ABI) is onset of delayed ischemic injury (DII) caused by complications during the days to weeks following the initial insult, resulting in reduced CBF and impaired oxygen delivery. The development of a safe, portable, and quantitative DCE-NIR method for measuring CBF in NCU patients is addressed by focusing on four main areas: designing a clinically compatible instrument, developing an appropriate analytical framework, creating a relevant ABI animal model, and validating the method against CT perfusion. In Chapter 2, depth-resolved continuous-wave NIR recovered values of CBF in a juvenile pig show strong correlation with CT perfusion CBF during mild ischemia and hyperemia (r=0.84, p\u3c0.001). In particular, subject-specific light propagation modeling reduces the variability caused by extracerebral layer contamination. In Chapter 3, time-resolved (TR) NIR improves the signal sensitivity to brain tissue, and a relative CBF index is be both sensitive and specific to flow changes in the brain. In particular, when compared with the change in CBF measured with CT perfusion during hypocapnia, the deconvolution-based index has an error of 0.8%, compared to 21.8% with the time-to-peak method. To enable measurement of absolute CBF, a method for characterizing the AIF is described in Chapter 4, and the theoretical basis for an advanced analytical framework—the kinetic deconvolution optical reconstruction (KDOR)—is provided in Chapter 5. Finally, a multichannel TR-NIR system is combined with KDOR to quantify CBF in an adult pig model of ischemia (Chapter 6). In this final study, measurements of CBF obtained with the DCE-NIR technique show strong agreement with CT perfusion measurements of CBF in mild and moderate ischemia (r=0.86, p\u3c0.001). The principle conclusion of this thesis is that the DCE-NIR method, combining multidistance TR instrumentation with the KDOR analytical framework, can recover CBF values that are in strong agreement with CT perfusion values of CBF. Ultimately, bedside CBF measurements could improve clinical management of ABI by detecting delayed ischemia before permanent brain damage occurs

    Responses and adaptation strategies of terrestrial ecosystems to climate change

    Get PDF
    Terrestrial ecosystems are likely to be affected by climate change, as climate change-induced shift of water and heat stresses patterns will have significant impacts on species composition, habitat distribution, and ecosystem functions, and thereby weaken the terrestrial carbon (C) sink and threaten global food security and biofuel production. This thesis investigates the responses of terrestrial ecosystems to climate change and is structured in four main chapters.;The first chapter of the thesis is directed towards the impacts of snow variation on ecosystem phenology. Variations in seasonal snowfall regulate regional and global climatic systems and vegetation growth by changing energy budgets of the lower atmosphere and land surface. We investigated the effects of snow on the start of growing season (SGS) of temperate vegetation in China. Across the entire temperate region in China, the winter snow depth increased at a rate of 0.15 cm•yr-1 (p=0.07) during the period 1982-1998, and decreased at a rate of 0.36 cm•yr-1 (p=0.09) during the period 1998-2005. Correspondingly, the SGS advanced at a rate of 0.68 d•yr-1 (p\u3c0.01) during 1982 to 1998, and delayed at a rate of 2.13 d•yr-1 (p=0.07) during 1998 to 2005, against a warming trend throughout the entire study period of 1982-2005. Spring air temperature strongly regulated the SGS of both deciduous broad-leaf and coniferous forests; whilst the winter snow had a greater impact on the SGS of grassland and shrubs. Snow depth variation combined with air temperature contributed to the variability in the SGS of grassland and shrubs, as snow acted as an insulator and modulated the underground thermal conditions. Additionally, differences were seen between the impacts of winter snow depth and spring snow depth on the SGS; as snow depths increased, the effect associated went from delaying SGS to advancing SGS. The observed thresholds for these effects were snow depths of 6.8 cm (winter) and 4.0 cm (spring). The results of this study suggest that the response of the vegetation\u27s SGS to seasonal snow change may be attributed to the coupling effects of air temperature and snow depth associated with the soil thermal conditions.;The second chapter further addresses snow impacts on terrestrial ecosystem with focus on regional carbon exchange between atmosphere and biosphere. Winter snow has been suggested to regulate terrestrial carbon (C) cycling by modifying micro-climate, but the impacts of snow cover change on the annual C budget at the large scale are poorly understood. Our aim is to quantify the C balance under changing snow depth. Here, we used site-based eddy covariance flux data to investigate the relationship between snow cover depth and ecosystem respiration (Reco) during winter. We then used the Biome-BGC model to estimate the effect of reductions in winter snow cover on C balance of Northern forests in non-permafrost region. According to site observations, winter net ecosystem C exchange (NEE) ranged from 0.028-1.53 gC•m-2•day-1, accounting for 44 +/- 123% of the annual C budget. Model simulation showed that over the past 30 years, snow driven change in winter C fluxes reduced non-growing season CO2 emissions, enhancing the annual C sink of northern forests. Over the entire study area, simulated winter ecosystem respiration (Reco) significantly decreased by 0.33 gC•m-2•day -1•yr-1 in response to decreasing snow cover depth, which accounts for approximately 25% of the simulated annual C sink trend from 1982 to 2009. Soil temperature was primarily controlled by snow cover rather than by air temperature as snow served as an insulator to prevent chilling impacts. A shallow snow cover has less insulation potential, causing colder soil temperatures and potentially lower respiration rates. Both eddy covariance analysis and model-simulated results showed that both Reco and NEE were significantly and positively correlated with variation in soil temperature controlled by variation in snow depth. Overall, our results highlight that a decrease in winter snow cover restrains global warming through emitting less C to the atmosphere.;The third chapter focused on assessing drought\u27s impact on global terrestrial ecosystems. Drought can affect the structure, composition and function of terrestrial ecosystems, yet the drought impacts and post-drought recovery potential of different land cover types have not been extensively studied at a global scale. Here, we evaluated drought impacts on gross primary productivity (GPP), evapotranspiration (ET), and water use efficiency (WUE) of different global terrestrial ecosystems, as well as the drought-resilience of each ecosystem type during the period of 2000 to 2011. We found the rainfall and soil moisture during drought period were dramatically lower than these in non-drought period, while air temperatures were higher than normal during drought period with amplitudes varied by land cover types. The length of recovery days (LRD) presented an evident gradient of high (\u3e 60 days) in mid- latitude region and low (\u3c 60 days) in low (tropical area) and high (boreal area) latitude regions. As average GPP increased, the LRD showed a significantly decreasing trend among different land covers (R 2=0.53, p\u3c0.0001). Moreover, the most dramatic reduction of the drought-induced GPP was found in the mid-latitude region of north Hemisphere (48% reduction), followed by the low-latitude region of south Hemisphere (13% reduction). In contrast, a slightly enhanced GPP (10%) was showed in the tropical region under drought impact. Additionally, the highest drought-induced reduction of ET was found in the Mediterranean area, followed by Africa. The water use efficiency, however, showed a pattern of decreasing in the north Hemisphere and increasing in the south Hemisphere.;The last chapter compared the differences of performance in trading water for carbon in planted forest and natural forest, with specific focus on China. Planted forests have been widely established in China as an essential approach to improving the ecological environment and mitigating climate change. Large-scale forest planting programs, however, are rarely examined in the context of tradeoffs between carbon sequestration and water yield between planted and natural forests. We reconstructed evapotranspiration (ET) and gross primary production (GPP) data based on remote-sensing and ground observational data, and investigated the differences between natural and planted forests, in order to evaluate the suitability of tree-planting activity in different climate regions where the afforestation and reforestation programs have been extensively implemented during the past three decades in China. While the differences changed with latitude (and region), we found that, on average, planted forests consumed 5.79% (29.13mm) more water but sequestered 1.05% (-12.02 gC m-2 yr -1) less carbon than naturally generated forests, while the amplitudes of discrepancies varied with latitude. It is suggested that the most suitable lands in China for afforestation should be located in the moist south subtropical region (SSTP), followed by the mid-subtropical region (MSTP), to attain a high carbon sequestration potential while maintain a relatively low impact on regional water balance. The high hydrological impact zone, including the north subtropical region (NSTP), warm temperate region (WTEM), and temperate region (TEM) should be cautiously evaluated for future afforestation due to water yield reductions associated with plantations

    SATELLITE MICROWAVE MEASUREMENT OF LAND SURFACE PHENOLOGY: CLARIFYING VEGETATION PHENOLOGY RESPONSE TO CLIMATIC DRIVERS AND EXTREME EVENTS

    Get PDF
    The seasonality of terrestrial vegetation controls feedbacks to the climate system including land-atmosphere water, energy and carbon (CO2) exchanges with cascading effects on regional-to-global weather and circulation patterns. Proper characterization of vegetation phenology is necessary to understand and quantify changes in the earthÆs ecosystems and biogeochemical cycles and is a key component in tracking ecological species response to climate change. The response of both functional and structural vegetation phenology to climatic drivers on a global scale is still poorly understood however, which has hindered the development of robust vegetation phenology models. In this dissertation I use satellite microwave vegetation optical depth (VOD) in conjunction with an array of satellite measures, Global Positioning System (GPS) reflectometry, field observations and flux tower data to 1) clarify vegetation phenology response to water, temperature and solar irradiance constraints, 2) demonstrate the asynchrony between changes in vegetation water content and biomass and changes in greenness and leaf area in relation to land cover type and climate constraints, 3) provide enhanced assessment of seasonal recovery of vegetation biomass following wildfire and 4) present a method to more accurately model tropical vegetation phenology. This research will establish VOD as a useful and informative parameter for regional-to-global vegetation phenology modeling, more accurately define the drivers of both structural and functional vegetation phenology, and help minimize errors in phenology simulations within earth system models. This dissertation also includes the development of Gross Primary Productivity (GPP) and Net Primary Productivity (NPP) vegetation health climate indicators as part of a NASA funded project entitled Development and Testing of Potential Indicators for the National Climate Assessment; Translating EOS datasets into National Ecosystem Biophysical Indicators

    Development and Extrapolation of a General Light Use Efficiency Model for the Gross Primary Production

    Get PDF
    The global carbon cycle is one of the large biogeochemical cycles spanning all living and non-living compartments of the Earth system. Against the background of accelerating global change, the scientific community is highly interested in analyzing and understanding the dynamics of the global carbon cycle and its complex feedback mechanism with the terrestrial biosphere. The international network FLUXNET was established to serve this aim with measurement towers around the globe. The overarching objective of this thesis is to exploit the powerful combination of carbon flux measurements and satellite remote sensing in order to develop a simple but robust model for the gross primary production (GPP) of vegetation stands. Measurement data from FLUXNET sites as well as remote sensing data from the NASA sensor MODIS are exploited in a data-based model development approach. The well-established concept of light use efficiency is chosen as modeling framework. As a result, a novel gross primary production model is established to quantify the carbon uptake of forests and grasslands across a broad range of climate zones. Furthermore, an extrapolation scheme is derived, with which the model parameters calibrated at FLUXNET sites can be regionalized to pave the way for spatially continuous model applications

    Synthesis of Satellite Microwave Observations for Monitoring Global Land-Atmosphere CO2 Exchange

    Get PDF
    This dissertation describes the estimation, error quantification, and incorporation of land surface information from microwave satellite remote sensing for modeling global ecosystem land-atmosphere net CO2 exchange. Retrieval algorithms were developed for estimating soil moisture, surface water, surface temperature, and vegetation phenology from microwave imagery timeseries. Soil moisture retrievals were merged with model-based soil moisture estimates and incorporated into a light-use efficiency model for vegetation productivity coupled to a soil decomposition model. Results, including state and uncertainty estimates, were evaluated with a global eddy covariance flux tower network and other independent global model- and remote-sensing based products

    Multichannel near infrared spectroscopy to monitor cerebral oxygenation in infants and children supported in extracorporeal membrane oxygenation (ECMO)

    Get PDF
    Extracorporeal membrane oxygenation is a life support system for infants and children with intractable cardiorespiratory failure. The risk for developing neurological injuries in this group of patients is considerable. The causes are multifactorial and are not yet fully understood. Induction of ECMO involves ligation of the major neck vessels – common carotid artery and internal jugular vein – which may cause lateralised cerebrovascular injury. Physiologic changes such as hypoxia, hypotension and hypercarbia associated with ECMO can disrupt cerebral autoregulation. Near infrared spectroscopy (NIRS) offers the advantage of continuous non invasive means of monitoring cerebral oxygenation at the bedside. To date, NIRS systems used clinically are single or dual channel systems and do not allow evaluation of the status of cerebral circulation in the extended cerebral regions. This work involves the development of a multichannel NIRS system for use in paediatric cardiothoracic intensive care with specific application on patients supported on extracorporeal circulation. A novel flexible neonatal cap was designed and constructed to accommodate an array of sources and detectors that provide measurements of multisite cerebral oxygenation from 12 channels. Multimodal data collection (systemic and ECMO circuit parameters) simultaneous with multichannel NIRS was established to allow monitoring of multisite cerebral oxygenation and haemodynamics. A novel method of analysis, wavelet cross-correlation, was generated to study the concordance between multisite oxyhaemoglobin concentration (HbO2) and mean arterial pressure (MAP) as a means to investigate regional variations in cerebral circulation and assess cerebral autoregualtion. Group data of 6 patients showed statistically significant differences in WCC between right and left hemispheres during sequential changes in ECMO circuit blood flow. WCC between HbO2 and MAP provides a useful method to investigate the dynamics of cerebral autoregulation during ECMO. Modest manipulations of ECMO flows are associated with regional changes in cerebral autoregulation which may potentially have an important bearing on clinical outcome

    HUMAN AND CLIMATE IMPACTS ON FLOODING VIA REMOTE SENSING, BIG DATA ANALYTICS, AND MODELING

    Get PDF
    Over the last 20 years, the amount of streamflow has greatly increased and spring snowmelt floods have occurred more frequently in the north-central U.S. In the Red River of the North Basin (RRB) overlying portions of North Dakota and Minnesota, six of the 13 major floods over the past 100 years have occurred since the late 1990s. Based on numerous previous studies as well as senior flood forecasters’ experiences, recent hydrological changes related to human modifications [e.g. artificial subsurface drainage (SSD) expansion] and climate change are potential causes of notable forecasting failures over the past decade. My dissertation focuses on the operational and scientific gaps in current forecasting models and observational data and provides insights and value to both the practitioner and the research community. First, the current flood forecasting model needs both the location and installation timing of SSD and SSD physics. SSD maps were developed using satellite “big” data and a machine learning technique. Next, using the maps with a land surface model, the impacts of SSD expansion on regional hydrological changes were quantified. In combination with model physics, the inherent uncertainty in the airborne gamma snow survey observations hinders the accurate flood forecasting model. The operational airborne gamma snow water equivalent (SWE) measurements were improved by updating antecedent surface moisture conditions using satellite observations on soil moisture. From a long-term perspective, flood forecasters and state governments need knowledge of historical changes in snowpack and snowmelt to help flood management and to develop strategies to adapt to climate changes. However, historical snowmelt trends have not been quantified in the north-central U.S. due to the limited historical snow data. To overcome this, the current available historical long-term SWE products were evaluated across diverse regions and conditions. Using the most reliable SWE product, a trend analysis quantified the magnitude of change extreme snowpack and melt events over the past 36 years. Collectively, this body of research demonstrates that human and climate impacts, as well as limited and noisy data, cause uncertainties in flood prediction in the great plains, but integrated approaches using remote sensing, big data analytics, and modeling can quantify the hydrological changes and reduce the uncertainties. This dissertation improves the practice of flood forecasting in Red River of the North Basin and advances research in hydrology and snow science

    Examining Ecosystem Drought Responses Using Remote Sensing and Flux Tower Observations

    Get PDF
    Indiana University-Purdue University Indianapolis (IUPUI)Water is fundamental for plant growth, and vegetation response to water availability influences water, carbon, and energy exchanges between land and atmosphere. Vegetation plays the most active role in water and carbon cycle of various ecosystems. Therefore, comprehensive evaluation of drought impact on vegetation productivity will play a critical role for better understanding the global water cycle under future climate conditions. In-situ meteorological measurements and the eddy covariance flux tower network, which provide meteorological data, and estimates of ecosystem productivity and respiration are remarkable tools to assess the impacts of drought on ecosystem carbon and water cycles. In regions with limited in-situ observations, remote sensing can be a very useful tool to monitor ecosystem drought status since it provides continuous observations of relevant variables linked to ecosystem function and the hydrologic cycle. However, the detailed understanding of ecosystem responses to drought is still lacking and it is challenging to quantify the impacts of drought on ecosystem carbon balance and several factors hinder our explicit understanding of the complex drought impacts. This dissertation addressed drought monitoring, ecosystem drought responses, trends of vegetation water constraint based on in-situ metrological observations, flux tower and multi-sensor remote sensing observations. This dissertation first developed a new integrated drought index applicable across diverse climate regions based on in-situ meteorological observations and multi-sensor remote sensing data, and another integrated drought index applicable across diverse climate regions only based on multi-sensor remote sensing data. The dissertation also evaluated the applicability of new satellite dataset (e.g., solar induced fluorescence, SIF) for responding to meteorological drought. Results show that satellite SIF data could have the potential to reflect meteorological drought, but the application should be limited to dry regions. The work in this dissertation also accessed changes in water constraint on global vegetation productivity, and quantified different drought dimensions on ecosystem productivity and respiration. Results indicate that a significant increase in vegetation water constraint over the last 30 years. The results highlighted the need for a more explicit consideration of the influence of water constraints on regional and global vegetation under a warming climate

    Land surface temperature and evapotranspiration estimation in the Amazon evergreen forests using remote sensing data

    Get PDF
    Amazonian tropical forests play a significant role in global water, carbon and energy cycles. Considering the relevance of this biome and the climate change projections which predict a hotter and drier climate for the region, the monitoring of the vegetation status of these forests becomes of significant importance. In this context, vegetation temperature and evapotranspiration (ET) can be considered as key variables. Vegetation temperature is directly linked with plant physiology. In addition, some studies have shown the existing relationship between this variable and the CO2 absorption capacity and biomass loss of these forests. Evapotranspiration resulting from the combined processes of transpiration and evaporation links the terrestrial water, carbon and surface energy exchanges of these forests. How this variable will response to the changing climate is critical to understand the stability of these forests. Satellite remote sensing is presented as a feasible means in order to provide accurate spatially-distributed estimates of these variables. Nevertheless, the use of satellite passive imagery for analysing this region still has some limitations being of special importance the proper cloud masking of the satellite data which becomes a difficult task due to the continuous cloud cover of the region. Under the light of the aforementioned issues, the present doctoral thesis is aimed at estimating the land surface temperature and evapotranspiration of the Amazonian tropical forests using remote sensing data. In addition, as cloud screening of satellite imagery is a critical step in the processing chain of the previous magnitudes and becomes of special importance for the study region this topic has also been included in this thesis. We have mainly focused on the use of data from the Moderate Resolution Imaging Spectroradiometer (MODIS) which is amongst major tools for studying this region. Regarding the cloud detection topic, the potential of supervised learning algorithms for cloud masking was studied in order to overcome the cloud contamination issue of the current satellite products. Models considered were: Gaussian NaĂŻve Bayes (GNB), Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Random Forests (RF), Support Vector Machine (SVM) and Multilayer Perceptron (MLP). These algorithms are able to provide a continuous measure of cloud masking uncertainty (i.e. a probability estimate of each pixel belonging to clear and cloudy class) and therefore can be used under the light of a probabilistic approach. Reference dataset (a priori knowledge) requirement was satisfied by considering the collocation of Cloud Profiling Radar (CPR) and Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) observations with MODIS sensor. Model performance was tested using three independent datasets: 1) collocated CPR/CALIOP and MODIS data, 2) MODIS manually classified images and 3) in-situ ground data. For the case of satellite image and in-situ testing, results were additionally compared to current operative MYD35 (version 6.1) and Multi-Angle Implementation of the Atmospheric Correction (MAIAC) cloud masking algorithms. These results showed that machine learning algorithms were able to improve MODIS operative cloud masking performance over the region. MYD35 and MAIAC tended to underestimate and overestimate the cloud cover, respectively. Amongst the models considered, LDA stood out as the best candidate because of its maximum accuracy (difference in Kappa coefficient of 0.293/0.155 (MYD35 /MAIAC respectively)) and minimum computational associated. Regarding the estimation of land surface temperature (LST), the aim of this study was to generate specific LST products for the Amazonian tropical forests. This goal was accomplished by using a tuned split-window (SW) equation. Validation of the LST products was obtained by direct comparison between LST estimates as derived from the algorithms and two types of different LST observations: in-situ LST (T-based validation) and LST derived from the R-based method. In addition, LST algorithms were validated using independent simulated data. In-situ LST was retrieved from two infrared radiometers (SI-100 and IR-120) and a CNR4 net radiometer, situated at Tambopata test site (12.832 S, 62.282 W) in the Peruvian Amazon. Apart from this, current satellite LST products were also validated and compared to the tuned split-window. Although we have mainly focus on MODIS LST products which derive from three different LST algorithms: split-window, day and night (DN) and Temperature Emissivity Separation (TES), we have also considered the inclusion of the Visible Infrared Imaging Radiometer Suite (VIIRS) sensor. In addition, a first assessment of the Sea and Land Surface Temperature Radiometer (SLSTR) is presented. Validation was performed separately for daytime and nighttime conditions. For MODIS sensor, current LST products showed Root Mean Square Errors (RMSE) in LST estimations between 2 K and 3K for daytime and 1 K and 2 K for nighttime. In the best case (with a restrictive cloud screening) RMSE errors decrease to values below 2K and around 1 K, respectively. The proposed LST showed RMSE values of 1K to approximately 2 K and 0.7-1.5 K (below 1.5 K and below 1 K in the best case) for daytime and nighttime conditions, thus improving current LST MODIS products. This is also in agreement with the R-based validation results, which show a RMSE reduction of 0.7 K to 1.7 K in comparison to MODIS LST products. For the case of VIIRS sensor daytime conditions, VIIRS-TES algorithm provides the best performance with a difference of 0.2 K to around 0.3 K in RMSE regarding the split window algorithm (in the best case it reduces to 0.2 K). All VIIRS LST products considered have RMSE values between 2 K and 3K. At nighttime, however VIIRS-TES is not able to outperform the SW algorithm. A difference of 0.7 K to 0.8 K in RMSE is obtained. Contrary to MODIS and the SW LST products, VIIRS-TES tends to overestimate in-situ LST values. Regarding SLSTR sensor, the L2 product provides a better agreement with in-situ observations than the proposed algorithm (daytime difference in RMSE around 0.6 K and up 0.07 K at nighttime). In the estimation of the ET, we focused on the evaluation of four commonly used remote-sensing based ET models. These were: i) Priestley-Taylor Jet Propulsion Laboratory (PT-JPL), ii) Penman-Monteith MODIS operative parametrization (PM-Mu), iii) Surface Energy Balance System (SEBS), and iv) Satellite Application Facility on Land Surface Analysis (LSASAF). These models were forced using remote-sensing data from MODIS and two ancillary meteorological data sources: i) in-situ data extracted from Large-Scale Biosphere-Atmosphere Experiment (LBA) stations (scenario I), and ii) three reanalysis datasets (scenario II), including Modern-Era Retrospective analysis for Research and Application (MERRA-2), European Centre for Medium-range Weather Forecasts (ECMWF) Re-Analysis-Interim (ERA-Interim), and Global Land Assimilation System (GLDAS-2.1). Performance of algorithms under the two scenarios was validated using in-situ eddy-covariance measurements. For scenario I, PT-JPL provided the best agreement with in-situ ET observations (RMSE = 0.55 mm/day, R = 0.88). Neglecting water canopy evaporation resulted in an underestimation of ET measurements for LSASAF. SEBS performance was similar to that of PT-JPL, nevertheless SEBS estimates were limited by the continuous cloud cover of the region. A physically-based ET gap-filling method was used in order to alleviate this issue. PM-Mu also with a similar performance to PT-JPL tended to overestimate in-situ ET observations. For scenario II, quality assessment of reanalysis input data demonstrated that MERRA-2, ERA-Interim and GLDAS-2.1 contain biases that impact model performance. In particular, biases in radiation inputs were found the main responsible of the observed biases in ET estimates. For the region, MERRA-2 tends to overestimate daily net radiation and incoming solar radiation. ERA-Interim tends to underestimate both variables, and GLDAS-2.1 tends to overestimate daily radiation while underestimating incoming solar radiation. Discrepancies amongst these inputs resulted in large absolute deviations in spatial patterns (deviations greater than 500 mm/year) and temporal patterns
    • …
    corecore