CN113837620A - Method for evaluating temperature extreme value change - Google Patents

Method for evaluating temperature extreme value change Download PDF

Info

Publication number
CN113837620A
CN113837620A CN202111132485.6A CN202111132485A CN113837620A CN 113837620 A CN113837620 A CN 113837620A CN 202111132485 A CN202111132485 A CN 202111132485A CN 113837620 A CN113837620 A CN 113837620A
Authority
CN
China
Prior art keywords
temperature
extreme value
data
temperature extreme
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111132485.6A
Other languages
Chinese (zh)
Other versions
CN113837620B (en
Inventor
翟媛媛
黄国和
周雄
吴莹辉
鲁晨
宋唐女
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN202111132485.6A priority Critical patent/CN113837620B/en
Publication of CN113837620A publication Critical patent/CN113837620A/en
Application granted granted Critical
Publication of CN113837620B publication Critical patent/CN113837620B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/10Devices for predicting weather conditions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Administration (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Game Theory and Decision Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Environmental Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

The invention provides a method for evaluating temperature extreme value change, which comprises the following steps of firstly, acquiring temperature extreme value index data, circulation field data and large-scale circulation mode index data; then, performing EOF decomposition on the data based on different time scales and different seasons, and further performing trend change analysis of the developed temperature extreme value according to least square regression and Mann-Kendall test; secondly, identifying factors of a single circulation field and large-scale circulation mode indexes which affect the temperature extreme value through Pearson correlation analysis; and finally, identifying factors which significantly contribute to the temperature extreme value based on a factor analysis method, and quantifying the influence of the interaction effect among the key factors on the temperature extreme value.

Description

Method for evaluating temperature extreme value change
Technical Field
The invention belongs to the field of climate research, and particularly relates to a method for evaluating temperature extreme value change.
Background
Over the past few decades, temperature and extreme temperature anomalies have had extremely serious negative impacts on society, economy, and people's life. The research on the extreme climate value and the influence mechanism thereof at home and abroad is attached with importance, and the understanding of the change mechanism of the extreme climate is very important for predicting the occurrence of the extreme climate event and taking measures to reduce the relevant influence. There are intricate and complex interaction effects among various elements (such as potential height, wind field, remote correlation factor, etc.) which affect the extreme temperature, but these effects are difficult to be characterized by the traditional functional form. Currently, most of the research in this area is only directed to one or more individual elements, and the influence of the interaction between different elements on the temperature extreme value is rarely considered.
Object of the Invention
The invention aims to solve the problems in the prior art and provide a method for evaluating the change of the extreme temperature value, thereby providing favorable support for the future extreme climate change.
Disclosure of Invention
The invention provides a method for evaluating temperature extreme value change, which comprises the following steps:
step a, respectively acquiring temperature extreme value data, circulation field data and large-scale circulation mode index data; the definition of the temperature extreme value data is formulated by a world weather organization WMO climate change detection and index expert group ETCCDI, and is obtained by a land climate extreme event index data set HadEX3 covering the land climate extreme event data set with the precision of 1.875 multiplied by 1.25 longitude and latitude from 1901 to 2018, and specific indexes comprise a daily maximum air temperature value TXx, a daily minimum air temperature value TNx, a daily maximum air temperature value TXn, a daily minimum air temperature value TNn, a cold night TN10p, a warm night TN90p, a cold day TX10p and a hot day TX90 p; the circulation field data comprises sea level pressure SLP, 500hPa wind field and 500hPa potential height of ERA5 reanalysis data; the large-scale circulation pattern data are indexes for measuring a natural variability pattern, are time sequence data with a scale of days and comprise index items of Erleno and southern billow ENSO, Pacific decade billow PDO, North Atlantic billow NAO, Arctic billow AO and Pacific-North Atlantic related PNA;
b, analyzing temperature extreme values from a time sequence and a space scale respectively according to different time dimensions of each year, years, summer and winter by using an empirical orthogonal function EOF method, and identifying areas with extremely high temperature and low temperature in summer in the northern hemisphere area; the specific processing procedure of the EOF method is as follows: setting the number of lattice points of the temperature extreme value as m and the time sequence length as n, firstly, processing the temperature extreme value data into a distance form to obtain a matrix X related to the temperature extreme value time and spacem×n=[xij]I is 1,2, …, m, j is 1,2, …, n, where the ith row represents the value of the temperature extremum x at location i and the jth column represents the value of the temperature extremum x at time j; then, for the matrix Xm×nPerforming orthogonal decomposition, wherein the feature vector of the orthogonal decomposition corresponds to a space function S capable of reflecting the space distribution characteristics of the variable, the main component of the orthogonal decomposition corresponds to a time function T capable of reflecting the change of the space mode along with the time, and the orthogonal decomposition is expressed as Xm×nSxt ═ sxt; the jth spatial field of the temperature extreme value is represented as linear superposition of m typical spatial fields according to different time coefficients, wherein each column of a spatial function S represents a typical spatial field, each row of a time function T represents a time coefficient, and when a research variable is a summer temperature extreme value of 30 years in the northern hemisphere area, a first feature vector obtained by orthogonal decomposition by using the EOF method is a feature field with the 30-year summer temperature extreme value most similar to a flat field, so that the extremely high-temperature and low-temperature area in summer in the northern hemisphere area is identified;
c, respectively calculating linear trends of the temperature extreme value time sequence according to different time dimensions of every year, years, summer and winter by a least square regression method, and evaluating the statistical significance of the calculated linear trends of the temperature extreme value time sequence on the significance level of 0.05 by a nonparametric Mann-Kendall test method;
d, determining the relationship between the temperature extreme value data and the annular flow field data as well as the large-scale annular flow mode index data by adopting a Pearson correlation analysis method;
and e, quantifying the statistical significance of the annular flow field data and the large-scale circulation mode index data based on a factor analysis method, respectively identifying factors which contribute to the temperature extreme values, and reflecting the interactive relationship among the factors, namely quantifying the relationship between a single annular flow field or a single large-scale circulation mode and the temperature extreme values, the influence of the interaction among a plurality of annular flow fields on the temperature extreme values, and the influence of the interaction of a plurality of large-scale circulation modes or positive and negative phases thereof on the temperature extreme values.
Drawings
Fig. 1 is a flowchart of a method for evaluating a change in an extreme temperature value according to an embodiment of the present invention.
Detailed Description
The following detailed description of embodiments of the invention refers to the accompanying drawings.
The invention provides a method for evaluating temperature extreme value change, which comprises the following steps:
step a, respectively acquiring temperature extreme value data, circulation field data and large-scale circulation mode index data; the definition of the temperature extreme value data is formulated by a world weather organization WMO climate change detection and index expert group ETCCDI, and is obtained by a land climate extreme event index data set HadEX3 covering the land climate extreme event data set with the precision of 1.875 multiplied by 1.25 longitude and latitude from 1901 to 2018, and specific indexes comprise a daily maximum air temperature value TXx, a daily minimum air temperature value TNx, a daily maximum air temperature value TXn, a daily minimum air temperature value TNn, a cold night TN10p, a warm night TN90p, a cold day TX10p and a hot day TX90 p; the circulation field data comprises sea level pressure SLP, 500hPa wind field and 500hPa potential height of ERA5 reanalysis data; the large-scale circulation pattern data are indexes for measuring a natural variability pattern, are time sequence data with a scale of days and comprise index items of Erleno and southern billow ENSO, Pacific decade billow PDO, North Atlantic billow NAO, Arctic billow AO and Pacific-North Atlantic related PNA;
b, utilizing an empirical orthogonal function EOF method to respectively perform different time dimensions of every year, years, summer and winterAnalyzing the temperature extreme value from the time sequence and the spatial scale, and identifying the regions with extremely high temperature and low temperature in the northern hemisphere region in summer; the specific processing procedure of the EOF method is as follows: setting the number of lattice points of the temperature extreme value as m and the time sequence length as n, firstly, processing the temperature extreme value data into a distance form to obtain a matrix X related to the temperature extreme value time and spacem×n=[xij]I is 1,2, …, m, j is 1,2, …, n, where the ith row represents the value of the temperature extremum x at location i and the jth column represents the value of the temperature extremum x at time j; then, for the matrix Xm×nPerforming orthogonal decomposition, wherein the feature vector of the orthogonal decomposition corresponds to a space function S capable of reflecting the space distribution characteristics of the variable, the main component of the orthogonal decomposition corresponds to a time function T capable of reflecting the change of the space mode along with the time, and the orthogonal decomposition is expressed as Xm×nSxt ═ sxt; the jth spatial field of the temperature extreme value is represented as linear superposition of m typical spatial fields according to different time coefficients, wherein each column of a spatial function S represents a typical spatial field, each row of a time function T represents a time coefficient, and when a research variable is a summer temperature extreme value of 30 years in the northern hemisphere area, a first feature vector obtained by orthogonal decomposition by using the EOF method is a feature field with the 30-year summer temperature extreme value most similar to a flat field, so that the extremely high-temperature and low-temperature area in summer in the northern hemisphere area is identified;
c, respectively calculating linear trends of the temperature extreme value time sequence according to different time dimensions of every year, years, summer and winter by a least square regression method, and evaluating the statistical significance of the calculated linear trends of the temperature extreme value time sequence on the significance level of 0.05 by a nonparametric Mann-Kendall test method;
d, determining the relationship between the temperature extreme value data and the annular flow field data as well as the large-scale annular flow mode index data by adopting a Pearson correlation analysis method;
and e, quantifying the statistical significance of the annular flow field data and the large-scale circulation mode index data based on a factor analysis method, respectively identifying factors which contribute to the temperature extreme values, and reflecting the interactive relationship among the factors, namely quantifying the relationship between a single annular flow field or a single large-scale circulation mode and the temperature extreme values, the influence of the interaction among a plurality of annular flow fields on the temperature extreme values, and the influence of the interaction of a plurality of large-scale circulation modes or positive and negative phases thereof on the temperature extreme values.
In order to make the objects, features and advantages of the present invention more comprehensible, the present invention is described in detail with reference to the accompanying drawings and examples.
Example (b): as shown in fig. 1, the present invention provides a method for evaluating the change in the temperature extremes. Firstly, acquiring extreme value index data related to temperature through a high-precision lattice point land climate extreme event index data set of HadEX 3; analyzing data again through ERA5 to obtain circulation field data such as 500hPa wind field, 500hPa potential altitude, sea level pressure and the like; large-scale circulation mode index data such as ENSO, PDO, PNA, AO, NAO and the like are obtained through the national weather administration. Secondly, performing EOF decomposition on the data, and further performing trend change analysis of the developed temperature extreme value according to minimum two-component regression and Mann-Kendall test based on different time scales and different seasons; extracting corresponding annular flow field and large-scale annular flow mode data according to the time of the temperature extreme value, and identifying factors of a single annular flow field and large-scale annular flow mode indexes which influence the temperature extreme value through Pearson correlation analysis; further, factors which significantly contribute to the temperature extreme value are identified based on a factor analysis method, and the influence of the interaction effect between the key factors on the temperature extreme value is quantified. The invention can simultaneously consider single and multiple circulation fields and large-scale circulation mode indexes, and quantize the significance of the influence of the interaction between the single and multiple circulation fields on the change of the temperature extreme value, thereby disclosing the mechanism causing the trend change of the temperature extreme value.
Compared with the prior art, the invention has the beneficial effects that: the method for analyzing the temperature extreme values and the influence mechanisms of the temperature extreme values in different time scales is provided, and the influence of interaction of a plurality of influence factors on the temperature extreme values is quantified simultaneously.

Claims (1)

1. A method for evaluating a change in an extreme temperature value, comprising the steps of:
step a, respectively acquiring temperature extreme value data, circulation field data and large-scale circulation mode index data; the definition of the temperature extreme value data is formulated by a world weather organization WMO climate change detection and index expert group ETCCDI, and is obtained by a land climate extreme event index data set HadEX3 covering the land climate extreme event data set with the precision of 1.875 multiplied by 1.25 longitude and latitude from 1901 to 2018, and specific indexes comprise a daily maximum air temperature value TXx, a daily minimum air temperature value TNx, a daily maximum air temperature value TXn, a daily minimum air temperature value TNn, a cold night TN10p, a warm night TN90p, a cold day TX10p and a hot day TX90 p; the circulation field data comprises sea level pressure SLP, 500hPa wind field and 500hPa potential height of ERA5 reanalysis data; the large-scale circulation pattern data are indexes for measuring a natural variability pattern, are time sequence data with a scale of days and comprise index items of Erleno and southern billow ENSO, Pacific decade billow PDO, North Atlantic billow NAO, Arctic billow AO and Pacific-North Atlantic related PNA;
b, analyzing temperature extreme values from a time sequence and a space scale respectively according to different time dimensions of each year, years, summer and winter by using an empirical orthogonal function EOF method, and identifying areas with extremely high temperature and low temperature in summer in the northern hemisphere area; the specific processing procedure of the EOF method is as follows: setting the number of lattice points of the temperature extreme value as m and the time sequence length as n, firstly, processing the temperature extreme value data into a distance form to obtain a matrix X related to the temperature extreme value time and spacem×n=[xij]I is 1,2, …, m, j is 1,2, …, n, where the ith row represents the value of the temperature extremum x at location i and the jth column represents the value of the temperature extremum x at time j; then, for the matrix Xm×nPerforming orthogonal decomposition, wherein the feature vector of the orthogonal decomposition corresponds to a space function S capable of reflecting the space distribution characteristics of the variable, the main component of the orthogonal decomposition corresponds to a time function T capable of reflecting the change of the space mode along with the time, and the orthogonal decomposition is expressed as Xm×nSxt ═ sxt; the jth spatial field of the temperature extremum is represented as m typical spatial fields according to different time coefficientsLinear superposition, wherein each column of the space function S represents a typical space field, each line of the time function T represents a time coefficient, and when the research variable is a summer temperature extreme value of the northern hemisphere area for 30 years, a first feature vector obtained by orthogonal decomposition by using the EOF method is a feature field with the 30-year summer temperature extreme value most similar to a flat field, so that the area with extremely high temperature and low temperature in summer in the northern hemisphere area is identified;
c, respectively calculating linear trends of the temperature extreme value time sequence according to different time dimensions of every year, years, summer and winter by a least square regression method, and evaluating the statistical significance of the calculated linear trends of the temperature extreme value time sequence on the significance level of 0.05 by a nonparametric Mann-Kendall test method;
d, determining the relationship between the temperature extreme value data and the annular flow field data as well as the large-scale annular flow mode index data by adopting a Pearson correlation analysis method;
and e, quantifying the statistical significance of the annular flow field data and the large-scale circulation mode index data based on a factor analysis method, respectively identifying factors which contribute to the temperature extreme values, and reflecting the interactive relationship among the factors, namely quantifying the relationship between a single annular flow field or a single large-scale circulation mode and the temperature extreme values, the influence of the interaction among a plurality of annular flow fields on the temperature extreme values, and the influence of the interaction of a plurality of large-scale circulation modes or positive and negative phases thereof on the temperature extreme values.
CN202111132485.6A 2021-09-26 2021-09-26 Method for evaluating temperature extremum change Active CN113837620B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111132485.6A CN113837620B (en) 2021-09-26 2021-09-26 Method for evaluating temperature extremum change

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111132485.6A CN113837620B (en) 2021-09-26 2021-09-26 Method for evaluating temperature extremum change

Publications (2)

Publication Number Publication Date
CN113837620A true CN113837620A (en) 2021-12-24
CN113837620B CN113837620B (en) 2024-02-02

Family

ID=78970317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111132485.6A Active CN113837620B (en) 2021-09-26 2021-09-26 Method for evaluating temperature extremum change

Country Status (1)

Country Link
CN (1) CN113837620B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114462247A (en) * 2022-02-14 2022-05-10 中国人民解放军61540部队 Method and system for identifying annual modal associations of surface salinity of North Pacific ocean

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011088473A2 (en) * 2010-01-18 2011-07-21 The Regents Of The University Of California System and method for identifying patterns in and/or predicting extreme climate events
WO2018081559A1 (en) * 2016-10-27 2018-05-03 Ohio University Air data probe
CN110058328A (en) * 2019-01-30 2019-07-26 沈阳区域气候中心 Summer Precipitation in Northeast China multi-mode combines NO emissions reduction prediction technique
CN111856621A (en) * 2020-07-17 2020-10-30 中国气象科学研究院 Integrated evolution SVD conversion method based on mode and observation data fusion

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011088473A2 (en) * 2010-01-18 2011-07-21 The Regents Of The University Of California System and method for identifying patterns in and/or predicting extreme climate events
WO2018081559A1 (en) * 2016-10-27 2018-05-03 Ohio University Air data probe
CN110058328A (en) * 2019-01-30 2019-07-26 沈阳区域气候中心 Summer Precipitation in Northeast China multi-mode combines NO emissions reduction prediction technique
CN111856621A (en) * 2020-07-17 2020-10-30 中国气象科学研究院 Integrated evolution SVD conversion method based on mode and observation data fusion

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DAN ZHANG等: "Spatiotemporal Variations of Extreme Precipitation Events in the Jinsha River Basin, Southwestern China", 《ADVANCES IN METEOROLOGY》, vol. 2020, pages 1 - 13 *
ZHAOFEI LIU等: "Evaluation of Extreme Cold and Drought over the Mongolian Plateau", 《WATER》, vol. 11, no. 74, pages 1 - 17 *
尹红等: "基于ETCCDI 指数2017 年中国极端温度和降水特征分析", 《气候变化研究进展》, vol. 15, no. 4, pages 363 - 373 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114462247A (en) * 2022-02-14 2022-05-10 中国人民解放军61540部队 Method and system for identifying annual modal associations of surface salinity of North Pacific ocean
CN114462247B (en) * 2022-02-14 2022-10-21 中国人民解放军61540部队 Method and system for identifying annual representative modality of sea surface salinity of North Pacific ocean

Also Published As

Publication number Publication date
CN113837620B (en) 2024-02-02

Similar Documents

Publication Publication Date Title
Zhang et al. Using species distribution modeling to delineate the botanical richness patterns and phytogeographical regions of China
Anderson-Frey et al. Investigation of near-storm environments for tornado events and warnings
MacKenzie et al. Quantifying environmental heterogeneity: habitat size necessary for successful development of cod Gadus morhua eggs in the Baltic Sea
Lobo et al. Spatial and environmental determinants of vascular plant species richness distribution in the Iberian Peninsula and Balearic Islands
Correa-Metrio et al. Pollen distribution along climatic and biogeographic gradients in northern Central America
Mackey Special paper: a spatial analysis of the environmental relations of rainforest structural types
Musa et al. A climate distribution model of malaria transmission in Sudan
CN113837620A (en) Method for evaluating temperature extreme value change
CN111737321A (en) Urban atmospheric pollution joint defense joint control area division method based on data mining
Lai et al. Reconciling debates on the controls on surface urban heat island intensity: Effects of scale and sampling
Srinivas Regionalization of precipitation in India–a review
Zampieri et al. Simple statistical approach for computing land cover types and potential natural vegetation
Li et al. Surface urban heat islands in 932 urban region agglomerations in China during the morning and before midnight: spatial-temporal changes, drivers, and simulation
Buechling et al. Predictive habitat modeling of rare plant species in Pacific Northwest forests
Koolagudi Long-range prediction of Indian summer monsoon rainfall using data mining and statistical approaches
Wolfe et al. Population genetics of Penstemon albomarginatus (Plantaginaceae), a rare Mojave Desert species of conservation concern
CN113190633B (en) Climate space partitioning method based on double distances
CN107863150B (en) Web-based haze weather crowd health exposure reaction relation analysis method
CN113610436A (en) Disaster-bearing body dynamic vulnerability assessment method and system
Xie et al. Regional forecasting of PM2. 5 concentrations: A novel model based on the empirical orthogonal function analysis and Nadaraya–Watson kernel regression estimator
CN115688544B (en) Microwave attenuation snowfield chromatography reconstruction method based on MLP neural network
CN116645011B (en) Quantitative index calculation method for evaluating building climate partition performance
Huang et al. The Factors Influencing China’s Population Distribution and Spatial Heterogeneity: Based on Multi-Source Remote Sensing Data
Chen et al. RBF neural network modeling based on PCA clustering analysis
CN116738159B (en) Global ionosphere space weather response extraction method based on complex network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant