Thursday, August 31, 2023

Manitoba: March Precipitation Analysis

Based on the three precipitation trend charts, there is a significant decrease in March precipitation, especially in the last 50 years (1975-2023). The most notable annual decrease in precipitation is observed in the northeastern region.

Wednesday, August 30, 2023

Manitoba: Feburary Precipitation Analysis

Between 1925 and 1974, only a small area in the southwest of Manitoba showed a significant increase in precipitation. However, from 1975 to 2023, there has been a notable decrease in precipitation in large portions of the northeastern and central areas. Throughout the period from 1925 to 2023, there has been a trend of decreasing precipitation near The Pas and in the northern areas of Winnipeg.

Manitoba: January Precipitation Analysis

Thursday, August 24, 2023

Manitoba: Analysis for the Month of July

Matlab: Render a single-band file and export it as an image example

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
clear;
clc;

% Source Directory
SouDir='input';
% Destination Directory
DesDir='output';

% Manitoba Boundary
manitshap='D:\MB.Temp.Prec.20230805\files\boundary.shp';
% All Files
files=dir(fullfile(SouDir, '*.tif'));

% Loop to Export figures
for ii=1:numel(files)
    filepath=fullfile(files(ii).folder, files(ii).name);
    [~, fname, ~]=fileparts(filepath);

    yymm=strsplit(fname, '_');
    yymm=yymm{2};
    yymm=[yymm(1:4), '-', yymm(5:6)];

    otpath=fullfile(DesDir, [fname, '.png']);

    [A, R]=readgeoraster(filepath, "OutputType","double");
    latlim=R.LatitudeLimits;
    lonlim=R.LongitudeLimits;
    worldmap(latlim, lonlim);
    setm(gca, 'mapprojection', 'mercator');
    geoshow(A, R,"DisplayType","surface");
    % setm(gca,"Frame","on","Grid","on")
    plabel off;
    setm(gca,"PLineLocation", 3);
    mlabel off;
    setm(gca,"MLineLocation", 3);
    cb=colorbar('east');
    clim([-1.346 25.911]);
    colorbar('off');
    colormap('jet');

    hold on;
    mani=readgeotable(manitshap);
    kk=geoshow(mani, 'FaceColor', 'none');
    geoshow(49.8954, -97.1385, "Marker","*");
    textm(49.8954, -97.1385,"  Winnipeg");



    textm((max(latlim)-min(latlim))*0.93+min(latlim), ...
        (max(lonlim)-min(lonlim))*0.02+min(lonlim), yymm, ...
        "Color", "black","FontSize", 30);


    %set(gcf,'outerposition',get(0,'screensize'));
    %export_fig(otpath, '-png', '-painters', '-r300')
    %pause(3);
    export_fig(otpath, '-png', '-r400','-q300');
    close all;


end

disp('****************************************************');

Tuesday, August 15, 2023

Manitoba: Analysis for the Month of May

For May temperatures in Manitoba, between 1915 and 1964, only a very small area on the northeastern edge showed a significant cooling trend. From 1915 to 2023, the area experiencing cooling shifted to a central location, while a warming trend appeared in the eastern region.

Monday, August 14, 2023

Manitoba: Analysis for the Month of April

Over the continuous span of 109 years, the temperature in April has been higher in the southwest and southern parts.
The April temperature trend is still divided into three periods. Unlike the previous three months, between 1915 and 1964, a very small part of the northeastern edge of Manitoba exhibited a significant cooling trend. From 1965 to 2023, the small cooling area shifted to the central part. Throughout the period from 1915 to 2023, a large region in the east still shows a warming trend.

Manitoba: Analysis for the Month of March

The temperature in March over the continuous span of 109 years has mostly presented a pattern of being higher in the southwest and lower in the northeast. Among this, the March temperature in 2023 was relatively low compared to the past ten years.
The March temperature trend is also divided into three periods. Between 1915 and 1964, there were no temperature trends in March that passed the significance level. Between 1965 and 2023, only the area near Winnipeg in the entire Manitoba had some significant warming trends, and the vast majority of areas had no noticeable temperature trends. Looking at the full period from 1915 to 2023, the warming area significantly expanded, covering almost all of Manitoba.

Saturday, August 12, 2023

Manitoba: Analysis for the Month of February

As can be seen from the following figures, the February temperature in Manitoba has generally maintained a pattern of being higher in the south and lower in the north over 109 years. Upon closer examination of February 2023, unusually low temperatures were observed in the northeastern region, a situation not seen in many years.
For the calculation of the February temperature's time trend, we have also divided it into three periods. In the first period from 1915 to 1964, there is absolutely no time trend that reaches a significant level. In the other two periods, although there are regions with warming that passes the significance level, the warming areas between the two periods do not overlap, and their warming magnitude and range are both smaller than in January.

Wednesday, August 9, 2023

Manitoba: Analysis for the Month of January

The weather observation data from the Manitoba meteorological station sites have been used as the input, and the local spatial distribution of temperature data for over a hundred continuous months has been obtained through ANUSPLIN interpolation. This is the foundation of our data analysis in this series.
This is the first part of the series, a January analysis, with a time range from 1915 to 2023. Temperatures are represented with deep red for higher temperatures and deep blue for lower temperatures. As shown in the following chart, over each decade, the frequency of warmer January temperatures has notably increased compared to cooler ones.
The temperature in January has been continuous for a total of 109 years, and our research divides them into three periods to calculate and identify their time series trends. From the calculated results, there is absolutely no significant time series trend in the study area between 1915 and 1964. In the other two periods, as shown in the following chart, the vast majority of areas show a significant warming trend. It was also found that the warming trend between 1965 and 2023 is noticeably higher than between 1915 and 2023.

Sunday, August 6, 2023

ANUSPLIN: Configuration of Interpolation Range

The original elevation has a spatial resolution of 0.0041666665°, with a total of 3124 columns and 2641 rows, as shown in the figure below:
Observing the above figure, we can see that the longitude range is from -102.0083363186 to -88.9916701726, and the latitude range is from 48.99791691515 to 60.00208314165. If the cmd file of the lapgrd process directly uses the mentioned range as the interpolation range, it will result in an error. The reason for this error is that:
1
2
-102.0083363186+3124*0.0041666665=-88.9917
48.99791691515+2641*0.0041666665=60.0021
So, the range in the cmd file of lapgrd should be represented as follows:
1
2
-102.0083363186 -88.9917 0.0041666665
48.99791691515 60.0021 0.0041666665
If the error message 'GETLIM - LIMITS EXCEED SURFACE COORDINATE LIMITS' is encountered, the root cause of the issue may be that the range set in the 'splina' command is smaller than the actual DEM (Digital Elevation Model) range. To resolve this, simply adjust the range in the 'splina' command to be slightly larger to cover the entire extent of the DEM. This should allow the interpolation process to proceed without encountering the error.






Share another successful case study:
Lapgrd:
1
2
0.29531831E+06 0.17793583E+07 30
0.24045270E+07 0.40906770E+07 30
When the spatial extent of the lapgrd's cmd file is difficult to determine, consider specifying the spatial range precisely in the splina's cmd file, and use the .sur file obtained from splina's computation to determine the spatial extent of lapgrd.