ABSTRACT Coverage, resolution, and accuracy in the spatial and temporal estimates of remotely sensed precipitation from space satellites, along with the number of instruments deployed to deliver these observations, are increasing. Of key interest in this study is the unsurpassed opportunity offered by their broad and continuous coverage to complement sparse, but more accurate, in situ rain gauge measurements for building climate resilience at the local, regional, and global scales. For many parts of the globe, this opportunity remains untapped. Australia is no exception and provides a unique challenge, given the small fraction of the continent that has rain gauges, the highly diverse climate due to its large size, and the apparent worsening of extreme weather events in both frequency and intensity. Notwithstanding this great impetus, a continent-wide record of multi-satellite-gauge fused precipitation data for Australia remains lacking. This missing data asset is a prerequisite for understanding the emerging spatiotemporal dynamics of precipitation, without which reliable forecasts would be difficult if not impossible. This study seeks to address this need. Here, we develop a method which can synergistically fuse precipitation data from different sources. More specifically, the aim of this study is to develop Precipitation Profiler-Observation Fusion and Estimation (PPrOFusE), a tool to deliver high-quality gauge and multi-satellite fused precipitation data. We test and apply this tool for Australia, but it is by no means limited in scope to this region. By design, PPrOFusE has a built-in capability to assess the strengths and weaknesses of each platform. In this case study, we fused data for a period of 22 years (2000–2022) using the rain gauge network data from the Australian Bureau of Meteorology (BOM) together with satellite data from the Japan Aerospace Exploration Agency’s (JAXA) Global Satellite Mapping of Precipitation (GSMaP) and National Oceanic and Atmospheric Administration’s (NOAA) Climate Prediction Center Morphing technique (CMORPH). Our proposed precipitation data fusion method consists of two steps. Step 1, the relationship among the three sources of data is modeled by multiple linear regression at each rain gauge location, returning the least squares estimates for the associated regression coefficient vector. Step 2, such regression coefficient vector estimates for all rain gauge locations are fitted by a spatial autoregression model, whereafter the multiple linear regression coefficient vectors for those locations void of rain gauges are predicted by spatial interpolation. Key findings are twofold. First, CMORPH is more accurate for most regions of Australia than GSMaP. Second, a clustering analysis of the fused precipitation over the last 20 years suggests two key trends on Australia’s changing climate, relative to BOM’s six major climate zones, from the previous century: (a) increased spatial variability to the north, consistent with meteorological expectations, amid a southwards expansion of the wet summer dominant zones across the continent; (b) the edge of the arid region shifts southwards and pushes out Mediterranean climate and winter dominant rainfall zones across southern Australia.